Synthesis and Characterization of Two Novel Naphthalenediimide/Zinc Phosphonate Crystalline Materials Precipitated from Different Solvents

Hybrid naphthalenediimide/zinc phosphonate materials (NDI/Zn) were prepared by mixing solutions of N,N'-bis(2-phosphonoethyl)-1,4,5,8-naphthalenediimide (PNDI) and zinc nitrate, resulting in the precipitation of the desired compounds. Samples precipitated from water and N,N-dimethylformamide (DMF) were produced. The obtained samples had the expected elemental composition, and the presence of naphthalenediimides (NDI) was ascertained by infrared and UV-visible spectroscopy. All the samples were crystalline, according to powder X-ray diffraction. Nitrogen adsorption isotherms showed the presence of porosity in the NDI/Zn samples. Mesopores with a diameter = 4.1 nm were present in the sample from DMF, with total pore volume reaching 0.13 cm3/g.

Hydrogen bonding patterns and C-H...π interactions in the structure of the antiparkinsonian drug (R)-rasagiline mesylate determined using laboratory and synchrotron X-ray powder diffraction data

The structure of (R)-rasagiline mesylate [(R)-RasH+·Mes-], an active pharmaceutical ingredient used to treat Parkinson's disease, is presented. The structure was determined from laboratory and synchrotron powder diffraction data, refined using the Rietveld method, and validated and optimized using dispersion-corrected DFT calculations. The unit-cell parameters obtained in both experiments are in good agreement and the refinement with both datasets converged to good agreement factors. The final parameters obtained from laboratory data were a = 5.4905 (8), b = 6.536 (2), c = 38.953 (3) Å, V = 1398.0 (4) Å3 and from synchrotron powder data were a = 5.487530 (10) Å, b = 6.528939 (12) Å, c = 38.94313 (9) Å, V = 1395.245 (5) Å3 with Z = 4 and space group P212121. Preferred orientation was properly accounted for using the synchrotron radiation data, leading to a March-Dollase parameter of 1.140 (1) instead of the 0.642 (1) value obtained from laboratory data. In the structure, (R)-RasH+ moieties form layers parallel to the ab plane connected by mesylate ions through N-H...O and C-H...O hydrogen bonds. These layers stack along the c axis and are further connected by C-H...π interactions. Hirshfeld surface analysis and fingerprint plot calculations indicate that the main interactions are: H...H (50.9%), H...C/C...H (27.1%) and H...O/O...H (21.1%).

Predicting crystal form stability under real-world conditions

The physicochemical properties of molecular crystals, such as solubility, stability, compactability, melting behaviour and bioavailability, depend on their crystal form1. In silico crystal form selection has recently come much closer to realization because of the development of accurate and affordable free-energy calculations2-4. Here we redefine the state of the art, primarily by improving the accuracy of free-energy calculations, constructing a reliable experimental benchmark for solid-solid free-energy differences, quantifying statistical errors for the computed free energies and placing both hydrate crystal structures of different stoichiometries and anhydrate crystal structures on the same energy landscape, with defined error bars, as a function of temperature and relative humidity. The calculated free energies have standard errors of 1-2 kJ mol-1 for industrially relevant compounds, and the method to place crystal structures with different hydrate stoichiometries on the same energy landscape can be extended to other multi-component systems, including solvates. These contributions reduce the gap between the needs of the experimentalist and the capabilities of modern computational tools, transforming crystal structure prediction into a more reliable and actionable procedure that can be used in combination with experimental evidence to direct crystal form selection and establish control5.

Crystal structure from X-ray powder diffraction data, DFT-D calculation, Hirshfeld surface analysis, and energy frameworks of (RS)-trichlormethiazide

The structure of the racemic form of the diuretic drug trichlorme­thia­zide was determined from laboratory X-ray powder diffraction data: the extended structure features an intricate combination of N—H⋯O hydrogen bonds and π–π and C—Cl⋯π inter­actions.

