Benchmark accuracy of predicted NMR observables for quadrupolar 14 N and 17 O nuclei in molecular crystals.

Nuclear quadrupole resonances for 14 N and 17 O nuclei are exquisitely sensitive to interactions with surrounding atoms. As a result, nitrogen and oxygen solid-state nuclear magnetic resonance (ssNMR) provides a powerful tool for investigating structure and dynamics in complex systems. First-principles calculations are increasingly used to facilitate spectral assignment and to evaluate and adjust crystal structures. Recent work combining the strengths of planewave density functional theory (DFT) calculations with a single molecule correction obtained using a higher level of theory has proven successful in improving the accuracy of predicted chemical shielding (CS) tensors and 17 O quadrupolar coupling constants ( C q ). Here we extend this work by examining the accuracy of predicted 14 N and 17 O electric field gradient (EFG) tensor components obtained using alternative planewave-corrections involving cluster and two-body fragment-based calculations. We benchmark the accuracy of CS and EFG tensor predictions for both nitrogen and oxygen using planewave, two-body fragment, and enhanced planewave-corrected techniques. Combining planewave and two-body fragment calculations reduces the error in predicted 17 O C q values by 35% relative to traditional planewave calculations. These enhanced planewave-correction methods improve the accuracy of 17 O and 14 N EFG tensor components by 15% relative to planewave DFT but yield minimal improvement relative to a simple molecular correction. However, in structural environments involving either high symmetry or strong intermolecular interactions, enhanced planewave-corrected methods provide a distinct advantage.

13C CPMAS NMR as an Alternative Method to Verify the Quality of Dietary Supplements Containing Curcumin.

Turmeric is a traditional Indian spice that has recently become very popular worldwide because it contains a powerful ingredient called curcumin, which has strong anti-inflammatory properties. Hence, dietary supplements containing extracts rich in curcumin have gained great popularity. The main problems related to curcumin-containing dietary supplements are poor water solubility and the fact that they are often faked by using synthetic curcumin instead of the plant extract. In this article, we propose the use of the 13C CPMAS NMR method to control the quality of dietary supplements. The analysis of 13C CPMAS NMR spectra supported by GIPAW computations allowed us to identify a polymorphic form present in dietary supplements (which affected the solubility of curcumin) and to point out a dietary supplement that could be faked by using synthetic curcumin. Further PXRD and HPLC investigations confirmed that the examined supplement contained synthetic curcumin instead of the genuine extract. Our method can be used for routine control, especially because the investigation is performed directly from the capsule/tablet content and does not require any special sample preparation.

Density Functional Theory and Density Functional Tight Binding Studies of Thiamine Hydrochloride Hydrates.

Thiamine hydrochloride (THCL), also known as vitamin B1, is an active pharmaceutical ingredient (API), present on the list of essential medicines developed by the WHO, which proves its importance for public health. THCL is highly hygroscopic and can occur in the form of hydrates with varying degrees of hydration, depending on the air humidity. Although experimental characterization of the THCL hydrates has been described in the literature, the questions raised in previously published works suggest that additional research and in-depth analysis of THCL dehydration behavior are still needed. Therefore, the main aim of this study was to characterize, by means of quantum chemical calculations, the behavior of thiamine hydrates and explain the previously obtained results, including changes in the NMR spectra, at the molecular level. To achieve this goal, a series of DFT (CASTEP) and DFTB (DFTB+) calculations under periodic boundary conditions have been performed, including molecular dynamics simulations and GIPAW NMR calculations. The obtained results explain the differences in the relative stability of the studied forms and changes in the spectra observed for the samples of various degrees of hydration. This work highlights the application of periodic DFT calculations in the analysis of various solid forms of APIs.

Solid-State Study of the Structure, Dynamics, and Thermal Processes of Safinamide Mesylate─A New Generation Drug for the Treatment of Neurodegenerative Diseases.

Safinamide mesylate (SM), the pure active pharmaceutical ingredient (API) recently used in Parkinson disease treatment, recrystallized employing water-ethanol mixture of solvents (vol/vol 1:9) gives a different crystallographic form compared to SM in Xadago tablets. Pure SM crystallizes as a hemihydrate in the monoclinic system with the P21 space group. Its crystal and molecular structure were determined by means of cryo X-ray crystallography at 100 K. SM in the Xadago tablet exists in anhydrous form in the orthorhombic crystallographic system with the P212121 space group. The water migration and thermal processes in the crystal lattice were monitored by solid-state NMR spectroscopy, differential scanning calorimetry, and thermogravimetric analysis. SM in Xadago in the high-humidity environment undergoes phase transformation to the P21 form which can be easily reversed just by heating up to 80 °C. For the commercial form of the API, there is also a reversible thermal transformation observed between Z' = 1 ↔ Z' = 3 crystallographic forms in the 0-20 °C temperature range. Analysis of molecular motion in the crystal lattice proves that the observed conformational polymorphism is forced by intramolecular dynamics. All above-mentioned processes were analyzed and described employing the NMR crystallography approach with the support of advanced theoretical calculations.

