NMR crystallography of molecular organics

An NMR crystallographic characterisation of solid (+)-usnic acid

We use a combination of one- and two-dimensional solid-state nuclear magnetic resonance (NMR) spectroscopy and density functional theory (DFT) calculations to obtain a full assignment of the 1H and 13C signals for solid (+)-usnic acid, which contains two molecules in the asymmetric unit. By combining through-space 1H-1H correlation data with computation it is possible to assign signals not just to the same molecules (relative assignment) but to assign the signals to specific crystallographic molecules (absolute assignment). Variable-temperature measurements reveal that there is some variation in many of the 13C chemical shifts with temperature, likely arising from varying populations of different tautomeric forms of the molecule. The NMR spectrum of crystalline (+)-usnic acid is then compared with that of ground Usnea dasopoga lichen (the source material of the usnic acid). The abundance of usnic acid is so great in the lichen that this natural product can be observed directly in the NMR spectrum without further purification. This natural sample of usnic acid appears to have the same crystalline form as that in the pure commercial sample.

1H and 13C chemical shift-structure effects in anhydrous β-caffeine and four caffeine-diacid cocrystals probed by solid-state NMR experiments and DFT calculations

By using density functional theory (DFT) calculations, we refined the H atom positions in the structures of β-caffeine (C), α-oxalic acid (OA; (COOH)2), α-(COOH)2·2H2O, β-malonic acid (MA), β-glutaric acid (GA), and I-maleic acid (ME), along with their corresponding cocrystals of 2 : 1 (2C-OA, 2C-MA) or 1 : 1 (C-GA, C-ME) stoichiometry. The corresponding 13C/1H chemical shifts obtained by gauge including projector augmented wave (GIPAW) calculations agreed overall very well with results from magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy experiments. Chemical-shift/structure trends of the precursors and cocrystals were examined, where good linear correlations resulted for all COO1H sites against the H⋯O and/or H⋯N H-bond distance, whereas a general correlation was neither found for the aliphatic/caffeine-stemming 1H sites nor any 13C chemical shift against either the intermolecular hydrogen- or tetrel-bond distance, except for the 13COOH sites of the 2C-OA, 2C-MA, and C-GA cocrystals, which are involved in a strong COOH⋯N bond with caffeine that is responsible for the main supramolecular stabilization of the cocrystal. We provide the first complete 13C NMR spectral assignment of the structurally disordered anhydrous β-caffeine polymorph. The results are discussed in relation to previous literature on the disordered α-caffeine polymorph and the ordered hydrated counterpart, along with recommendations for NMR experimentation that will secure sufficient 13C signal-resolution for reliable resonance/site assignments.

Understanding Cu(i) local environments in MOFs via63/65Cu NMR spectroscopy

The field of metal-organic frameworks (MOFs) includes a vast number of hybrid organic and inorganic porous materials with wide-ranging applications. In particular, the Cu(i) ion exhibits rich coordination chemistry in MOFs and can exist in two-, three-, and four-coordinate environments, which gives rise to many structural motifs and potential applications. Direct characterization of the structurally and chemically important Cu(i) local environments is essential for understanding the sources of specific MOF properties. For the first time, 63/65Cu solid-state NMR has been used to investigate a variety of Cu(i) sites and local coordination geometries in Cu MOFs. This approach is a sensitive probe of the local Cu environment, particularly when combined with density functional theory calculations. A wide range of structurally-dependent 63/65Cu NMR parameters have been observed, including 65Cu quadrupolar coupling constants ranging from 18.8 to 74.8 MHz. Using the data from this and prior studies, a correlation between Cu quadrupolar coupling constants, Cu coordination number, and local Cu coordination geometry has been established. Links between DFT-calculated and experimental Cu NMR parameters are also presented. Several case studies illustrate the feasibility of 63/65Cu NMR for investigating and resolving inequivalent Cu sites, monitoring MOF phase changes, interrogating the Cu oxidation number, and characterizing the product of a MOF chemical reaction involving Cu(ii) reduction to Cu(i). A convenient avenue to acquire accurate 65Cu NMR spectra and NMR parameters from Cu(i) MOFs at a widely accessible magnetic field of 9.4 T is described, with a demonstrated practical application for tracking Cu(i) coordination evolution during MOF anion exchange. This work showcases the power of 63/65Cu solid-state NMR spectroscopy and DFT calculations for molecular-level characterization of Cu(i) centers in MOFs, along with the potential of this protocol for investigating a wide variety of MOF structural changes and processes important for practical applications. This approach has broad applications for examining Cu(i) centers in other weight-dilute systems.

NMR crystallography of amino acids

The development of NMR crystallography methods requires a reliable database of chemical shifts measured for systems with known crystal structure. We measured and assigned carbon and hydrogen chemical shifts of twenty solid natural amino acids of known polymorphic structure, meticulously determined using powder X-ray diffraction. We then correlated the experimental data with DFT-calculated isotropic shieldings. The small size of the unit cell of most amino acids allowed for advanced computations using various families of DFT functionals, including generalized gradient approximation (GGA), meta-GGA and hybrid DFT functionals. We tested several combinations of functionals for geometry optimizations and NMR calculations. For carbon shieldings, the widely used GGA functional PBE performed very well, although an improvement could be achieved by adding shielding corrections calculated for isolated molecules using a hybrid functional. For hydrogen nuclei, we observed the best performance for NMR calculations carried out with structures optimized at the hybrid DFT level. The high fidelity of the calculations made it possible to assign additional signals that could not be assigned based on experiments alone, for example signals of two non-equivalent molecules in the unit cell of some of the amino acids.

Systematic QM/MM Study for Predicting 31P NMR Chemical Shifts of Adenosine Nucleotides in Solution and Stages of ATP Hydrolysis in a Protein Environment

