Characterizing Slight Structural Disorder in Solids by Combined Solid-State NMR and First Principles Calculations

Role of Crystal Disorder and Mechanoactivation in Solid-State Stability of Pharmaceuticals

Common energy-intensive processes applied in oral solid dosage development, such as milling, sieving, blending, compaction, etc. generate particles with surface and bulk crystal disorder. An intriguing aspect of the generated crystal disorder is its evolution and repercussion on the physical- and chemical stabilities of drugs. In this review, we firstly examine the existing literature on crystal disorder and its implications on solid-state stability of pharmaceuticals. Secondly, we discuss the key aspects related to the generation and evolution of crystal disorder, dynamics of the disordered/amorphous phase, analytical techniques to measure/quantify them, and approaches to model the disordering propensity from first principles. The main objective of this compilation is to provide special impetus to predict or model the chemical degradation(s) resulting from processing-induced manifestation in bulk solid manufacturing. Finally, a generic workflow is proposed that can be useful to investigate the relevance of crystal disorder on the degradation of pharmaceuticals during stability studies. The present review will cater to the requirements for developing physically- and chemically stable drugs, thereby enabling early and rational decision-making during candidate screening and in assessing degradation risks associated with formulations and processing.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

Local-structure effects on 31P NMR chemical shift tensors in solid state

The effect of the local structure on the ³¹P NMR chemical shift tensor (CST) has been studied experimentally and simulated theoretically using the density functional theory gauge-independent-atomic-orbital approach. It has been shown that the dominating impact comes from a small number of noncovalent interactions between the phosphorus-containing group under question and the atoms of adjacent molecules. These interactions can be unambiguously identified using the Bader analysis of the electronic density. A robust and computationally effective approach designed to attribute a given experimental ³¹P CST to a certain local morphology has been elaborated. This approach can be useful in studies of surfaces, complex molecular systems, and amorphous materials.

Disorder in molecular crystals justified with the help of statistical mechanics: a case of two enantiomer solid solutions

An elegant statistical mechanics approach has been exploited in combination with accurate quantum chemical calculations to justify the disorder in two previously reported racemic solids.

Rationalization of the Color Properties of Fluorescein in the Solid State: A Combined Computational and Experimental Study

Fluorescein is known to exist in three tautomeric forms defined as quinoid, zwitterionic, and lactoid. In the solid state, the quinoid and zwitterionic forms give rise to red and yellow materials, respectively. The lactoid form has not been crystallized pure, although its cocrystal and solvate forms exhibit colors ranging from yellow to green. An explanation for the observed colors of the crystals is found using a combination of UV/Vis spectroscopy and plane-wave DFT calculations. The role of cocrystal coformers in modifying crystal color is also established. Several new crystal structures are determined using a combination of X-ray and electron diffraction, solid-state NMR spectroscopy, and crystal structure prediction (CSP). The protocol presented herein may be used to predict color properties of materials prior to their synthesis.

Combining solid-state NMR spectroscopy with first-principles calculations - a guide to NMR crystallography

Recent advances in the application of first-principles calculations of NMR parameters to periodic systems have resulted in widespread interest in their use to support experimental measurement. Such calculations often play an important role in the emerging field of "NMR crystallography", where NMR spectroscopy is combined with techniques such as diffraction, to aid structure determination. Here, we discuss the current state-of-the-art for combining experiment and calculation in NMR spectroscopy, considering the basic theory behind the computational approaches and their practical application. We consider the issues associated with geometry optimisation and how the effects of temperature may be included in the calculation. The automated prediction of structural candidates and the treatment of disordered and dynamic solids are discussed. Finally, we consider the areas where further development is needed in this field and its potential future impact.

Structural Properties, Order-Disorder Phenomena, and Phase Stability of Orotic Acid Crystal Forms

