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.

Influence of microsolvation on vibrational circular dichroism spectra in dimethyl sulfoxide solvent: A Bottom-Up approach using Quantum cluster growth

Chiroptical spectroscopic measurements serve as routine methods to assign the absolute configuration of chiral compounds and interpret their conformational behavior in solution. One common challenge is the use of strongly hydrogen-bonding solvents, which can significantly bias the conformational ensemble and affect the vibrational circular dichroism (VCD) active bands in solution. One such solvent is dimethyl sulfoxide (DMSO)-an excellent solvent for stubborn compounds-that must be explicitly considered in VCD analysis. Explicit consideration of solvent remains a critical challenge in chiroptical spectroscopy due to the need to explore solute-solvent conformational space and the computational expense in modeling these clusters. Interested in the recent development of the Quantum Cluster Growth (QCG) program by the Grimme lab, we set out to model and interpret previously reported VCD spectra for several molecules using their efficient program. Our purposes are two-fold: (1) to investigate the applicability of the QCG program to the problem of reproducing VCD spectra in DMSO solvent and (2) to identify limitations in using this approach. We find that we can conveniently model and analyze the VCD spectra of investigated molecules in DMSO. However, the final set of conformers used for VCD calculations are functional dependent and different sets of conformers can provide satisfactory quantitative agreement between experimental and predicted VCD spectra. We hope that this study provides guidance for future chiroptical studies in the challenging DMSO solvent.

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.

Modelling solute-solvent interactions in VCD spectra analysis with the micro-solvation approach

Vibrational circular dichroism (VCD) spectroscopy has become an important part of the (stereo-)chemists' toolbox as a reliable method for the determination of absolute configurations. Being the chiroptical version of infrared spectroscopy, it has also been recognized as being very sensitive to conformational changes and intermolecular interactions. This sensitivity originates from the fact that the VCD spectra of individual conformers are often more different than their IR spectra, so that changes in conformational distributions or band positions and intensities become more pronounced. What is an advantage for studies focussing on intermolecular interactions can, however, quickly turn into a major obstacle during AC determinations: solute-solvent interactions can have a strong influence on spectral signatures and they must be accurately treated when simulating VCD and IR spectra. In this perspective, we showcase selected examples which exhibit particularly pronounced solvent effects. It is demonstrated that it is typically sufficient to model solute-solvent interactions by placing single solvent molecules near hydrogen bonding sites of the solute and subsequently use the optimized structures for spectra simulations. This micro-solvation approach works reasonably well for medium-sized, not too conformationally flexible molecules. We thus also discuss its limitations and outline the next steps that method development needs to take in order to further improve the workflows for VCD spectra predictions.

Chirality transfer observed in Raman optical activity spectra

Chirality transfer (also called induced chirality) is a phenomenon present in chiroptical spectra that manifests itself as a new band or bands of an achiral molecule interacting with a chiral one. In the Raman optical activity (ROA) spectroscopy, the bands of achiral solvents have been recently observed, but the latest papers have shown that they corresponded to the new ECD-Raman (eCP-Raman) effect. Here, we show an unambiguous example of chirality transfer observed in the ROA spectra. The spectra registered for the (1:1) mixtures of achiral benzonitrile with the enantiomers of 2,2,2-trifluoro-1-phenylethanol, 1-phenylethanol, and 1-phenylethylamine exhibited the v(CN) vibration band at about 2230 cm^-1. The ROA measurements were repeated several times to ensure the reliability of the phenomenon. Calculations revealed the CN···HO or CN···HNH hydrogen bond formation accompanied by the π···π or CH···π interactions. The interaction strength was shown to be an important factor for the pronouncement of the ROA chirality transfer effect.

7-Azaindole breaks carboxylic acid dimers and simplifies VCD spectra analyses of natural products

The determination of absolute configurations of carboxylic acids by vibrational circular dichroism (VCD) spectroscopy is often complicated by self-aggregation and the subsequent need to compute the spectra of the aggregates. We show that 7-azaindole effectively breaks up these aggregates by stronger complemental hydrogen bonding to the COOH moiety, enabling drastic simplification and acceleration of VCD spectra calculations.

