Computation and NMR crystallography of terbutaline sulfate

Removal of bromophenol blue from polluted water using a novel azo-functionalized magnetic nano-adsorbent

Water pollution from organic dyes poses a serious danger to the environment. In the present work, we report a novel adsorbent (ADFS) based on azo-dye-functionalized superparamagnetic iron oxide nanoparticles (SPIONs) for the removal of the anionic dye bromophenol blue (BPB) from contaminated water. The fabricated SPIONs, azo dye, and ADFS adsorbent were characterized with FTIR and UV-vis absorption spectroscopy, 1HNMR spectroscopy, mass spectrometry, SEM imaging, dynamic light scattering (DLS), zeta potential measurements, vibrating sample magnetometry, thermogravimetric analysis, differential thermal analysis, and X-ray diffraction analysis. DLS measurements showed a particle size of 46.1 and 176.5 nm for the SPIONs and the ADFS, respectively. The adsorbent exhibited an adsorption capacity of 7.43 mg g-1 and followed the pseudo-second-order kinetics model (r2 = 0.9981). The ADFS could efficiently remove BPB from water after stirring for 120 minutes at room temperature and pH 2. The adsorption process was proved to occur via physisorption, as revealed by the Freundlich isotherm (n = 1.82 and KF = 11.5). Thermodynamic studies implied that the adsorption is spontaneous (-8.03 ≤ ΔG ≤ -0.58 kJ mol-1) and enthalpy-driven might take place via van der Waals interactions and/or hydrogen bonding (ΔH = -82.19 kJ mol-1 and ΔS = -0.24 kJ mol-1 K-1).

Zwitterionic or Not? Fast and Reliable Structure Determination by Combining Crystal Structure Prediction and Solid-State NMR

When it comes to crystal structure determination, computational approaches such as Crystal Structure Prediction (CSP) have gained more and more attention since they offer some insight on how atoms and molecules are packed in the solid state, starting from only very basic information without diffraction data. Furthermore, it is well known that the coupling of CSP with solid-state NMR (SSNMR) greatly enhances the performance and the accuracy of the predictive method, leading to the so-called CSP-NMR crystallography (CSP-NMRX). In this paper, we present the successful application of CSP-NMRX to determine the crystal structure of three structural isomers of pyridine dicarboxylic acid, namely quinolinic, dipicolinic and dinicotinic acids, which can be in a zwitterionic form, or not, in the solid state. In a first step, mono- and bidimensional SSNMR spectra, i.e., ¹H Magic-Angle Spinning (MAS), ¹³C and ¹⁵N Cross Polarisation Magic-Angle Spinning (CPMAS), ¹H Double Quantum (DQ) MAS, ¹H-¹³C HETeronuclear CORrelation (HETCOR), were used to determine the correct molecular structure (i.e., zwitterionic or not) and the local molecular arrangement; at the end, the RMSEs between experimental and computed ¹H and ¹³C chemical shifts allowed the selection of the correct predicted structure for each system. Interestingly, while quinolinic and dipicolinic acids are zwitterionic and non-zwitterionic, respectively, in the solid state, dinicotinic acid exhibits in its crystal structure a "zwitterionic-non-zwitterionic continuum state" in which the proton is shared between the carboxylic moiety and the pyridinic nitrogen. Very refined SSNMR experiments were carried out, i.e., ¹⁴N-¹H Phase-Modulated (PM) pulse and Rotational-Echo Saturation-Pulse Double-Resonance (RESPDOR), to provide an accurate N-H distance value confirming the hybrid nature of the molecule. The CSP-NMRX method showed a remarkable match between the selected structures and the experimental ones. The correct molecular input provided by SSNMR reduced the number of CSP calculations to be performed, leading to different predicted structures, while RMSEs provided an independent parameter with respect to the computed energy for the selection of the best candidate.

Exploring the Potential of Multinuclear Solid-State 1 H, 13 C, and 35 Cl Magnetic Resonance To Characterize Static and Dynamic Disorder in Pharmaceutical Hydrochlorides

