Solid-State NMR Characterization of the Molecular Conformation in Disordered Methyl α-l-Rhamnofuranoside

Modeling Small Structural and Environmental Differences in Solids with 15 N NMR Chemical Shift Tensors

The ability to theoretically predict accurate NMR chemical shifts in solids is increasingly important due to the role such shifts play in selecting among proposed model structures. Herein, two theoretical methods are evaluated for their ability to assign ¹⁵ N shifts from guanosine dihydrate to one of the two independent molecules present in the lattice. The NMR data consist of ¹⁵ N shift tensors from 10 resonances. Analysis using periodic boundary or fragment methods consider a benchmark dataset to estimate errors and predict uncertainties of 5.6 and 6.2 ppm, respectively. Despite this high accuracy, only one of the five sites were confidently assigned to a specific molecule of the asymmetric unit. This limitation is not due to negligible differences in experimental data, as most sites exhibit differences of >6.0 ppm between pairs of resonances representing a given position. Instead, the theoretical methods are insufficiently accurate to make assignments at most positions.

Unveiling the Molecular Structure of Antimalarial Drugs Chloroquine and Hydroxychloroquine in Solution through Analysis of 1H NMR Chemical Shifts

Chloroquine (CQ) and hydroxychloroquine (HCQ) have been standard antimalarial drugs since the early 1950s, and very recently, the possibility of their use for the treatment of COVID-19 patients has been considered. To understand the drug mode of action at the submicroscopic level (atoms and molecules), molecular modeling studies with the aid of computational chemistry methods have been of great help. A fundamental step in such theoretical investigations is the knowledge of the predominant drug molecular structure in solution, which is the real environment for the interaction with biological targets. Our strategy to access this valuable information is to perform density functional theory (DFT) calculations of ¹H NMR chemical shifts for several plausible molecular conformers and then find the best match with experimental NMR profile in solution (since it is extremely sensitive to conformational changes). Through this procedure, after optimizing 30 trial distinct molecular structures (ωB97x-D/6-31G(d,p)-PCM level of calculation), which may be considered representative conformations, we concluded that the global minimum (named M24), stabilized by an intramolecular N-H hydrogen bond, is not likely to be observed in water, chloroform, and dimethyl sulfoxide (DMSO) solution. Among fully optimized conformations (named M1 to M30, and MD1 and MD2), we found M12 (having no intramolecular H-bond) as the most probable structure of CQ and HCQ in water solution, which is a good approximate starting geometry in drug-receptor interaction simulations. On the other hand, the preferred CQ and HCQ structure in chloroform (and CQ in DMSO-d₆) solution was assigned as M8, showing the solvent effects on conformational preferences. We believe that the analysis of ¹H NMR data in solution can establish the connection between the macro level (experimental) and the sub-micro level (theoretical), which is not so apparent to us and appears to be more appropriate than the thermodynamic stability criterion in conformational analysis studies.

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.

CO Guest Interactions in SDB-Based Metal-Organic Frameworks: A Solid-State Nuclear Magnetic Resonance Investigation