Absolute Configuration of the Polyketide Natural Product (-)-Enterocin

The absolute configuration of the polyketide natural product (-)-enterocin was established by two independent approaches. In the first approach, synthetic enterocin with a defined configuration was compared to the natural product. While identical in all scalar properties, the compound displayed an opposite specific rotation and a different chiral HPLC retention time when compared with (-)-enterocin. In a second approach, the vibrational circular dichroism (VCD) of the natural product was measured and shown to be opposite to the calculated VCD of its enantiomer.

Dynamic Stereochemistry of a Biphenyl-Bisprolineamide Model Catalyst and its Imidazolidinone Intermediates

In this study, we characterize the dynamic stereochemistry of a biphenyl-2,2'-bis(proline amide) catalyst in chloroform and DMSO as representative weakly and strongly hydrogen bonding solvents. Using vibrational circular dichroism (VCD) spectroscopy and density functional theory (DFT) based spectra calculations, we show that the preferred axial stereochemistry of the catalyst is determined by solute-solvent interactions. Explicitly considering solvation with DMSO molecules is found to be essential to correctly predict the conformational preferences of the catalyst. Furthermore, we investigate the stereochemistry of the corresponding enamines and imidazolidinones that are formed upon reaction with isovaleraldehyde. The enamines are found to rapidly convert to endo-imidazolidinones and the thermodynamically favored exo-imidazolidinones are formed only slowly. The present study demonstrates that the stereochemistry of these imidazolidinones can be deduced directly from the VCD spectra analysis without any further detailed analysis of NMR spectra. Hence, we herein exemplify the use of VCD spectroscopy for an in situ characterization of intermediates relevant in asymmetric catalysts.

VCD spectroscopy reveals conformational changes of chiral crown ethers upon complexation of potassium and ammonium cations

Two chiral derivatives of 18-crown-6, namely the host molecules 2,3-diphenyl- and 2-phenyl-18c6, serve as model systems to investigate whether VCD spectroscopy can be used to monitor conformational changes occurring upon complexation of guests. Host-guest complexes of both crown ethers were prepared by addition of KNO₃. The more bulky 2,3-diphenyl-18c6 is found to undergo major conformational changes upon encapsulation of K⁺, which are revealed as characteristic changes of the VCD spectral signatures. In contrast, while 2-phenyl-18c6 also incorporates K⁺ into the macrocycle, strong conformational changes are not occurring and thus spectral changes are negligible. With an octyl ammonium cation as guest molecule, 2,3-diphenyl-18c6 shows the same conformational and spectral changes that were observed for K⁺-complexes. In addition, the asymmetric NH₃-deformation modes are found to gain VCD intensity through an induced VCD process. An analysis of the vibrational spectra enables a differentiation of VCD active and inactive guest modes: There appears to be a correlation between the symmetry of the vibrational mode and the induced VCD intensity. While this finding makes the host-guest complexes interesting systems for future theoretical studies on the origin of induced VCD signatures, the observations described in this study demonstrate that VCD spectroscopy is indeed a suitable technique for the characterization of supramolecular host-guest complexes.

Insights in the vibrational optical activity spectra of the antibiotic vancomycin in DMSO

Vibrational circular dichroism (VCD) and Raman optical activity (ROA) are two spectroscopic techniques that are sensitive towards the conformational behaviour of molecules, and are often complementary herein. In this work we pursue the determination of the conformational ensemble of the antibiotic glycopeptide vancomycin in DMSO through comparison of experimental and computational spectra, both for VCD and ROA. ROA is found to be highly suitable for the task, identifying an ensemble that strongly resembles the NMR conformation. In the case of VCD, however, a too high sensitivity of the intensities with respect to minor conformational changes hampers a reliable conformational analysis. Whence attempting to improve the match between the VCD experiment and calculations by any means - e.g., by inducing minor conformational changes or including solvent effects in the calculations - we show that there is the risk of going down the rabbit hole. In conclusion, this work contributes to the broader understanding of where, when and how VCD and ROA can be deployed as techniques for conformational analysis.

Solvation of the Boc-Val-Phe-nPr peptide characterized by VCD spectroscopy and DFT calculations