Crystallographic disorder, whether static or dynamic, can be detrimental to the physical and chemical stability, ease of crystallization and dissolution rate of an active pharmaceutical ingredient. Disorder can result in a loss of manufacturing control leading to batch-to-batch variability and can lengthen the process of structural characterization. The range of NMR active nuclei makes solid-state NMR a unique technique for gaining nucleus-specific information about crystallographic disorder. Here, we explore the use of high-field ³⁵ Cl solid-state NMR at 23.5 T to characterize both static and dynamic crystallographic disorder: specifically, dynamic disorder occurring in duloxetine hydrochloride (1), static disorder in promethazine hydrochloride (2), and trifluoperazine dihydrochloride (3). In all structures, the presence of crystallographic disorder was confirmed by ¹³ C cross-polarization magic-angle spinning (CPMAS) NMR and supported by GIPAW-DFT calculations, and in the case of 3, ¹ H solid-state NMR provided additional confirmation. Applying ³⁵ Cl solid-state NMR to these compounds, we show that higher magnetic fields are beneficial for resolving the crystallographic disorder in 1 and 3, while broad spectral features were observed in 2 even at higher fields. Combining the data obtained from ¹ H, ¹³ C, and ³⁵ Cl NMR, we show that 3 exhibits a unique case of disorder involving the ⁺ N-H hydrogen positions of the piperazinium ring, driving the chloride anions to occupy three distinct sites.

Combining X-ray and NMR Crystallography to Explore the Crystallographic Disorder in Salbutamol Oxalate

Salbutamol is an active pharmaceutical ingredient commonly used to treat respiratory distress and is listed by the World Health Organization as an essential medicine. Here, we establish the crystal structure of its oxalate form, salbutamol oxalate, and explore the nature of its crystallographic disorder by combined X-ray crystallography and ¹³C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR. The *C-OH chiral center of salbutamol (note that the crystal structures are a racemic mixture of the two enantiomers of salbutamol) is disordered over two positions, and the tert-butyl group is rotating rapidly, as revealed by ¹³C solid-state NMR. The impact of crystallization conditions on the disorder was investigated, finding variations in the occupancy ratio of the *C-OH chiral center between single crystals and a consistency across samples in the bulk powder. Overall, this work highlights the contrast between investigating crystallographic disorder by X-ray diffraction and solid-state NMR experiment, and gauge-including projector-augmented-wave (GIPAW) density functional theory (DFT) calculations, with their combined use, yielding an improved understanding of the nature of the crystallographic disorder between the local (i.e., as viewed by NMR) and longer-range periodic (i.e., as viewed by diffraction) scale.

Solid-State NMR-Driven Crystal Structure Prediction of Molecular Crystals: The Case of Mebendazole

Among all possible NMR crystallography approaches for crystal-structure determination, crystal structure prediction - NMR crystallography (CSP-NMRX) has recently turned out to be a powerful method. In the latter, the original procedure exploited solid-state NMR (SSNMR) information during the final steps of the prediction. In particular, it used the comparison of computed and experimental chemical shifts for the selection of the correct crystal packing. Still, the prediction procedure, generally carried out with DFT methods, may require important computational resources and be quite time-consuming, especially if there are no available constraints to use at the initial stage. Herein, the successful application of this combined prediction method, which exploits NMR information also in the input step to reduce the search space of the predictive algorithm, is presented. Herein, this method was applied on mebendazole, which is characterized by desmotropism. The use of SSNMR data as constraints for the selection of the right tautomer and the determination of the number of independent molecules in the unit cell led to a considerably faster process, reducing the number of calculations to be performed. In this way, the crystal packing was successfully predicted for the three known phases of mebendazole. To evaluate the quality of the predicted structures, these were compared to the experimental ones. The crystal structure of phase B of mebendazole, in particular, was determined de novo by powder diffraction and is presented for the first time in this paper.

A toolbox for improving the workflow of NMR crystallography

