Elucidation of Structure and Dynamics in Solid Octafluoronaphthalene from Combined NMR, Diffraction, And Molecular Dynamics Studies

Unravelling Guest Dynamics in Crystalline Molecular Organics Using 2H Solid-State NMR and Molecular Dynamics Simulation

2H solid-state NMR and atomistic molecular dynamics (MD) simulations are used to understand the disorder of guest solvent molecules in two cocrystal solvates of the pharmaceutical furosemide. Traditional approaches to interpreting the NMR data fail to provide a coherent model of molecular behavior and indeed give misleading kinetic data. In contrast, the direct prediction of the NMR properties from MD simulation trajectories allows the NMR data to be correctly interpreted in terms of combined jump-type and libration-type motions. Time-independent component analysis of the MD trajectories provides additional insights, particularly for motions that are invisible to NMR. This allows a coherent picture of the dynamics of molecules restricted in molecular-sized cavities to be determined.

Exploring the Potential of Multinuclear Solid-State 1 H, 13 C, and 35 Cl Magnetic Resonance To Characterize Static and Dynamic Disorder in Pharmaceutical Hydrochlorides

Crystallographic disorder, whether static or dynamic, can be detrimental to the physical and chemical stability, ease of crystallization and dissolution rate of an active pharmaceutical ingredient. Disorder can result in a loss of manufacturing control leading to batch-to-batch variability and can lengthen the process of structural characterization. The range of NMR active nuclei makes solid-state NMR a unique technique for gaining nucleus-specific information about crystallographic disorder. Here, we explore the use of high-field ³⁵ Cl solid-state NMR at 23.5 T to characterize both static and dynamic crystallographic disorder: specifically, dynamic disorder occurring in duloxetine hydrochloride (1), static disorder in promethazine hydrochloride (2), and trifluoperazine dihydrochloride (3). In all structures, the presence of crystallographic disorder was confirmed by ¹³ C cross-polarization magic-angle spinning (CPMAS) NMR and supported by GIPAW-DFT calculations, and in the case of 3, ¹ H solid-state NMR provided additional confirmation. Applying ³⁵ Cl solid-state NMR to these compounds, we show that higher magnetic fields are beneficial for resolving the crystallographic disorder in 1 and 3, while broad spectral features were observed in 2 even at higher fields. Combining the data obtained from ¹ H, ¹³ C, and ³⁵ Cl NMR, we show that 3 exhibits a unique case of disorder involving the ⁺ N-H hydrogen positions of the piperazinium ring, driving the chloride anions to occupy three distinct sites.

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.

Collective dynamics of molecular rotors in periodic mesoporous organosilica: a combined solid-state 2H-NMR and molecular dynamics simulation study

Molecular rotors offer a platform to realize controlled dynamics and modulate the functions of solids. The motional mechanisms in arrays of rotors have not been explored in depth. Crystal-like porous organosilicas, comprising p-phenylene rotators pivoted onto a siloxane scaffold, were modelled using molecular dynamics (MD) simulations. Long simulations, on a microsecond scale, allowed to follow the reorientation statistics of rotor collections and single out group configurations and frequency distributions as a function of temperature. The motions observed in the MD simulations support a multiple-site model for rotor reorientations. Computed motional frequencies revealed a complex rotatory phenomenon combining an ultra-fast libration motion (oscillation up to 30°) with a slow and fast 180° flip reorientation. Adopting a multiple-site model provides a more accurate simulation of the ²H-NMR spectra and a rationalization of their temperature dependence. In particular, rotators endowed with distinct rates could be explained by the presence of slower rings locked in a T-shaped conformation.

Correlated motion and mechanical gearing in amphidynamic crystalline molecular machines

In this review we highlight the recent efforts towards the development of molecular gears with an emphasis on building molecular gears in the solid state and the role that molecular gearing and correlated motions may play in the function of crystalline molecular machines. We discuss current molecular and crystal engineering strategies, challenges associated with engineering correlated motion in crystals, and outline experimental and theoretical tools to explore gearing dynamics while highlighting key advances made to date.

Nanoscale rotational dynamics of four independent rotators confined in crowded crystalline layers

We report a study where Car-Parrinello molecular dynamics simulations and variable-temperature (30-300 K) ¹H spin-lattice relaxation time experiments nicely complement each other to characterize the dynamics within a set of four crystalline 1,4-diethynylbicyclo[2.2.2]octane (BCO) rotors assembled in the metal-organic rotor, {Li⁺₄(^-CO₂-Ph-BCO-py)₄(H₂O)₈}·2DMF. The remarkable finding of this work is that, despite the individual rotational barriers of four rotors being indiscernible and superimposed in a broad relaxation process, we were able to unravel a strongly interrelated series of rotational motions involving disrotatory and conrotatory motions in pairs as well as rotational steps of single rotators, all three processes with similar, sizeable rotational barriers of 6 kcal mol^-1. It is noteworthy that DFT molecular dynamics simulations and variable-temperature (30-300 K) proton spin-lattice relaxation time experiments deliver the same high value for the rotational barriers stressing the potential of the combined use of the two techniques in understanding rotational motion at the nanoscale.

Perspective: Current advances in solid-state NMR spectroscopy