The conformational preferences of peptides are strongly determined by hydrogen bonding interactions. Intermolecular solute-solvent interactions compete with intramolecular interactions, which typically stabilize the secondary structure of the peptide. The analysis of vibrational circular dichroism (VCD) spectra can give insights into solvation-induced changes in the conformational distribution of small peptides. Here we describe the VCD spectroscopic characterization of the model peptide Boc-Val-Phe-nPr in chloroform as representative for a weakly interacting solvent and dimethyl sulfoxide (DMSO-d₆) as a strongly hydrogen bonding solvent. We show that the conformational preferences of the peptide in chloroform are well-described by the computationally predicted distribution of the isolated molecule assuming only implicit solvation effects through a continuum solvation model. In order to simulate the spectra recorded in DMSO-d₆, solvation was accounted for explicitly by computed microsolvated structures containing one to three solvent molecules. A good match of the computed spectra with the experimental data is obtained by this method. Comparing the conformational distributions in deuterated chloroform-d₁ and DMSO-d₆, structures with intramolecular hydrogen bonds such as the (δ,δ)-conformer family contribute to the conformational distribution only when there is no strong interaction with the solvent. This is in contrast to the results for the related Boc-Pro-Phe-nPr studied before, for which the intramolecular interaction was found to persist in DMSO-d₆. Furthermore, we discuss the influence of hydrogen bonding to different numbers of solvent molecules on the spectral signatures and show that the structure of the peptide in DMSO-d₆ is best described as a mixture of twofold-solvated (δ,β)- and threefold-solvated (β,β)-conformers.

Fourier Transform Infrared Spectroscopy and Vibrational Circular Dichroism Assisted Elucidation of the Solution-State Supramolecular Speciation in Racemic and Enantiopure Ketoprofen

The molecular structure and solution-state molecular interactions in the popular non-steroidal anti-inflammatory drug, ketoprofen, are extensively studied with the aim of gaining a better understanding of the chemical behavior of its solution state and its connection to its nucleation pathway and crystallization outcome. Using as reference solid-state X-ray structures of enantiomeric and racemic forms of ketoprofen, a set of self-assembly models underpinned by density functional theory calculations has been considered for the analysis of spectroscopic data, infrared (IR) and vibrational circular dichroism (VCD), obtained for solutions of the samples as a function of composition and solvent. From our results it can be concluded that, contrary to the general belief for generic carboxylic acids, there are no cyclic dimeric structures of ketoprofen present in solution, but rather linear arrays made up of two (in high polar or diluted media) or more units (in low polar or low dilution media). This observation is in line with the idea that the weak contacts (other than H-bonding) would hold the key to molecular self-assembly, in agreement with recent studies on other aromatic carboxylic acids.

VCD spectroscopy distinguishes the enamine and iminium ion of a 1,1’-binaphthyl azepine

We present a VCD spectroscopic characterization of a chiral 1,1’-binaphthyl azepine catalyst and show that the VCD spectra of an in-situ generated enamine and an ex-situ prepared iminium ion are...

Vibrational CD study on the solution phase structures of the MacMillan catalyst and its corresponding iminium ion

We demonstrate that VCD spectroscopy can reveal insights into the conformational preferences of the iminium ion obtained from MacMillan's imidazolidinone catalyst. For both the isolated and in situ generated iminium ion, the comparison of experimental and computed VCD spectra directly confirms that conformer 2b ("Houk-conformer") is the dominant structure in solution. This conclusion is reached without any in-depth interpretation of the spectroscopic data, just by visual comparison of the spectral signatures. For the parent catalyst 1 and its salts 1·HCl and 1·HClO₄, we report a comprehensive analysis of the conformational preferences in two solvents. VCD spectroscopy is subsequently shown to be able to reveal small conformational changes induced by solute-solvent and solute-anion interactions.

Double Hydrogen Bonding Dimerization Propensity of Aqueous Hydroxy Acids Investigated Using Vibrational Optical Activity