NMR crystallography is a powerful tool with applications in structural characterization and crystal structure verification, to name two. However, applying this tool presents several challenges, especially for industrial users, in terms of consistency, workflow, time consumption, and the requirement for a high level of understanding of experimental solid-state NMR and GIPAW-DFT calculations. Here, we have developed a series of fully parameterized scripts for use in Materials Studio and TopSpin, based on the .magres file format, with a focus on organic molecules (e.g. pharmaceuticals), improving efficiency, robustness, and workflow. We separate these tools into three major categories: performing the DFT calculations, extracting & visualizing the results, and crystallographic modelling. These scripts will rapidly submit fully parameterized CASTEP jobs, extract data from the calculations, assist in visualizing the results, and expedite the process of structural modelling. Accompanied with these tools is a description on their functionality, documentation on how to get started and use the scripts, and links to video tutorials for guiding new users. Through the use of these tools, we hope to facilitate NMR crystallography and to harmonize the process across users.

Bayesian probabilistic assignment of chemical shifts in organic solids

A prerequisite for NMR studies of organic materials is assigning each experimental chemical shift to a set of geometrically equivalent nuclei. Obtaining the assignment experimentally can be challenging and typically requires time-consuming multidimensional correlation experiments. An alternative solution for determining the assignment involves statistical analysis of experimental chemical shift databases, but no such database exists for molecular solids. Here, by combining the Cambridge Structural Database with a machine learning model of chemical shifts, we construct a statistical basis for probabilistic chemical shift assignment of organic crystals by calculating shifts for more than 200,000 compounds, enabling the probabilistic assignment of organic crystals directly from their two-dimensional chemical structure. The approach is demonstrated with the ¹³C and ¹H assignment of 11 molecular solids with experimental shifts and benchmarked on 100 crystals using predicted shifts. The correct assignment was found among the two most probable assignments in more than 80% of cases.

Fast and Accurate Electric Field Gradient Calculations in Molecular Solids With Density Functional Theory

Modern approaches for calculating electric field gradient (EFF) tensors in molecular solids rely upon plane-wave calculations employing periodic boundary conditions (PBC). In practice, models employing PBCs are limited to generalized gradient approximation (GGA) density functionals. Hybrid density functionals applied in the context of gauge-including atomic orbital (GIAO) calculations have been shown to substantially improve the accuracy of predicted NMR parameters. Here we propose an efficient method that effectively combines the benefits of both periodic calculations and single-molecule techniques for predicting electric field gradient tensors in molecular solids. Periodic calculations using plane-wave basis sets were used to model the crystalline environment. We then introduce a molecular correction to the periodic result obtained from a single-molecule calculation performed with a hybrid density functional. Single-molecule calculations performed using hybrid density functionals were found to significantly improve the agreement of predicted 17O quadrupolar coupling constants (C q ) with experiment. We demonstrate a 31% reduction in the RMS error for the predicted 17O C q values relative to standard plane-wave methods using a carefully constructed test set comprised of 22 oxygen-containing molecular crystals. We show comparable improvements in accuracy using five different hybrid density functionals and find predicted C q values to be relatively insensitive to the choice of basis set used in the single molecule calculation. Finally, the utility of high-accuracy 17O C q predictions is demonstrated by examining the disordered 4-Nitrobenzaldehyde crystal structure.

Accurate fragment-based 51-V chemical shift predictions in molecular crystals