Orotic acid (OTA) is reported to exist in the anhydrous (AH), monohydrate (Hy1), and dimethyl sulfoxide monosolvate (S_DMSO) forms. In this study we investigate the (de)hydration/desolvation behavior, aiming at an understanding of the elusive structural features of anhydrous OTA by a combination of experimental and computational techniques, namely, thermal analytical methods, gravimetric moisture (de)sorption studies, water activity measurements, X-ray powder diffraction, spectroscopy (vibrational, solid-state NMR), crystal energy landscape, and chemical shift calculations. The Hy1 is a highly stable hydrate, which dissociates above 135 °C and loses only a small part of the water when stored over desiccants (25 °C) for more than one year. In Hy1, orotic acid and water molecules are linked by strong hydrogen bonds in nearly perfectly planar arranged stacked layers. The layers are spaced by 3.1 Å and not linked via hydrogen bonds. Upon dehydration the X-ray powder diffraction and solid-state NMR peaks become broader, indicating some disorder in the anhydrous form. The Hy1 stacking reflection (122) is maintained, suggesting that the OTA molecules are still arranged in stacked layers in the dehydration product. Desolvation of S_DMSO, a nonlayer structure, results in the same AH phase as observed upon dehydrating Hy1. Depending on the desolvation conditions, different levels of order–disorder of layers present in anhydrous OTA are observed, which is also suggested by the computed low energy crystal structures. These structures provide models for stacking faults as intergrowth of different layers is possible. The variability in anhydrate crystals is of practical concern as it affects the moisture dependent stability of AH with respect to hydration.

DFT calculations in the assignment of solid-state NMR and crystal structure elucidation of a lanthanum(iii) complex with dithiocarbamate and phenanthroline

Structure of a novel rare-earth lanthanum(iii) complex resolved by a combination of DFT modelling, NMR spectroscopy, and single crystal XRD.

Testing the limits of NMR crystallography: the case of caffeine–citric acid hydrate

The effects of geometry optimisation on the ability to predict linewidths due to disorder and crystal packing energies is investigated on a previously unreported caffeine citric acid cocrystal system.

Conjugated microporous polymers with dimensionality-controlled heterostructures for green energy devices

Dimensionality for conjugated micro-porous polymers (CMP-nD, n = 0, 1, 2) is proven to be of great importance for tailoring their photophysical properties. Moreover, CMP-nD can further be converted into boron and nitrogen co-doped porous carbons (nDBN, n = 0, 1, 2) with maintained 0D, 1D, and 2D nano-structures and highly efficient catalytic performance.

Fine refinement of solid state structure of racemic form of phospho-tyrosine employing NMR Crystallography approach

We present step by step facets important in NMR Crystallography strategy employing O-phospho-dl-tyrosine as model sample. The significance of three major techniques being components of this approach: solid state NMR (SS NMR), X-ray diffraction of powdered sample (PXRD) and theoretical calculations (Gauge Invariant Projector Augmented Wave; GIPAW) is discussed. Each experimental technique provides different set of structural constraints. From the PXRD measurement the size of the unit cell, space group and roughly refined molecular structure are established. SS NMR provides information about content of crystallographic asymmetric unit, local geometry, molecular motion in the crystal lattice and hydrogen bonding pattern. GIPAW calculations are employed for validation of quality of elucidation and fine refinement of structure. Crystal and molecular structure of O-phospho-dl-tyrosine solved by NMR Crystallography is deposited at Cambridge Crystallographic Data Center under number CCDC 1005924.

Fine refinement of solid-state molecular structures of Leu- and Met-enkephalins by NMR crystallography