Lactic and malic acids are key substances in a number of biochemical processes in living cells and are also utilized in industry. Vibrational spectroscopy represents an efficient and sensitive way to study their structure and interactions. Since water is the natural environment, proper understanding of their molecular dynamics in aqueous solutions is of critical importance. To this end, we employed Raman spectroscopy and Raman optical activity (ROA) to study the conformation of l-lactic and l-malic acids in water (while varying pH, temperature, and concentration), with special emphasis on their double hydrogen bonding dimerization propensity. Raman and ROA experimental data were supported by extensive theoretical calculations of the vibrational properties and by additional experiments (IR absorption, vibrational circular dichroism, and NMR). Conformational behavior of the acids in water was described by molecular dynamics simulations. Reliability of the results was verified by calculating the vibrational properties of populated conformers and by comparing thus obtained spectral features with the experimental data. Calculations estimated the incidence of H-bonded dimers in water to be low in lactic acid and comparable to monomers in malic acid. The "hybrid" approach presented here reveals limitations of relying on the experimental spectra alone to study dimer formation.

Induced VCD and conformational chirality in host-guest complexes of a chiral ammonium salt with crown ethers

The hydrogen bonded complexes of the chiral ammonium salt α-methylbenzyl ammonium chloride (MBA-H⁺Cl^-) and the achiral crown ethers 18c6 and 15c5 serve as model systems to investigate the effect of host-guest complex formation on the conformational preferences of the macrocycles. We demonstrate that the intermolecular interactions result in new VCD signatures, that can be assigned to vibrational modes of the crown ethers. Based on a detailed conformational analysis, we investigate the origin of these signatures and discuss induced VCD (iVCD) and conformational chirality as possible sources of VCD intensity. The macrocycle in the MBA-H⁺/18c6 complex prefers either an achiral D_3d-symmetric conformation, which gives rise to iVCD, or chiral conformations, that feature individual contributions to the VCD spectrum. For the MBA-H⁺/15c5 complex, the contributions of the macrocycle to the VCD signatures are less pronounced and found to arise solely from conformational chirality. Therefore, analysis of the VCD signatures confirms that the small chiral guest molecule is able to affect the conformational preferences of a macrocyclic host. The study thus demonstrates the suitability of VCD spectroscopy for the characterization of analogous supramolecular host-guest complexes.

Correlating solvation with conformational pathways of proteins in alcohol-water mixtures: a THz spectroscopic insight

Water, being an active participant in most of the biophysical processes, is important to trace how protein solvation changes as its conformation evolves in the presence of solutes or co-solvents. In this study, we investigate how the secondary structures of two diverse proteins - lysozyme and β-lactoglobulin - change in the aqueous mixtures of two alcohols - ethanol and 2,2,2-trifluoroethanol (TFE) using circular dichroism measurements. We observe that these alcohols change the secondary structures of these proteins and the changes are protein-specific. Subsequently, we measure the collective solvation dynamics of these two proteins both in the absence and in the presence of alcohols by measuring the frequency-dependent absorption coefficient (α(ν)) in the THz (0.1-1.2 THz) frequency domain. The alcohol-water mixtures exhibit a non-ideal behaviour with the highest absorption difference (Δα) obtained at X_alcohol = 0.2. The protein solvation in the presence of the alcohols shows an oscillating behaviour in which Δα_protein changes with X_alcohol. Such an oscillatory behaviour of protein solvation results from a delicate interplay between the protein-water, protein-alcohol and water-alcohol associations. We attempt to correlate the various structural conformations of the proteins with the associated solvation.

How important are the intermolecular hydrogen bonding interactions in methanol solvent for interpreting the chiroptical properties?

Two crispine A analogs and tetrahydrofuro[2,3-b]furan-3,3a(6aH)-diol, endowed with hydroxyl groups that can participate in intramolecular hydrogen bonding, have been synthesized and experimental vibrational circular dichroism (VCD) spectra and optical rotatory dispersion (ORD) data have been measured in CD₃OD/CH₃OH solvents. The absolute configurations (ACs) of these compounds have been determined using their synthetic schemes, supplemented wherever possible with X-ray diffraction data. The ACs are also analyzed with quantum chemical (QC) calculations of VCD and ORD utilizing implicit solvation as well as explicit solvation models, with the later employing classical molecular dynamics (MD) simulations. It is found that VCD calculations with implicit solvation model are adequate for determining the ACs, despite propensity of studied compounds for intermolecular hydrogen bonding between solute and solvent molecules. This observation is important because time-consuming MD simulations may not be necessary in the type of situations studied here. Additionally, it is found that the QC predicted VCD spectra provided enough diastereomer discrimination for determining the correct AC of studied compounds independently. The same observation did not apply to ORD.