Nuclear magnetic resonance (NMR) spectroscopy plays a crucial role in determining molecular structure for complex biological and pharmaceutical compounds. NMR investigations are increasingly reliant on computation for mapping spectral features to chemical structures. Here we benchmark the accuracy of fragment-based ⁵¹V chemical shielding tensor calculations using a training set comprised of 10 biologically and pharmaceutically relevant oxovanadium complexes. Using our self-consistent reproduction of the Madelung potential (SCRMP) electrostatic embedding model, we demonstrate comparable performance between fragment methods and computationally demanding cluster-based techniques. Specifically, fragment methods employing hybrid density functionals are capable of reproducing the experimental ⁵¹V isotropic chemical shifts with a training set rms error of ~9 ​ppm, representing a 20% improvement over traditional plane wave techniques. We provide training set-derived linear regression models for mapping the absolute shieldings obtained from computation to the experimentally determined chemical shifts using four common density functionals; PBE0, B3LYP, PBE, and BLYP. Finally, we establish the utility of fragment methods and the reported regression parameters examining four oxovanadium structures excluded from the training set including the tetracoordinate oxovanadium silicate [Formula: see text] , VO¹⁵NGlySalbz which contains redox-active ligands, and the solid-state form of the common ⁵¹V NMR reference compound VOCl₃.

Structural Refinement of Carbimazole by NMR Crystallography

The characterization of the three-dimensional structure of solids is of major importance, especially in the pharmaceutical field. In the present work, NMR crystallography methods are applied with the aim to refine the crystal structure of carbimazole, an active pharmaceutical ingredient used for the treatment of hyperthyroidism and Grave’s disease. Starting from previously reported X-ray diffraction data, two refined structures were obtained by geometry optimization methods. Experimental ¹H and ¹³C isotropic chemical shift measured by the suitable ¹H and ¹³C high-resolution solid state NMR techniques were compared with DFT-GIPAW calculated values, allowing the quality of the obtained structure to be experimentally checked. The refined structure was further validated through the analysis of ¹H-¹H and ¹H-¹³C 2D NMR correlation experiments. The final structure differs from that previously obtained from X-ray diffraction data mostly for the position of hydrogen atoms.

Structure of Silk I (Bombyx mori Silk Fibroin before Spinning) -Type II β-Turn, Not α-Helix

Recently, considerable attention has been paid to Bombyx mori silk fibroin by a range of scientists from polymer chemists to biomaterial researchers because it has excellent physical properties, such as strength, toughness, and biocompatibility. These appealing physical properties originate from the silk fibroin structure, and therefore, structural determinations of silk fibroin before (silk I) and after (silk II) spinning are a key to make wider applications of silk. There are discrepancies about the silk I structural model, i.e., one is type II β-turn structure determined using many solid-state and solution NMR spectroscopies together with selectively stable isotope-labeled model peptides, but another is α-helix or partially α-helix structure speculated using IR and Raman methods. In this review, firstly, the process that led to type II β-turn structure by the authors was introduced in detail. Then the problems in speculating silk I structure by IR and Raman methods were pointed out together with the problem in the assignment of the amide I band in the spectra. It has been emphasized that the conformational analyses of proteins and peptides from IR and Raman studies are not straightforward and should be very careful when the proteins contain β-turn structure using many experimental data by Vass et al. In conclusion, the author emphasized here that silk I structure should be type II β-turn, not α-helix.

Parametrizing the Spatial Dependence of 1H NMR Chemical Shifts in π-Stacked Molecular Fragments

Most recently a renewed interest in several areas has arisen in factors governing the 1H NMR chemical shift (1H CS) of protons in aromatic systems. Therefore, it is important to describe how 1H CS values are affected by π-stacking intermolecular interactions. The parametrization of radial and angular dependences of the 1H CS is proposed, which is based on conventional gauge-independent atomic orbital (GIAO) calculations of explicit molecular fragments. Such a parametrization is exemplified for a benzene dimer with intermonomer vertical and horizontal distances which are in the range of values often found in crystals of organic compounds. Results obtained by the GIAO calculations combined with B3LYP and MP2 methods were compared, and revealed qualitatively the same trends in the 1H CS data. The parametrization was found to be quantitatively correct for the T-shaped benzene dimers, and its limitations were discussed. Parametrized 1H CS surfaces should become useful for providing additional restraints in the search of site-specific information through an analysis of structurally induced 1H CS changes.

Centerband-Only Detection of Exchange NMR with Natural-Abundance Correction Reveals an Expanded Unit Cell in Phenylalanine Crystals

