Modifying conformational distribution of chiral tetrahydro-2-furoic acid through its interaction with water: a rotational spectroscopic and theoretical investigation

Conformational distributions of tetrahydro-2-turoic acid in water at different pH values by their IR and vibrational circular dichroism spectra

Infrared (IR) and vibrational circular dichroism (VCD) spectra of tetrahydro-2-furoic acid (THFA) in aqueous solutions under several different pH conditions were recorded. To interpret the IR and VCD spectra of THFA obtained in highly acidic and basic aqueous solutions, extensive and systematic conformational searches were conducted to acquire the low-energy minima for both the neutral and deprotonated forms of THFA species, as well as their hydrated clusters. This was accomplished by using the conformer-rotamer ensemble sampling tool (CREST) with an implicit solvation model for water. The CREST candidates were further optimized at the B3LYP-D3BJ/def2-TZVP level of theory. The simulated VCD spectra of the neutral THFA conformers in the polarizable continuum model (PCM) of water alone exhibit little agreement with the experimental data under highly acidic conditions. Applying the clusters-in-a-liquid solvation model by considering the monohydrate THFA conformers in the PCM of water, significantly improved agreement with the experimental data. Similarly, the deprotonated THFA species solvated with one to four explicit water molecules in the PCM of water were considered. While the IR and VCD spectra of the deprotonated THFA monohydrate conformers offer the best agreement with the experimental data, other larger hydrated clusters, particularly the dihydrates, also contribute to the spectra. Through the synergistic combined experimental and theoretical approach used in the study, comprehensive conformational distributions of the predominant THFA species across various pH conditions were extracted.

Assignment of the absolute configuration of molecules that are chiral by virtue of deuterium substitution using chiral tag molecular rotational resonance spectroscopy

Chiral tag molecular rotational resonance (MRR) spectroscopy is used to assign the absolute configuration of molecules that are chiral by virtue of deuterium substitution. Interest in the improved performance of deuterated active pharmaceutical ingredients has led to the development of precision deuteration reactions. These reactions often generate enantioisotopomer reaction products that pose challenges for chiral analysis. Chiral tag rotational spectroscopy uses noncovalent derivatization of the enantioisotopomer to create the diastereomers of the 1:1 molecular complexes of the analyte and a small, chiral molecule. Assignment of the absolute configuration requires high-confidence determinations of the structures of these weakly bound complexes. A general search method, CREST, is used to identify candidate geometries. Subsequent geometry optimization using dispersion corrected density functional theory gives equilibrium geometries with sufficient accuracy to identify the isomers of the chiral tag complexes produced in the pulsed jet expansion used to introduce the sample into the MRR spectrometer. Rotational constant scaling based on the fact that the diastereomers have the same equilibrium geometry gives accurate predictions allowing identification of the homochiral and heterochiral tag complexes and, therefore, assignment of absolute configuration. The method is successfully applied to three oxygenated substrates from enantioselective Cu-catalyzed alkene transfer hydrodeuteration reaction chemistry.

Exploring the distinct conformational preferences of allyl ethyl ether and allyl ethyl sulfide using rotational spectroscopy and computational chemistry

The conformational energy landscapes of allyl ethyl ether (AEE) and allyl ethyl sulfide (AES) were investigated using Fourier transform microwave spectroscopy in the frequency range of 5-23 GHz aided by density functional theory B3LYP-D3(BJ)/aug-cc-pVTZ calculations. The latter predicted highly competitive equilibria for both species, including 14 unique conformers of AEE and 12 for the sulfur analog AES within 14 kJ mol-1. The experimental rotational spectrum of AEE was dominated by transitions arising from its three lowest energy conformers, which differ in the arrangement of the allyl side chain, while in AES, transitions due to the two most stable forms, distinct in the orientation of the ethyl group, were observed. Splitting patterns attributed to methyl internal rotation were analyzed for AEE conformers I and II, and the corresponding V3 barriers were determined to be 12.172(55) and 12.373(32) kJ mol-1, respectively. The experimental ground state geometries of both AEE and AES were derived using the observed rotational spectra of the 13C and 34S isotopic species and are highly dependent on the electronic properties of the linking chalcogen (oxygen vs sulfur). The observed structures are consistent with a decrease in hybridization in the bridging atom from oxygen to sulfur. The molecular-level phenomena that drive the conformational preferences are rationalized through natural bond orbital and non-covalent interaction analyses. These show that interactions involving the lone pairs on the chalcogen atom with the organic side chains favor distinct geometries and energy orderings for the conformers of AEE and AES.

