An analytical derivative procedure for the calculation of vibrational Raman optical activity spectra

Understanding the surrounding effects on Raman optical activity signatures of a chiral cage system: Cryptophane-PP-111

Cryptophane molecules are cage-like structures consisting in two hemispheres, each made of three benzene rings. These hemispheres are bound together with three O(CH2)nOlinkers of various lengths giving rise to a plethora of cryptophane derivatives. Moreover, they are able to encapsulate neutral guests: CH2Cl2, CHCl3, …; and charged species: Cs+, Tl+, …. Finally, they exhibit chiroptical properties thanks to the anti arrangement of the linkers between the hemispheres. This work focuses on the Raman optical activity (ROA) signatures of Cryptophane-111 (n=1 for each linker). More specifically, we aim at simulating accurately its ROA spectra with and without a xenon atom inside its cavity. Experimental data (Buffeteau et al., 2017) have already demonstrated the effect of the encapsulation in the low-wavenumbers region. To generate the initial structures, we rely on the novel Conformer-Rotamer Ensemble Sampling Tool (CREST) program, developed by S. Grimme and co-workers. This is required due to the flexibility provided by the linkers. The CREST algorithm seems promising and has already been used to sample the potential energy surface (PES) of target systems before the simulation of their vibrational spectroscopies (Eikås et al., 2022). We observe large similarities between the two sets of conformers (one with and one without Xe encapsulated), demonstrating the robustness of the CREST algorithm. For corresponding structures, the presence of xenon pushed the two hemispheres slightly further apart. After optimization at the DFT level, only one unique conformer has a Boltzmann population ratio greater than 1%, pointing out the relative rigidity of the cage. Based on this unique conformer, our simulations are in good agreement with the experimental data. Regarding xenon encapsulation, the (experimental and theoretical) ROA signatures at low wavenumbers are impacted: slight shifts in wavenumbers are observed as well as a decrease in relative ROA intensity for bands around 150 cm-1. The wavenumber shifts were very well reproduced by our simulations, but the experimental decrease in the ROA intensity was unfortunately not reproduced.

Intense chiral signal from α-helical poly-L-alanine observed in low-frequency Raman optical activity

Raman optical activity (ROA) spectral features reliably indicate the structure of peptides and proteins, but the signal is often weak. However, we observed significantly enhanced low-frequency bands for α-helical poly-L-alanine (PLA) in solution. The biggest ROA signal at ∼100 cm^-1 is about 10 times stronger than higher-frequency bands described previously, which facilitates the detection. The low-frequency bands of PLA were compared to those of α-helical proteins. For PLA, density functional simulations well reproduced the experimental spectra and revealed that about 12 alanine residues within two turns of the α-helix generate the strong ROA band. Averaging based on molecular dynamics (MD) provided an even more realistic spectrum compared to the static model. The low-frequency bands could be largely related to a collective motion of the α-helical backbone, partially modulated by the solvent. Helical and intermolecular vibrational coordinates have been introduced and the helical unwinding modes were assigned to the strongest ROA signal at 101-128 cm^-1. Further analysis indicated that the helically arranged amide and methyl groups are important for the strong chiral signal of PLA, while the local chiral centers C_αH contribute in a minor way only. The strong low-frequency ROA can thus provide precious information about the motions of the peptide backbone and facilitate future protein studies.

Measurement and Theory of Resonance Raman Optical Activity for Gases, Liquids, and Aggregates. What It Tells about Molecules

Resonance Raman Optical Activity (RROA) appeared as a natural extension of the nonresonance branch. It combines the structural sensitivity of chiroptical spectroscopy with the signal enhancement coming from the resonance of molecular electronic transitions with the excitation laser light. However, the idea has been hampered by many technical and theoretical problems that are being clarified only in recent years. We provide the theoretical basis and several examples documenting the problems, achievements, and potential of RROA, in particular in biomolecular studies.

Theory and algorithms for chiroptical properties and spectroscopies of aqueous systems

Chiroptical properties and spectroscopies are valuable tools to study chiral molecules and assign absolute configurations. The spectra that result from chiroptical measurements may be very rich and complex, and hide much of their information content. For this reason, the interplay between experiments and calculations is especially useful, provided that all relevant physico-chemical interactions that are present in the experimental sample are accurately modelled. The inherent difficulty associated to the calculation of chiral signals of systems in aqueous solutions requires the development of specific tools, able to account for the peculiarities of water-solute interactions, and especially its ability to form hydrogen bonds. In this perspective we discuss a multiscale approach, which we have developed and challenged to model the most used chiroptical techniques.