The NMR pulse sequence CODEX (centerband-only detection of exchange) is a widely used method to report on the number of magnetically inequivalent spins that exchange magnetization via spin diffusion. For crystals, this rules out certain symmetries, and the rate of equilibration is sensitive to distances. Here we show that for 13 C CODEX, consideration of natural abundance spins is necessary for crystals of high complexity, demonstrated here with the amino acid phenylalanine. The NMR data rule out the C2 space group that was originally reported for phenylalanine, and are only consistent with a larger unit cell containing eight magnetically inequivalent molecules. Such an expanded cell was recently described based on single crystal data. The large unit cell dictates the use of long spin diffusion times of more than 200 seconds, in order to equilibrate over the entire unit cell volume of 1622 Å3 .

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.

Resolving alternative organic crystal structures using density functional theory and NMR chemical shifts

Alternative (‘repeat’) determinations of organic crystal structures deposited in the Cambridge Structural Database are analysed to characterise the nature and magnitude of the differences between structure solutions obtained by diffraction methods. Of the 3132 structure pairs considered, over 20% exhibited local structural differences exceeding 0.25 Å. In most cases (about 83%), structural optimisation using density functional theory (DFT) resolved the differences. Many of the cases where distinct and chemically significant structural differences remained after optimisation involved differently positioned hydroxyl groups, with obvious implications for the correct description of hydrogen bonding. ¹H and ¹³C chemical shifts from solid-state NMR experiments are proposed as an independent methodology in cases where DFT optimisation fails to resolve discrepancies.

DFT optimisation often resolves conflicting crystal structure determinations, with NMR shifts helping in cases where optimisation diverges to different structures.

NMR Crystallography at Fast Magic-Angle Spinning Frequencies: Application of Novel Recoupling Methods

Chemical characterisation of active pharmaceutical compounds can be challenging, especially when these molecules exhibit tautomeric or desmotropic behaviour. The complexity can increase manyfold if these molecules are not susceptible to crystallisation. Solid-state NMR has been employed effectively for characterising such molecules. However, characterisation of a molecule is just a first step in identifying the differences in the crystalline structure. 1 H solid-state Nuclear Magnetic Resonance (ssNMR) studies on these molecules at fast magic-angle-spinning frequencies can provide a wealth of information and may be used along with ab initio calculations to predict the crystal structure in the absence of X-ray crystallographic studies. In this work, we attempted to use solid-state NMR to measure 1 H - 1 H distances that can be used as restraints for crystal structure calculations. We performed studies on the desmotropic forms of albendazole.

Furosemide's one little hydrogen atom: NMR crystallography structure verification of powdered molecular organics

The potential of NMR crystallography to verify molecular crystal structures deposited in structural databases is evaluated, with two structures of the pharmaceutical furosemide serving as examples. While the structures differ in the placement of one H atom, using this approach, we verify one of the structures in the Cambridge Structural Database using quantitative tools, while establishing that the other structure does not meet the verification criteria.

Ab initio random structure searching of organic molecular solids: assessment and validation against experimental data

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.

Exploring Systematic Discrepancies in DFT Calculations of Chlorine Nuclear Quadrupole Couplings

Previous studies have revealed significant discrepancies between density functional theory (DFT)-calculated and experimental nuclear quadrupolar coupling constants (C_Q) for chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effects of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the quadrupolar coupling calculations, and that functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data.

Structure and physicochemical characterization of a naproxen-picolinamide cocrystal

