A combined first principles computational and solid-state NMR study of a molecular crystal: flurbiprofen

On the Use of Solomon Echoes in 27 Al NMR Studies of Complex Aluminium Hydrides

The quadrupole coupling constant CQ and the asymmetry parameter η have been determined for two complex aluminium hydrides from 27 Al NMR spectra recorded for stationary samples by using the Solomon echo sequence. The thus obtained data for KAlH4 (CQ =(1.30±0.02) MHz, η=(0.64±0.02)) and NaAlH4 (CQ =(3.11±0.02) MHz, η<0.01) agree very well with data previously determined from MAS NMR spectra. The accuracy with which these parameters can be determined from static spectra turned out to be at least as good as via the MAS approach. The experimentally determined parameters (δiso , CQ and η) are compared with those obtained from DFT-GIPAW (density functional theory - gauge-including projected augmented wave) calculations. Except for the quadrupole coupling constant for KAlH4 , which is overestimated in the GIPAW calculations by about 30 %, the agreement is excellent. Advantages of the application of the Solomon echo sequence for the measurement of less stable materials or for in situ studies are discussed.

Intramolecular interactions (O-H∙∙∙O, C-H∙∙∙N, N-H∙∙∙π) in isomers of neutral, cation, and anion dopamine molecules: A DFT study on the influence of solvents (water and ethanol)

CONTEXT

Dopamine (DA) is one of the most important neurotransmitters associated with numerous neural disorders. This investigation reports the intramolecular interactions present in the isomers of neutral (DA⁰), anionic (DA^-), and cationic (DA⁺) dopamine isomers in gas, water, and ethanol mediums. Neutral and anion isomers have O-H∙∙∙O, C-H∙∙∙N intramolecular hydrogen bonds and N-H∙∙∙π interactions. All the interactions are electrostatic in nature. Isomers of cation dopamine show no intramolecular interactions in the solvent. Natural charges from natural bond orbital (NBO) analysis show that O-H∙∙∙O bonds and the N-H∙∙∙π interactions are the most and least polar, respectively. ¹H NMR study reveals the inverse linear correlation between shielding constant and electron density in a solvent medium. HOMO-LUMO energy gap indicates higher stability for neutral and cationic forms of dopamine isomers in water and ethanol medium.

METHODS

We have optimized all the structural forms of dopamine molecule using the Becke three hybrid exchange and Lee-Yang-Parr correlation functional with Grimme's dispersion correction, B3LYP-D3(BJ), and aug-cc-pVTZ basis set using the Gaussian16 software. Vibrational frequency analysis with no imaginary frequencies confirms the nature of global minima. The solvent studies (water and ethanol) were carried out using the SCRF keyword and the polarisable continuum model (PCM) of Miertus and Tomasi. NBO analysis and NMR studies were also performed for all conformers. Topology analysis was explored using the software Multiwfn.

GIPAW Pseudopotentials of d Elements for Solid-State NMR

Computational methods are increasingly used to support interpreting, assigning and predicting the solid-state nuclear resonance magnetic spectra of materials. Currently, density functional theory is seen to achieve a good balance between efficiency and accuracy in solid-state chemistry. To be specific, density functional theory allows the assignment of signals in nuclear resonance magnetic spectra to specific sites and can help identify overlapped or missing signals from experimental nuclear resonance magnetic spectra. To avoid the difficulties correlated to all-electron calculations, a gauge including the projected augmented wave method was introduced to calculate nuclear resonance magnetic parameters with great success in organic crystals in the last decades. Thus, we developed a gauge including projected augmented pseudopotentials of 21 d elements and tested them on, respectively, oxides or nitrides (semiconductors), calculating chemical shift and quadrupolar coupling constant. This work can be considered the first step to improving the ab initio prediction of nuclear magnetic resonance parameters, and leaves open the possibility for inorganic compounds to constitute an alternative standard compound, with respect to tetramethylsilane, to calculate the chemical shift. Furthermore, this work represents the possibility to obtain results from first-principles calculations, to train a machine-learning model to solve or refine structures using predicted nuclear magnetic resonance spectra.

