Solid-state vibrational circular dichroism for pharmaceutical applications: Polymorphs and cocrystal of sofosbuvir

X-ray diffraction is a commonly used technique in the pharmaceutical industry for the determination of the atomic and molecular structure of crystals. However, it is costly, sometimes time-consuming, and it requires a considerable degree of expertise. Vibrational circular dichroism (VCD) spectroscopy resolves these limitations, while also exhibiting substantial sensitivity to subtle modifications in the conformation and molecular packaging in the solid state. This study showcases VCD's ability to differentiate between various crystal structures of the same molecule (polymorphs, cocrystals). We examined the most effective approach for producing high-quality spectra and unveiled the intricate link between structure and spectrum via quantum-chemical computations. We rigorously assessed, using alanine as a model compound, multiple experimental conditions on the resulting VCD spectra, with the aim of proposing an optimal and efficient procedure. The proposed approach, which yields reliable, reproducible, and artifact-free results with maximal signal-to-noise ratio, was then validated using a set comprising of three amino acids (serine, alanine, tyrosine), one hydroxy acid (tartaric acid), and a monosaccharide (ribose) to mimic active pharmaceutical components. Finally, the optimized approach was applied to distinguish three polymorphs of the antiviral drug sofosbuvir and its cocrystal with piperazine. Our results indicate that solid-state VCD is a prompt, cost-effective, and easy-to-use technique to identify crystal structures, demonstrating potential for application in pharmaceuticals. We also adapted the cluster and transfer approach to calculate the spectral properties of molecules in a periodic crystal environment. Our findings demonstrate that this approach reliably produces solid-state VCD spectra of model compounds. Although for large molecules with many atoms per unit cell, such as sofosbuvir, this approach has to be simplified and provides only a qualitative match, spectral calculations, and energy analysis helped us to decipher the observed differences in the experimental spectra of sofosbuvir.

Vibrational Circular Dichroism Spectroscopy with a Classical Polarizable Force Field: Alanine in the Gas and Condensed Phases

Polarizable force fields are an essential component for the chemically accurate modeling of complex molecular systems with a significant degree of fluxionality, beyond harmonic or perturbative approximations. In this contribution we examine the performance of such an approach for the vibrational spectroscopy of the alanine amino acid, in the gas and condensed phases, from the Fourier transform of appropriate time correlation functions generated along molecular dynamics (MD) trajectories. While the infrared (IR) spectrum only requires the electric dipole moment, the vibrational circular dichroism (VCD) spectrum further requires knowledge of the magnetic dipole moment, for which we provide relevant expressions to be used with polarizable force fields. The AMOEBA force field was employed here to model alanine in the neutral and zwitterionic isolated forms, solvated by water or nitrogen, and as a crystal. Within this framework, comparison of the electric and magnetic dipole moments to those obtained with nuclear velocity perturbation theory based on density-functional theory for the same MD trajectories are found to agree well with one another. The statistical convergence of the IR and VCD spectra is examined and found to be more demanding in the latter case. Comparisons with experimental frequencies are also provided for the condensed phases.

Differentiation of solvatomorphs of active pharmaceutical ingredients (API) by solid-state vibrational circular dichroism (VCD)

Here, we present the new application of solid-state Vibrational Circular Dichroism (VCD) spectroscopy to differentiate several dutasteride (DS) solvatomorphs - the model active pharmaceutical ingredient (API). Several crystalline DS hydrochloride hydrates solvated with methanol, ethanol, acetonitrile, acetone, and acetic acid were prepared. In contrast to almost identical IR spectra, the VCD ones were very sensitive to changes in the sample composition. We marked significant differences in the shape of VCD spectra of studied DS solvatomorphs, DS hydrates, and DS polymorphic forms. Our findings, supported by DFT calculations, show that VCD spectroscopy has the pronounced ability to distinguish their crystal arrangements. We believe that this contribution will extend the use of VCD in the pharmaceutical industry for developing and designing new chiral drug products for the identification, description, and in-depth probing of several pharmaceutical solvatomorphs in the future.