The effect of protein backbone hydration on the amide vibrations in Raman and Raman optical activity spectra

Raman and specifically Raman optical activity (ROA) spectroscopy are very sensitive to the solution structure and conformation of biomolecules. Because of this strong conformational sensitivity, density functional theory (DFT) calculations are often used to get a better understanding of the experimentally observed spectral patterns. While e.g. for carbohydrate structure the water molecules that surround the solute have been demonstrated to be of vital importance to get accurate modelled ROA spectra, the effect of explicit water molecules on the calculated ROA patterns of peptides and proteins is less well studied. Here, the effect of protein backbone hydration was studied using DFT calculations of HCO-(l-Ala)5-NH2 in specific secondary structure conformations with different treatments of the solvation. The effect of the explicit water molecules on the calculated spectra mainly arises from the formation of hydrogen bonds with the amide C[double bond, length as m-dash]O and N-H groups. Hydrogen bonding of water with the C[double bond, length as m-dash]O group determines the shape and position of the amide I band. The C[double bond, length as m-dash]O bond length increases upon formation of C[double bond, length as m-dash]OH2O hydrogen bonds. The effect of the explicit water molecules on the amide III vibrations arises from hydrogen bonding of the solvent with both the C[double bond, length as m-dash]O and N-H group, but their contributions to this spectral region differ: geometrically, the formation of a C[double bond, length as m-dash]OH2O bond decreases the C-N bond length, while upon forming a N-HH2O hydrogen bond, the N-H bond length increases.

Localized molecular orbitals for calculation and analysis of vibrational Raman optical activity

First calculations of vibrational Raman optical activity based on localized molecular orbitals are presented, which pave the way for novel insight into spectroscopic signatures of chiral systems.

Ramachandran mapping of peptide conformation using a large database of computed Raman and Raman optical activity spectra

A novel ROA database is reported that assigns peptide structures in detail by pattern recognition of the experimental spectrum.

Calculation of Raman optical activity spectra for vibrational analysis

This review provides the necessary knowledge to accurately model ROA spectra of solvated systems and interpret their vibrational characteristics.

Where does the Raman optical activity of [Rh(en)3]3+ come from? Insight from a combined experimental and theoretical approach

Raman optical activity (ROA) spectra are measured and calculated for Δ- and Λ-tris-(ethylenediamine)rhodium(iii) chloride in aqueous solution.

Three types of induced tryptophan optical activity compared in model dipeptides: theory and experiment

The tryptophan (Trp) aromatic residue in chiral matrices often exhibits a large optical activity and thus provides valuable structural information. However, it can also obscure spectral contributions from other peptide parts. To better understand the induced chirality, electronic circular dichroism (ECD), vibrational circular dichroism (VCD), and Raman optical activity (ROA) spectra of Trp-containing cyclic dipeptides c-(Trp-X) (where X = Gly, Ala, Trp, Leu, nLeu, and Pro) are analyzed on the basis of experimental spectra and density functional theory (DFT) computations. The results provide valuable insight into the molecular conformational and spectroscopic behavior of Trp. Whereas the ECD is dominated by Trp π-π* transitions, VCD is dominated by the amide modes, well separated from minor Trp contributions. The ROA signal is the most complex. However, an ROA marker band at 1554 cm(-1) indicates the local χ(2) angle value in this residue, in accordance with previous theoretical predictions. The spectra and computations also indicate that the peptide ring is nonplanar, with a shallow potential so that the nonplanarity is primarily induced by the side chains. Dispersion-corrected DFT calculations provide better results than plain DFT, but comparison with experiment suggests that they overestimate the stability of the folded conformers. Molecular dynamics simulations and NMR results also confirm a limited accuracy of the dispersion-DFT model in nonaqueous solvents. Combination of chiral spectroscopies with theoretical analysis thus significantly enhances the information that can be obtained from the induced chirality of the Trp aromatic residue.

Molecular structure determination using chiroptical spectroscopy: where we may go wrong?