Pseudocyclic Form of 4-Hydroxypyrrolidine-2-carboxanilide Podands with Trioxyethylene Chain: Modeling, Conformational Search, and NMR Analysis

The 4-hydroxypyrrolidine-2-carboxanilide podand salt demonstrates catalytic activity in asymmetric Biginelli reaction. The systematic search for prevalent conformational state of the cation was carried out by computer simulations in combination with one- and two-dimensional NMR experiments. For that purpose, we proposed a novel algorithm for the generation and selection of conformers based on molecular dynamics and clustering in the space of principal components. The search had found an important trend of the podand to form a pseudocyclic structure with a horseshoe-shaped conformation of the oligooxyethylene fragment. This conformation is stabilized by different types of intramolecular hydrogen bonds between the acidic and basic centers of the two 4-hydroxypyrrolidine-2-carboxanilide residuals (branches). The proposed approach had made it possible to identify the major structural factors, providing a correlation between the calculated and experimental chemical shifts of hydrogen atoms in the ¹H NMR spectra of the protonated podand.

Automated fragmentation quantum mechanical calculation of 13C and 1H chemical shifts in molecular crystals

In this work, the automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) approach was applied to calculate the ¹³C and ¹H nuclear magnetic resonance (NMR) chemical shifts in molecular crystals. Two benchmark sets of molecular crystals were selected to calculate the NMR chemical shifts. Systematic investigation was conducted to examine the convergence of AF-QM/MM calculations and the impact of various density functionals with different basis sets on the NMR chemical shift prediction. The result demonstrates that the calculated NMR chemical shifts are close to convergence when the distance threshold for the QM region is larger than 3.5 Å. For ¹³C chemical shift calculations, the mPW1PW91 functional is the best density functional among the functionals chosen in this study (namely, B3LYP, B3PW91, M06-2X, M06-L, mPW1PW91, OB98, and OPBE), while the OB98 functional is more suitable for the ¹H NMR chemical shift prediction of molecular crystals. Moreover, with the B3LYP functional, at least a triple-ζ basis set should be utilized to accurately reproduce the experimental ¹³C and ¹H chemical shifts. The employment of diffuse basis functions will further improve the accuracy for ¹³C chemical shift calculations, but not for the ¹H chemical shift prediction. We further proposed a fragmentation scheme of dividing the central molecule into smaller fragments. By comparing with the results of the fragmentation scheme using the entire central molecule as the core region, the AF-QM/MM calculations with the fragmented central molecule can not only achieve accurate results but also reduce the computational cost. Therefore, the AF-QM/MM approach is capable of predicting the ¹³C and ¹H NMR chemical shifts for molecular crystals accurately and effectively, and could be utilized for dealing with more complex periodic systems such as macromolecular polymers and biomacromolecules. The AF-QM/MM program for molecular crystals is available at https://github.com/shiman1995/NMR.

Accurate location of hydrogen atoms in hydrogen bonds of tizoxanide from the combination of experimental and theoretical models

To obtain detailed information about the position of hydrogen atoms in hydrogen bonds, HBs, of crystalline organic molecular compounds is not an easy task. In this work we propose a combination of ssNMR experimental data with theoretical procedures to get such information. Furthermore, the combination of experimental and theoretical models provides us with well-defined grounds to analyse the strength of π-stacking interactions between layers of hydrogen bonded molecules. Two different theoretical models were considered, both approaches being quite different. The first one is a solid-state model, so that the periodicity of a crystalline system underlies calculations of the electronic energy, the electronic density and NMR parameters. The other one is a molecular model in which molecules are taken as isolated monomers, dimers and tetramers. These two models were applied to the tizoxanide, TIZ, molecular crystal though it can widely be applied to any other molecular crystal. By the application of the quantum molecular model it was possible to learn about the way the intermolecular HBs affect the position of hydrogen atoms that belong to HBs in TIZ. This molecule has two intermolecular HBs that stabilize the structure of a basic dimer, but it also has an intramolecular HB in each monomer whose position should be optimized together with the other ones. We found that by doing this it is possible to obtain reliable results of calculations of NMR spectroscopic parameters. Working with the solid-state model we found that any local variation of the TIZ crystalline structure is correlated with the variation of the values of the NMR parameters of each nucleus. The excellent agreement between experimental and calculated chemical shifts leads to the conclusion that the N10-H10 bond distance should be (1.00 ± 0.02) Å.