NMR (nuclear magnetic resonance) spectroscopy allows for important atomistic insights into the structure and dynamics of biological macromolecules; however, reliable assignments of experimental spectra are often difficult. Herein, quantum mechanical/molecular mechanical (QM/MM) calculations can provide crucial support. A major problem for the simulations is that experimental NMR signals are time-averaged over much longer time scales, and since computed chemical shifts are highly sensitive to local changes in the electronic and structural environment, sufficiently large averages over representative structural ensembles are essential. This entails high computational demands for reliable simulations. For NMR measurements in biological systems, a nucleus of major interest is 31P since it is both highly present (e.g., in nucleic acids) and easily observable. The focus of our present study is to develop a robust and computationally cost-efficient framework for simulating 31P NMR chemical shifts of nucleotides. We apply this scheme to study the different stages of the ATP hydrolysis reaction catalyzed by p97. Our methodology is based on MM molecular dynamics (MM-MD) sampling, followed by QM/MM structure optimizations and NMR calculations. Overall, our study is one of the most comprehensive QM-based 31P studies in a protein environment and the first to provide computed NMR chemical shifts for multiple nucleotide states in a protein environment. This study sheds light on a process that is challenging to probe experimentally and aims to bridge the gap between measured and calculated NMR spectroscopic properties.

Axial-Equatorial Halide Ordering in Layered Hybrid Perovskites from Isotropic-Anisotropic 207 Pb NMR

Bandgap-tuneable mixed-halide 3D perovskites are of interest for multi-junction solar cells, but suffer from photoinduced spatial halide segregation. Mixed-halide 2D perovskites are more resistant to halide segregation and are promising coatings for 3D perovskite solar cells. The properties of mixed-halide compositions depend on the local halide distribution, which is challenging to study at the level of single octahedra. In particular, it has been suggested that there is a preference for occupation of the distinct axial and equatorial halide sites in mixed-halide 2D perovskites. 207 Pb NMR can be used to probe the atomic-scale structure of lead-halide materials, but although the isotropic 207 Pb shift is sensitive to halide stoichiometry, it cannot distinguish configurational isomers. Here, we use 2D isotropic-anisotropic correlation 207 Pb NMR and relativistic DFT calculations to distinguish the [PbX6 ] configurations in mixed iodide-bromide 3D FAPb(Br1-x Ix )3 perovskites and 2D BA2 Pb(Br1-x Ix )4 perovskites based on formamidinium (FA+ ) and butylammonium (BA+ ), respectively. We find that iodide preferentially occupies the axial site in BA-based 2D perovskites, which may explain the suppressed halide mobility.

Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data

The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.

Carbon-13 chemical shift tensor measurements for nitrogen-dense compounds

This paper reports the principal values of the 13 C chemical shift tensors for five nitrogen-dense compounds (i.e., cytosine, uracil, imidazole, guanidine hydrochloride, and aminoguanidine hydrochloride). Although these are all fundamentally important compounds, the majority do not have 13 C chemical shift tensors reported in the literature. The chemical shift tensors are obtained from 1 H→13 C cross-polarization magic-angle spinning (CP/MAS) experiments that were conducted at a high field of 18.8 T to suppress the effects of 14 N-13 C residual dipolar coupling. Quantum chemical calculations using density functional theory are used to obtain the 13 C magnetic shielding tensors for these compounds. The best agreement with experiment arises from calculations using the hybrid functional PBE0 or the double-hybrid functional PBE0-DH, along with the triple-zeta basis sets TZ2P or pc-3, respectively, and intermolecular effects modeled using large clusters of molecules with electrostatic embedding through the COSMO approach. These measurements are part of an ongoing effort to expand the catalog of accurate 13 C chemical shift tensor measurements, with the aim of creating a database that may be useful for benchmarking the accuracy of quantum chemical calculations, developing nuclear magnetic resonance (NMR) crystallography protocols, or aiding in applications involving machine learning or data mining. This work was conducted at the National High Magnetic Field Laboratory as part of a 2-week school for introducing undergraduate students to practical laboratory experience that will prepare them for scientific careers or postgraduate studies.

Can simple 'molecular' corrections outperform projector augmented-wave density functional theory in the prediction of 35 Cl electric field gradient tensor parameters for chlorine-containing crystalline systems?

Many-body expansion (MBE) fragment approaches have been applied to accurately compute nuclear magnetic resonance (NMR) parameters in crystalline systems. Recent examples demonstrate that electric field gradient (EFG) tensor parameters can be accurately calculated for 14 N and 17 O. A key additional development is the simple molecular correction (SMC) approach, which uses two one-body fragment (i.e., isolated molecule) calculations to adjust NMR parameter values established using 'benchmark' projector augmented-wave (PAW) density functional theory (DFT) values. Here, we apply a SMC using the hybrid PBE0 exchange-correlation (XC) functional to see if this can improve the accuracy of calculated 35 Cl EFG tensor parameters. We selected eight organic and two inorganic crystal structures and considered 15 chlorine sites. We find that this SMC improves the accuracy of computed values for both the 35 Cl quadrupolar coupling constant (CQ ) and the asymmetry parameter ( η Q ) by approximately 30% compared with benchmark PAW DFT values. We also assessed a SMC that offers local improvements not only in terms of the quality of the XC functional but simultaneously in the quality of the description of relativistic effects via the inclusion of spin-orbit effects. As the inorganic systems considered contain heavy atoms bonded to the chlorine atoms, we find further improvements in the accuracy of calculated 35 Cl EFG tensor parameters when both a hybrid functional and spin-orbit effects are included in the SMC. On the contrary, for chlorine-containing organics, the inclusion of spin-orbit relativistic effects using a SMC does not improve the accuracy of computed 35 Cl EFG tensor parameters.

Opportunities and Challenges in Applying Solid-State NMR Spectroscopy in Organic Mechanochemistry

