Atomic-level structure determination of amorphous molecular solids by NMR

Spiers Memorial Lecture: NMR crystallography

Chemical function is directly related to the spatial arrangement of atoms. Consequently, the determination of atomic-level three-dimensional structures has transformed molecular and materials science over the past 60 years. In this context, solid-state NMR has emerged to become the method of choice for atomic-level characterization of complex materials in powder form. In the following we present an overview of current methods for chemical shift driven NMR crystallography, illustrated with applications to complex materials.

An Updated Structure of Oxybutynin Hydrochloride

Oxybutynin (Ditropan), a widely distributed muscarinic antagonist for treating the overactive bladder, has been awaiting a definitive crystal structure for ≈50 years due to the sample and technique limitations. Past reports used powder X-ray diffraction (PXRD) to shed light on the possible packing of the molecule however their model showed some inconsistencies when compared with the 2D chemical structure. These are largely attributed to X-ray-induced photoreduction. Here microcrystal electron diffraction (MicroED) is used to successfully unveil the experimental 3D structure of oxybutynin hydrochloride showing marked improvement over the reported PXRD structure. Using the improved model, molecular docking is applied to investigate the binding mechanism between M3 muscarinic receptor (M3R) and (R)-oxybutynin, revealing essential contacts/residues and conformational changes within the protein pocket. A possible universal conformation is proposed for M3R antagonists, which is valuable for future drug development and optimization. This study underscores the immense potential of MicroED as a complementary technique for elucidating unknown pharmaceutical structures, as well as for protein-drug interactions.

The essential synergy of MD simulation and NMR in understanding amorphous drug forms

Molecular dynamics (MD) simulations and chemical shifts from machine learning are used to predict 15N, 13C and 1H chemical shifts for the amorphous form of the drug irbesartan. The local environments are observed to be highly dynamic well below the glass transition, and averaging over the dynamics is essential to understanding the observed NMR shifts. Predicted linewidths are about 2 ppm narrower than observed experimentally, which is hypothesised to largely result from susceptibility effects. Previously observed differences in the 13C shifts associated with the two tetrazole tautomers can be rationalised in terms of differing conformational dynamics associated with the presence of an intramolecular interaction in one tautomer. 1H shifts associated with hydrogen bonding can also be rationalised in terms of differing average frequencies of transient hydrogen bonding interactions.

Crystal structure validation of verinurad via proton-detected ultra-fast MAS NMR and machine learning

The recent development of ultra-fast magic-angle spinning (MAS) (>100 kHz) provides new opportunities for structural characterization in solids. Here, we use NMR crystallography to validate the structure of verinurad, a microcrystalline active pharmaceutical ingredient. To do this, we take advantage of 1H resolution improvement at ultra-fast MAS and use solely 1H-detected experiments and machine learning methods to assign all the experimental proton and carbon chemical shifts. This framework provides a new tool for elucidating chemical information from crystalline samples with limited sample volume and yields remarkably faster acquisition times compared to 13C-detected experiments, without the need to employ dynamic nuclear polarization.

Atomic-level structure of the amorphous drug atuliflapon via NMR crystallography

We determine the complete atomic-level structure of the amorphous form of the drug atuliflapon, a 5-lipooxygenase activating protein (FLAP) inhibitor, via chemical-shift-driven NMR crystallography. The ensemble of preferred structures allows us to identify a number of specific conformations and interactions that stabilize the amorphous structure. These include preferred hydrogen-bonding motifs with water and with other drug molecules, as well as conformations of the cyclohexane and pyrazole rings that stabilize structure by indirectly allowing for optimization of hydrogen bonding.

The interplay of density functional selection and crystal structure for accurate NMR chemical shift predictions