The structure of racemic (RS)-trichlorme­thia­zide [systematic name: (RS)-6-chloro-3-(di­chloro­meth­yl)-1,1-dioxo-3,4-di­hydro-2H-1λ⁶,2,4-benzo­thia­di­azine-7-sulfonamide], C₈H₈Cl₃N₃O₄S₂ (RS-TCMZ), a diuretic drug used in the treatment of oedema and hypertension, was determined from laboratory X-ray powder diffraction data using DASH [David et al. (2006 ▸). J. Appl. Cryst.39, 910–915.], refined by the Rietveld method with TOPAS-Academic [Coelho (2018 ▸). J. Appl. Cryst.51, 210–218], and optimized using DFT-D calculations. The extended structure consists of head-to-tail dimers connected by π–π inter­actions which, in turn, are connected by C—Cl⋯π inter­actions. They form chains propagating along [101], further connected by N—H⋯O hydrogen bonds to produce layers parallel to the ac plane that stack along the b-axis direction, connected by additional N—H⋯O hydrogen bonds. The Hirshfeld surface analysis indicates a major contribution of H⋯O and H⋯Cl inter­actions (32.2 and 21.7%, respectively). Energy framework calculations confirm the major contribution of electrostatic inter­actions (E_elec) to the total energy (E_tot). A comparison with the structure of S-TCMZ is also presented.

Structural analysis of metastable pharmaceutical loratadine form II, by 3D electron diffraction and DFT+D energy minimisation

Coupling 3D electron diffraction and density functional theory provided the metastable pharmaceutical crystal structure within nanometre range, under ambient conditions.

A jumping crystal predicted with molecular dynamics and analysed with TLS refinement against powder diffraction data

By running a temperature series of molecular dynamics (MD) simulations starting from the known low-temperature phase, the experimentally observed phase transition in a 'jumping crystal' was captured, thereby providing a prediction of the unknown crystal structure of the high-temperature phase and clarifying the phase-transition mechanism. The phase transition is accompanied by a discontinuity in two of the unit-cell parameters. The structure of the high-temperature phase is very similar to that of the low-temperature phase. The anisotropic displacement parameters calculated from the MD simulations readily identified libration as the driving force behind the phase transition. Both the predicted crystal structure and the phase-transition mechanism were verified experimentally using TLS (translation, libration, screw) refinement against X-ray powder diffraction data.

Four interpenetrating hydrogen-bonded three-dimensional networks in divanillin