In contrast to the rapid and revolutionary impact of solution-state Nuclear Magnetic Resonance (NMR) on modern chemistry, the field of solid-state NMR has matured more slowly. This reflects the major technical challenges of much reduced spectral resolution and sensitivity in solid-state as compared to solution-state spectra, as well as the relative complexity of the solid state. In this perspective, we outline the technique developments that have pushed resolution to intrinsic limits and the approaches, including ongoing major developments in the field of Dynamic Nuclear Polarisation, that have enhanced spectral sensitivity. The information on local structure and dynamics that can be obtained using these gains in sensitivity and resolution is illustrated with a diverse range of examples from large biomolecules to energy materials and pharmaceuticals and from both ordered and highly disordered materials. We discuss how parallel developments in quantum chemical calculation, particularly density functional theory, have enabled experimental data to be translated directly into information on local structure and dynamics, giving rise to the developing field of "NMR crystallography."

A heuristic approach to evaluate peri interactions versus intermolecular interactions in an overcrowded naphthalene

Octachloronaphthalene (OCN), a serious environmental pollutant, has been investigated by charge density analysis to unravel several unexplored factors responsible for steric overcrowding. The topological features of the enigmatic peri interactions contributing to steric overcrowding are qualified and quantified from experimental and theoretical charge-density studies. A new facet in the fundamental understanding of peri interactions is revealed by NCI (non-covalent interaction) analysis. The potential role of these interactions in deforming the molecular geometry and subsequent effect on aromaticity are substantiated from NICS (Nuclear Independent Chemical Shift) and QTAIM (Quantum Theory of Atoms in Molecules) calculations. The eye-catching dissimilarity in the out-of-plane twisting of OCN renders the molecule in an asymmetric geometry in the crystalline phase compared with symmetric geometry in the optimized solvated phase. This is uniquely characterized by their molecular electrostatic potential (MESP), respectively, and is explained in terms of conflict between two opposing forces - peri interactions, and symbiotic intermolecular Cl⋯Cl and Cl⋯π contacts.

Characterising the role of water in sildenafil citrate by NMR crystallography

A combination of solid-state NMR techniques, including ¹³C/¹H correlation, ²H magic-angle spinning NMR and first principles calculation are employed to characterise the role of water in different hydration states of sildenafil citrate.

Perfluoro-1,1′-biphenyl and perfluoronaphthalene and their derivatives as π-acceptors for anions

Perfluoro-1,1′-biphenyl and perfluoronaphthalene are identified experimentally as π-acceptors for anions by studying the parent systems as well as cationic derivatives.

Well-tempered metadynamics as a tool for characterizing multi-component, crystalline molecular machines

The well-tempered, smoothly converging form of the metadynamics algorithm has been implemented in classical molecular dynamics simulations and used to obtain an estimate of the free energy surface explored by the molecular rotations in the plastic crystal, octafluoronaphthalene. The biased simulations explore the full energy surface extremely efficiently, more than 4 orders of magnitude faster than unbiased molecular dynamics runs. The metadynamics collective variables used have also been expanded to include the simultaneous orientations of three neighboring octafluoronaphthalene molecules. Analysis of the resultant three-dimensional free energy surface, which is sampled to a very high degree despite its significant complexity, demonstrates that there are strong correlations between the molecular orientations. Although this correlated motion is of limited applicability in terms of exploiting dynamical motion in octafluoronaphthalene, the approach used is extremely well suited to the investigation of the function of crystalline molecular machines.

Dynamic nuclear polarization NMR spectroscopy of microcrystalline solids

Dynamic nuclear polarization (DNP) solid-state NMR has been applied to powdered microcrystalline solids to obtain sensitivity enhancements on the order of 100. Glucose, sulfathiazole, and paracetamol were impregnated with bis-nitroxide biradical (bis-cyclohexyl-TEMPO-bisketal, bCTbK) solutions of organic solvents. The organic solvents were carefully chosen to be nonsolvents for the compounds, so that DNP-enhanced solid-state NMR spectra of the unaltered solids could be acquired. A theoretical model is presented that illustrates that for externally doped organic solids characterized by long spin-lattice relaxation times (T(1)((1)H) > 200 s), (1)H-(1)H spin diffusion can relay enhanced polarization over micrometer length scales yielding substantial DNP enhancements (ε). ε on the order of 60 are obtained for microcrystalline glucose and sulfathiazole at 9.4 T and with temperatures of ca. 105 K. The large gain in sensitivity enables the rapid acquisition of (13)C-(13)C correlation spectra at natural isotopic abundance. It is anticipated that this will be a general method for enhancing the sensitivity of solid-state NMR experiments of organic solids.

Vibrational averaging of the chemical shift in crystalline α-glycine

Averaging of the chemical shift over the molecular motion improves the simulated data and provides additional information about the temperature dependence and system dynamics. However, crystal modeling is difficult due to the limited precision of the plane-wave density functional theory (DFT) methods and approximate vibrational schemes. On the glycine example, we investigate how the averaging can be achieved within the periodic boundary conditions at the DFT level. The nuclear motion is modeled with the vibrational configuration interaction, with other simplified quantum anharmonic schemes, and the classical Born-Oppenheimer molecular dynamics (BOMD). The results confirm a large vibrational contribution to the isotropic shielding values. Both the first and second derivatives of the shielding were found important for the quantum averaging. The first derivatives influence the shielding mostly due to the anharmonic character of the CH and NH stretching modes, whereas second derivatives produce most vibrational corrections associated with the lower-frequency vibrational modes. Temperature excitations of the lowest-frequency vibrational states and the expansion of the crystal cell both determine the temperature dependence of nuclear magnetic resonance parameters. The vibrational quantum approach as well as classical BOMD schemes provided temperature dependencies of the chemical shifts that are consistent with the previous experimental data.