Chiroptical spectroscopy is being widely used for determining the three-dimensional molecular structures (i.e., absolute configurations and conformations) of chiral molecules. The general procedure used with any of the chiroptical spectroscopic methods is to analyze the experimental data using corresponding quantum chemical predictions. Such analysis involves multiple steps, including consideration of conformations, solvent effects, electronic transitions, stereoisomers, and experimental artifacts, each of which possesses certain limitations. These limitations, when not recognized or properly taken into account, may lead to incorrect conclusions. This review emphasizes on selected examples that illustrate the potential limitations in utilizing the chiroptical spectroscopic methods. The examples used include hibiscus acid dimethylester, hibiscus acid disodium salt, 3,3'-diphenyl-[2,2'-binaphthalene]-1,1'-diol, tartaric acid esters, and 6,6'-dibromo-[1,1'-binaphthalene]-2,2'-diol.

Towards modelling the vibrational signatures of functionalized surfaces: carboxylic acids on H-Si(111) surfaces

In this work, we investigate the adsorption process of two carboxylic acids (stearic and undecylenic) on a H-Si(111) surface via the calculation of structural and energy changes as well as the simulation of their IR and Raman spectra. The two molecules adsorb differently at the surface since the stearic acid simply physisorbs while the undecylenic acid undergoes a chemical reaction with the hydrogen atoms of the surface. This difference is observed in the change of geometry during the adsorption. Indeed, the chemisorption of the undecylenic acid has a bigger impact on the structure than the physisorption of the stearic acid. Consistently, the former is also characterized by a larger value of adsorption energy and a smaller value of the tilting angle with respect to the normal plane. For both the IR and Raman signatures, the spectra of both molecules adsorbed at the surface are in a first approximation the superposition of the spectra of the Si cluster and of the carboxylic acid considered individually. The main deviation from this simple observation is the peak of the stretching Si-H (ν(Si-H)) mode, which is split into two peaks upon adsorption. As expected, the splitting is bigger for the chemisorption than the physisorption. The modes corresponding to atomic displacements close to the adsorption site display a frequency upshift by a dozen wavenumbers. One can also see the disappearance of the peaks associated with the C=C double bond when the undecylenic acid chemisorbs at the surface. The Raman and IR spectra are complementary and one can observe here that the most active Raman modes are generally IR inactive. Two exceptions to this are the two ν(Si-H) modes which are active in both spectroscopies. Finally, we compare our simulated spectra with some experimental measurements and we find an overall good agreement.

Theoretical investigation of Raman optical activity signatures of Tröger's base

The Raman and VROA spectra of (S,S)-Tröger's base are simulated. We mainly discuss the peaks in the 1140-1400 cm(-1) wavenumber range where an intense VROA signature is found. In this range, nearly all of the Raman-active bands belong to the irreducible representation A (C(2) point group), whereas no such observation is made for the VROA spectrum. The vibrational normal modes associated with the peaks in this range mainly consist of wagging and twisting motions of the hydrogen atoms. From the atomic contribution patterns (ACPs) and the group coupling matrices (GCMs), one finds that the VROA backward-scattering intensities mainly arise from hydrogen and carbon atoms in the vicinity of the two chiral nitrogen atoms. The VROA signatures in the 1140-1400 cm(-1) range are therefore a fingerprint of the local chirality around the two chiral nitrogen centers.

Theoretical approaches to the calculation of Raman optical activity spectra

In this article, we will give a brief account of the different approaches that have been presented in the literature for calculating Raman optical activity (ROA) spectra by ab initio methods. We will also outline the general structure of a self-consistent-field-based approach for analytic calculations of ROA spectra, including also contributions from London orbitals. The use of London orbitals ensures that the relevant ROA parameters are gauge origin independent. We will also give an outlook on the future of ab initio calculations of Raman optical activity spectra.

Vibrational Raman optical activity of 1-phenylethanol and 1-phenylethylamine: revisiting old friends

The samples used for the first observations of vibrational Raman optical activity (ROA) in 1972, namely both enantiomers of 1-phenylethanol and 1-phenylethylamine, have been revisited using a modern commercial ROA instrument together with state-of-the-art ab initio calculations. The simulated ROA spectra reveal for the first time the vibrational origins of the first reported ROA signals, which comprised similar couplets in the alcohol and amine in the spectral range approximately 280-400 cm(-1). The results demonstrate how easy and routine ROA measurements have become, and how current ab initio quantum-chemical calculations are capable of simulating experimental ROA spectra quite closely provided sufficient averaging over accessible conformations is included. Assignment of absolute configuration is, inter alia, completely secure from results of this quality. Anharmonic corrections provided small improvements in the simulated Raman and ROA spectra. The importance of conformational averaging emphasized by this and previous related work provides the underlying theoretical background to ROA studies of dynamic aspects of chiral molecular and biomolecular structure and behavior.