Navigating the Complex Solid Form Landscape of the Quercetin Flavonoid Molecule

Quercetin, a naturally occurring bioflavonoid substance widely used in the nutraceutical and food industries, exists in various solid forms that can have different physicochemical properties, thus impacting this compound's performance in various applications. In this work, we will clarify the complex solid-form landscape of this molecule. Two elusive isostructural solvates of quercetin were obtained from ethanol and methanol. The obtained crystals were characterized experimentally, but the crystallographic structure could not be solved due to their high instability. Nevertheless, the desolvated structure resulting from a high-temperature treatment (or prolonged storage at ambient conditions) of both these two labile crystals was characterized and solved via powder X-ray diffraction and solid-state nuclear magnetic resonance (SSNMR). This anhydrous crystal structure was compared with another anhydrous quercetin form obtained in our previous work, indicating that, at least, two different anhydrous polymorphs of quercetin exist. Navigating the solid-form landscape of quercetin is essential to ensure accurate control of the functional properties of food, nutraceutical, or pharmaceutical products containing crystal forms of this substance.

The Raman Active Vibrations of Flavone and Quercetin: The Impact of Conformers and Hydrogen Bonding on Fingerprint Modes

The detection and analysis of flavonoids by Raman spectroscopy are of interest in many fields, including medicinal chemistry, food science, and astrobiology. Spectral interpretation would benefit from better identification of the fingerprint vibrational peaks of different flavonoids and how they are affected by intermolecular interactions. The Raman spectra of two flavonoids, flavone and quercetin, were investigated through comparisons between spectra recorded from pure powders and spectra calculated with time dependent density functional theory (TDDFT). For both flavone and quercetin, 17 peaks were assigned to specific molecular vibrations. Both flavonoids were found to have a split peak between 1250-1350 cm^-1 that is not predicted by TDDFT calculations on isolated molecules. In each case, it is shown that the addition of hydrogen bonded molecules arranged based on crystal structures reproduces the split peaks. These peaks were due to a stretching vibration of the bond between benzopyrone and phenyl rings and represent a characteristic spectral feature of flavonoids. Spectra of pollen grains from Quercus virginiana were also recorded and exhibit several peaks that correspond to the quercetin spectrum.

Hydrogen-Mediated Noncovalent Interactions in Solids: What Can NMR Crystallography Tell About?

Hydrogen atoms play a crucial role in the aggregation of organic (bio)molecules through diverse number of noncovalent interactions that they mediate, such as electrostatic in proton transfer systems, hydrogen bonding, and CH-π interactions, to mention only the most prominent. To identify and adequately describe such low-energy interactions, increasingly sensitive methods have been developed over time, among which quantum chemical computations have witnessed impressive advances in recent years. For reaching the present state-of-the-art, computations had to rely on a pool of relevant experimental data, needed at least for validation, if not also for other purposes. In the case of molecular crystals, the best illustration for the synergy between computations and experiment is given by the so-called NMR crystallography approach. Originally designed to increase the confidence level in crystal structure determination of organic compounds from powders, NMR crystallography is able now to offer also a wealth of information regarding the noncovalent interactions that drive molecules to pack in a given crystalline pattern or another. This is particularly true for the noncovalent interactions which depend on the exact location of labile hydrogen atoms in the system: in such cases, NMR crystallography represents a valuable characterization tool, in some cases complementing even the standard single-crystal X-ray diffraction technique. A concise introduction in the field is made in this mini-review, which is aimed at providing a comprehensive picture with respect to the current accuracy level reached by NMR crystallography in the characterization of hydrogen-mediated noncovalent interactions in organic solids. Different types of practical applications are illustrated with the example of molecular crystals studied by our research group, but references to other representative developments reported in the literature are also made. By summarizing the major concepts and methodological progresses, the present work is also intended to be a guide to the practical potential of this relatively recent analytical tool for the scientists working in areas where crystal engineering represents the main approach for rational design of novel materials.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

Crystal structure determination of an elusive methanol solvate – hydrate of catechin using crystal structure prediction and NMR crystallography

A useful short-cut was developed to limit the number of molecular conformations that need to be regarded in crystal structure prediction calculations, which led to the crystal structure determination of new methanol solvate – hydrate of catechin.

An NMR crystallography investigation of furosemide

This paper presents an NMR crystallography study of three polymorphs of furosemide. Experimental magic-angle spinning (MAS) solid-state NMR spectra are reported for form I of furosemide, and these are assigned using density-functional theory (DFT)-based gauge-including projector augmented wave (GIPAW) calculations. Focusing on the three known polymorphs, we examine the changes to the NMR parameters due to crystal packing effects. We use a recently developed formalism to visualise which regions are responsible for the chemical shielding of particular sites and hence understand the variation in NMR parameters between the three polymorphs.