Metal-organic frameworks (MOFs) are promising materials for greener carbon monoxide (CO) capture and separation processes. SDB-based (SDB = 4,4^'-sulfonyldibenzoate) MOFs are particularly attractive due to their remarkable gas adsorption capacity under humid conditions. However, to the best of our knowledge, their CO adsorption abilities have yet to be investigated. In this report, CO-loaded PbSDB and CdSDB were characterized using variable-temperature (VT) ¹³C solid-state nuclear magnetic resonance (SSNMR) spectroscopy. These MOFs readily captured CO, with the adsorbed CO exhibiting dynamics as indicated by the temperature-dependent changes in the SSNMR spectra. Spectral simulations revealed that the CO simultaneously undergoes a localized wobbling about the adsorption site and a nonlocalized hopping between adjacent adsorption sites. The wobbling and hopping angles were also found to be temperature-dependent. From the appearance of the VT spectra and the extracted motional data, the CO adsorption mechanism was concluded to be analogous to that of CO₂. To gain a better understanding on the gas adsorption properties of these MOFs and their CO capture abilities, we subsequently compared the motional data to those reported for CO₂ in SDB-based MOFs and CO in MOF-74, respectively. A significant contrast in adsorption strength was observed in both cases because of the different physical properties of the guests (i.e., CO vs CO₂) and the MOF frameworks (i.e., SDB-based MOFs vs MOFs with open metal sites). Our results demonstrate that SSNMR spectroscopy can be employed to probe variations in binding behavior.

Accurate 13-C and 15-N molecular crystal chemical shielding tensors from fragment-based electronic structure theory

Standard nuclear magnetic resonance (NMR) spectroscopy experiments measure isotropic chemical shifts, but measuring the chemical shielding anisotropy (CSA) tensor can provide additional insights into solid state chemical structures. Interpreting the principal components of these tensors is facilitated by first-principles chemical shielding tensor predictions. Here, the ability to predict molecular crystal CSA tensor components for ¹³C and ¹⁵N nuclei with fragment-based electronic structure techniques is explored. Similar to what has been found previously for isotropic chemical shifts, the benchmarking demonstrates that fragment-based techniques can accurately reproduce CSA tensor components. The use of hybrid density functionals like PBE0 or B3LYP provide higher accuracy than generalized gradient approximation functionals like PBE. Unlike for planewave density functional techniques, hybrid density functionals can be employed routinely with modest computational cost in fragment approaches. Finally, good consistency between the regression parameters used to map either isotropic shieldings or CSA tensor components is demonstrated, providing further evidence for the quality of the models and highlighting that models trained for isotropic shifts can also be applied to CSA tensor components.

Does Z' equal 1 or 2? Enhanced powder NMR crystallography verification of a disordered room temperature crystal structure of a p38 inhibitor for chronic obstructive pulmonary disease

The crystal structure of the Form A polymorph of N-cyclopropyl-3-fluoro-4-methyl-5-[3-[[1-[2-[2-(methylamino)ethoxy]phenyl]cyclopropyl]amino]-2-oxo-pyrazin-1-yl]benzamide (i.e., AZD7624), determined using single-crystal X-ray diffraction (scXRD) at 100 K, contains two molecules in the asymmetric unit (Z' = 2) and has regions of local static disorder. This substance has been in phase IIa drug development trials for the treatment of chronic obstructive pulmonary disease, a disease which affects over 300 million people and contributes to nearly 3 million deaths annually. While attempting to verify the crystal structure using nuclear magnetic resonance crystallography (NMRX), we measured ¹³C solid-state NMR (SSNMR) spectra at 295 K that appeared consistent with Z' = 1 rather than Z' = 2. To understand this surprising observation, we used multinuclear SSNMR (¹H, ¹³C, ¹⁵N), gauge-including projector augmented-wave density functional theory (GIPAW DFT) calculations, crystal structure prediction (CSP), and powder XRD (pXRD) to determine the room temperature crystal structure. Due to the large size of AZD7624 (ca. 500 amu, 54 distinct ¹³C environments for Z' = 2), static disorder at 100 K, and (as we show) dynamic disorder at ambient temperatures, NMR spectral assignment was a challenge. We introduce a method to enhance confidence in NMR assignments by comparing experimental ¹³C isotropic chemical shifts against site-specific DFT-calculated shift distributions established using CSP-generated crystal structures. The assignment and room temperature NMRX structure determination process also included measurements of ¹³C shift tensors and the observation of residual dipolar coupling between ¹³C and ¹⁴N. CSP generated ca. 90 reasonable candidate structures (Z' = 1 and Z' = 2), which when coupled with GIPAW DFT results, room temperature pXRD, and the assigned SSNMR data, establish Z' = 2 at room temperature. We find that the polymorphic Form A of AZD7624 is maintained at room temperature, although dynamic disorder is present on the NMR timescale. Of the CSP-generated structures, 2 are found to be fully consistent with the SSNMR and pXRD data; within this pair, they are found to be structurally very similar (RMSD₁₆ = 0.30 Å). We establish that the CSP structure in best agreement with the NMR data possesses the highest degree of structural similarity with the scXRD-determined structure (RMSD₁₆ = 0.17 Å), and has the lowest DFT-calculated energy amongst all CSP-generated structures with Z' = 2.

Dynamic Disorder and Electronic Structures of Electron-Precise Dianionic Diboranes: Insights from Solid-State Multinuclear Magnetic Resonance Spectroscopy

The J(¹¹B,¹¹B) coupling constants of various salts of the electron-precise hexacyanodiborane(6) dianion, [B₂(CN)₆]^2-, were obtained using ¹¹B double-quantum-filtered (DQF) J-resolved solid-state nuclear magnetic resonance (SSNMR) spectroscopy. Our results show that the magnitude of the DQF J splitting is influenced by both the crystallographic symmetry of the system and the presence of dynamics. The splittings are amplified by a factor of 3 as compared to the corresponding theoretical J coupling constants for cases where (1) there is an absence of dynamics but the boron pairs are crystallographically equivalent or (2) the boron pairs are crystallographically inequivalent but are rendered magnetically equivalent on the time scale of the experiment due to dynamic disorder, which was identified by ¹¹B and ¹³C SSNMR experiments. Consequently, molecular motions need to be taken into consideration when interpreting the results of DQF J-resolved experiments, and conversely, these experiments may be used to identify dynamic disorder. Variable-temperature NMR data support the notion of three different motional processes with correlation times ranging from 10² to 10⁶ s^-1 over the temperature range of 248-306 K. When molecular motion and crystallographic symmetry are both accounted for, the J(¹¹B,¹¹B) coupling constants for various [B₂(CN)₆]^2- salts were measured to range from 29.4 to 35.8 Hz, and their electronic origins were determined using natural localized molecular orbital and natural bond orbital analyses. The coupling constants were found to strongly correlate to the hybridization states of the boron orbitals that form the B-B bonds and to the strength of the B-B bonds. This study provides a novel tool to study dynamics in ordered and disordered solids and provides new perspectives on electron-precise dianionic diboranes featuring two-center-two-electron bonds in the context of related compounds featuring multiply and singly bonded boron spin pairs.

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.

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.