The important role of non-covalent interactions for the vibrational circular dichroism of lactic acid in aqueous solution

Vibrational circular dichroism spectra of proline in water at different pH values

Recording VCD spectra of aqueous solution poses a particular challenge as water is a strong infrared absorber. Likewise, the computational analysis of VCD spectra by means of DFT-based spectral calculations requires the consideration of explicit solvent molecules, thus posing an even greater challenge. Several studies suggested that by modeling the solvent environment with a few water molecules in a micro-solvation approach would be sufficient to describe experimental spectra. For example, using proline at different pH values, we herein show that a change in the relative spatial orientation of a single water molecule in five-fold solvated structures strongly affects the computed VCD spectral signatures and that Boltzmann-weighted spectra do not correctly reproduce the experiment. We thus explored an approach based on molecular dynamics and subsequent DFT-calculations, in which we considered 30 water molecules (about 1.5 solvation shells). Once again, it was found that the Boltzmann-weighted spectra obtained on the basis of several hundred structures did not correctly reproduce experimental signatures, and a simple averaging scheme resulted in well-matching spectra with comparable bandwidths. The rationale behind the procedure was that sampling the configurational space of the solvent molecules is as equally important as the conformational sampling of the solute. For conformationally more flexible molecules, it is assumed that a much larger set of structures will have to be computed in order to properly sample the conformational space of both solute and solvent.

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.

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.

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.

Infrared and vibrational circular dichroism spectra of methyl β-D-glucopyranose in water: The application of the quantum cluster growth and clusters-in-a-liquid solvation models

The infrared (IR) and vibrational circular dichroism (VCD) spectra of methyl β-D-glucopyranose in water were measured. Both implicit and explicit solvation models were utilized to explain the observed spectra. The vast body of existing experimental and theoretical data suggested that about eight explicit water molecules are needed to account for the solvent effects, supported by the current Quantum Cluster Growth (QCG) analysis. Extensive manual and systematic conformational searches of the molecular target and its water clusters were carried out by using a recently developed conformational searching tool, conformer-rotamer ensemble sampling tool (CREST), and the microsolvation model in the associated QCG code. The Boltzmann averaged IR and VCD spectra of the methyl β-D-glucopyranose-(water)n (n = 8) conformers in the PCM of water provide better agreement with the experimental ones than those with n = 0, 1, and 2. The explicit solvation with eight water molecules was shown to greatly modify the conformational preference of methyl β-D-glucopyranose from its monomeric form. Further analyses show that the result is consistent with the existence of long-lived methyl β-D-glucopyranose monohydrates with the additional explicit water effects being accounted for with the quantum mechanical treatment of the other seven close-by water molecules in the PCM of water.

Isolated and solvated chiroptical behavior in conformationally flexible butanamines

The nonresonant optical activity of two highly flexible aliphatic amines, (2R)-3-methyl-2-butanamine (R-MBA) and (2R)-(3,3)-dimethyl-2-butanamine (R-DMBA), has been probed under isolated and solvated conditions to examine the roles of conformational isomerism and to explore the influence of extrinsic perturbations. The optical rotatory dispersion (ORD) measured in six solvents presented uniformly negative rotatory powers over the 320-590 nm region, with the long-wavelength magnitude of chiroptical response growing nearly monotonically as the dielectric constant of the surroundings diminished. The intrinsic specific optical rotation,α λ T (in deg dm-1 [g/mL]-1 ), extracted for ambient vapor-phase samples of R-MBA [-11.031(98) and -2.29 (11)] and R-DMBA [-9.434 (72) and -1.350 (48)] at 355 and 633 nm were best reproduced by counterintuitive solvents of high polarity (yet low polarizability) like acetonitrile and methanol. Attempts to interpret observed spectral signatures quantitatively relied on the linear-response frameworks of density-functional theory (B3LYP, cam-B3LYP, and dispersion-corrected analogs) and coupled-cluster theory (CCSD), with variants of the polarizable continuum model (PCM) deployed to account for the effects of implicit solvation. Building on the identification of several low-lying equilibrium geometries (nine for R-MBA and three for R-DMBA), ensemble-averaged ORD profiles were calculated at T = 300 K by means of the independent-conformer ansatz, which enabled response properties predicted for the optimized structure of each isomer to be combined through Boltzmann-weighted population fractions derived from corresponding relative internal-energy or free-energy values, the latter of which stemmed from composite CBS-APNO and G4 analyses. Although reasonable accord between theory and experiment was realized for the isolated (vapor-phase) species, the solution-phase results were less satisfactory and tended to degrade progressively as the solvent polarity increased. These trends were attributed to solvent-mediated changes in structural parameters and energy metrics for the transition states that separate and putatively isolate the equilibrium conformations supported by the ground electronic potential-energy surface, with the resulting displacement of barrier locations and/or decrease of barrier heights compromising the underlying premise of the independent-conformer ansatz.