Combining heteronuclear correlation NMR with spin-diffusion to detect relayed Cl-H-H and N-H-H proximities in molecular solids.

Analysis of short-to-intermediate range intermolecular interactions offers a great way of characterizing the solid-state organization of small molecules and materials. This can be achieved by two-dimensional (2D) homo- and heteronuclear correlation NMR spectroscopy, for example, by carrying out experiments at high magnetic fields in conjunction with fast magic-angle spinning (MAS) techniques. But, detecting 2D peaks for heteronuclear dipolar coupled spin pairs separated by greater than 3 Å is not always straightforward, particularly when low-gamma quadrupolar nuclei are involved. Here, we present a 2D correlation NMR experiment that combines the advantages of heteronuclear-multiple quantum coherence (HMQC) and proton-based spin-diffusion (SD) pulse sequences using radio-frequency-driven-recouping (RFDR) to probe inter and intramolecular 1H-X (X = 14N, 35Cl) interactions. This experiment can be used to acquire 2D 1H{X}-HMQC filtered 1H-1H correlation as well as 2D 1H-X HMQC spectra. Powder forms of dopamine·HCl and l-histidine·HCl·H2O are characterized at high fields (21.1 T and 18.8 T) with fast MAS (60 kHz) using the 2D HMQC-SD-RFDR approach. Solid-state NMR results are complemented with NMR crystallography analyses using the gauge-including projector augmented wave (GIPAW) approach. For histidine·HCl·H2O, 2D peaks associated with 14N-1H-1H and 35Cl-1H-1H distances of up to 4.4 and 3.9 Å have been detected. This is further corroborated by the observation of 2D peaks corresponding to 14N-1H-1H and 35Cl-1H-1H distances of up to 4.2 and 3.7 Å in dopamine·HCl, indicating the suitability of the HMQC-SD-RFDR experiments for detecting medium-range proximities in molecular solids.

Improving the accuracy of GIPAW chemical shielding calculations with cluster and fragment corrections.

Ab initio methods for predicting NMR parameters in the solid state are an essential tool for assigning experimental spectra and play an increasingly important role in structural characterizations. Recently, a molecular correction (MC) technique has been developed which combines the strengths of plane-wave methods (GIPAW) with single molecule calculations employing Gaussian basis sets. The GIPAW + MC method relies on a periodic calculation performed at a lower level of theory to model the crystalline environment. The GIPAW result is then corrected using a single molecule calculation performed at a higher level of theory. The success of the GIPAW + MC method in predicting a range of NMR parameters is a result of the highly local character of the tensors underlying the NMR observable. However, in applications involving strong intermolecular interactions we find that expanding the region treated at the higher level of theory more accurately captures local many-body contributions to the N15 NMR chemical shielding (CS) tensor. We propose alternative corrections to GIPAW which capture interactions between adjacent molecules at a higher level of theory using either fragment or cluster-based calculations. Benchmark calculations performed on N15 and C13 data sets show that these advanced GIPAW-corrected calculations improve the accuracy of chemical shielding tensor predictions relative to existing methods. Specifically, cluster-based N15 corrections show a 24% and 17% reduction in RMS error relative to GIPAW and GIPAW + MC calculations, respectively. Comparing the benchmark data sets using multiple computational models demonstrates that N15 CS tensor calculations are significantly more sensitive to intermolecular interactions relative to C13. However, fragment and cluster-based corrections that include direct hydrogen bond partners are sufficient for optimizing the accuracy of GIPAW-corrected methods. Finally, GIPAW-corrected methods are applied to the particularly challenging NMR spectral assignment of guanosine dihydrate which contains two guanosine molecules in the asymmetric unit.

New salts of teriflunomide (TFM) - Single crystal X-ray and solid state NMR investigation.