The crystal structure of divainillin (systematic name: 6,6'-dihydroxy-5,5'-dimethoxy-[1,1'-biphenyl]-3,3'-dicarbaldehyde), C₁₆H₁₄O₆, was determined from laboratory powder X-ray diffraction data using the software EXPO2013 (direct methods) and WinPSSP (direct-space approach). Divanillin molecules crystallize in the orthorhombic space group Pba2 (No. 32), with two molecules per unit cell (Z' = 1/2). Each divanillin molecule, with twofold symmetry, is linked through strong alcohol-aldehyde hydrogen bonds to four equivalent molecules, defining a three-dimensional hydrogen-bonding network, with rings made up of six divanillin units (a diamond-like arrangement). Each molecule is also connected through π-π interactions to a translation-equivalent molecule along c. Four consecutive molecules stacked along [001] belong to four different three-dimensional hydrogen-bonding networks defining a quadruple array of interpenetrating networks. This complex hydrogen-bonding array is proposed as an explanation for the aging process experienced by divanillin powders.

Crystal structure analysis of a star-shaped triazine compound: a combination of single-crystal three-dimensional electron diffraction and powder X-ray diffraction

The solid-state structure of star-shaped 2,4,6-tris{(E)-2-[4-(dimethylamino)-phenyl]ethenyl}-1,3,5-triazine is determined from a powder sample by exploiting the respective strengths of single-crystal three-dimensional electron diffraction and powder X-ray diffraction data. The unit-cell parameters were determined from single crystal electron diffraction data. Using this information, the powder X-ray diffraction data were indexed, and the crystal structure was determined from the powder diffraction profile. The compound crystallizes in a noncentrosymmetric space group, P2₁2₁2₁. The molecular conformation in the crystal structure was used to calculate the molecular dipole moment of 3.22 Debye, which enables the material to show nonlinear optical effects.

How many ritonavir cases are there still out there?

The number of dormant ritonavir cases is estimated based on 41 commercial pharmaceutical crystal structure prediction studies.

Validation of missed space-group symmetry in X-ray powder diffraction structures with dispersion-corrected density functional theory

More than 600 molecular crystal structures with correct, incorrect and uncertain space-group symmetry were energy-minimized with dispersion-corrected density functional theory (DFT-D, PBE-D3). For the purpose of determining the correct space-group symmetry the required tolerance on the atomic coordinates of all non-H atoms is established to be 0.2 Å. For 98.5% of 200 molecular crystal structures published with missed symmetry, the correct space group is identified; there are no false positives. Very small, very symmetrical molecules can end up in artificially high space groups upon energy minimization, although this is easily detected through visual inspection. If the space group of a crystal structure determined from powder diffraction data is ambiguous, energy minimization with DFT-D provides a fast and reliable method to select the correct space group.

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.

Crystallographic and Dynamic Aspects of Solid-State NMR Calibration Compounds: Towards ab Initio NMR Crystallography

The excellent results of dispersion-corrected density functional theory (DFT-D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT-D calculations is a target, especially for the field of molecular NMR crystallography. Four (13) C ss-NMR calibration compounds are investigated by single-crystal X-ray diffraction, molecular dynamics and DFT-D calculations. The crystal structure of 3-methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated (13) C chemical shifts of these compounds are in excellent agreement with experiment, with a root-mean-square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT-D chemical shift calculation improves the accuracy of calculated chemical shifts.

Crystallographic and Dynamic Aspects of Solid-State NMR Calibration Compounds: Towards ab Initio NMR Crystallography

The excellent results of dispersion-corrected density functional theory (DFT-D) calculations for static systems have been well established over the past decade. The introduction of dynamics into DFT-D calculations is a target, especially for the field of molecular NMR crystallography. Four (13) C ss-NMR calibration compounds are investigated by single-crystal X-ray diffraction, molecular dynamics and DFT-D calculations. The crystal structure of 3-methylglutaric acid is reported. The rotator phases of adamantane and hexamethylbenzene at room temperature are successfully reproduced in the molecular dynamics simulations. The calculated (13) C chemical shifts of these compounds are in excellent agreement with experiment, with a root-mean-square deviation of 2.0 ppm. It is confirmed that a combination of classical molecular dynamics and DFT-D chemical shift calculation improves the accuracy of calculated chemical shifts.

Local structure in the disordered solid solution ofcis- andtrans-perinones

Thecis- andtrans-isomers of the polycyclic aromatic compound perinone, C26H12N4O2, form a solid solution (Vat Red 14). This solid solution is isotypic to the crystal structures ofcis-perinone (Pigment Red 194) andtrans-perinone (Pigment Orange 34) and exhibits a combined positional and orientational disorder: In the crystal, each molecular position is occupied by either acis- ortrans-perinone molecule, both of which have two possible molecular orientations. The structure ofcis-perinone exhibits a twofold orientational disorder, whereas the structure oftrans-perinone is ordered. The crystal structure of the solid solution was determined by single-crystal X-ray analysis. Extensive lattice-energy minimizations with force-field and DFT-D methods were carried out on combinatorially complete sets of ordered models. For the disordered systems, local structures were calculated, including preferred local arrangements, ordering lengths, and probabilities for the arrangement of neighbouring molecules. The superposition of the atomic positions of all energetically favourable calculated models corresponds well with the experimentally determined crystal structures, explaining not only the atomic positions, but also the site occupancies and anisotropic displacement parameters.

One hydrogen bond does not a separation make, or does it? Resolution of amines by diacetoneketogulonic acid

Diacetoneketogulonic acid was used to separate primary amines from their racemic modifications and the selectivity of the acid was rationalized by lattice energy calculations and analyzing the weak interactions around the captured amines.

Structure of Pigment Yellow 181 dimethylsulfoxide N-methyl-2-pyrrolidone (1:1:1) solvate from XRPD + DFT-D

With only a 2.6 Å resolution laboratory powder diffraction pattern of the θ phase of Pigment Yellow 181 (P.Y. 181) available, crystal-structure solution and Rietveld refinement proved challenging; especially when the crystal structure was shown to be a triclinic dimethylsulfoxide N-methyl-2-pyrrolidone (1:1:1) solvate. The crystal structure, which in principle has 28 possible degrees of freedom, was determined in three stages by a combination of simulated annealing, partial Rietveld refinement with dummy atoms replacing the solvent molecules and further simulated annealing. The θ phase not being of commercial interest, additional experiments were not economically feasible and additional dispersion-corrected density functional theory (DFT-D) calculations were employed to confirm the correctness of the crystal structure. After the correctness of the structure had been ascertained, the bond lengths and valence angles from the DFT-D minimized crystal structure were fed back into the Rietveld refinement as geometrical restraints (`polymorph-dependent restraints') to further improve the details of the crystal structure; the positions of the H atoms were also taken from the DFT-D calculations. The final crystal structure is a layered structure with an elaborate network of hydrogen bonds.

Validation of molecular crystal structures from powder diffraction data with dispersion-corrected density functional theory (DFT-D)

In 2010 we energy-minimized 225 high-quality single-crystal (SX) structures with dispersion-corrected density functional theory (DFT-D) to establish a quantitative benchmark. For the current paper, 215 organic crystal structures determined from X-ray powder diffraction (XRPD) data and published in an IUCr journal were energy-minimized with DFT-D and compared to the SX benchmark. The on average slightly less accurate atomic coordinates of XRPD structures do lead to systematically higher root mean square Cartesian displacement (RMSCD) values upon energy minimization than for SX structures, but the RMSCD value is still a good indicator for the detection of structures that deserve a closer look. The upper RMSCD limit for a correct structure must be increased from 0.25 Å for SX structures to 0.35 Å for XRPD structures; the grey area must be extended from 0.30 to 0.40 Å. Based on the energy minimizations, three structures are re-refined to give more precise atomic coordinates. For six structures our calculations provide the missing positions for the H atoms, for five structures they provide corrected positions for some H atoms. Seven crystal structures showed a minor error for a non-H atom. For five structures the energy minimizations suggest a higher space-group symmetry. For the 225 SX structures, the only deviations observed upon energy minimization were three minor H-atom related issues. Preferred orientation is the most important cause of problems. A preferred-orientation correction is the only correction where the experimental data are modified to fit the model. We conclude that molecular crystal structures determined from powder diffraction data that are published in IUCr journals are of high quality, with less than 4% containing an error in a non-H atom.

In situ monitoring of hydration and dehydration in pharmaceutical solids

This talk will describe a multi-technique study of hydration and dehydration in pharmaceutical solids, using sodium naproxen as a case study. The aim of the work is to establish molecular-level structural understanding of the chemical changes that take place in the solid state as a function of temperature and relative humidity. Dynamic vapour sorption (DVS) analysis on the anhydrous compound carried out as a function of temperature provides a preliminary overview of the solid-form landscape and identifies static conditions to obtain four different hydrate forms [1]. Differences in the sorption and desorption cycles indicates the existence of a polymorphic dihydrate, and the two polymorphs show significant differences in their dehydration behaviour. Crystal structures are established for all phases in the system using either single-crystal or powder X-ray diffraction data, supplemented by dispersion-corrected density functional theory (DFT-D) calculations. The hydration and dehydration processes are monitored by powder X-ray diffraction (PXRD), as a function of relative humidity and temperature, and by variable-temperature solid-state 13C and 23Na NMR. Synchrotron PXRD is applied to the two dihydrate polymorphs to monitor the dehydration processes in approximately real time. The kinetic and structural details of dehydration are established by applying parametric Rietveld refinement [2] to the synchrotron data. This approach adds a structural picture to the kinetic processes. The PXRD studies indicate an essentially continuous dehydration pathway from one of the dihydrate polymorphs to the monohydrate, but a stepped dehydration pathway for the other dihydrate polymorph. The different mechanisms are linked to different degrees of structural similarity, and in particular to the existence of topotactic or non-topotactic transformations between the dihydrate polymorphs and the unique monohydrate and anhydrate phases.

Parametric Rietveld refinement applied to dehydration of sodium naproxen

The identification of the hydration and dehydration transformations of drugs is vital to establish stable pharmaceutical components. Our aim is to monitor the kinetics of hydration and dehydration processes in pharmaceutical solids, and to relate these to the molecular-level crystal structures. One of the primary tools to achieve this is parametric Rietveld refinement. The dehydration of two dihydrate polymorphs of the non-steroidal anti-inflammatory drug (NSAID) sodium naproxen was monitored using synchrotron powder X-ray diffraction measurements at Beamline I711, MAX IV Laboratory, Lund University. Diffraction patterns were measured in the range 300-400 K at 1 K increments. Both polymorphs dehydrate to form identical monohydrate then anhydrous phases. Independent Rietveld refinements were initially performed for each unique phase in order to establish initial values for the parametric refinement. The refinements were performed using TOPAS-Academic [1]. The structures are molecular and both dihydrate polymorphs display pseudosymmetry, thereby requiring an extensive set of restraints. One of the key advantages of the parametric Rietveld refinement is the possibility to introduce algebraic equations that describe the evolution of various parameters [2]. The kinetics of the dehydration processes were monitored using a sigmoid function applied to the scale factors of the various phases (see Figure). The evolution of the unit-cell parameters and atomic displacement parameters were also treated parametrically, and the influence of using multiple coordinate sets (one model per temperature) or a single common coordinate set for each phase was examined. One dihydrate polymorph shows a smooth and continuous transition to the monohydrate, whereas the other polymorph shows an abrupt transition. These differences are linked to the existence of topotactic or non-topotactic chemical transformations between the dihydrate polymorphs and the monohydrate phase.

Interpreting Disordered PXRD Structures using Molecular Dynamics

A central topic in the formulation of solid medicinal products is the identification of a suitable solid form of an active compound to obtain optimal physicochemical properties. To this end, disorder may be important for relevant crystal properties like stability. For example, disorder may account for more than 10% of the crystal volume. A rational approach to solid-form selection is typically based on structural information at atomic resolution. In practice, pharmaceutical compounds are not always well-behaved and especially in the study of polymorphs or compounds with flexible groups it can be challenging to obtain crystals suitable for single-crystal X-ray diffraction. Powder X-ray diffraction (PXRD) is a popular alternative, but it generally requires supplementary information like molecular connectivity in simple cases or computational models to solve larger structures. Computational modeling has come a long way and accurate and reliable structures of pharmaceutically relevant compounds can indeed be obtained using laboratory PXRD measurements and quantum-mechanical calculations [1]. The major limitation of quantum mechanical calculations, however, is that they do not consider time nor temperature but only static structures at zero temperature. Thus, these methods cannot model phenomena related to disorder. The molecular dynamics (MD) method can add temperature as well as time and spatial resolution to a model and has in recent years developed to be a scalable, reliable and increasingly available technique. As more and more groups from academia as well as industry employ MD in their work, the development will increase to gain momentum in the coming years. We use MD in a high-performance setting to study crystal properties that are relevant for pharmaceutical research. Using a combination of models from first principles and MD we are able to study highly disordered structures and polymorphs on the basis of PXRD data.

Computational Pharmaceutical Materials Science: Beyond Static Structures

The five Crystal-Structure Prediction (CSP) Blind Tests have shown that molecular-mechanics force fields are not accurate enough for crystal structure prediction[1]. The first--and only--method to successfully predict all four target crystal structures of one of the CSP Blind Tests was dispersion-corrected Density Functional Theory (DFT-D), and this is what we use for our work. However, quantum-mechanical methods (such as DFT-D), are too slow to allow simulations that include the effects of time and temperature, certainly for the size of molecules that are common in pharmaceutical industry. Including the effects of time and temperature therefore still requires molecular dynamics (MD) with less accurate force fields. In order to combine the accuracy of the successful DFT-D method with the speed of a force field to enable molecular dynamics, our group uses Tailor-Made Force Fields (TMFFs) as described by Neumann[2]. In Neumann's TMFF approach, the force field for each chemical compound of interest is parameterised from scratch against reference data from DFT-D calculations; in other words, the TMFF is fitted to mimic the DFT-D energy potential. Parameterising a dedicated force field for each individual compound requires an investment of several weeks, but has the advantage that the resulting force field is more accurate than a transferable force field. Combining crystal-structure prediction with DFT-D followed by molecular dynamics with a tailor-made force field allows us to calculate e.g. the temperature-dependent unit-cell expansion of each predicted polymorph, as well as possible temperature-dependent disorder. This is relevant for example when comparing the calculated X-ray powder diffraction patterns of the predicted crystal structures against experimental data.

Crystal structure prediction of indomethacin

Indomethacin is a non-steroidal anti-inflammatory and antipyretic agent. Because different packing arrangements of the same drug can greatly affect drug properties such as colours, solubility, stability, melting point, dissolution rate and so forth, it is important to predict its polymorphs. The computational prediction of the stable form will reduce undesirable risks in both clinical trials and manufacturing. Reported polymorphs of indomethacin include α, β, γ, δ, ε, η and ζ [1], of which only the thermodynamically stable form γ and the metastable form α are determined. Density functional theory with dispersion-correction (DFT-D) has been used extensively to study molecular crystal structures[2]. It gives better results with a compromise between the computational cost and accuracy towards the reproduction of molecular crystal structures. In the fourth blind test of crystal structure prediction in 2007, the DFT-D method gave a very successful result that predicted all four structures correctly. Rather than using transferable force fields, a dedicated tailor-made force field (TMFF) parameterised by DFT-D calculations[3] is used for every chemical compound. The force field is used to generate a set of crystal structures and delimit a candidate window for energy ranking. The powder diffraction patterns of predicted polymorphs are calculated to compare with experimental data.

On the correlation between hydrogen bonding and melting points in the inositols

Inositol, 1,2,3,4,5,6-hexahydroxycyclohexane, exists in nine stereoisomers with different crystal structures and melting points. In a previous paper on the relationship between the melting points of the inositols and the hydrogen-bonding patterns in their crystal structures [Simperler et al. (2006 ▶). CrystEngComm 8, 589], it was noted that although all inositol crystal structures known at that time contained 12 hydrogen bonds per molecule, their melting points span a large range of about 170 °C. Our preliminary investigations suggested that the highest melting point must be corrected for the effect of molecular symmetry, and that the three lowest melting points may need to be revised. This prompted a full investigation, with additional experiments on six of the nine inositols. Thirteen new phases were discovered; for all of these their crystal structures were examined. The crystal structures of eight ordered phases could be determined, of which seven were obtained from laboratory X-ray powder diffraction data. Five additional phases turned out to be rotator phases and only their unit cells could be determined. Two previously unknown melting points were measured, as well as most enthalpies of melting. Several previously reported melting points were shown to be solid-to-solid phase transitions or decomposition points. Our experiments have revealed a complex picture of phases, rotator phases and phase transitions, in which a simple correlation between melting points and hydrogen-bonding patterns is not feasible.

Structures of cefradine dihydrate and cefaclor dihydrate from DFT-D calculations

The crystal structure of cefradine dihydrate, C16H19N3O4S·2H2O, is considered in the pharmaceutical sciences to be the epitome of an isolated-site hydrate. The structure from single-crystal X-ray data was described in 1976, but atomic coordinates were not published. The atomic coordinates are determined here by combining the information available from the published single-crystal data with a dispersion-corrected density functional theory (DFT-D) method that has been validated to reproduce molecular crystal structures very accurately. Additional proof for the correctness of the structure comes from comparison with cefaclor dihydrate, C15H14ClN3O4S·2H2O, which is isomorphous and for which more complete single-crystal data are available. H-atom positions have not previously been published for either compound. The DFT-D calculations confirm that both cefradine and cefaclor are present in the zwitterionic form in the two dihydrate structures. A potential ambiguity concerning the orientation of the cyclohexadienyl ring in cefradine dihydrate is also clarified, and on the basis of the calculated energies it is shown that disorder should not be expected at room temperature. The DFT-D methods can be applied to recover full structural data in cases where only partial information is available, and where it may not be possible or desirable to obtain new experimental data.

Reinterpretation of the monohydrate of clarithromycin from X-ray powder diffraction data as a trihydrate

Noguchi, Fujiki, Iwao, Miura & Itai [Acta Cryst. (2012), E68, o667-o668] recently reported the crystal structure of clarithromycin monohydrate from synchrotron X-ray powder diffraction data. Voids in the crystal structure suggested the possible presence of two more water molecules. After successful location of the two additional water molecules, the Rietveld refinement still showed minor problems. These were resolved by noticing that one of the chiral centres in the molecule had been inverted. The corrected crystal structure of clarithromycin trihydrate, refined against the original data, is now reported. Dispersion-corrected density functional theory calculations were used to check the final crystal structure and to position the H atoms.

Complementing high-throughput X-ray powder diffraction data with quantum-chemical calculations: Application to piroxicam form III

High-throughput crystallisation and characterisation platforms provide an efficient means to carry out solid-form screening during the pre-formulation phase. To determine the crystal structures of identified new solid phases, however, usually requires independent crystallisation trials to produce single crystals or bulk samples of sufficient quantity to carry out high-quality X-ray diffraction measurements. This process could be made more efficient by a robust procedure for crystal structure determination directly from high-throughput X-ray powder diffraction (XRPD) data. Quantum-chemical calculations based on dispersion-corrected density functional theory (DFT-D) have now become feasible for typical small organic molecules used as active pharmaceutical ingredients. We demonstrate how these calculations can be applied to complement high-throughput XRPD data by determining the crystal structure of piroxicam form III. These combined experimental/quantum-chemical methods can provide access to reliable structural information in the course of an intensive experimentally based solid-form screening activity or in other circumstances wherein single crystals might never be viable, for example, for polymorphs obtained only during high-energy processing such as spray drying or milling.

Ab-initio crystal structure analysis and refinement approaches of oligo p-benzamides based on electron diffraction data

Ab-initio crystal structure analysis of organic materials from electron diffraction data is presented. The data were collected using the automated electron diffraction tomography (ADT) technique. The structure solution and refinement route is first validated on the basis of the known crystal structure of tri-p-benzamide. The same procedure is then applied to solve the previously unknown crystal structure of tetra-p-benzamide. In the crystal structure of tetra-p-benzamide, an unusual hydrogen-bonding scheme is realised; the hydrogen-bonding scheme is, however, in perfect agreement with solid-state NMR data.

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.

Validation of experimental molecular crystal structures with dispersion-corrected density functional theory calculations

The accuracy of a dispersion-corrected density functional theory method is validated against 241 experimental organic crystal structures from Acta Cryst. Section E.

This paper describes the validation of a dispersion-corrected density functional theory (d-DFT) method for the purpose of assessing the correctness of experimental organic crystal structures and enhancing the information content of purely experimental data. 241 experimental organic crystal structures from the August 2008 issue of Acta Cryst. Section E were energy-minimized in full, including unit-cell parameters. The differences between the experimental and the minimized crystal structures were subjected to statistical analysis. The r.m.s. Cartesian displacement excluding H atoms upon energy minimization with flexible unit-cell parameters is selected as a pertinent indicator of the correctness of a crystal structure. All 241 experimental crystal structures are reproduced very well: the average r.m.s. Cartesian displacement for the 241 crystal structures, including 16 disordered structures, is only 0.095 Å (0.084 Å for the 225 ordered structures). R.m.s. Cartesian displacements above 0.25 Å either indicate incorrect experimental crystal structures or reveal interesting structural features such as exceptionally large temperature effects, incorrectly modelled disorder or symmetry breaking H atoms. After validation, the method is applied to nine examples that are known to be ambiguous or subtly incorrect.