Ab initio chemical shift prediction plays a central role in nuclear magnetic resonance (NMR) crystallography, and the accuracy with which chemical shifts can be predicted relative to experiment impacts the confidence with which structures can be assigned. For organic crystals, periodic density functional theory calculations with the gauge-including projector augmented wave (GIPAW) approximation and the PBE functional are widely used at present. Many previous studies have examined how using more advanced density functionals can increase the accuracy of predicted chemical shifts relative to experiment, but nearly all of those studies employed crystal structures that were optimized with generalized-gradient approximation (GGA) functionals. Here, we investigate how the accuracy of the predicted chemical shifts in organic crystals is affected by replacing GGA-level PBE-D3(BJ) crystal geometries with more accurate hybrid functional PBE0-D3(BJ) ones. Based on benchmark data sets containing 132 13C and 35 15N chemical shifts, plus case studies on testosterone, acetaminophen, and phenobarbital, we find that switching from GGA-level geometries and chemical shifts to hybrid-functional ones reduces 13C and 15N chemical shift errors by ∼40-60% versus experiment. However, most of the improvement stems from the use of the hybrid functional for the chemical shift calculations, rather than from the refined geometries. In addition, even with the improved geometries, we find that double-hybrid functionals still do not systematically increase chemical shift agreement with experiment beyond what hybrid functionals provide. In the end, these results suggest that the combination of GGA-level crystal structures and hybrid-functional chemical shifts represents a particularly cost-effective combination for NMR crystallography in organic systems.

An Updated Structure of Oxybutynin Hydrochloride

Oxybutynin (Ditropan), a widely distributed muscarinic antagonist for treating the overactive bladder, has been awaiting a definitive crystal structure for nearly 50 years due to the sample and technique limitations. Past reports used powder X-ray diffraction (PCRD) to shed light on the possible packing of the molecule however a 3D structure remained elusive. Here we used Microcrystal Electron Diffraction (MicroED) to successfully unveil the 3D structure of oxybutynin hydrochloride. We identify several inconsistencies between the reported PXRD analyses and the experimental structure. Using the improved model, molecular docking was applied to investigate the binding mechanism between M3 muscarinic receptor (M3R) and (R)-oxybutynin, revealing essential contacts/residues and conformational changes within the protein pocket. A possible universal conformation was proposed for M3R antagonists, which is valuable for future drug development and optimization. This study underscores the immense potential of MicroED as a complementary technique for elucidating the unknown pharmaceutical crystal structures, as well as for the protein-drug interactions.

Salt-induced Fmoc-tripeptide supramolecular hydrogels: a combined experimental and computational study of the self-assembly

Delving into the mechanism behind the molecular interactions at the atomic level of short-sequence peptides plays a key role in the development of nanomaterials with specific structure-property-function relationships from a bottom-up perspective. Due to their poor water solubility, the self-assembly of Fmoc-bearing peptides is usually induced by dissolution in an organic solvent, followed by a dilution step in water, pH changes, and/or a heating-cooling process. Herein, we report a straightforward methodology for the gelation of Fmoc-FFpY (F: phenylalanine; Y: tyrosine; and p: PO42-), a negatively charged tripeptide, in NaCl solution. The electrostatic interactions between Fmoc-FFpY and Na+ ions give rise to different nanofibrillar hydrogels with rheological properties and nanofiber sizes modulated by the NaCl concentration in pure aqueous media. Initiated by the electrostatic interactions between the peptide phosphate groups and the Na+ ions, the peptide self-assembly is stabilized thanks to hydrogen bonds between the peptide backbones and the π-π stacking of aromatic Fmoc and phenyl units. The hydrogels showed self-healing and thermo-responsive properties for potential biomedical applications. Molecular dynamics simulations from systems devoid of prior training not only confirm the aggregation of peptides at a critical salt concentration and the different interactions involved, but also corroborate the secondary structure of the hydrogels at the microsecond timescale. It is worth highlighting the remarkable achievement of reproducing the morphological behavior of the hydrogels using atomistic simulations. To our knowledge, this study is the first to report such a correspondence.

Optimal sensitivity for 1H detected relayed DNP of organic solids at fast MAS