Exploring the Influence of the Phosphorus-Heteroatom Substitution in Nicotine on Its Electronic and Vibrational Spectroscopic Properties

The influence of phosphorus substitution of nitrogen in heterocyclic compounds on the vibrational spectroscopy as well as frontier molecular orbitals are analyzed. Nicotine with two nitrogen atoms in its structure is taken as the sample system to be studied computationally. By replacing the nitrogen atom in one or both rings of this molecule with phosphorus, three nicotine derivatives are created. The vibrational circular dichroism and infrared spectra of these four molecules in their monomer state, as well as the assemblies up to trimers are determined. The aforementioned spectra are calculated using static quantum chemical calculations employing a cluster-weighted approach. The calculated gas phase spectra of nicotine are compared to their respective experimental spectra. It is observed that the nicotine derivatives with phosphorus in the methylpyrrolidine ring have considerably different gas phase and bulk phase vibrational circular dichroism spectra when compared to nicotine. The phosphorus substitution reduces the energy gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital as well as altering the polarizability and reactivity of the investigated molecules.

VCD spectra of chiral naphthalene-1-carboxamides in the solid-state

The VCD spectra of chiral 2,3-dihydro-1H-benzo[de]isoquinolin-1-one (8-substituted naphthalene-1-carboxamide, BIQ) were studied in KBr pellets. The X-ray diffractometry revealed that the Me, Ph, and pClPh BIQs crystalize in the monoclinic P21, while nBu, pMePh, and oMeOPh BIQs in the orthorhombic P212121 space group. Only the Me-BIQ crystal exhibits the presence of cyclic amide dimers, while the others contain chains of the amid group hydrogen bonds. For all BIQs, except pMePh, the most intense IR band in the 1750-1550 cm-1 region is located at ca. 1680 cm-1 and is accompanied by two weak ones at ca. 1618 and 1590 cm-1. For the pMePh derivative, four almost equally intense IR bands at 1662, 1639, 1614, and 1588 cm-1 are observed. This region of the IR spectra of BIQs, but pMePh, is well reproduced by calculations based on BIQ monomers. On the other hand, the complex IR pattern of pMePh is computationally reproduced when larger crystal fragments, like octamers, are considered. Registration of the VCD spectra enabled recognizing the complexity of IR contours at ca. 1680 cm-1 by the corresponding VCD motives. For (i) Me, Ph and pClPh (R)-enantiomers, two (+)(-) bands were distinguished and for (ii) nBu and pMePh ones, one VCD band with right-side asymmetry was found. For (iii) oMeOPh the VCD pattern cannot be unambiguously assigned. Thus, the VCD spectra in the ν(C=O) range diverse the studied compounds. Among the set of molecules, pMePh has exceptional crystal geometry. Therefore, its most intense ν(C=O) band position and shape can be connected with the geometry of the hydrogen bonds, interactions, and crystal packing. Interpretation of the VCD spectra is based on linear and packed BIQ octamers. This cluster model can reproduce the main features of the solid-state VCD of BIQs.

Matrix effects in MI-VCD spectra of two chiral oxiranes and their potential microscopic origin

Combining vibrational circular dichroism (VCD) spectroscopy with the matrix isolation (MI) technique opens up interesting possibilities to study chiral molecules. MI involves the isolation of guest species in inert solid matrices at cryogenic temperatures. Hence, MI-VCD measures are solid-state VCD measurements, and as such, can suffer from mostly birefringance-related artefacts in the same way as common solid-state VCD measurements. In this contribution, we demonstrate that the sample preparation condition have tremendous impact on the quality and reliability of the recorded MI-VCD spectra. While MI-IR spectra are basically blind to these artefacts, the variation of deposition temperatures and host flow rates seem to control whether high quality MI-VCD spectra are obtained or if depolarization effects lead to completely obscured spectra. For two selected examples, styrene oxide (SO) and 1-phenyl propylene oxide (PPO), we discuss how the various experimental conditions may lead to the aforementioned effects and give a microscopic description of their origin.