Quantitative description of 1H SQ and DQ coherences for the hydroxyl disorder within hydrous ringwoodite

Proton-containing point defects in solid materials are important for a variety of properties ranging from ionic transport over thermal conductivity up to compressibility. Ultrafast magic-angle spinning techniques nowadays offer high-resolution solid-state NMR spectra, even for 1H, and thus open up possibilities to study the underlying defect chemistry. Nevertheless, disorder within such defects again leads to heavy spectral overlap of 1H resonances, which prevents quantitative analysis of defect concentrations, if several defect types are present. Here, we present a strategy to overcome this limitation by simulating the 1H lineshape as well as 1H-1H double-quantum buildup curves, which we then validate against the experimental data in a joint cost function. To mimic the local structural disorder, we use molecular dynamics simulations at the DFT level. It turned out to be advantageous for the joint refinement to put the computational effort into the structural optimisation to derive accurate proton positions and to use empirical correlations for the relation between isotropic and anisotropic 1H chemical shifts and structural elements. The expressiveness of this approach is demonstrated on ringwoodite's (γ-Mg2SiO4) OH defect chemistry containing four different defect types in octahedral and tetrahedral voids with both pure Mg and mixed Si and Mg cation environments. Still, we determine the ratio for each defect type with an accuracy of about 5% as a result of the minimization of the joint cost function. We expect that our approach is generally applicable for local proton disorder and might prove to be a valuable alternative to the established AIRSS and Monte Carlo methods, respectively.

Closer to the polydopamine structure: new insights from a combined 13C/1H/2H solid-state NMR study on deuterated samples

¹³C/¹H/²H ss-NMR on deuterated samples provide strong experimental evidence for the most probable monomer connectivity, π–π stacking, and the water dynamics in polydopamine.

Hydrogen Atomic Positions of O-H···O Hydrogen Bonds in Solution and in the Solid State: The Synergy of Quantum Chemical Calculations with ¹H-NMR Chemical Shifts and X-ray Diffraction Methods

The exact knowledge of hydrogen atomic positions of O-H···O hydrogen bonds in solution and in the solid state has been a major challenge in structural and physical organic chemistry. The objective of this review article is to summarize recent developments in the refinement of labile hydrogen positions with the use of: (i) density functional theory (DFT) calculations after a structure has been determined by X-ray from single crystals or from powders; (ii) ¹H-NMR chemical shifts as constraints in DFT calculations, and (iii) use of root-mean-square deviation between experimentally determined and DFT calculated ¹H-NMR chemical shifts considering the great sensitivity of ¹H-NMR shielding to hydrogen bonding properties.

Preparation of efficient quercetin delivery system on Zn-modified mesoporous SBA-15 silica carrier

Mesoporous silica material type SBA-15 was modified with different amounts of Zn (2 and 4wt.%) by incipient wetness impregnation method in ethanol. The parent, Zn-modified and quercetin loaded samples, were characterized by XRD, N₂ physisorption, TEM, thermal gravimetric analysis, UV-vis and FT-IR spectroscopies and in vitro release of quercetin at pH5.5 which is typical of dermal formulations. By this loading method anhydrous quercetin was formed on the silica carrier It was found that the different hydrate forms of quercetin (dihydrate, monohydrate, anhydrite) significantly influence the physico-chemical properties of the delivery system. It was found that hydrate forms of quercetin can be differentiated by XRD and by FT-IR spectroscopic methods. Thus, by evaluating the interaction of the drug with the silica carrier the changes due to its hydration state always have to be taken into account. Formation of Zn-quercetin complex was evidenced on zinc modified SBA-15 silica by FT-IR spectroscopy. High quercetin loading capacity (over 40wt.%) could be achieved on the parent and Zn-containing SBA-15 samples. The in-vitro release process at pH=5.5 showed slower quercetin release from Zn-modified SBA-15 samples compared to the parent one. Additionally, the comparative cytotoxic experiments evidenced that quercetin encapsulated in Zn-modified silica carriers has superior antineoplastic potential against HUT-29 cells compared to free drug. Zn-modified SBA-15 silica particles could be promising carriers for dermal delivery of quercetin.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Highly sensitive solid forms discrimination on the whole tablet of the active ingredients in quercetin dietary supplements by NMR crystallography approaches