New salts of teriflunomide TFM (drug approved for Multiple Sclerosis treatment) with inorganic counterions: lithium (TFM_Li), sodium (TFM_Na), potassium (TFM_K), rubidium (TFM_Rb), caesium (TFM_Cs) and ammonium (TFM_NH4) were prepared and investigated employing solid state NMR Spectroscopy, Powder X-ray Diffraction PXRD and Single Crystal X-ray Diffraction (SC XRD). Crystal and molecular structures of three salts: TFM_Na (CCDC: 2173257), TFM_Cs (CCDC: 2165288) and TFM_NH4 (CCDC: 2165281) were determined and deposited. Compared to the native TFM, for all crystalline salt structures, a conformational change of the teriflunomide molecule involving about 180-degree rotation of the end group, forming an intramolecular hydrogen bond N-H⋯O is observed. By applying a complementary multi-technique approach, employing 1D and 2D solid state MAS NMR techniques, single and powder X-ray diffraction measurements, as well as the DFT-based GIPAW calculations of NMR chemical shifts for TFM_Na and TFM_Cs allowed to propose structural features of TFM_Li for which it was not possible to obtain adequate material for single crystal X-Ray measurement.

13C CP MAS NMR and DFT Studies of 6-Chromanyl Ethereal Derivatives.

Vitamin E consists of a group of compounds including α- ß- γ- and δ-tocopherols and α- ß- γ- and δ-tocotrienols, containing the chroman-6-ol system. The recognition of the structural and dynamic properties of this system, present in all vitamers, seems to be important for the full explanation of the mechanism of the biological activity of vitamin E. This paper presents results of the structural analysis of the chosen 6-chromanyl ethereal derivatives using experimental (13 C NMR-in solution and solid state, as well as variable temperature experiments; single crystal X-ray diffraction) and theoretical (DFT) methods. For one of the studied compounds, 2,2,5,7,8-pentamethyl-6-((tetrahydro-2H-pyran-2-yl)oxy) chroman, the splitting of some signals was observed in the 13C dynamic NMR spectra. This observation was explained by the application of a conformational analysis and subsequent DFT optimization, followed by the calculation of NMR properties.

A structural analysis of 2,5-diaryl-4H-2,4-dihydro-3H-1,2,4-triazol-3-ones: NMR in the solid state, X-ray crystallography, and GIPAW calculations.

The 1 H, 13 C, 15 N, and 19 F nuclear magnetic resonance (NMR) spectra of 11 2,5-diaryl-2,4-dihydro-3H-1,2,4-triazol-3-ones have been acquired in DMSO-d6 solution and the 13 C, 15 N, and 19 F NMR spectra have also been acquired in the solid state (solid-state nuclear magnetic resonance [SSNMR] and magic angle spinning [MAS]). The X-ray structures of Compounds 3, 5, and 6 have been determined by X-ray diffraction. Theoretical calculations at the gauge-independent atomic orbital (GIAO)/B3LYP/6-311++G(d,p) level have provided a set of 321 chemical shifts that were compared with 310 experimental values in DMSO-d6 . To obtain good agreements, some effects need to be included. The SSNMR chemical shifts have been compared with gauge-including projector-augmented wave (GIPAW) calculations and with the heavy atom-light atom (HALA) effects.

Looking into the dynamics of molecular crystals of ibuprofen and terephthalic acid using 17 O and 2 H nuclear magnetic resonance analyses.

Oxygen-17 and deuterium are two quadrupolar nuclei that are of interest for studying the structure and dynamics of materials by solid-state nuclear magnetic resonance (NMR). Here, 17 O and 2 H NMR analyses of crystalline ibuprofen and terephthalic acid are reported. First, improved 17 O-labelling protocols of these molecules are described using mechanochemistry. Then, dynamics occurring around the carboxylic groups of ibuprofen are studied considering variable temperature 17 O and 2 H NMR data, as well as computational modelling (including molecular dynamics simulations). More specifically, motions related to the concerted double proton jump and the 180° flip of the H-bonded (-COOH)2 unit in the crystal structure were looked into, and it was found that the merging of the C=O and C-OH 17 O resonances at high temperatures cannot be explained by the sole presence of one of these motions. Lastly, preliminary experiments were performed with a 2 H-17 O diplexer connected to the probe. Such configurations can allow, among others, 2 H and 17 O NMR spectra to be recorded at different temperatures without needing to tune or to change probe configurations. Overall, this work offers a few leads which could be of use in future studies of other materials using 17 O and 2 H NMR.

A Review on Combination of Ab Initio Molecular Dynamics and NMR Parameters Calculations.