Vibrational Circular Dichroism Spectroscopy of Chiral Molecular Crystals: Insights from Theory

Chirality is a curious phenomenon that appears in various forms. While the concept of molecular (RS-)chirality is ubiquitous in chemistry, there are also more intricate forms of structural chirality. One of them is the enantiomorphism of crystals, especially molecular crystals, that describes the lack of mirror symmetry in the unit cell. Its relation to molecular chirality is not obvious, but still an open question, which can be addressed with chiroptical tools. Vibrational circular dichroism (VCD) denotes chiral infrared (IR) spectroscopy that is susceptible to both, the molecular as well as the intermolecular space by means of vibrational transitions. When carried out in the solid state, VCD delivers a very rich set of non-local contributions that are determined by crystal packing and collective motion. Since its discovery in the 1970s, VCD has become the method of choice for the determination of absolute configurations, but its applicability reaches beyond towards the study of different crystal forms and polymorphism. This brief review summarises the theoretical concepts of crystal chirality and how computations of solid-state VCD can shed light into the intimate connection of chiral structure and vibrational optical activity.

How Crystal Symmetry Dictates Non-Local Vibrational Circular Dichroism in the Solid State

Solid-State Vibrational Circular Dichroism (VCD) can be used to determine the absolute structure of chiral crystals, but its interpretation remains a challenge in modern spectroscopy. In this work, we investigate the effect of a twofold screw axis on the solid-state VCD spectrum in a combined experimental and theoretical analysis of P2₁ crystals of (S)-(+)-1-indanol. Even though the space group is achiral, a single proper symmetry operation has an important impact on the VCD spectrum, which reflects the supramolecular chirality of the crystal. Distinguishing between contributions originating from molecular chirality and from chiral crystal packing, we find that while IR absorption hardly depends on the symmetry of the space group, the situation is different for VCD, where completely new non-local patterns emerge. Understanding the two underlying mechanisms, namely gauge transport and direct coupling, will help to use VCD to distinguish polymorphic forms.

How change in chirality prevents β-amyloid type interaction in a protonated cyclic dipeptide dimer

The protonated dimers of the diketopiperazine dipeptide cyclo (LPhe-LHis) and cyclo (LPhe-DHis) are studied by laser spectroscopy combined with mass spectrometry to shed light on the influence of stereochemistry on the clustering propensity of cyclic dipeptides. The marked spectroscopic differences experimentally observed in the hydride stretch region are well accounted for by the results of DFT calculations. Both diastereomeric protonated dimers involve a strong ionic hydrogen bond from the protonated imidazole ring of one monomer to the neutral imidazole nitrogen of the other. While this strong interaction is accompanied by a single NH⋯O hydrogen bond between the amide functions of the two moieties for the protonated dimer of cyclo (LPhe-DHis), that of cyclo (LPhe-LHis) involves two NH⋯O interactions, forming the motif of an antiparallel β sheet. Therefore, a change in chirality of the residue prevents the formation of the β sheet pattern observed in the amyloid type aggregation. These results emphasize the peculiar role of the histidine residue in peptide structure and interaction.

Unambiguous Stereochemical Assignment of Cyclo(Phe-Pro), Cyclo(Leu-Pro), and Cyclo(Val-Pro) by Electronic Circular Dichroic Spectroscopy

2,5-diketopiperazines (DKPs) are cyclic dipeptides ubiquitously found in nature. In particular, cyclo(Phe-Pro), cyclo(Leu-Pro), and cyclo(Val-Pro) are frequently detected in many microbial cultures. Each of these DKPs has four possible stereoisomers due to the presence of two chirality centers. However, absolute configurations of natural DKPs are often ambiguous due to the lack of a simple, sensitive, and reproducible method for stereochemical assignment. This is an important problem because stereochemistry is a key determinant of biological activity. Here, we report a synthetic DKP library containing all stereoisomers of cyclo(Phe-Pro), cyclo(Leu-Pro), and cyclo(Val-Pro). The library was subjected to spectroscopic characterization using mass spectrometry, NMR, and electronic circular dichroism (ECD). It turned out that ECD can clearly differentiate DKP stereoisomers. Thus, our ECD dataset can serve as a reference for unambiguous stereochemical assignment of cyclo(Phe-Pro), cyclo(Leu-Pro), and cyclo(Val-Pro) samples from natural sources. The DKP library was also subjected to a biological screening using assays for E. coli growth and biofilm formation, which revealed distinct biological effects of cyclo(D-Phe-L-Pro).