Conformational Distributions of Phenyl β-D-Glucopyranoside and Gastrodin in Solution by Vibrational Optical Activity and Theoretical Calculations

The conformational landscapes of two highly flexible monosaccharide derivatives, namely phenyl β-D-glucopyranoside (ph-β-glu) and 4-(hydroxymethyl)phenyl β-D-glucopyranoside, also commonly known as gastrodin, were explored using a combined experimental and theoretical approach. For the infrared, Raman, and the associated vibrational optical activity (VOA), i.e., vibrational circular dichroism and Raman optical activity, experiments of these two compounds in DMSO and in water were carried out. Extensive and systematic conformational searches were performed using a recently developed conformational searching tool called CREST (conformer-rotamer ensemble sampling tool) in the two solvents. Fourteen and twenty-four low-energy conformers were identified at the DFT level for ph-β-glu and gastrodin, respectively. The spectral simulations of individual conformers were done at the B3LYP-D3BJ/def2-TZVPD level with the polarizable continuum model of the solvents. The VOA spectral features exhibit much higher specificity to conformational differences than their parent infrared and Raman. The excellent agreements achieved between the experimental and simulated VOA spectra allow for the extraction of experimental conformational distributions of these two carbohydrates in solution directly. The experimental percentage abundances based on the hydroxymethyl (at the pyranose ring) conformations G+, G-, and T for ph-β-glu were obtained to be 15%, 75%, and 10% in DMSO and 53%, 40%, and 7% in water, respectively, in comparison to the previously reported gas phase values of 68%, 25%, and 7%, highlighting the important role of solvents in conformational preferences. The corresponding experimental distributions for gastrodin are 56%, 22%, and 22% in DMSO and 70%, 21%, and 9% in water.

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.

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.

The solid state VCD of a novel N-acylhydrazone trifluoroacetate

A novel N-acylhydrazone with pharmaceutical importance was subject of structural and IR/VCD investigations in the solid state. In the crystal structure, dimers of anion-cation pairs are stabilized by H-bonding and ionic interactions. Some less common interaction types, like C=N···C-NH₃⁺ (σ-hole) interactions, hydrazone-aromatic interactions and dispersive contacts of the CF₃ groups are also present in the crystal. Satisfactory reproduction of the solid state IR and VCD spectra required that quantum-chemical calculations be done on a tetramer (four cation-anion pairs) cut out from the crystal structure, exhibiting key intermolecular interactions. Ten DFT functionals were assessed as to the agreement between the calculated and experimental spectra. Various approaches to scaling of the calculated frequencies were applied. The best results were yielded with individual (optimized) frequency scaling factors (FSFs) and band half-widths at half maximum-(HWHM) for four separate spectral subregions. The best matching between the experimental and theoretical spectra (according to SimIR, SimVCD and SimVDF indices) was found for the B3PW91 functional, however, a few other functionals follow closely in the ranking. Based on the quantum chemical calculations, spectral assignments have been made.