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

Co-crystal structures of furosemide:urea and carbamazepine:indomethacin determined from powder x-ray diffraction data

Co-crystallization is a promising approach to improving both the solubility and the dissolution rate of active pharmaceutical ingredients. Crystal structure determination from powder diffraction data plays an important role in determining co-crystal structures, especially those generated by mechanochemical means. Here, two new structures of pharmaceutical interest are reported: a 1:1 co‑crystal of furosemide with urea formed by liquid-assisted grinding and a second polymorphic form of a 1:1 co‑crystal of carbamazepine with indomethacin, formed by solvent evaporation. Energy minimization using dispersion-corrected density functional theory was used in finalizing both structures. In the case of carbamazepine:indomethacin, this energy minimization step was essential in obtaining a satisfactory final Rietveld refinement

Spray drying of fenofibrate loaded nanostructured lipid carriers

AbstractThe conversion of aqueous dispersion of nanostructured lipid carriers (NLCs) into dry powder by spray drying could be a useful approach to render NLCs with better physical chemical stability than the aqueous dispersion. In this study, aqueous NLC dispersion containing fenofibrate was converted into dry, easily reconstitutable powder using spray drying. A central composite face centered design (CCFD) was used to investigate the influence of the ratio of lipid to protectant (mannitol and trehalose) and crystallinity of spray-dried powder on the particle size, yield and residual moisture content of the dried powder. A linear relationship (R2 = 0.9915) was established between the crystalline content of the spray-dried powders against the ratio of mannitol to trehalose from 3:7 to 10:0 (w/w). Spray drying of NLC aqueous dispersion using a mannitol and trehalose mixture resulted in an increase in particle size of the NLCs after reconstitution in water as compared to that in the initial aqueous dispersion. The decrease in crystallinity of the dry powder by reducing the ratio of mannitol to trehalose could improve the reconstitution of the NLCs in water. However the yield and residual moisture content of dry powder decreased with an increase in the ratio of mannitol to trehalose. Lipid nanoparticles were able to retain the drug incorporation and the prolonged drug release profile after spray drying. The experimental model was robust, and suggested that spray drying is a viable technique for the conversion of NLCs into dry powder

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[Simperler, A., Watt, S. W., Bonnet, P. A., Jones, W. & Motherwell, W. D. S. (2006). CrystEngComm, 8, 589-600.]). 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

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

Computational polymorph screening reveals late-appearing and poorly-soluble form of rotigotine

The active pharmaceutical ingredient rotigotine—a dopamine agonist for the treatment of Parkinson’s and restless leg diseases—was known to exist in only one polymorphic form since 1985. In 2008, the appearance of a thermodynamically more stable and significantly less soluble polymorph led to a massive batch recall followed by economic and public health implications. Here, we carry out state-of-the-art computational crystal structure prediction, revealing the late-appearing polymorph without using any prior information. In addition, we predict a third crystalline form of rotigotine having thermodynamic stability between forms I and II. We provide quantitative description of the relative stability and solubility of the rotigotine polymorphs. Our study offers new insights into a challenging polymorphic system and highlights the robustness of contemporary computational crystal structure prediction during pharmaceutical development

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

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