Solvation and the secondary structure of a proline-containing dipeptide: insights from VCD spectroscopy

In this study we investigate the IR and VCD spectra of the diastereomeric dipeptide Boc-Pro-Phe-(n-propyl) 1 in chloroform-d₁ (CDCl₃) and the strongly hydrogen bonding solvent dimethylsulfoxide-d₆ (DMSO-d₆). From comparison of the experimental spectra, the amide II spectral region is identified as marker signature for the stereochemistry of the dipeptide: the homochiral LL-1 features a (+/-)-pattern in the amide II region of the VCD spectrum, while the amide II signature of the diastereomer LD-1 is inverted. Computational analysis of the IR and VCD spectra of LL-1 reveals that the experimentally observed amide II signature is characteristic for a β_I-turn structure of the peptide. Likewise, the inverted pattern found for LD-1 arises from a β_II-turn structure of the dipeptide. Following a micro-solvation approach, the experimental spectra recorded in DMSO-d₆ are computationally well reproduced by considering only a single solvent molecule in a hydrogen bond with N-H groups. Considering a second solvent molecule, which would lead to a cleavage of intramolecular hydrogen bonds in 1, is found to give a significantly worse match with the experiment. Hence, the detailed computational analysis of the spectra of LL- and LD-1 recorded in DMSO-d₆ confirms that the intramolecular hydrogen bonding pattern, that stabilizes the β-turns and other conformations of LL- and LD-1 in apolar solvents, remains intact. Our findings also show that it is essential to consider solvation explicitly in the analysis of the IR and VCD spectra of dipeptides in strongly hydrogen bonding solvents. As the solute-solvent interactions affect both conformational preferences and spectral signatures, it is also demonstrated that this inclusion of solvent molecules cannot be circumvented by applying fitting procedures to non-solvated structures.

Solvent Effects and Aggregation Phenomena Studied by Vibrational Optical Activity and Molecular Dynamics: The Case of Pantolactone

Raman and Raman optical activity (ROA), IR, and vibrational circular dichroism (VCD) spectra of (R)- and (S)-pantolactone have been recorded in three solvents. ROA has been employed on water and DMSO solutions, VCD on DMSO and CCl₄ solutions. In the last solvent, monomer-dimer equilibrium is present. Due to the low conformational flexibility of the isolated molecule and to the possibility of aggregation, this compound has been used here to test different protocols for computation of the spectroscopic responses taking into account solvent effects. Molecular dynamics (MD) simulations have been carried out together with statistical clustering methods based on collective variables to extract the structures needed to calculate the spectra. Quantum mechanical DFT calculations based on PCM are compared with approaches based on different representations of the solvent shell (MM or QM level). Appropriate treatment of the solvent permits obtaining of good band-shapes, with the added advantage that the MD analysis allows one to take into account flexibility of dimeric structures justifying the broadness of observed bands and the absence of intense VCD couplets in the carbonyl and OH stretching regions.

More complex, less complicated? Explicit solvation of hydroxyl groups for the analysis of VCD spectra

With increasing size of the molecules, hydrogen bonding induced solvent effects on the IR and VCD spectra become more negligible.

A combined Raman optical activity and vibrational circular dichroism study on artemisinin-type products

Artemisinin and two of its derivatives, dihydroartemisinin and artesunate, front line drugs against malaria, were studied using Raman optical activity (ROA), vibrational circular dichroism (VCD) experiments and density functional theory (DFT) calculations.

How many solvent molecules are required to solvate chiral 1,2-diols with hydrogen bonding solvents? A VCD spectroscopic study

A VCD spectroscopic analysis of selected model systems for solute–solvent interactions of chiral diols with hydrogen bonding solvents DMSO and ACN.

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.