In recent years it is shown that mechanochemical strategies can be beneficial in directed conversions of organic compounds. Finding new reactions proved difficult, and due to the lack of mechanistic understanding of mechanochemical reaction events, respective efforts have mostly remained empirical. Spectroscopic techniques are crucial in shedding light on these questions. In this overview, the opportunities and challenges of solid-state nuclear magnetic resonance (NMR) spectroscopy in the field of organic mechanochemistry are discussed. After a brief discussion of the basics of high-resolution solid-state NMR under magic-angle spinning (MAS) conditions, seven opportunities for solid-state NMR in the field of organic mechanochemistry are presented, ranging from ex situ approaches to structurally elucidated reaction products obtained by milling to the potential and limitations of in situ solid-state NMR approaches. Particular strengths of solid-state NMR, for instance in differentiating polymorphs, in NMR-crystallographic structure-determination protocols, or in detecting weak noncovalent interactions in molecular-recognition events employing proton-detected solid-state NMR experiments at fast MAS frequencies, are discussed.

Synthesis of 3,3'-dihydroxy-2,2'-diindan-1,1'-dione derivatives for tautomeric organic semiconductors exhibiting intramolecular double proton transfer

To investigate potential applications of the 3,3'-dihydroxy-2,2'-biindan-1,1'-dione (BIT) structure as an organic semiconductor with intramolecular hydrogen bonds, a new synthetic route under mild conditions is developed based on the addition reaction of 1,3-dione to ninhydrin and the subsequent hydrogenation of the hydroxyl group. This route affords several new BIT derivatives, including asymmetrically substituted structures that are difficult to access by conventional high-temperature synthesis. The BIT derivatives exhibit rapid tautomerization by intramolecular double proton transfer in solution. The tautomerizations are also observed in the solid state by variable temperature measurements of X-ray diffractometry and magic angle spinning 13C solid-state NMR. Possible interplay between the double proton transfer and the charge transport is suggested by quantum chemical calculations. The monoalkylated BIT derivative with a lamellar packing structure suitable for lateral charge transport in films shows a hole mobility of up to 0.012 cm2 V-1 s-1 with a weak temperature dependence in an organic field effect transistor.

A focus on detection of polymorphs by dynamic nuclear polarization solid-state nuclear magnetic resonance spectroscopy