This paper presents a methodology that allows the fine refinement of the crystal and molecular structure for compounds for which the data deposited in the crystallographic bases are of poor quality. Such species belong to the group of samples with molecular disorder. In the Cambridge Crystallographic Data Center (CCDC), there are approximately 22,000 deposited structures with an R-factor over 10. The powerful methodology we present employs crystal data for Leu-enkephalin (two crystallographic forms) with R-factor values of 14.0 and 8.9 and for Met-enkephalin (one form) with an R-factor of 10.5. NMR crystallography was employed in testing the X-ray data and the quality of the structure refinement. The GIPAW (gauge invariant projector augmented wave) method was used to optimize the coordinates of the enkephalins and to compute NMR parameters. As we reveal, this complementary approach makes it possible to generate a reasonable set of new coordinates that better correlate to real samples. This methodology is general and can be employed in the study of each compound possessing magnetically active nuclei.

Freezing of Molecular Motions Probed by Cryogenic Magic Angle Spinning NMR

Cryogenic magic angle spinning makes it possible to obtain the NMR spectra of solids at temperatures low enough to freeze out most molecular motions. We have applied cryogenic magic angle spinning NMR to a crystalline small-molecule solid (ibuprofen sodium salt), which displays a variety of molecular dynamics. Magic angle (13)C NMR spectra are shown for a wide range of temperatures, including in the cryogenic regime down to 20 K. The hydrophobic and hydrophilic regions of the molecular structure display different behavior in the cryogenic regime, with the hydrophilic region remaining well-structured, while the hydrophobic region exhibits a broad frozen conformational distribution.

Applications of NMR crystallography to problems in biomineralization: refinement of the crystal structure and 31P solid-state NMR spectral assignment of octacalcium phosphate

By combining X-ray crystallography, first-principles density functional theory calculations, and solid-state nuclear magnetic resonance spectroscopy, we have refined the crystal structure of octacalcium phosphate (OCP), reassigned its (31)P NMR spectrum, and identified an extended hydrogen-bonding network that we propose is critical to the structural stability of OCP. Analogous water networks may be related to the critical role of the hydration state in determining the mechanical properties of bone, as OCP has long been proposed as a precursor phase in bone mineral formation. The approach that we have taken in this paper is broadly applicable to the characterization of crystalline materials in general, but particularly to those incorporating hydrogen that cannot be fully characterized using diffraction techniques.

One-phase crystal disorder in pharmaceutical solids and its implication for solid-state stability

Solid-state disorders of active pharmaceutical ingredients have been characterized by means of X-ray diffraction techniques and solid-state nuclear magnetic resonance spectroscopy. The results determined that the pleuromutilin-derivative, I, displays a unique continuous conformational disorder while retaining its long-range crystalline structure. The propionic acid (PA) version of this compound displayed partial crystalline order and site disorder of PA, depending on the quantity of PA incorporated in the structure. Thus, I is a unique example of one-phase crystalline-amorphous model. Physical and chemical stability data was acquired on these disordered systems and discussed in relation with the characterized disorder present in the crystal systems. Analysis of the results showed that in contrast to phase-separated amorphous, restrained disorders do not influence the stability.

The PAW/GIPAW approach for computing NMR parameters: a new dimension added to NMR study of solids

In 2001, Mauri and Pickard introduced the gauge including projected augmented wave (GIPAW) method that enabled for the first time the calculation of all-electron NMR parameters in solids, i.e. accounting for periodic boundary conditions. The GIPAW method roots in the plane wave pseudopotential formalism of the density functional theory (DFT), and avoids the use of the cluster approximation. This method has undoubtedly revitalized the interest in quantum chemical calculations in the solid-state NMR community. It has quickly evolved and improved so that the calculation of the key components of NMR interactions, namely the shielding and electric field gradient tensors, has now become a routine for most of the common nuclei studied in NMR. Availability of reliable implementations in several software packages (CASTEP, Quantum Espresso, PARATEC) make its usage more and more increasingly popular, maybe indispensable in near future for all material NMR studies. The majority of nuclei of the periodic table have already been investigated by GIPAW, and because of its high accuracy it is quickly becoming an essential tool for interpreting and understanding experimental NMR spectra, providing reliable assignments of the observed resonances to crystallographic sites or enabling a priori prediction of NMR data. The continuous increase of computing power makes ever larger (and thus more realistic) systems amenable to first-principles analysis. In the near future perspectives, as the incorporation of dynamical effects and/or disorder are still at their early developments, these areas will certainly be the prime target.