Naproxen (NPX) is a nonsteroidal anti-inflammatory drug with pain- and fever-relieving properties, currently marketed in the sodium salt form to overcome solubility problems; however, alternative solutions for improving its solubility across all pH values are desirable. NPX is suitable for cocrystal formation, with hydrogen-bonding possibilities via the COOH group. The crystal structure is presented of a 1:1 cocrystal of NPX with picolinamide as a coformer [systematic name: (S)-2-(6-methoxynaphthalen-2-yl)propanoic acid-pyridine-2-carboxamide (1/1), C14H14O3·C6H6N2O]. The pharmaceutically relevant physical properties were investigated and the intrinsic dissolution rate was found to be essentially the same as that of commercial naproxen. An NMR crystallography approach was used to investigate the H-atom positions in the two crystallographically unique COOH-CONH hydrogen-bonded dimers. 1H solid-state NMR distinguished the two carboxyl protons, despite the very similar crystallographic environments. The nature of the hydrogen bonding was confirmed by solid-state NMR and density functional theory calculations.

13C and 19F solid-state NMR and X-ray crystallographic study of halogen-bonded frameworks featuring nitrogen-containing heterocycles

Halogen bonding is a noncovalent interaction between the electrophilic region of a halogen (σ-hole) and an electron donor. We report a crystallographic and structural analysis of halogen-bonded compounds by applying a combined X-ray diffraction (XRD) and solid-state nuclear magnetic resonance (SSNMR) approach. Single-crystal XRD was first used to characterize the halogen-bonded cocrystals formed between two fluorinated halogen-bond donors (1,4-diiodotetrafluorobenzene and 1,3,5-trifluoro-2,4,6-triiodobenzene) and several nitrogen-containing heterocycles (acridine, 1,10-phenanthroline, 2,3,5,6-tetramethylpyrazine, and hexamethylenetetramine). New structures are reported for the following three cocrystals, all in the P2₁/c space group: acridine-1,3,5-trifluoro-2,4,6-triiodobenzene (1/1), C₆F₃I₃·C₁₃H₉N, 1,10-phenanthroline-1,3,5-trifluoro-2,4,6-triiodobenzene (1/1), C₆F₃I₃·C₁₂H₈N₂, and 2,3,5,6-tetramethylpyrazine-1,3,5-trifluoro-2,4,6-triiodobenzene (1/1), C₆F₃I₃·C₈H₁₂N₂. ¹³C and ¹⁹F solid-state magic-angle spinning (MAS) NMR is shown to be a convenient method to characterize the structural features of the halogen-bond donor and acceptor, with chemical shifts attributable to cocrystal formation observed in the spectra of both nuclides. Cross polarization (CP) from ¹⁹F to ¹³C results in improved spectral sensitivity in characterizing the perfluorinated halogen-bond donor when compared to conventional ¹H CP. Gauge-including projector-augmented wave density functional theory (GIPAW DFT) calculations of magnetic shielding constants, along with optimization of the XRD structures, provide a final set of structures in best agreement with the experimental ¹³C and ¹⁹F chemical shifts. Data for carbons bonded to iodine remain outliers due to well-known relativistic effects.

Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions

Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state ¹³C and ¹⁵N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.

Glycerin-Induced Conformational Changes in Bombyx mori Silk Fibroin Film Monitored by (13)C CP/MAS NMR and ¹H DQMAS NMR

In order to improve the stiff and brittle characteristics of pure Bombyx mori (B. mori) silk fibroin (SF) film in the dry state, glycerin (Glyc) has been used as a plasticizer. However, there have been very limited studies on the structural characterization of the Glyc-blended SF film. In this study, (13)C Cross Polarization/Magic Angle Spinning nuclear magnetic resonance (CP/MAS NMR) was used to monitor the conformational changes in the films by changing the Glyc concentration. The presence of only 5 wt % Glyc in the film induced a significant conformational change in SF where Silk I* (repeated type II β-turn and no α-helix) newly appeared. Upon further increase in Glyc concentration, the percentage of Silk I* increased linearly up to 9 wt % Glyc and then tended to be almost constant (30%). This value (30%) was the same as the fraction of Ala residue within the Silk I* form out of all Ala residues of SF present in B. mori mature silkworm. The ¹H DQMAS NMR spectra of Glyc-blended SF films confirmed the appearance of Silk I* in the Glyc-blended SF film. A structural model of Glyc-SF complex including the Silk I* form was proposed with the guidance of the Molecular Dynamics (MD) simulation using ¹H-¹H distance constraints obtained from the ¹H Double-Quantum Magic Angle Spinning (DQMAS) NMR spectra.