Identifying aspirin polymorphs from combined DFT-based crystal structure prediction and solid-state NMR

A combined experimental and computational approach was used to distinguish between different polymorphs of the pharmaceutical drug aspirin. This method involves the use of ab initio random structure searching (AIRSS), a density functional theory (DFT)-based crystal structure prediction method for the high-accuracy prediction of polymorphic structures, with DFT calculations of nuclear magnetic resonance (NMR) parameters and solid-state NMR experiments at natural abundance. AIRSS was used to predict the crystal structures of form-I and form-II of aspirin. The root-mean-square deviation between experimental and calculated ¹ H chemical shifts was used to identify form-I as the polymorph present in the experimental sample, the selection being successful despite the large similarities between the molecular environments in the crystals of the two polymorphs.

Investigating discrepancies between experimental solid-state NMR and GIPAW calculation: NC-N 13C and OH⋯O 1H chemical shifts in pyridinium fumarates and their cocrystals

An NMR crystallography analysis is presented for four solid-state structures of pyridine fumarates and their cocrystals, using crystal structures deposited in the Cambridge Crystallographic Data Centre, CCDC. Experimental one-dimensional one-pulse ¹H and ¹³C cross-polarisation (CP) magic-angle spinning (MAS) nuclear magnetic resonance (NMR) and two-dimensional ¹⁴N-¹H heteronuclear multiple-quantum coherence MAS NMR spectra are compared with gauge-including projector augmented wave (GIPAW) calculations of the ¹H and ¹³C chemical shifts and the ¹⁴N shifts that additionally depend on the quadrupolar interaction. Considering the high ppm (>10 ​ppm) ¹H resonances, while there is good agreement (within 0.4 ​ppm) between experiment and GIPAW calculation for the hydrogen-bonded NH moieties, the hydrogen-bonded fumaric acid OH resonances are 1.2-1.9 ​ppm higher in GIPAW calculation as compared to experiment. For the cocrystals of a salt and a salt formed by 2-amino-5-methylpyridinium and 2-amino-6-methylpyridinium ions, a large discrepancy of 4.2 and 5.9 ​ppm between experiment and GIPAW calculation is observed for the quaternary ring carbon ¹³C resonance that is directly bonded to two nitrogens (in the ring and in the amino group). By comparison, there is excellent agreement (within 0.2 ​ppm) for the quaternary ring carbon ¹³C resonance directly bonded to the ring nitrogen for the salt and cocrystal of a salt formed by 2,6-lutidinium and 2,5-lutidinium, respectively.

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.

31P MAS NMR and DFT study of crystalline phosphate matrices

We present the study of the phosphorus local environment by using ³¹P MAS NMR in a series of seven double monophosphates M^IIM^IV(PO₄)₂ (M^II and M^IV being divalent and tetravalent cations, respectively) of yavapaiite and low-yavapaiite type crystal structures. Solid-state and cluster DFT calculations were found to be efficient for predicting the ³¹P isotropic chemical shift and chemical shift anisotropy. To achieve this performance, however, a proper computational optimisation of the experimental structural data was required. From the three optimisation methods tested, the full optimisation provided the best reference structure for the calculation of the NMR parameters of the studied phosphates. Also, a better prediction of the chemical shifts was possible by using a correction to the GIPAW calculated shielding.

Hydrogen bonding networks in gabapentin protic pharmaceutical salts: NMR and in silico studies