Similarly to synthetic drugs, the exact crystalline form of active ingredients in solid formulations of dietary supplements may directly influence the dissolution rate, bioavailability, and stability of the final product, but this information is usually not provided by manufacturers. Working on the examples of two commercial quercetin dietary supplements a quick, reliable, and sensitive method is introduced for quercetin solid forms discrimination directly on the marketed products, without the need for prior sample preparation. It exploits the complementarity between solid-state Nuclear Magnetic Resonance (ss-NMR) and Powder X-Ray Diffraction (PXRD), which proved essential for performing a complete and accurate solid-state characterization of the two commercial products, and for obtaining new insights into the complex quercetin solid-forms landscape. The method can be readily generalized also to other dietary supplements based on bio-flavonoids/polyphenols.

NMR crystallography of ezetimibe co-crystals

An efficient, simplified protocol for solvent-drop assisted co-crystal preparation of ezetimibe (a drug for the treatment of primary hypercholesterolemia) with both imidazole and l-proline has been derived. The structures of the white powders were successfully solved via "NMR crystallography" combining solid-state NMR, powder X-ray diffraction and DFT chemical shift computations. Detailed insights into the likely crystallization mechanism were obtained from competition experiments, where efficient co-crystallization was feasible using ezetimibe monohydrate as precursor indicating that the crystal water acts as "molecular catalyst". It was also found that co-crystallization of imidazole is favored over l-proline, thus suggesting a clear preference of neutral hydrogen bonds compared to charge-assisted motifs.

Can the conformation of flexible hydroxyl groups be constrained by simple NMR crystallography approaches? The case of the quercetin solid forms

Hydrogen atoms in systems with many flexible hydroxyl side-groups are difficult to be exactly located from experimental X-Ray diffraction and/or solid-state NMR data, thus often leading to wrong conclusions with respect to the hydrogen bonding network established in crystal lattice. A simple computational method is proposed in the present work to tackle this problem, which may be readily incorporated in conventional NMR crystallography protocols. The method is based on ranking all possible conformations of the flexible hydroxyls according to their lattice energy in crystalline environments. Its effectiveness is investigated on two distinct solid forms of quercetin, for which only two out of the five hydroxyl side-groups can be well constrained from experimental/theoretical data. For this purpose, first-principle quantum-mechanical computations were combined with calculations at the molecular mechanics (MM) level of theory, and previous ss-NMR and X-Ray diffraction data. To assess accuracy in ranking the identified conformers, tests have been performed first on quercetin dihydrate, for which an X-Ray single-crystal structure is available. The possibility of applying this method in a real NMR crystallography context has been investigated finally on anhydrous quercetin, for which only powder X-Ray crystal structure has been reported so far.

Distinguishing tautomerism in the crystal structure of (Z)-N-(5-ethyl-2,3-dihydro-1,3,4-thiadiazol-2-ylidene)-4-methylbenzenesulfonamide using DFT-D calculations and (13)C solid-state NMR

The crystal structure of (Z)-N-(5-ethyl-2,3-di­hydro-1,3,4-thia­diazol-2-yl­idene)-4-methyl­benzene­sulfonamide contains an imine tautomer, rather than the previously reported amine tautomer. The tautomers can be distinguished using dispersion-corrected density functional theory calculations and by comparison of calculated and measured ¹³C solid-state NMR spectra.

The crystal structure of the title compound, C₁₁H₁₃N₃O₂S₂, has been determined previously on the basis of refinement against laboratory powder X-ray diffraction (PXRD) data, supported by comparison of measured and calculated ¹³C solid-state NMR spectra [Hangan et al. (2010 ▶). Acta Cryst. B66, 615–621]. The mol­ecule is tautomeric, and was reported as an amine tautomer [systematic name: N-(5-ethyl-1,3,4-thia­diazol-2-yl)-p-toluene­sulfonamide], rather than the correct imine tautomer. The protonation site on the mol­ecule’s 1,3,4-thia­diazole ring is indicated by the inter­molecular contacts in the crystal structure: N—H⋯O hydrogen bonds are established at the correct site, while the alternative protonation site does not establish any notable inter­molecular inter­actions. The two tautomers provide essentially identical Rietveld fits to laboratory PXRD data, and therefore they cannot be directly distinguished in this way. However, the correct tautomer can be distinguished from the incorrect one by previously reported qu­anti­tative criteria based on the extent of structural distortion on optimization of the crystal structure using dispersion-corrected density functional theory (DFT-D) calculations. Calculation of the ¹³C SS-NMR spectrum based on the correct imine tautomer also provides considerably better agreement with the measured ¹³C SS-NMR spectrum.

Polymorphs and co-crystals of haloprogin: an antifungal agent

Haloprogin is a widely used antifungal agent. Here we report the first polymorphs and halogen-bonded co-crystals ever described.