New perspectives in the PAW/GIPAW approach: J(P-O-Si) coupling constants, antisymmetric parts of shift tensors and NQR predictions

In 2001, Pickard and Mauri implemented the gauge including projected augmented wave (GIPAW) protocol for first-principles calculations of NMR parameters using periodic boundary conditions (chemical shift anisotropy and electric field gradient tensors). In this paper, three potentially interesting perspectives in connection with PAW/GIPAW in solid-state NMR and pure nuclear quadrupole resonance (NQR) are presented: (i) the calculation of J coupling tensors in inorganic solids; (ii) the calculation of the antisymmetric part of chemical shift tensors and (iii) the prediction of (14)N and (35)Cl pure NQR resonances including dynamics. We believe that these topics should open new insights in the combination of GIPAW, NMR/NQR crystallography, temperature effects and dynamics. Points (i), (ii) and (iii) will be illustrated by selected examples: (i) chemical shift tensors and heteronuclear (2)J(P-O-Si) coupling constants in the case of silicophosphates and calcium phosphates [Si(5)O(PO(4))(6), SiP(2)O(7) polymorphs and α-Ca(PO(3))(2)]; (ii) antisymmetric chemical shift tensors in cyclopropene derivatives, C(3)X(4) (X = H, Cl, F) and (iii) (14)N and (35)Cl NQR predictions in the case of RDX (C(3)H(6)N(6)O(6)), β-HMX (C(4)H(8)N(8)O(8)), α-NTO (C(2)H(2)N(4)O(3)) and AlOPCl(6). RDX, β-HMX and α-NTO are explosive compounds.

Powder crystallography by combined crystal structure prediction and high-resolution 1H solid-state NMR spectroscopy

A fast method for crystal structure determination using crystal structure prediction and solid-state (1)H NMR is presented. This technique does not need any prior knowledge except the chemical formula; resonance assignment is not necessary. Starting from an ensemble of predicted crystal structures for powdered thymol, comparison between experimental and calculated (1)H solid-state isotropic NMR chemical shifts is sufficient to determine which predicted structure corresponds to the powder under study. The same approach using proton-proton spin-diffusion data is successful and can be used for cross-validation.

Magnetic resonance tensors in uracil: calculation of 13C, 15N, 17O NMR chemical shifts, 17O and 14N electric field gradients and measurement of 13C and 15N chemical shifts

The experimental (13)C NMR chemical shift components of uracil in the solid state are reported for the first time (to our knowledge), as well as newer data for the (15)N nuclei. These experimental values are supported by extensive calculated data of the (13)C, (15)N and (17)O chemical shielding and (17)O and (14)N electric field gradient (EFG) tensors. In the crystal, uracil forms a number of strong and weak hydrogen bonds, and the effect of these on the (13)C and (15)N chemical shift tensors is studied. This powerful combination of the structural methods and theoretical calculations gives a very detailed view of the strong and weak hydrogen bond formation by this molecule. Good calculated results for the optimized cluster in most cases (except for the EFG values of the (14)N3 and (17)O4 nuclei) certify the accuracy of our optimized coordinates for the hydrogen nuclei. Our reported RMSD values for the calculated chemical shielding and EFG tensors are smaller than those reported previously. In the optimized cluster the 6-311+G** basis set is the optimal one in the chemical shielding and EFG calculations, except for the EFG calculations of the oxygen nuclei, in which the 6-31+G** basis set is the optimal one. The optimal method for the chemical shielding and EFG calculations of the oxygen and nitrogen nuclei is the PW91PW91 method, while for the chemical shielding calculations of the (13)C nuclei the B3LYP method gives the best results.