Dynamic nuclear polarization (DNP) combined with high magnetic fields and fast magic angle spinning (MAS) has opened up a new avenue for the application of exceptionally sensitive 1H NMR detection schemes to study protonated solids. Recently, it has been shown that DNP experiments at fast MAS rates lead to slower spin diffusion and hence reduced DNP enhancements for impregnated materials. However, DNP enhancements alone do not determine the overall sensitivity of a NMR experiment. Here we measure the overall sensitivity of one-dimensional 1H detected relayed DNP experiments as a function of the MAS rate in the 20-60 kHz regime using 0.7 mm diameter rotors at 21.2 T. Although faster MAS rates are detrimental for the DNP enhancement on the target material, due to slower spin diffusion, we find that with increasing spinning rates the gain in sensitivity due to 1H line-narrowing and the folding-in of sideband intensity compensates a large part of the loss of overall hyperpolarization. We find that sensitivity depends on the atomic site in the molecule, and is maximised at between 40 and 50 kHz MAS for the sample of L-histidine.HCl·H2O studied here. There is a 10-20 % difference in sensitivity between the optimum MAS rate and the fastest rate currently accessible (60 kHz).

Rational Design of Dinitroxide Polarizing Agents for Dynamic Nuclear Polarization to Enhance Overall NMR Sensitivity

We evaluate the overall sensitivity gains provided by a series of eighteen nitroxide biradicals for dynamic nuclear polarization (DNP) solid-state NMR at 9.4 T and 100 K, including eight new biradicals. We find that in the best performing group the factors contributing to the overall sensitivity gains, namely the DNP enhancement, the build-up time, and the contribution factor, often compete with each other leading to very similar overall sensitivity across a range of biradicals. NaphPol and HydroPol are found to provide the best overall sensitivity factors, in organic and aqueous solvents respectively. One of the new biradicals, AMUPolCbm, provides high sensitivity for all three solvent formulations measured here, and can be considered to be a "universal" polarizing agent.

Unraveling the Structure of Meclizine Dihydrochloride with MicroED

Meclizine (Antivert, Bonine) is a first-generation H1 antihistamine used in the treatment of motion sickness and vertigo. Despite its wide medical use for over 70 years, its crystal structure and the details of protein-drug interactions remained unknown. Single-crystal X-ray diffraction (SC-XRD) is previously unsuccessful for meclizine. Today, microcrystal electron diffraction (MicroED) enables the analysis of nano- or micro-sized crystals that are merely a billionth the size needed for SC-XRD directly from seemingly amorphous powder. In this study, MicroED to determine the 3D crystal structure of meclizine dihydrochloride is used. Two racemic enantiomers (R/S) are found in the unit cell, which is packed as repetitive double layers in the crystal lattice. The packing is made of multiple strong N-H-Cl- hydrogen bonding interactions and weak interactions like C-H-Cl- and pi-stacking. Molecular docking reveals the binding mechanism of meclizine to the histamine H1 receptor. A comparison of the docking complexes between histamine H1 receptor and meclizine or levocetirizine (a second-generation antihistamine) shows the conserved binding sites. This research illustrates the combined use of MicroED and molecular docking in unraveling elusive drug structures and protein-drug interactions for precision drug design and optimization.

Probing Molecular Packing of Amorphous Pharmaceutical Solids Using X-ray Atomic Pair Distribution Function and Solid-State NMR

The structural investigation of amorphous pharmaceuticals is of paramount importance in comprehending their physicochemical stability. However, it has remained a relatively underexplored realm primarily due to the limited availability of high-resolution analytical tools. In this study, we utilized the combined power of X-ray pair distribution functions (PDFs) and solid-state nuclear magnetic resonance (ssNMR) techniques to probe the molecular packing of amorphous posaconazole and its amorphous solid dispersion at the molecular level. Leveraging synchrotron X-ray PDF data and employing the empirical potential structure refinement (EPSR) methodology, we unraveled the existence of a rigid conformation and discerned short-range intermolecular C-F contacts within amorphous posaconazole. Encouragingly, our ssNMR 19F-13C distance measurements offered corroborative evidence supporting these findings. Furthermore, employing principal component analysis on the X-ray PDF and ssNMR data sets enabled us to gain invaluable insights into the chemical nature of the intermolecular interactions governing the drug-polymer interplay. These outcomes not only furnish crucial structural insights facilitating the comprehension of the underlying mechanisms governing the physicochemical stability but also underscore the efficacy of synergistically harnessing X-ray PDF and ssNMR techniques, complemented by robust modeling strategies, to achieve a high-resolution exploration of amorphous structures.