Multinuclear Solid-State Magnetic Resonance as a Sensitive Probe of Structural Changes upon the Occurrence of Halogen Bonding in Co-crystals

Although the understanding of intermolecular interactions, such as hydrogen bonding, is relatively well-developed, many additional weak interactions work both in tandem and competitively to stabilize a given crystal structure. Due to a wide array of potential applications, a substantial effort has been invested in understanding the halogen bond. Here, we explore the utility of multinuclear (13C, 14/15N, 19F, and 127I) solid-state magnetic resonance experiments in characterizing the electronic and structural changes which take place upon the formation of five halogen-bonded co-crystalline product materials. Single-crystal X-ray diffraction (XRD) structures of three novel co-crystals which exhibit a 1:1 stoichiometry between decamethonium diiodide (i.e., [(CH3)3N+(CH 2)10N+(CH3)3][2 I -]) and different para-dihalogen-substituted benzene moieties (i.e., p-C6X2Y4, X=Br, I; Y=H, F) are presented. 13C and 15N NMR experiments carried out on these and related systems validate sample purity, but also serve as indirect probes of the formation of a halogen bond in the co-crystal complexes in the solid state. Long-range changes in the electronic environment, which manifest through changes in the electric field gradient (EFG) tensor, are quantitatively measured using 14N NMR spectroscopy, with a systematic decrease in the 14N quadrupolar coupling constant (CQ) observed upon halogen bond formation. Attempts at 127I solid-state NMR spectroscopy experiments are presented and variable-temperature 19F NMR experiments are used to distinguish between dynamic and static disorder in selected product materials, which could not be conclusively established using solely XRD. Quantum chemical calculations using the gauge-including projector augmented-wave (GIPAW) or relativistic zeroth-order regular approximation (ZORA) density functional theory (DFT) approaches complement the experimental NMR measurements and provide theoretical corroboration for the changes in NMR parameters observed upon the formation of a halogen bond

Elucidating an amorphous form stabilization mechanism of tenapanor hydrochloride: crystal structure analysis using Xray diffraction, NMR crystallography and molecular modelling

By the combined use of powder and single crystal X-ray diffraction, solid-state NMR, and molecular modelling, the crystal structures of two systems containing the unusually large tenapanor drug molecule have been determined: the free form, ANHY and a dihydrochloride salt form, 2HCl. Dynamic nuclear polarization (DNP) assisted solid-state NMR (SSNMR) crystallography investigations were found essential for the final assignment, and were used to validate the crystal structure of ANHY. From the structural informatics analysis of ANHY and 2HCl, conformational ring differences in one part of the molecule were observed which influences the relative orientation of a methyl group on a ring nitrogen and thereby impacts the crystallizability of the dihydrochloride salt. From quantum chemistry calculations, the dynamics between different ring conformations in tenapanor is predicted to be fast. Addition of HCl to tenapanor results in general in a mixture of protonated ring conformers and hence a statistical mix of diastereoisomers which builds up the amorphous form, a-2HCl. This was qualitatively verified by 13C CP/MAS NMR investigations of the amorphous form. Thus, to form any significant amount of the crystalline material 2HCl, which originates from the minor (i.e., energetically less stable) ring conformations, one needs to involve nitrogen deprotonation to allow exchange between minor and major conformations of ANHY in solution. Thus, by controlling the solution pH value to well below the pKa of ANHY, the equilibrium between ANHY and 2HCl can be controlled and by this mechanism the crystallization of 2HCl can be avoided and the amorphous form of the dichloride salt can therefore be stabilized

Furosemide’s One Little Hydrogen Atom: NMR Crystallography Structure Verification of Powdered Molecular Organics

The potential of NMR crystallography to verify molecular crystal structures deposited in structural databases is evaluated, with two structures of the pharmaceutical furosemide serving as examples. While the structures differ in the placement of one H atom, using this approach, we verify one of the structures in the Cambridge Structural Database using quantitative tools, while establishing that the other structure does not meet the verification criteria