Hydrogen bonds (HBs) play a key role in the supramolecular arrangement of crystalline solids and, although they have been extensively studied, the influence of their strength and geometry on crystal packing remains poorly understood. Here we describe the crystal structures of two novel protic gabapentin (GBP) pharmaceutical salts prepared with the coformers methanesulfonic acid (GBP:METHA) and ethanesulfonic acid (GBP:ETHA). This study encompasses experimental and computational electronic structure analyses of ¹ H NMR chemical shifts (CSs), upon in silico HB cleavage. GBP:METHA and GBP:ETHA crystal packing comprise two main structural domains: an ionic layer (characterized by the presence of charge-assisted ⁺ NH_GBP ⋯O^-_METHA/ETHA HB interactions) and a neutral layer generated in a different way for each salt, mainly due to the presence of bifurcated HB interactions. A comprehensive study of HB networks is presented for GBP:METHA, by isolating molecular fragments involved in distinct HB types (NH⋯O, OH⋯O, and CH⋯O) obtained from in silico disassembling of an optimized three-dimensional packing structure. Formation of HB leads to calculated ¹ H NMR CS changes from 0.4 to ~5.8 ppm. This study further attempts to assess how ¹ H NMR CS of protons engaged in certain HB are affected when other nearby HB, involving bifurcated or geminal/vicinal hydrogen atoms, are removed.

14N NQR spectroscopy reveals the proton position in N–H⋯N bonds: a case study with proton sponges

The position of the proton in intramolecular N–H⋯N hydrogen bonds has been determined to a high accuracy with ¹⁴N Nuclear Quadrupole Resonance (NQR) spectroscopy.

An XRD and NMR crystallographic investigation of the structure of 2,6-lutidinium hydrogen fumarate

A crystallographic study highlighting the benefits of a combined XRD and NMR approach in investigating both stability and variation within an organic multicomponent crystal.

Perspective: Current advances in solid-state NMR spectroscopy

In contrast to the rapid and revolutionary impact of solution-state Nuclear Magnetic Resonance (NMR) on modern chemistry, the field of solid-state NMR has matured more slowly. This reflects the major technical challenges of much reduced spectral resolution and sensitivity in solid-state as compared to solution-state spectra, as well as the relative complexity of the solid state. In this perspective, we outline the technique developments that have pushed resolution to intrinsic limits and the approaches, including ongoing major developments in the field of Dynamic Nuclear Polarisation, that have enhanced spectral sensitivity. The information on local structure and dynamics that can be obtained using these gains in sensitivity and resolution is illustrated with a diverse range of examples from large biomolecules to energy materials and pharmaceuticals and from both ordered and highly disordered materials. We discuss how parallel developments in quantum chemical calculation, particularly density functional theory, have enabled experimental data to be translated directly into information on local structure and dynamics, giving rise to the developing field of "NMR crystallography."

Recent Advances in Solid-State Nuclear Magnetic Resonance Spectroscopy

The sensitivity of nuclear magnetic resonance (NMR) spectroscopy to the local atomic-scale environment offers great potential for the characterization of a diverse range of solid materials. Despite offering more information than its solution-state counterpart, solid-state NMR has not yet achieved a similar level of recognition, owing to the anisotropic interactions that broaden the spectral lines and hinder the extraction of structural information. Here, we describe the methods available to improve the resolution of solid-state NMR spectra and the continuing research in this area. We also highlight areas of exciting new and future development, including recent interest in combining experiment with theoretical calculations, the rise of a range of polarization transfer techniques that provide significant sensitivity enhancements, and the progress of in situ measurements. We demonstrate the detailed information available when studying dynamic and disordered solids and discuss the future applications of solid-state NMR spectroscopy across the chemical sciences.

Refinement of labile hydrogen positions based on DFT calculations of 1H NMR chemical shifts: comparison with X-ray and neutron diffraction methods

