Importance of Nuclear Quantum Effects for Molecular Cocrystals with Short Hydrogen Bonds

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.

Spectroscopically Visualizing the Evolution of Hydrogen-Bonding Interactions

Understanding chemical bond variations is the soul of chemistry as it is essential for any chemical process. The evolution of hydrogen bonds is one of the most fundamental and emblematic events during proton transfer; however, its experimental visualization remains a formidable challenge because of the transient timescales. Herein, by subtly regulating the proton-donating ability of distinct proton donors (zeolites or tungstophosphoric acid), a series of different hydrogen-bonding configurations were precisely manipulated. Then, an advanced two-dimensional (2D) heteronuclear correlation nuclear magnetic resonance (NMR) spectroscopic technique was utilized to simultaneously monitor the electronic properties of proton donors and acceptors (2-13C-acetone or trimethylphosphine oxide) through chemical shifts. Parabolic 1H-13C NMR relationships combined with single-well and double-well potential energy surfaces derived from theoretical simulations quantitatively identified the hydrogen bond types and allowed the evolution of hydrogen bonds to be visualized in diverse acid-base interaction complexes during proton transfer. Our findings provide a new perspective to reveal the nature and evolution of hydrogen bonds and confirm the superiority of 2D NMR techniques in identifying the subtle distinctions of various hydrogen-bonding configurations.

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.

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.