Powder crystallography of pharmaceutical materials by combined crystal structure prediction and solid-state 1H NMR spectroscopy

A protocol for the ab initio crystal structure determination of powdered solids at natural isotopic abundance by combining solid-state NMR spectroscopy, crystal structure prediction, and DFT chemical shift calculations was evaluated to determine the crystal structures of four small drug molecules: cocaine, flutamide, flufenamic acid, and theophylline. For cocaine, flutamide and flufenamic acid, we find that the assigned 1H isotropic chemical shifts provide sufficient discrimination to determine the correct structures from a set of predicted structures using the root-mean-square deviation (rmsd) between experimentally determined and calculated chemical shifts. In most cases unassigned shifts could not be used to determine the structures. This method requires no prior knowledge of the crystal structure, and was used to determine the correct crystal structure to within an atomic rmsd of less than 0.12 Å with respect to the known reference structure. For theophylline, the NMR spectra are too simple to allow for unambiguous structure selection

Exploring Systematic Discrepancies in DFT Calculations of Chlorine Nuclear Quadrupole Couplings

Previous studies have revealed significant discrepancies between density functional theory (DFT)-calculated and experimental nuclear quadrupolar coupling constants (CQ) for chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effects of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the quadrupolar coupling calculations, and that functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data

Polymorphs of Theophylline Characterized by DNP Enhanced Solid-State NMR

We show how dynamic nuclear polarization (DNP) enhanced solid-state NMR spectroscopy can be used to characterize polymorphs and solvates of organic solids. We applied DNP to three polymorphs and one hydrated form of the asthma drug molecule theophylline. For some forms of theophylline, sample grinding and impregnation with the radical-containing solution, which are necessary to prepare the samples for DNP, were found to induce polymorphic transitions or desolvation between some forms. We present protocols for sample preparation for solid-state magic-angle spinning (MAS) DNP experiments that avoid the polymorphic phase transitions in theophylline. These protocols include cryogrinding, grinding under inert atmosphere, and the appropriate choice of the impregnating liquid. By applying these procedures, we subsequently demonstrate that two-dimensional correlation experiments, such as H-1-C-13 and H-1-N-15 HETCOR or C-13-C-13 INADEQUATE, can be obtained at natural isotopic abundance in reasonable times, thus enabling more advanced structural characterization of polymorphs

Weak Halogen Bonding in Solid Haloanilinium Halides Probed Directly via Chlorine-35, Bromine-81, and Iodine-127 NMR Spectroscopy

A series of monohaloanilinium halides exhibiting weak halogen bonding (XB) has been prepared and characterized by ³⁵Cl, ⁸¹Br, and ¹²⁷I solid-state nuclear magnetic resonance (SSNMR) spectroscopy in magnetic fields of up to 21.1 T. The quadrupolar and chemical shift (CS) tensor parameters for halide ions (Cl^–, Br^–, I^–) which act as electron density donors in the halogen bonds of these compounds are measured to provide insight into the possible relationship between halogen bonding and NMR observables. The NMR data for certain series of related compounds are strongly indicative of when such compounds pack in the same space group, thus providing practical structural information. Careful interpretation of the NMR data in the context of novel and previously available X-ray crystallographic data, and new gauge-including projector-augmented-wave density functional theory (GIPAW DFT) calculations has revealed several notable trends. When a series of related compounds pack in the same space group, it has been possible to interpret trends in the NMR data in terms of the strength of the halogen bond. For example, in isostructural series, the halide quadrupolar coupling constant was found to increase as the halogen bond weakens. In the case of a series of haloanilinium bromides, the ⁸¹Br isotropic chemical shift and CS tensor span both decrease as the bromide–halogen XB is weakened. These trends were reproduced using both GIPAW DFT and cluster-model calculations of the bromide ion magnetic shielding tensor. Such trends are particularly exciting given the well-known role that NMR has played historically in the characterization of hydrogen bonding