Numerous gas phase electron diffraction, ultra-fast electron diffraction, X-ray and neutron diffraction experiments on β-dicarbonyl compounds exhibiting enol-enol tautomeric equilibrium, with emphasis on acetylacetone and dibenzoylmethane, have so far been reported with conflicting results on the structural details of the O-HO intramolecular hydrogen bond and resulted in alternative hypotheses on the intramolecular hydrogen bond potential function either a double minimum potential corresponding to two tautomeric forms in equilibrium or a single symmetrical one. We demonstrate herein, firstly, that the DFT calculated OH ¹H NMR chemical shifts of acetylacetone and dibenzoylmethane exhibit a strong linear dependence on the computed OO hydrogen bond length of ∼-50 ppm Å^-1 and as a function of the O-HO bond angle of ∼1 ppm per degree, upon the transfer of the hydrogen atom from the ground state toward the transition state. Secondly, the refinement of labile hydrogen atomic positions in intramolecular hydrogen bonds based on the root-mean-square deviation between experimentally determined and DFT calculated ¹H NMR chemical shifts in solution can provide high resolution structures of O-H and O(H)O bond lengths and O-HO bond angles with an accuracy of ∼10^-2 Å and ∼0.5°, respectively. Thirdly, the calculated ¹H NMR chemical shifts in solution of the two ground state tautomers in equilibrium of acetylacetone and dibenzoylmethane are in excellent agreement with the experimental value, even for moderate basis sets for energy minimization. In contrast, the single symmetrical structure in a strongly delocalized system is a transition state with calculated ¹H NMR chemical shifts which strongly deviate from the experimental value. Fourth, the DFT calculated ground state O-H bond lengths of acetylacetone and dibenzoylmethane are in quantitative agreement with the literature data which take into account the effect of quantum nuclear motion. The DFT structural results are critically discussed with respect to the state-of-the-art variable temperature X-ray and neutron diffraction methods.

1H CSA parameters by ultrafast MAS NMR: Measurement and applications to structure refinement

A ¹H anisotropic-isotropic chemical shift correlation experiment which employs symmetry-based recoupling sequences to reintroduce the chemical shift anisotropy in ν₁ and ultrafast MAS to resolve ¹H sites in ν₂ is described. This experiment is used to measure ¹H shift parameters for L-ascorbic acid, a compound with a relatively complex hydrogen-bonding network in the solid. The ¹H CSAs of hydrogen-bonded sites with resolved isotropic shifts can be extracted directly from the recoupled lineshapes. In combination with DFT calculations, hydrogen positions in crystal structures obtained from X-ray and neutron diffraction are refined by comparison with simulations of the full two-dimensional NMR spectrum. The improved resolution afforded by the second dimension allows even unresolved hydrogen-bonded sites ¹H to be assigned and their shift parameters to be obtained.

NMR and IR Investigations of Strong Intramolecular Hydrogen Bonds

For the purpose of this review, strong hydrogen bonds have been defined on the basis of experimental data, such as OH stretching wavenumbers, νOH, and OH chemical shifts, δOH (in the latter case, after correction for ring current effects). Limits for O-H···Y systems are taken as 2800 > νOH > 1800 cm-1, and 19 ppm > δOH > 15 ppm. Recent results as well as an account of theoretical advances are presented for a series of important classes of compounds such as β-diketone enols, β-thioxoketone enols, Mannich bases, proton sponges, quinoline N-oxides and diacid anions. The O···O distance has long been used as a parameter for hydrogen bond strength in O-H···O systems. On a broad scale, a correlation between OH stretching wavenumbers and O···O distances is observed, as demonstrated experimentally as well as theoretically, but for substituted β-diketone enols this correlation is relatively weak.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

A 13C solid-state NMR investigation of four cocrystals of caffeine and theophylline

We report an analysis of the 13C solid-state NMR chemical shift data in a series of four cocrystals involving two active pharmaceutical ingredient (API) mimics (caffeine and theophylline) and two diacid coformers (malonic acid and glutaric acid). Within this controlled set, we make comparisons of the isotropic chemical shifts and the principal values of the chemical shift tensor. The dispersion at 14.1 T (600 MHz 1H) shows crystallographic splittings in some of the resonances in the magic angle spinning spectra. By comparing the isotropic chemical shifts of individual C atoms across the four cocrystals, we are able to identify pronounced effects on the local electronic structure at some sites. We perform a similar analysis of the principal values of the chemical shift tensors for the anisotropic C atoms (most of the ring C atoms for the API mimics and the carbonyl C atoms of the diacid coformers) and link them to differences in the known crystal structures. We discuss the future prospects for extending this type of study to incorporate the full chemical shift tensor, including its orientation in the crystal frame of reference.

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.