Rotational Spectroscopy of 2-Furoic Acid and Its Dimer: Conformational Distribution and Double Proton Tunneling

Structural and tunneling properties of the 2-furoic acid (FA) monomer and dimer were investigated using rotational spectroscopy and DFT calculations. CREST, a conformational ensemble space exploration tool, was used to identify all possible low energy conformations of the FA monomer and dimer, followed by the DFT geometry optimization and harmonic frequency calculations. Broadband rotational spectra in the 2-6 GHz and in the 8-12 GHz regions were recorded in a supersonic jet expansion. The monomeric FA was found to exist dominantly as three different conformers: I , II , and III in a jet, with I and II taking on the trans -COOH configuration while III having the cis -COOH configuration. For the FA dimer, only the I - II conformer was observed experimentally, whereas the symmetric I - I and II - II conformers were not observed because of their zero dipole moments. The analysis of the splittings in the rotational transitions of I - II allowed one to extract the tunneling splitting to be 1056.0(12) MHz. The barrier height was determined to be ~442 cm -1 using the scaled potential energy scans at several different levels of theory.

Raman Optical Activity of N-Acetyl-L-Cysteine in water and in methanol: the "clusters-in-a-liquid" model and ab initio molecular dynamics simulations

Raman and Raman Optical Activity (ROA) spectra of N-acetyl-L-cysteine (NALC), a flexible chiral molecule, were measured in water and in methanol to evaluate the solvent effects. Two different solvation approaches, i.e. the DFT based clusters-in-a-liquid solvent model and the ab initio molecular dynamics (AIMD) simulations, were applied to simulate the Raman and ROA spectra. Systematic conformational searches were carried out using a recently developed conformational searching tool, CREST, with the inclusion of polarizable continuum model of water and of methanol. The CREST candidates of NALC and the NALC-solvent complexes were re-optimized and their Raman and ROA simulations were done at the B3LYP-D3BJ/def2-TZVP and the B3LYP-aug-cc-pVDZ//cc-pVTZ levels. Also, AIMD simulations , which includes some anharmonic effects and all intermolecular interactions in solution, were performed. By empirically weighting the computed Raman and ROA spectra of each conformer, good agreements with the experimental data were achieved with both approaches, while AIMD offered some improvements in the carbonyl and in the low wavenumber regions over the static DFT approach. The pros and cons of these two different approaches for accounting the solvent effects on Raman and ROA of this flexible chiral system will also be discussed.

Assignment-free chirality detection in unknown samples via microwave three-wave mixing

Straightforward identification of chiral molecules in multi-component mixtures of unknown composition is extremely challenging. Current spectrometric and chromatographic methods cannot unambiguously identify components while the state of the art spectroscopic methods are limited by the difficult and time-consuming task of spectral assignment. Here, we introduce a highly sensitive generalized version of microwave three-wave mixing that uses broad-spectrum fields to detect chiral molecules in enantiomeric excess without any prior chemical knowledge of the sample. This method does not require spectral assignment as a necessary step to extract information out of a spectrum. We demonstrate our method by recording three-wave mixing spectra of multi-component samples that provide direct evidence of enantiomeric excess. Our method opens up new capabilities in ultrasensitive phase-coherent spectroscopic detection that can be applied for chiral detection in real-life mixtures, raw products of chemical reactions and difficult to assign novel exotic species.

The structure of isolated thalidomide as reference for its chirality-dependent biological activity: a laser-ablation rotational study

Thalidomide is a drug that presents two enantiomers with markedly different pharmacological and toxicological activities. It is sadly famous due to its teratogenic effects mostly caused by the preferential docking of the (S)-enantiomer to the target protein cereblon (CRBN). To compare the structure of the bound CRBN thalidomide enantiomers with that of the isolated molecule, the rotational spectrum of laser-ablated thalidomide has been studied by chirp-pulsed Fourier transform microwave spectroscopy in supersonic jets complemented by theoretical computations. A new setup of the laser ablation nozzle used is presented. Two stable equatorial and axial conformers of thalidomide have been predicted corresponding to the two possible bent conformations exhibited by the glutarimide moiety. Only the most stable equatorial conformer has been detected. The comparison of its structure with those of the (S)- and (R)-enantiomers bound to CBRN shows that the bound (S) species is only slightly distorted. On the contrary, the bound (R)-enantiomer exhibits a highly distorted structure which affects the degree of puckering of the glutarimide ring and especially to the orientation of the phtalimide and glutarimide subunits. This is consistent with a less stable (R)-enantiomer and the known preference of (S)-thalidomide to bind CRBN, which starts the process leading to teratogenic effects.