Competition between inter and intramolecular hydrogen bond evidenced by vibrational circular dichroism spectroscopy: The case of (1S,2R)-(-)-cis-1-amino-2-indanol

The infrared (IR) absorption and vibrational circular dichroism (VCD) spectra of an intramolecularly hydrogen-bonded chiral amino-alcohol, (1S,2R)-(-)-cis-1-amino-2-indanol, are studied in DMSO-d₆ . The spectra are simulated at the density functional theory (DFT) level within the frame of the cluster-in-the-liquid model. Both IR and VCD spectra show a clear signature of the formation of intermolecular hydrogen bonds at the detriment of the intramolecular OH … N interaction present in the isolated molecule. Two solvent molecules are necessary to reproduce the experimental spectra. Whereas the first DMSO molecule captures the main spectral modifications due to hydrogen bond formation between the solute and the solvent, the second DMSO molecule is necessary for a good description of the Boltzmann contribution of the different complexes, based on their Gibbs free energy.

Computation of Solid-State Vibrational Circular Dichroism in the Periodic Gauge

We introduce a new theoretical formalism to compute solid-state vibrational circular dichroism (VCD) spectra from molecular dynamics simulations. Having solved the origin-dependence problem of the periodic magnetic gauge, we present IR and VCD spectra of (1S,2S)-trans-1,2-cyclohexanediol obtained from first-principles molecular dynamics calculations and nuclear velocity perturbation theory, along with the experimental results. Because the structure model imposes periodic boundary conditions, the common origin of the rotational strength has hitherto been ill-defined and was approximated by means of averaging multiple origins. The new formalism reconnects the periodic model with the finite physical system and restores gauge freedom. It nevertheless fully accounts for nonlocal spatial couplings from the gauge transport term. We show that even for small simulation cells the rich nature of solid-state VCD spectra found in experiments can be reproduced to a very satisfactory level.

Efficient incorporation and protection of lansoprazole in cyclodextrin metal-organic frameworks

Lansoprazole (LPZ) is an acid pump inhibitor, which readily degrades upon acidic or basic conditions and under heating. We investigated here LPZ stability upon incorporation in particles made of cyclodextrin metal-organic frameworks (CD-MOFs). LPZ loaded CD-MOFs were successfully synthesized, reaching high LPZ payloads of 23.2 ± 2.1 wt%, which correspond to a molar ratio of 1:1 between LPZ and γ-CD. The homogeneity of LPZ loaded CD-MOFs in terms of component distribution was confirmed by elemental mapping by STEM-EDX. Both CTAB, the surfactant used in the CD-MOFs synthesis, and LPZ compete for their inclusion in the CD cavities. CTAB allowed obtaining regular cubic particles of around 5 µm with 15 wt% residual CTAB amounts. When LPZ was incorporated, the residual CTAB amount was less than 0.1 wt%, suggesting a higher affinity of LPZ for the CDs than CTAB. These findings were confirmed by molecular simulations. Vibrational circular dichroism studies confirmed the LPZ incorporation inside the CDs. Solid-state NMR showed that LPZ was located in the CDs and that it remained intact even after three years storage. Remarkably, the CD-MOFs matrix protected the drug upon thermal decomposition. This study highlights the interest of CD-MOFs for the incorporation and protection of LPZ.

Assessing cluster models of solvation for the description of vibrational circular dichroism spectra: synergy between static and dynamic approaches

Solvation effects are essential for defining the shape of vibrational circular dichroism (VCD) spectra.