Structure elucidation of a complex CO2-based organic framework material by NMR crystallography

A three-dimensional structural model of a complex CO2-based organic framework made from high molecular weight, self-assembled, flexible and multi-functional oligomeric constituents has been determined de novo by solid-state NMR including DNP-enhanced experiments. The complete assignment of the 15N, 13C and 1H resonances was obtained from a series of two-dimensional through space and through bond correlation experiments. MM-QM calculations were used to generate different model structures for the material which were then evaluated by comparing multiple experimental and calculated NMR parameters. Both NMR and powder X-ray diffraction were evaluated as tools to determine the packing by crystal modelling, and at the level of structural modelling used here PXRD was found not to be a useful complement. The structure determined reveals a highly optimised H-bonding network that explains the unusual selectivity of the self-assembly process which generates the material. The NMR crystallography approach used here should be applicable for the structure determination of other complex solid materials.

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.

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.

Combining Nuclear Magnetic Resonance Spectroscopy and Density Functional Theory Calculations to Characterize Carvedilol Polymorphs

The experiments of carvedilol form II, form III, and hydrate by (13)C and (15)N cross-polarization magic-angle spinning (CP MAS) are reported. The GIPAW (gauge-including projector-augmented wave) method from DFT (density functional theory) calculations was used to simulate (13)C and (15)N chemical shifts. A very good agreement was found for the comparison between the global results of experimental and calculated nuclear magnetic resonance (NMR) chemical shifts for carvedilol polymorphs. This work aims a comprehensive understanding of carvedilol crystalline forms employing solution and solid-state NMR as well as DFT calculations.

Silica capsules enclosing P123 triblock copolymer micelles for flurbiprofen storage and release

Flurbiprofen was incorporated in 200-400 nm silica capsules filled with Pluronic P123 polyethylene oxide-polypropylene oxide-polyethylene oxide triblock copolymer. The assembly process of the capsules and the molecular organization of drug molecule, surfactant and silica were investigated using SAXS, TGA, SEM, DLS, DSC, ¹³C single-pulse, CPMAS and ¹H-¹H two-dimensional NMR. Flurbiprofen molecules are molecularly dispersed inside polypropylene cores of P123 surfactant micelles occluded in a 20-30 nm thick silica shell. Flurbiprofen molecules in polypropylene cores of P123 micelles exhibit high mobility and are easily released after introduction in simulated gastrointestinal fluid and the solubility limit is reached within minutes. Release rates are favored at high pH due to acid dissociation of the carboxylic acid group of the flurbiprofen molecule. The molecular environment of flurbiprofen in these P123 filled silica capsules is different from ordered mesoporous silica materials synthesized using P123 as sacrificial template having the drug molecules adsorbed on the silica surface of pore walls. These findings uncover a new concept of storage and release of hydrophobic bioactive molecules.

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.

The conformational behaviour of naproxen and flurbiprofen in solution by NMR spectroscopy

The conformational equilibrium of common anti-inflammatory drugs has been studied experimentally in solution by NMR in weakly ordered PBLG phases.

Experimental and solid-state computational study of structural and dynamic properties in the equilibrium form of temazepam

Structural properties and rotational dynamics of methyl groups in the most stable form of temazepam were investigated by means of (13)C CP MAS NMR, quasielastic neutron scattering (QENS), and (1)H NMR spin-lattice relaxation methods. The QENS and (1)H NMR studies reveal the inequivalency of methyl groups, delivering their activation parameters. The structural properties of the system were explored in frame of periodic density functional theory (DFT) computations, giving insight into the reorientational barriers and providing understanding of the solid-state NMR results. The theoretical computations are shedding light on the intermolecular interactions along their relation with particular asymmetric structural units.

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.