This review focuses on a combination of ab initio molecular dynamics (aiMD) and NMR parameters calculations using quantum mechanical methods. The advantages of such an approach in comparison to the commonly applied computations for the structures optimized at 0 K are presented. This article was designed as a convenient overview of the applied parameters such as the aiMD type, DFT functional, time step, or total simulation time, as well as examples of previously studied systems. From the analysis of the published works describing the applications of such combinations, it was concluded that including fast, small-amplitude motions through aiMD has a noticeable effect on the accuracy of NMR parameters calculations.

New Pharmaceutical Salts of Trazodone.

New pharmaceutically acceptable salts of trazodone (trazodone hydrogen bromide and trazodone 1-hydroxy-2-naphthonic acid) for the treatment of central nervous system disorders are synthesized and described. Although trazodone salts are poorly crystalline, single-crystal X-ray diffraction data for trazodone 1-hydroxy-2-naphthonic acid were collected and analyzed as well as compared to the previously described crystal structure of commercially available trazodone hydrochloride. The powder samples of all new salts were characterized by Fourier transform infrared spectroscopy, X-ray diffraction and 13C solid-state nuclear magnetic resonance spectroscopy. Spectroscopic studies were supported by gauge including projector augmented wave (GIPAW) calculations of carbon chemical shielding constants. The main goal of our research was to find salts with better physicochemical properties and to make an attempt to associate them with both the anion structure and the most prominent interactions exhibited by the protonated trazodone cation. The dissolution profiles of trazodone from tablets prepared from various salts with lactose monohydrate were investigated. The studies revealed that salts with simple anions show a fast release of the drug while the presence of more complex anion, more strongly interacting with the cation, effects a slow-release profile of the active substance and can be used for the preparation of the tables with a delay or prolonged mode of action.

Structure Solution of Nano-Crystalline Small Molecules Using MicroED and Solid-State NMR Dipolar-Based Experiments.

Three-dimensional electron diffraction crystallography (microED) can solve structures of sub-micrometer crystals, which are too small for single crystal X-ray crystallography. However, R factors for the microED-based structures are generally high because of dynamic scattering. That means R factor may not be reliable provided that kinetic analysis is used. Consequently, there remains ambiguity to locate hydrogens and to assign nuclei with close atomic numbers, like carbon, nitrogen, and oxygen. Herein, we employed microED and ssNMR dipolar-based experiments together with spin dynamics numerical simulations. The NMR dipolar-based experiments were 1H-14N phase-modulated rotational-echo saturation-pulse double-resonance (PM-S-RESPDOR) and 1H-1H selective recoupling of proton (SERP) experiments. The former examined the dephasing effect of a specific 1H resonance under multiple 1H-14N dipolar couplings. The latter examined the selective polarization transfer between a 1H-1H pair. The structure was solved by microED and then validated by evaluating the agreement between experimental and calculated dipolar-based NMR results. As the measurements were performed on 1H and 14N, the method can be employed for natural abundance samples. Furthermore, the whole validation procedure was conducted at 293 K unlike widely used chemical shift calculation at 0 K using the GIPAW method. This combined method was demonstrated on monoclinic l-histidine.

A GIPAW versus GIAO-ZORA-SO study of 13C and 15N CPMAS NMR chemical shifts of aromatic and heterocyclic bromo derivatives.

Theoretical simulation of NMR parameters in compounds bearing heavy atoms generally requires the application of relativistic corrections. We report herein the theoretical characterization of 13C and 15N CPMAS NMR of known bromo-derivative crystals by using both the GIPAW and the combined GIAO-ZORA-SO approximation methods. Several statistical analyses were performed to compare both approaches, with non-relativistic GIPAW method being more useful to predict the 13C and 15N chemical shifts. The problem of applying GIPAW to crystal structures showing static or dynamic crystalline disorder of the special class resulting in half-protons will be discussed in detail.

5-amino-2-methylpyridinium hydrogen fumarate: An XRD and NMR crystallography analysis.