Benchmark fragment-based (1)H, (13)C, (15)N and (17)O chemical shift predictions in molecular crystals

The performance of fragment-based ab initio(1)H, (13)C, (15)N and (17)O chemical shift predictions is assessed against experimental NMR chemical shift data in four benchmark sets of molecular crystals. Employing a variety of commonly used density functionals (PBE0, B3LYP, TPSSh, OPBE, PBE, TPSS), we explore the relative performance of cluster, two-body fragment, and combined cluster/fragment models. The hybrid density functionals (PBE0, B3LYP and TPSSh) generally out-perform their generalized gradient approximation (GGA)-based counterparts. (1)H, (13)C, (15)N, and (17)O isotropic chemical shifts can be predicted with root-mean-square errors of 0.3, 1.5, 4.2, and 9.8 ppm, respectively, using a computationally inexpensive electrostatically embedded two-body PBE0 fragment model. Oxygen chemical shieldings prove particularly sensitive to local many-body effects, and using a combined cluster/fragment model instead of the simple two-body fragment model decreases the root-mean-square errors to 7.6 ppm. These fragment-based model errors compare favorably with GIPAW PBE ones of 0.4, 2.2, 5.4, and 7.2 ppm for the same (1)H, (13)C, (15)N, and (17)O test sets. Using these benchmark calculations, a set of recommended linear regression parameters for mapping between calculated chemical shieldings and observed chemical shifts are provided and their robustness assessed using statistical cross-validation. We demonstrate the utility of these approaches and the reported scaling parameters on applications to 9-tert-butyl anthracene, several histidine co-crystals, benzoic acid and the C-nitrosoarene SnCl2(CH3)2(NODMA)2.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Characterising the role of water in sildenafil citrate by NMR crystallography

A combination of solid-state NMR techniques, including ¹³C/¹H correlation, ²H magic-angle spinning NMR and first principles calculation are employed to characterise the role of water in different hydration states of sildenafil citrate.

Polymorphs of Theophylline Characterized by DNP Enhanced Solid-State NMR

We show how dynamic nuclear polarization (DNP) enhanced solid-state NMR spectroscopy can be used to characterize polymorphs and solvates of organic solids. We applied DNP to three polymorphs and one hydrated form of the asthma drug molecule theophylline. For some forms of theophylline, sample grinding and impregnation with the radical-containing solution, which are necessary to prepare the samples for DNP, were found to induce polymorphic transitions or desolvation between some forms. We present protocols for sample preparation for solid-state magic-angle spinning (MAS) DNP experiments that avoid the polymorphic phase transitions in theophylline. These protocols include cryogrinding, grinding under inert atmosphere, and the appropriate choice of the impregnating liquid. By applying these procedures, we subsequently demonstrate that two-dimensional correlation experiments, such as (1)H-(13)C and (1)H-(15)N HETCOR or (13)C-(13)C INADEQUATE, can be obtained at natural isotopic abundance in reasonable times, thus enabling more advanced structural characterization of polymorphs.

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine

An NMR crystallography study of the hemihydrate of 2', 3'-O-isopropylidineguanosine (Gace) is presented, together with powder X-ray diffraction and thermogravimetric analysis. (1)H double-quantum and (14)N-(1)H HMQC spectra recorded at 850MHz and 75kHz MAS (using a JEOL 1mm probe) are presented together with a (1)H-(13)C refocused INEPT spectrum recorded at 500MHz and 12.5kHz MAS using eDUMBO-122(1)H homonuclear decoupling. NMR chemical shieldings are calculated using the GIPAW (gauge-including projector augmented wave) method; good two-dimensional agreement between calculation and experiment is observed for (13)C and (1)H chemical shifts for directly bonded CH and CH3 peaks. There are two Gace molecules in the asymmetric unit cell: differences in specific (1)H chemical shifts are rationalised in terms of the strength of CH-π and intermolecular hydrogen bonding interactions.