Solid-state nuclear magnetic resonance (ssNMR) spectroscopy has found increasing application as a method for quantification and structure determination of solid forms (polymorphs) of organic solids and active pharmaceutical ingredients (APIs). However, ssNMR spectroscopy suffers from low sensitivity and resolution, making it challenging to detect dilute solid forms that may be present after recrystallization or reaction with co-formers. Cousin et al. (S. F. Cousin et al., Chem. Sci., 2023, https://doi.org/10.1039/D3SC02063K) have demonstrated that dynamic nuclear polarization (DNP) enhanced 13C cross-polarization (CP) saturation recovery experiments can be used to detect dilute polymorphic forms that are present within a mixture of solid forms. Enhancement of the NMR signal by DNP and differences in signal build-up rates for different polymorphs provide the sensitivity and contrast needed to resolve NMR signals from minor polymorphic forms. This method demonstrated by Cousin et al. should aid the discovery of solid drug forms.

Exploiting solid-state dynamic nuclear polarization NMR spectroscopy to establish the spatial distribution of polymorphic phases in a solid material

Solid-state DNP NMR can enhance the ability to detect minor amounts of solid phases within heterogenous materials. Here we demonstrate that NMR contrast based on the transport of DNP-enhanced polarization can be exploited in the challenging case of early detection of a small amount of a minor polymorphic phase within a major polymorph, and we show that this approach can yield quantitative information on the spatial distribution of the two polymorphs. We focus on the detection of a minor amount (<4%) of polymorph III of m-aminobenzoic acid within a powder sample of polymorph I at natural isotopic abundance. Based on proposed models of the spatial distribution of the two polymorphs, simulations of 1H spin diffusion allow NMR data to be calculated for each model as a function of particle size and the relative amounts of the polymorphs. A comparison between simulated and experimental NMR data allows the model(s) best representing the spatial distribution of the polymorphs in the system to be established.

Review of Applications of Density Functional Theory (DFT) Quantum Mechanical Calculations to Study the High-Pressure Polymorphs of Organic Crystalline Materials

Since its inception, chemistry has been predominated by the use of temperature to generate or change materials, but applications of pressure of more than a few tens of atmospheres for such purposes have been rarely observed. However, pressure is a very effective thermodynamic variable that is increasingly used to generate new materials or alter the properties of existing ones. As computational approaches designed to simulate the solid state are normally tuned using structural data at ambient pressure, applying them to high-pressure issues is a highly challenging test of their validity from a computational standpoint. However, the use of quantum chemical calculations, typically at the level of density functional theory (DFT), has repeatedly been shown to be a great tool that can be used to both predict properties that can be later confirmed by experimenters and to explain, at the molecular level, the observations of high-pressure experiments. This article's main goal is to compile, analyze, and synthesize the findings of works addressing the use of DFT in the context of molecular crystals subjected to high-pressure conditions in order to give a general overview of the possibilities offered by these state-of-the-art calculations.

Atomic-level structure determination of amorphous molecular solids by NMR

Structure determination of amorphous materials remains challenging, owing to the disorder inherent to these materials. Nuclear magnetic resonance (NMR) powder crystallography is a powerful method to determine the structure of molecular solids, but disorder leads to a high degree of overlap between measured signals, and prevents the unambiguous identification of a single modeled periodic structure as representative of the whole material. Here, we determine the atomic-level ensemble structure of the amorphous form of the drug AZD4625 by combining solid-state NMR experiments with molecular dynamics (MD) simulations and machine-learned chemical shifts. By considering the combined shifts of all 1H and 13C atomic sites in the molecule, we determine the structure of the amorphous form by identifying an ensemble of local molecular environments that are in agreement with experiment. We then extract and analyze preferred conformations and intermolecular interactions in the amorphous sample in terms of the stabilization of the amorphous form of the drug.

Discovering the Solid-State Secrets of Lorlatinib by NMR Crystallography: To Hydrogen Bond or not to Hydrogen Bond

Lorlatinib is an active pharmaceutical ingredient (API) used in the treatment of lung cancer. Here, an NMR crystallography analysis is presented whereby the single-crystal X-ray diffraction structure (CSD: 2205098) determination is complemented by multinuclear (1H, 13C, 14/15N, 19F) magic-angle spinning (MAS) solid-state NMR and gauge-including projector augmented wave (GIPAW) calculation of NMR chemical shifts. Lorlatinib crystallises in the P21 space group, with two distinct molecules in the asymmetric unit cell, Z' = 2. Three of the four NH2 hydrogen atoms form intermolecular hydrogen bonds, N30-H…N15 between the two distinct molecules and N30-H…O2 between two equivalent molecules. This is reflected in one of the NH21H chemical shifts being significantly lower, 4.0 ppm compared to 7.0 ppm. Two-dimensional 1H-13C, 14N-1H and 1H (double-quantum, DQ)-1H (single-quantum, SQ) MAS NMR spectra are presented. The 1H resonances are assigned and specific HH proximities corresponding to the observed DQ peaks are identified. The resolution enhancement at a 1H Larmor frequency of 1 GHz as compared to 500 or 600 MHz is demonstrated.

Optically Enhanced Solid-State 1H NMR Spectroscopy

Low sensitivity is the primary limitation to extending nuclear magnetic resonance (NMR) techniques to more advanced chemical and structural studies. Photochemically induced dynamic nuclear polarization (photo-CIDNP) is an NMR hyperpolarization technique where light is used to excite a suitable donor-acceptor system, creating a spin-correlated radical pair whose evolution drives nuclear hyperpolarization. Systems that exhibit photo-CIDNP in solids are not common, and this effect has, up to now, only been observed for ¹³C and ¹⁵N nuclei. However, the low gyromagnetic ratio and natural abundance of these nuclei trap the local hyperpolarization in the vicinity of the chromophore and limit the utility for bulk hyperpolarization. Here, we report the first example of optically enhanced solid-state ¹H NMR spectroscopy in the high-field regime. This is achieved via photo-CIDNP of a donor-chromophore-acceptor molecule in a frozen solution at 0.3 T and 85 K, where spontaneous spin diffusion among the abundant strongly coupled ¹H nuclei relays polarization through the whole sample, yielding a 16-fold bulk ¹H signal enhancement under continuous laser irradiation at 450 nm. These findings enable a new strategy for hyperpolarized NMR beyond the current limits of conventional microwave-driven DNP.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

DELTA50: A Highly Accurate Database of Experimental 1H and 13C NMR Chemical Shifts Applied to DFT Benchmarking

Density functional theory (DFT) benchmark studies of ¹H and ¹³C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured ¹H and ¹³C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for ¹H and ¹³C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for ¹H and 0.5 to 2.9 for ¹³C. Maximum deviations were less than 0.5 and 6.5 ppm for ¹H and ¹³C, respectively.

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

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

Monitoring the influence of additives on the crystallization processes of glycine with dynamic nuclear polarization solid-state NMR

Crystallization is fundamental in many domains, and the investigation of the sequence of solid phases produced as a function of crystallization time is thus key to understand and control crystallization processes. Here, we used a solid-state nuclear magnetic resonance strategy to monitor the crystallization process of glycine, which is a model compound in polymorphism, under the influence of crystallizing additives, such as methanol or sodium chloride. More specifically, our strategy is based on a combination of low-temperatures and dynamic nuclear polarization (DNP) to trap and detect transient crystallizing forms, which may be present only in low quantities. Interestingly, our results show that these additives yield valuable DNP signal enhancements even in the absence of glycerol within the crystallizing solution.

The hydrogen bond continuum in solid isonicotinic acid

The understanding and correct description of intermolecular hydrogen bonds are crucial in the field of multicomponent pharmaceutical solids, such as salts and cocrystals. Solid isonicotinic acid can serve as a suitable model for the development of methods that can accurately characterize these hydrogen bonds. Experimental solid-state NMR has revealed a remarkable temperature dependence and deuterium-isotope-induced changes of the chemical shifts of the atoms involved in the intermolecular hydrogen bond; these NMR data are related to changes of the average position of the hydrogen atom. These changes of NMR parameters were interpreted using periodic DFT path-integral molecular dynamics (PIMD) simulations. The small size of the unit cell of isonicotinic acid allowed for PIMD simulations with the computationally demanding hybrid DFT functional. Calculations of NMR parameters based on the hybrid-functional PIMD simulations are in excellent agreement with experiment. It is thus demonstrated that an accurate characterization of intermolecular hydrogen bonds can be achieved by a combination of NMR experiments and advanced computations.

A Machine Learning Model of Chemical Shifts for Chemically and Structurally Diverse Molecular Solids

Nuclear magnetic resonance (NMR) chemical shifts are a direct probe of local atomic environments and can be used to determine the structure of solid materials. However, the substantial computational cost required to predict accurate chemical shifts is a key bottleneck for NMR crystallography. We recently introduced ShiftML, a machine-learning model of chemical shifts in molecular solids, trained on minimum-energy geometries of materials composed of C, H, N, O, and S that provides rapid chemical shift predictions with density functional theory (DFT) accuracy. Here, we extend the capabilities of ShiftML to predict chemical shifts for both finite temperature structures and more chemically diverse compounds, while retaining the same speed and accuracy. For a benchmark set of 13 molecular solids, we find a root-mean-squared error of 0.47 ppm with respect to experiment for ¹H shift predictions (compared to 0.35 ppm for explicit DFT calculations), while reducing the computational cost by over four orders of magnitude.

Narrowing down the conformational space with solid-state NMR in crystal structure prediction of linezolid cocrystals

Many solids crystallize as microcrystalline powders, thus precluding the application of single crystal X-Ray diffraction in structural elucidation. In such cases, a joint use of high-resolution solid-state NMR and crystal structure prediction (CSP) calculations can be successful. However, for molecules showing significant conformational freedom, the CSP-NMR protocol can meet serious obstacles, including ambiguities in NMR signal assignment and too wide conformational search space to be covered by computational methods in reasonable time. Here, we demonstrate a possible way of avoiding these obstacles and making as much use of the two methods as possible in difficult circumstances. In a simple case, our experiments led to crystal structure elucidation of a cocrystal of linezolid (LIN), a wide-range antibiotic, with 2,3-dihydroxybenzoic acid, while a significantly more challenging case of a cocrystal of LIN with 2,4-dihydroxybenzoic acid led to the identification of the most probable conformations of LIN inside the crystal. Having four rotatable bonds, some of which can assume many discreet values, LIN molecule poses a challenge in establishing its conformation in a solid phase. In our work, a set of 27 conformations were used in CSP calculations to yield model crystal structures to be examined against experimental solid-state NMR data, leading to a reliable identification of the most probable molecular arrangements.

Shared Hydrogen Atom Location and Chemical Composition in Picolinic Acid and Pyridoxal Hydrochloride Derivatives Determined by X-ray Crystallography

Three new single-crystal structures were isolated for picolinic acid (2), the trifluoroacetate salt of picolinic acid (1), and pyridoxal hydrochloride (3). These compounds displayed unconventional crystallographic features that must be considered when structural refinements are carried out. Thus, the generated Fourier differences map obtained with the diffraction data collected at 100 K was crucial to visualize electron densities, which were balanced by either one hydrogen atom or a hydrogen atom with an occupancy factor of 1/2 located between either two carboxylate moieties, two phenolic oxygen atoms, or two pyridinic nitrogen atoms. Moreover, NMR studies were conducted to analyze the bulk chemical composition of single crystals of 2-pyridinecarboxylic acid obtained from the gem-diol/hemiacetal forms and the polymerization products after the treatment of 2-pyridinecarboxaldehyde with TFA:H₂O (1) or a diluted Cu(NO₃)₂ solution (2). The quantitative yield of the pyridoxal hydrochloride crystalline material (3) obtained from a diluted CuCl₂ solution was exhaustively characterized by solid-state NMR methods. These methods allowed the resolution of the signals corresponding to the protons of the hydroxyl moiety of the intramolecular hemiacetal group and the phenolic hydrogen. Theoretical calculations using DFT methods were done to complement the atomic location of the hydrogen atoms obtained from the X-ray analysis.

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

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

19F Solid-state NMR characterization of pharmaceutical solids

Solid-state NMR has been increasingly recognized as a high-resolution and versatile spectroscopic tool to characterize drug substances and products. However, the analysis of pharmaceutical materials is often carried out at natural isotopic abundance and a relatively low drug loading in multi-component systems and therefore suffers from challenges of low sensitivity. The fact that fluorinated therapeutics are well represented in pipeline drugs and commercial products offers an excellent opportunity to utilize fluorine as a molecular probe for pharmaceutical analysis. We aim to review recent advancements of ¹⁹F magic angle spinning NMR methods in modern drug research and development. Applications to polymorph screening at the micromolar level, structural elucidation, and investigation of molecular interactions at the Ångström to submicron resolution in drug delivery, stability, and quality will be discussed.

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 ¹H-X (X = ¹⁴N, ³⁵Cl) interactions. This experiment can be used to acquire 2D ¹H{X}-HMQC filtered ¹H-¹H correlation as well as 2D ¹H-X HMQC spectra. Powder forms of dopamine·HCl and l-histidine·HCl·H₂O 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·H₂O, 2D peaks associated with ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H distances of up to 4.4 and 3.9 Å have been detected. This is further corroborated by the observation of 2D peaks corresponding to ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H 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.

Molecular and Electronic Structure of Isolated Platinum Sites Enabled by the Expedient Measurement of 195Pt Chemical Shift Anisotropy

Techniques that can characterize the molecular structures of dilute surface species are required to facilitate the rational synthesis and improvement of Pt-based heterogeneous catalysts. ¹⁹⁵Pt solid-state NMR spectroscopy could be an ideal tool for this task because ¹⁹⁵Pt isotropic chemical shifts and chemical shift anisotropy (CSA) are highly sensitive probes of the local chemical environment and electronic structure. However, the characterization of Pt surface-sites is complicated by the typical low Pt loadings that are between 0.2 and 5 wt% and broadening of ¹⁹⁵Pt solid-state NMR spectra by CSA. Here, we introduce a set of solid-state NMR methods that exploit fast MAS and indirect detection using a sensitive spy nucleus (¹H or ³¹P) to enable the rapid acquisition of ¹⁹⁵Pt MAS NMR spectra. We demonstrate that high-resolution wideline ¹⁹⁵Pt MAS NMR spectra can be acquired in minutes to a few hours for a series of molecular and single-site Pt species grafted on silica with Pt loading of only 3-5 wt%. Low-power, long-duration, sideband-selective excitation, and saturation pulses are incorporated into t₁-noise eliminated dipolar heteronuclear multiple quantum coherence, perfect echo resonance echo saturation pulse double resonance, or J-resolved pulse sequences. The complete ¹⁹⁵Pt MAS NMR spectrum is then reconstructed by recording a series of 1D NMR spectra where the offset of the ¹⁹⁵Pt pulses is varied in increments of the MAS frequency. Analysis of the ¹⁹⁵Pt MAS NMR spectra yields the ¹⁹⁵Pt chemical shift tensor parameters. Zeroth order approximation density functional theory calculations accurately predict ¹⁹⁵Pt CS tensor parameters. Simple and predictive orbital models relate the CS tensor parameters to the Pt electronic structure and coordination environment. The methodology developed here paves the way for the detailed structural and electronic analysis of dilute platinum surface-sites.

Solid-state multinuclear magnetic resonance and X-ray crystallographic investigation of the phosphorus...iodine halogen bond in a bis(dicyclohexylphenylphosphine)(1,6-diiodoperfluorohexane) cocrystal

Halogen bonding to phosphorus atoms remains uncommon, with relatively few examples reported in the literature. Here, the preparation and investigation of the cocrystal bis(dicyclohexylphenylphosphine)(1,6-diiodoperfluorohexane) by X-ray crystallography and solid-state multinuclear magnetic resonance spectroscopy is described. The crystal structure features two crystallographically unique C-I...P halogen bonds [d_I...P = 3.090 (5) Å, 3.264 (5) Å] and crystallographic disorder of one of the 1,6-diiodoperfluorohexane molecules. The first of these is the shortest and most linear I...P halogen bond reported to date. ¹³C, ¹⁹F, and ³¹P magic angle spinning solid-state NMR spectra are reported. A ³¹P chemical shift change of -7.0 p.p.m. in the cocrystal relative to pure dicyclohexylphenylphosphine, consistent with halogen bond formation, is noted. This work establishes iodoperfluoroalkanes as viable halogen bond donors when paired with phosphorus acceptors, and also shows that dicyclohexylphenylphosphine can act as a practical halogen bond acceptor.

A structural investigation of organic battery anode materials by NMR crystallography

Conjugated alkali metal dicarboxylates have recently received attention for applications as organic anode materials in lithium- and sodium-ion batteries. In order to understand and optimise these materials, it is important to be able to characterise both the long-range and local aspects of the crystal structure, which may change during battery cycling. Furthermore, some materials can display polymorphism or hydration behaviour. NMR crystallography, which combines long-range crystallographic information from diffraction with local information from solid-state NMR via interpretation aided by DFT calculations, is one such approach, but this has not yet been widely applied to conjugated dicarboxylates. In this work, we evaluate the application of NMR crystallography for a set of model lithium and sodium dicarboxylate salts. We investigate the effect of different DFT geometry optimisation strategies and find that the calculated NMR parameters are not systematically affected by the choice of optimisation method, although the inclusion of dispersion correction schemes is important to accurately reproduce the experimental unit cell parameters. We also observe hydration behaviour for two of the sodium salts and provide insight into the structure of an as-yet uncharacterised structure of sodium naphthalenedicarboxylate. This highlights the importance of sample preparation and characterisation for organic sodium-ion battery anode materials in particular.

Characterization of crystalline and amorphous forms of irbesartan by multi-nuclear solid-state NMR

Irbesartan (IRB) is an antihypertensive drug which exhibits the rare phenomenon of desmotropy; its 1H- and 2H- tetrazole tautomers can be isolated as distinct crystalline forms. The crystalline forms of IRB are poorly soluble, hence the amorphous form is potentially of interest for its faster dissolution rate. The tautomeric form and the nature of hydrogen bonding in amorphous IRB are unknown. In this study, crystalline form A and amorphous form of irbesartan were studied using ¹³C, ¹⁵N and ¹H solid-state NMR. Variable-temperature ¹³C SSMNR studies showed alkyl chain disorder in the crystalline form of IRB, which may explain the conflicting literature crystal structures of form A (the marketed form). ¹⁵N NMR indicates that the amorphous material contains an approximately 2:1 ratio of 1H- and 2H-tetrazole tautomers. Static ¹H SSNMR and relaxation time measurements confirmed different molecular mobilities of the samples and provided molecular-level insight into the nature of the glass transition. SSNMR is shown to be a powerful technique to investigate the solid state of disordered active pharmaceutical ingredients.

Resolving alternative structure determinations of indapamide using 13C solid-state NMR

The conflict between alternative crystal structures in the Cambridge Structural Database for the diuretic drug indapamide hemihydrate (IND) has been resolved with the aid of ¹³C solid-state NMR. IND is seen to contain multiple distinct molecules in the asymmetric unit (Z' = 4) rather than exhibiting disorder in the orientation of sulfonamide groups. The NMR crystallographic approach is a more effective tool for distinguishing between alternative structures than naïve judgements of quality based on crystallographic refinement agreement factors.

Resolving Atomic-Scale Interactions in Nonfullerene Acceptor Organic Solar Cells with Solid-State NMR Spectroscopy, Crystallographic Modelling, and Molecular Dynamics Simulations

Fused-ring core nonfullerene acceptors (NFAs), designated "Y-series," have enabled high-performance organic solar cells (OSCs) achieving over 18% power conversion efficiency (PCE). Since the introduction of these NFAs, much effort has been expended to understand the reasons for their exceptional performance. While several studies have identified key optoelectronic properties that govern high PCEs, little is known about the molecular level origins of large variations in performance, spanning from 5% to 18% PCE, for example, in the case of PM6:Y6 OSCs. Here, a combined solid-state NMR, crystallography, and molecular modeling approach to elucidate the atomic-scale interactions in Y6 crystals, thin films, and PM6:Y6 bulk heterojunction (BHJ) blends is introduced. It is shown that the Y6 morphologies in BHJ blends are not governed by the morphology in neat films or single crystals. Notably, PM6:Y6 blends processed from different solvents self-assemble into different structures and morphologies, whereby the relative orientations of the sidechains and end groups of the Y6 molecules to their fused-ring cores play a crucial role in determining the resulting morphology and overall performance of the solar cells. The molecular-level understanding of BHJs enabled by this approach will guide the engineering of next-generation NFAs for stable and efficient OSCs.

Atomic-resolution chemical characterization of (2x)72-kDa tryptophan synthase via four- and five-dimensional 1H-detected solid-state NMR

NMR chemical shifts provide detailed information on the chemical properties of molecules, thereby complementing structural data from techniques like X-ray crystallography and electron microscopy. Detailed analysis of protein NMR data, however, often hinges on comprehensive, site-specific assignment of backbone resonances, which becomes a bottleneck for molecular weights beyond 40 to 45 kDa. Here, we show that assignments for the (2x)72-kDa protein tryptophan synthase (665 amino acids per asymmetric unit) can be achieved via higher-dimensional, proton-detected, solid-state NMR using a single, 1-mg, uniformly labeled, microcrystalline sample. This framework grants access to atom-specific characterization of chemical properties and relaxation for the backbone and side chains, including those residues important for the catalytic turnover. Combined with first-principles calculations, the chemical shifts in the β-subunit active site suggest a connection between active-site chemistry, the electrostatic environment, and catalytically important dynamics of the portal to the β-subunit from solution.

Enantiotropy of Simvastatin as a Result of Weakened Interactions in the Crystal Lattice: Entropy-Driven Double Transitions and the Transient Modulated Phase as Seen by Solid-State NMR Spectroscopy

In crystalline molecular solids, in the absence of strong intermolecular interactions, entropy-driven processes play a key role in the formation of dynamically modulated transient phases. Specifically, in crystalline simvastatin, the observed fully reversible enantiotropic behavior is associated with multiple order–disorder transitions: upon cooling, the dynamically disordered high-temperature polymorphic Form I is transformed to the completely ordered low-temperature polymorphic Form III via the intermediate (transient) modulated phase II. This behavior is associated with a significant reduction in the kinetic energy of the rotating and flipping ester substituents, as well as a decrease in structural ordering into two distinct positions. In transient phase II, the conventional three-dimensional structure is modulated by periodic distortions caused by cooperative conformation exchange of the ester substituent between the two states, which is enabled by weakened hydrogen bonding. Based on solid-state NMR data analysis, the mechanism of the enantiotropic phase transition and the presence of the transient modulated phase are documented.

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate Cβ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

Multi-Component Crystals Containing Urea: Mechanochemical Synthesis and Characterization by 35Cl Solid-State NMR Spectroscopy and DFT Calculations

Mechanochemical synthesis provides new pathways for the rational design of multi-component crystals (MCCs) involving anionic or cationic components, which offer molecular-level architectures unavailable to MCCs comprised of strictly neutral components....

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

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

Submerged Eutectic-Assisted, Solvent-Free Mechanochemical Formation of a Propranolol Salt and Its Other Multicomponent Solids

Salt preparation via a solid-state reaction offers a solution to challenges posed by current pharmaceutical research, which include combining development of novel forms of active pharmaceutical ingredients with greener, sustainable synthesis. This work investigated in detail the mechanism of salt formation between propranolol (PRO) and capric acid (CAP) and explored the solid eutectic phases comprising this salt, propranolol caprate (PRC). The salt structure was solved by X-ray diffraction, and the properties in the crystalline and supercooled states were fully characterised using thermal analysis, nuclear magnetic resonance, Fourier-transform infrared spectroscopy and broadband dielectric spectroscopy (BDS). PRC forms via a submerged eutectic phase composed of PRO and CAP, below room temperature, by mechanochemistry without an extra input of energy. Two other solid eutectic phases are composed of PRC and either CAP or PRO, at 0.28 and 0.82 mol fraction of PRO, respectively. BDS indicated that the supercooled PRC has ionic character, whereas the supercooled PRC-PRO eutectic had predominantly non-ionic properties despite comprising the salt. In conclusion, knowledge of the mechanism of formation of multicomponent systems can help in designing more sustainable pharmaceutical processes.

A toolbox for improving the workflow of NMR crystallography

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

Modeling the Structure of Crystalline Alamethicin and Its NMR Chemical Shift Tensors

Alamethicin (ALM) is an antimicrobial peptide that is frequently employed in studies of the mechanism of action of pore-forming molecules. Advanced techniques of solid-state NMR spectroscopy (SSNMR) are important in these studies, as they are capable of describing the alignment of helical peptides, such as ALM, in lipid bilayers. Here, it is demonstrated how an analysis of the SSNMR measurements can benefit from fully periodic calculations, which employ the plane-wave density-functional theory (PW DFT) of the solid-phase geometry and related spectral parameters of ALM. The PW DFT calculations are used to obtain the structure of desolvated crystalline ALM and predict the NMR chemical shift tensors (CSTs) of its nuclei. A variation in the CSTs of the amidic nitrogens and carbonyl carbons along the ALM backbone is evaluated and included in simulations of the orientation-dependent anisotropic 15N and 13C chemical shift components. In this way, the influence of the site-specific structural effects on the experimentally determined orientation of ALM is shown in models of cell membranes.

Crystal Engineering of Supramolecular 1,4-Benzene Bisamides by Side-Chain Modification - Towards Tuneable Anisotropic Morphologies and Surfaces

Benzene bisamides are promising building blocks for supramolecular nano‐objects. Their functionality depends on morphology and surface properties. However, a direct link between surface properties and molecular structure itself is missing for this material class. Here, we investigate this interplay for two series of 1,4‐benzene bisamides with symmetric and asymmetric peripheral substitution. We elucidated the crystal structures, determined the nano‐object morphologies and derived the wetting behaviour of the preferentially exposed surfaces. The crystal structures were solved by combining single‐crystal and powder X‐ray diffraction, solid‐state NMR spectroscopy and computational modelling. Bulky side groups, here t‐butyl groups, serve as a structure‐directing motif into a packing pattern, which favours the formation of thin platelets. The use of slim peripheral groups on both sides, in our case linear perfluorinated, alkyl chains, self‐assemble the benzene bisamides into a second packing pattern which leads to ribbon‐like nano‐objects. For both packing types, the preferentially exposed surfaces consist of the ends of the peripheral groups. Asymmetric substitution with bulky and slim groups leads to an ordered alternating arrangement of the groups exposed to the surface. This allows the hydrophobicity of the surfaces to be gradually altered. We thus identified two leitmotifs for molecular packings of benzene bisamides providing the missing link between the molecular structure, the anisotropic morphologies and adjustable surface properties of the supramolecular nano‐objects.

Refined Structures of O-Phospho-l-serine and Its Calcium Salt by New Multinuclear Solid-State NMR Crystallography Methods

O-phospho-l-serine (Pser) and its Ca salt, Ca[O-phospho-l-serine]·H2O (CaPser), play important roles for bone mineralization and were recently also proposed to account for the markedly improved bone-adhesive properties of Pser-doped calcium phosphate-based cements for biomedical implants. However, the hitherto few proposed structural models of Pser and CaPser were obtained by X-ray diffraction, thereby leaving the proton positions poorly defined. Herein, we refine the Pser and CaPser structures by density functional theory (DFT) calculations and contrast them with direct interatomic-distance constraints from two-dimensional (2D) nuclear magnetic resonance (NMR) correlation experimentation at fast magic-angle spinning (MAS), encompassing double-quantum-single-quantum (2Q-1Q) 1H NMR along with heteronuclear 13C{1H} and 31P{1H} correlation NMR experiments. The Pser and CaPser structures before and after refinements by DFT were validated against sets of NMR-derived effective 1H-1H, 1H-31P, and 1H-13C distances, which confirmed the improved accuracy of the refined structures. Each distance set was derived from one sole 2D NMR experiment applied to a powder without isotopic enrichment. The distances were extracted without invoking numerical spin-dynamics simulations or approximate phenomenological models. We highlight the advantages and limitations of the new distance-extraction procedure. Isotropic 1H, 13C, and 31P chemical shifts obtained by DFT calculations using the gauge including projector augmented wave (GIPAW) method agreed very well with the experimental results. We discuss the isotropic and anisotropic 13C and 31P chemical-shift parameters in relation to the previous literature, where most data on CaPser are reported herein for the first time.

Drug-Polymer Interactions in Acetaminophen/Hydroxypropylmethylcellulose Acetyl Succinate Amorphous Solid Dispersions Revealed by Multidimensional Multinuclear Solid-State NMR Spectroscopy

The bioavailability of insoluble crystalline active pharmaceutical ingredients (APIs) can be enhanced by formulation as amorphous solid dispersions (ASDs). One of the key factors of ASD stabilization is the formation of drug-polymer interactions at the molecular level. Here, we used a range of multidimensional and multinuclear nuclear magnetic resonance (NMR) experiments to identify these interactions in amorphous acetaminophen (paracetamol)/hydroxypropylmethylcellulose acetyl succinate (HPMC-AS) ASDs at various drug loadings. At low drug loading (<20 wt %), we showed that 1H-13C through-space heteronuclear correlation experiments identify proximity between aromatic protons in acetaminophen with cellulose backbone protons in HPMC-AS. We also show that 14N-1H heteronuclear multiple quantum coherence (HMQC) experiments are a powerful approach in probing spatial interactions in amorphous materials and establish the presence of hydrogen bonds (H-bond) between the amide nitrogen of acetaminophen with the cellulose ring methyl protons in these ASDs. In contrast, at higher drug loading (40 wt %), no acetaminophen/HPMC-AS spatial proximity was identified and domains of recrystallization of amorphous acetaminophen into its crystalline form I, the most thermodynamically stable polymorph, and form II are identified. These results provide atomic scale understanding of the interactions in the acetaminophen/HPMC-AS ASD occurring via H-bond interactions.

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 ¹H-¹⁴N phase-modulated rotational-echo saturation-pulse double-resonance (PM-S-RESPDOR) and ¹H-¹H selective recoupling of proton (SERP) experiments. The former examined the dephasing effect of a specific ¹H resonance under multiple ¹H-¹⁴N dipolar couplings. The latter examined the selective polarization transfer between a ¹H-¹H 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 ¹H and ¹⁴N, 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.

Structural Refinement of Carbimazole by NMR Crystallography

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

Analyzing Discrepancies in Chemical-Shift Predictions of Solid Pyridinium Fumarates

Highly accurate chemical-shift predictions in molecular solids are behind the success and rapid development of NMR crystallography. However, unusually large errors of predicted hydrogen and carbon chemical shifts are sometimes reported. An understanding of these deviations is crucial for the reliability of NMR crystallography. Here, recently reported large deviations of predicted hydrogen and carbon chemical shifts of a series of solid pyridinium fumarates are thoroughly analyzed. The influence of the geometry optimization protocol and of the computational level of NMR calculations on the accuracy of predicted chemical shifts is investigated. Periodic calculations with GGA, meta-GGA and hybrid functionals are employed. Furthermore, molecular corrections at the coupled-cluster singles-and-doubles (CCSD) level are calculated. The effect of nuclear delocalization on the structure and NMR shielding is also investigated. The geometry optimization with a computationally demanding hybrid functional leads to a substantial improvement in proton chemical-shift predictions.

Paramagnetic solid-state NMR assignment and novel chemical conversion of the aldehyde group to dihydrogen ortho ester and hemiacetal moieties in copper(ii)- and cobalt(ii)-pyridinecarboxaldehyde complexes

The complex chemical functionalization of aldehyde moieties in Cu(ii)- and Co(ii)-pyridinecarboxaldehyde complexes was studied. X-ray studies demonstrated that the aldehyde group (RCHO) of the four pyridine molecules is converted to dihydrogen ortho ester (RC(OCH₃)(OH)₂) and hemiacetal (RCH(OH)(OCH₃)) moieties in both 4-pyridinecarboxaldehyde copper and cobalt complexes. In contrast, the aldehyde group is retained when the 3-pyridinecarboxaldehyde ligand is complexed with cobalt. In the different copper complexes, similar paramagnetic ¹H resonance lines were obtained in the solid state; however, the connectivity with the carbon structure and the ¹H vicinities were done with 2D ¹H–¹³C HETCOR, ¹H–¹H SQ/DQ and proton spin diffusion (PSD) experiments. The strong paramagnetic effect exerted by the cobalt center prevented the observation of ¹³C NMR signals and chemical information could only be obtained from X-ray experiments. 2D PSD experiments in the solid state were useful for the proton assignments in both Cu(ii) complexes. The combination of X-ray crystallography experiments with DFT calculations together with the experimental results obtained from EPR and solid-state NMR allowed the assignment of NMR signals in pyridinecarboxaldehyde ligands coordinated with copper ions. In cases where the crystallographic information was not available, as in the case of the 3-pyridinecarboxaldehyde Cu(ii) complex, the combination of these techniques allowed not only the assignment of NMR signals but also the study of the functionalization of the substituent group.

The complex chemical functionalization of the aldehyde group was elucidated in copper and cobalt complexes for 4- and 3-pyridinecarboxaldehyde ligands.