Single-crystal X-ray diffraction structures of the 5-amino-2-methylpyridinium hydrogen fumarate salt have been solved at 150 and 300 K (CCDC 1952142 and 1952143). A base-acid-base-acid ring is formed through pyridinium-carboxylate and amine-carboxylate hydrogen bonds that hold together chains formed from hydrogen-bonded hydrogen fumarate ions. 1 H and 13 C chemical shifts as well as 14 N shifts that additionally depend on the quadrupolar interaction are determined by experimental magic angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) and gauge-including projector-augmented wave (GIPAW) calculation. Two-dimensional homonuclear 1 H-1 H double-quantum (DQ) MAS and heteronuclear 1 H-13 C and 14 N-1 H spectra are presented. Only small differences of up to 0.1 and 0.6 ppm for 1 H and 13 C are observed between GIPAW calculations starting with the two structures solved at 150 and 300 K (after geometry optimisation of atomic positions, but not unit cell parameters). A comparison of GIPAW-calculated 1 H chemical shifts for isolated molecules and the full crystal structures is indicative of hydrogen bonding strength.

Monitoring the Site-Specific Solid-State NMR Data in Oligopeptides.

Reliable values of the solid-state NMR (SSNMR) parameters together with precise structural data specific for a given amino acid site in an oligopeptide are needed for the proper interpretation of measurements aiming at an understanding of oligopeptides' function. The periodic density functional theory (DFT)-based computations of geometries and SSNMR chemical shielding tensors (CSTs) of solids are shown to be accurate enough to support the SSNMR investigations of suitably chosen models of oriented samples of oligopeptides. This finding is based on a thorough comparison between the DFT and experimental data for a set of tripeptides with both 13Cα and 15Namid CSTs available from the single-crystal SSNMR measurements and covering the three most common secondary structural elements of polypeptides. Thus, the ground is laid for a quantitative description of local spectral parameters of crystalline oligopeptides, as demonstrated for the backbone 15Namid nuclei of samarosporin I, which is a pentadecapeptide (composed of five classical and ten nonproteinogenic amino acids) featuring a strong antimicrobial activity.

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.

Polymorphic Forms of Valinomycin Investigated by NMR Crystallography.

A dodecadepsipeptide valinomycin (VLM) has been most recently reported to be a potential anti-coronavirus drug that could be efficiently produced on a large scale. It is thus of importance to study solid-phase forms of VLM in order to be able to ensure its polymorphic purity in drug formulations. The previously available solid-state NMR (SSNMR) data are combined with the plane-wave DFT computations in the NMR crystallography framework. Structural/spectroscopical predictions (the PBE functional/GIPAW method) are obtained to characterize four polymorphs of VLM. Interactions which confer a conformational stability to VLM molecules in these crystalline forms are described in detail. The way how various structural factors affect the values of SSNMR parameters is thoroughly analyzed, and several SSNMR markers of the respective VLM polymorphs are identified. The markers are connected to hydrogen bonding effects upon the corresponding (13C/15N/1H) isotropic chemical shifts of (CO, Namid, Hamid, Hα) VLM backbone nuclei. These results are expected to be crucial for polymorph control of VLM and in probing its interactions in dosage forms.

Computational and Experimental 1H-NMR Study of Hydrated Mg-Based Minerals.

Magnesium oxide (MgO) can convert to different magnesium-containing compounds depending on exposure and environmental conditions. Many MgO-based phases contain hydrated species allowing 1H-nuclear magnetic resonance (NMR) spectroscopy to be used in the characterization and quantification of proton-containing phases; however, surprisingly limited examples have been reported. Here, 1H-magic angle spinning (MAS) NMR spectra of select Mg-based minerals are presented and assigned. These experimental results are combined with computational NMR density functional theory (DFT) periodic calculations to calibrate the predicted chemical shielding results. This correlation is then used to predict the NMR shielding for a series of different MgO hydroxide, magnesium chloride hydrate, magnesium perchlorate, and magnesium cement compounds to aid in the future assignment of 1H-NMR spectra for complex Mg phases.

Local Structure of Glassy Lithium Phosphorus Oxynitride Thin Films: A Combined Experimental and Ab†Initio Approach.

Lithium phosphorus oxynitride (LiPON) is an amorphous solid-state lithium ion conductor displaying exemplary cyclability against lithium metal anodes. There is no definitive explanation for this stability due to the limited understanding of the structure of LiPON. Herein, we provide a structural model of RF-sputtered LiPON. Information about the short-range structure results from 1D and 2D solid-state NMR experiments. These results are compared with first principles chemical shielding calculations of Li-P-O/N crystals and ab initio molecular dynamics-generated amorphous LiPON models to unequivocally identify the glassy structure as primarily isolated phosphate monomers with N incorporated in both apical and as bridging sites in phosphate dimers. Structural results suggest LiPON's stability is a result of its glassy character. Free-standing LiPON films are produced that exhibit a high degree of flexibility, highlighting the unique mechanical properties of glassy materials.