Can the conformation of flexible hydroxyl groups be constrained by simple NMR crystallography approaches? The case of the quercetin solid forms

Hydrogen atoms in systems with many flexible hydroxyl side-groups are difficult to be exactly located from experimental X-Ray diffraction and/or solid-state NMR data, thus often leading to wrong conclusions with respect to the hydrogen bonding network established in crystal lattice. A simple computational method is proposed in the present work to tackle this problem, which may be readily incorporated in conventional NMR crystallography protocols. The method is based on ranking all possible conformations of the flexible hydroxyls according to their lattice energy in crystalline environments. Its effectiveness is investigated on two distinct solid forms of quercetin, for which only two out of the five hydroxyl side-groups can be well constrained from experimental/theoretical data. For this purpose, first-principle quantum-mechanical computations were combined with calculations at the molecular mechanics (MM) level of theory, and previous ss-NMR and X-Ray diffraction data. To assess accuracy in ranking the identified conformers, tests have been performed first on quercetin dihydrate, for which an X-Ray single-crystal structure is available. The possibility of applying this method in a real NMR crystallography context has been investigated finally on anhydrous quercetin, for which only powder X-Ray crystal structure has been reported so far.

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.

De novo determination of the crystal structure of a large drug molecule by crystal structure prediction-based powder NMR crystallography

The crystal structure of form 4 of the drug 4-[4-(2-adamantylcarbamoyl)-5-tert-butyl-pyrazol-1-yl]benzoic acid is determined using a protocol for NMR powder crystallography at natural isotopic abundance combining solid-state (1)H NMR spectroscopy, crystal structure prediction, and density functional theory chemical shift calculations. This is the first example of NMR crystal structure determination for a molecular compound of previously unknown structure, and at 422 g/mol this is the largest compound to which this method has been applied so far.

Exploiting the Synergy of Powder X-ray Diffraction and Solid-State NMR Spectroscopy in Structure Determination of Organic Molecular Solids

We report a strategy for structure determination of organic materials in which complete solid-state nuclear magnetic resonance (NMR) spectral data is utilized within the context of structure determination from powder X-ray diffraction (XRD) data. Following determination of the crystal structure from powder XRD data, first-principles density functional theory-based techniques within the GIPAW approach are exploited to calculate the solid-state NMR data for the structure, followed by careful scrutiny of the agreement with experimental solid-state NMR data. The successful application of this approach is demonstrated by structure determination of the 1:1 cocrystal of indomethacin and nicotinamide. The ¹H and ¹³C chemical shifts calculated for the crystal structure determined from the powder XRD data are in excellent agreement with those measured experimentally, notably including the two-dimensional correlation of ¹H and ¹³C chemical shifts for directly bonded ¹³C-¹H moieties. The key feature of this combined approach is that the quality of the structure determined is assessed both against experimental powder XRD data and against experimental solid-state NMR data, thus providing a very robust validation of the veracity of the structure.

Powder crystallography of pharmaceutical materials by combined crystal structure prediction and solid-state 1H NMR spectroscopy

A protocol for the ab initio crystal structure determination of powdered solids at natural isotopic abundance by combining solid-state NMR spectroscopy, crystal structure prediction, and DFT chemical shift calculations was evaluated to determine the crystal structures of four small drug molecules: cocaine, flutamide, flufenamic acid, and theophylline. For cocaine, flutamide and flufenamic acid, we find that the assigned (1)H isotropic chemical shifts provide sufficient discrimination to determine the correct structures from a set of predicted structures using the root-mean-square deviation (rmsd) between experimentally determined and calculated chemical shifts. In most cases unassigned shifts could not be used to determine the structures. This method requires no prior knowledge of the crystal structure, and was used to determine the correct crystal structure to within an atomic rmsd of less than 0.12 Å with respect to the known reference structure. For theophylline, the NMR spectra are too simple to allow for unambiguous structure selection.