Theoretical approaches to the calculation of Raman optical activity spectra

Conformations and hydration of halopropionic acids studied by molecular dynamics and Raman optical activity

Chiral 2-halopropionic acids and their derivatives were synthesized and their properties studied computationally using Raman and Raman optical activity (ROA) spectroscopy. For neat acids present as liquids small amount of water led to significant changes in the spectra, resulting even to flipping of some ROA band signs. We find this interesting for the role water plays in interpretation of vibrational optical activity spectra of biomolecules. Analysis of the results shows that when the water is present, it can change ROA band signs due to the changes in acidobasic equilibrium. Corresponding esters without acidic hydrogens do not exhibit such effects.

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.

Photochemical synthesis of pink silver and its use for monitoring peptide nitration via surface enhanced Raman spectroscopy (SERS)

Oxidative stress may cause extended tyrosine posttranslational modifications of peptides and proteins. The 3-nitro-L-tyrosine (Nit), which is typically formed, affects protein behavior during neurodegenerative processes, such as Alzheimer's and Parkinson's diseases. Such metabolic products may be conveniently detected at very low concentrations by surface enhanced Raman spectroscopy (SERS). Previously, we have explored the SERS detection of the Nit NO₂ bending vibrational bands in a presence of hydrogen chloride (Niederhafner et al., Amino Acids 53:517-532, 2021, ibid). In this article, we describe performance of a new SERS substrate, "pink silver", synthesized photochemically. It provides SERS even without the HCl induction, and the acid further decreases the detection limit about 9 times. Strong SERS bands were observed in the asymmetric (1550-1475 cm^-1) and symmetric (1360-1290 cm^-1) NO stretching in the NO₂ group. The bending vibration was relatively weak, but appeared stronger when HCl was added. The band assignments were supported by density functional theory modeling.

Low-frequency Raman Optical Activity Provides Insight into Structure of Chiral Liquids

Vibrational frequencies of modes involving intermolecular motions in liquids are relatively small, in the Raman scattering close to the excitation frequency, and the bands may merge into a diverging uninterpretable...

Inspecting the structural characteristics of chiral drug penicillamine under different pH conditions using Raman optical activity spectroscopy and DFT calculations

The investigation of the structural characteristics of chiral drugs in physiological environments is a challenging research topic, which may lead to a better understanding of how the drugs work. Raman optical activity (ROA) spectroscopy in combination with density functional theory (DFT) calculations was exploited to inspect the structural changes in penicillamine under different acid-base states in aqueous solutions. The B3LYP/aug-cc-PVDZ method was employed and the implicit solvation model density (SMD) was considered for describing the solvation effect in H₂O. The conformations of penicillamine varied with pH, but penicillamine was liable to stabilize in the form of the P_C conformation (the sulfur atom is in a trans orientation with respect to carboxylate) in most cases for both D- and L-isomers. The relationship between the conformations of penicillamine and the ROA peaks, as well as peak assignments, were comprehensively studied and elucidated. In the fingerprint region, two ROA couplets and one ROA triplet with different patterns were recognized. The intensity, sign and frequency of the corresponding peaks also changed with varying pH. Deuteration was carried out to identify the vibrational modes, and the ROA peaks of the deuterated amino group in particular are sensitive to change in the ambient environment. The results are expected not only to serve as a reference for the interpretation of the ROA spectra of penicillamine and other chiral drugs with analogous structures but also to evaluate the structural changes of chiral molecules in physiological environments, which will form the basis of further exploration of the effects of structural characteristics on the pharmacological and toxicological properties of chiral drugs.

Advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase

The purpose of this review is to demonstrate advances, challenges and perspectives of quantum chemical approaches in molecular spectroscopy of the condensed phase. Molecular spectroscopy, particularly vibrational spectroscopy and electronic spectroscopy, has been used extensively for a wide range of areas of chemical sciences and materials science as well as nano- and biosciences because it provides valuable information about structure, functions, and reactions of molecules. In the meantime, quantum chemical approaches play crucial roles in the spectral analysis. They also yield important knowledge about molecular and electronic structures as well as electronic transitions. The combination of spectroscopic approaches and quantum chemical calculations is a powerful tool for science, in general. Thus, our article, which treats various spectroscopy and quantum chemical approaches, should have strong implications in the wider scientific community. This review covers a wide area of molecular spectroscopy from far-ultraviolet (FUV, 120-200 nm) to far-infrared (FIR, 400-10 cm-1)/terahertz and Raman spectroscopy. As quantum chemical approaches, we introduce several anharmonic approaches such as vibrational self-consistent field (VSCF) and the combination of periodic harmonic calculations with anharmonic corrections based on finite models, grid-based techniques like the Numerov approach, the Cartesian coordinate tensor transfer (CCT) method, Symmetry-Adapted Cluster Configuration-Interaction (SAC-CI), and the ZINDO (Semi-empirical calculations at Zerner's Intermediate Neglect of Differential Overlap). One can use anharmonic approaches and grid-based approaches for both infrared (IR) and near-infrared (NIR) spectroscopy, while CCT methods are employed for Raman, Raman optical activity (ROA), FIR/terahertz and low-frequency Raman spectroscopy. Therefore, this review overviews cross relations between molecular spectroscopy and quantum chemical approaches, and provides various kinds of close-reality advanced spectral simulation for condensed phases.

α-Synuclein conformations followed by vibrational optical activity. Simulation and understanding of the spectra

For α-synuclein novel structural markers were identified in vibrational optical activity spectra and supported by theoretical modeling.

Monitoring peptide tyrosine nitration by spectroscopic methods

Oxidative stress can lead to various derivatives of the tyrosine residue in peptides and proteins. A typical product is 3-nitro-L-tyrosine residue (Nit), which can affect protein behavior during neurodegenerative processes, such as those associated with Alzheimer's and Parkinson's diseases. Surface enhanced Raman spectroscopy (SERS) is a technique with potential for detecting peptides and their metabolic products at very low concentrations. To explore the applicability to Nit, we use SERS to monitor tyrosine nitration in Met-Enkephalin, rev-Prion protein, and α-synuclein models. Useful nitration indicators were the intensity ratio of two tyrosine marker bands at 825 and 870 cm^-1 and a bending vibration of the nitro group. During the SERS measurement, a conversion of nitrotyrosine to azobenzene containing peptides was observed. The interpretation of the spectra has been based on density functional theory (DFT) simulations. The CAM-B3LYP and ωB97XD functionals were found to be most suitable for modeling the measured data. The secondary structure of the α-synuclein models was monitored by electronic and vibrational circular dichroism (ECD and VCD) spectroscopies and modeled by molecular dynamics (MD) simulations. The results suggest that the nitration in these peptides has a limited effect on the secondary structure, but may trigger their aggregation.

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.

Vibrational optical activity for structural characterization of natural products

This review presents the recent progress towards elucidating the structures of chiral natural products and applications using vibrational optical activity (VOA) spectroscopy.

Stereostructural Elucidation of Glucose Phosphorylation by Raman Optical Activity

Phosphorylation of glucose is the prime step in sugar metabolism and energy storage. Two key glucose phosphates are involved, that is, glucose 6-phosphate (G6P) and α-glucose 1-phosphate (αG1P). The chiral conformation of glucose, G6P, and αG1P plays an essential role in enzyme-mediated conversions. However, few techniques were able to give a direct view of the conformational changes from glucose to G6P and αG1P. Here, Raman optical activity (ROA) was used to elucidate the stereochemical evolution of glucose upon phosphorylation. ROA was found to be extremely sensitive to different phosphorylation sites. A characteristic ROA marker of (+)980 cm^-1, originated from the phosphate group symmetric stretching vibration, is observed for αG1P with phosphorylation at chiral C1, while no corresponding ROA signal for G6P with phosphorylation at achiral C6 is observed. Phosphorylation-induced gauch-gauch (gg)/gauch-trans (gt) rotamer distribution changes can be sensitively probed by the sign of the ROA band around 1460 cm^-1. A positive ROA band at 1465 cm^-1 of glucose corresponds to a higher gt ratio, while a negative band at 1455 cm^-1 of G6P suggests a dominant gg population, and the disappearance of this ROA band for αG1P indicates a nearly balanced gg/gt distribution. Meanwhile, the phosphorylation at C6 and C1 could cause dramatic reduction of the conformational flexibility of the adjacent C4-OH and C2-OH, respectively. These stereochemical changes revealed by ROA spectra offer a structural basis on the understanding of sugar phosphorylation from the perspective of chirality.

Solution Structure of Mannobioses Unravelled by Means of Raman Optical Activity

Structural analysis of carbohydrates is a complicated endeavour, due to the complexity and diversity of the samples at hand. Herein, we apply a combined computational and experimental approach, employing molecular dynamics (MD) and density functional theory (DFT) calculations together with NMR and Raman optical activity (ROA) measurements, in the structural study of three mannobiose disaccharides, consisting of two mannoses with varying glycosidic linkages. The disaccharide structures make up the scaffold of high mannose glycans and are therefore important targets for structural analysis. Based on the MD population analysis and NMR, the major conformers of each mannobiose were identified and used as input for DFT analysis. By systematically varying the solvent models used to describe water interacting with the molecules and applying overlap integral analysis to the resulting calculational ROA spectra, we found that a full quantum mechanical/molecular mechanical approach is required for an optimal calculation of the ROA parameters. Subsequent normal mode analysis of the predicted vibrational modes was attempted in order to identify possible marker bands for glycosidic linkages. However, the normal mode vibrations of the mannobioses are completely delocalised, presumably due to conformational flexibility in these compounds, rendering the identification of isolated marker bands unfeasible.

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.

Establishing the link between fibril formation and Raman optical activity spectra of insulin

Molecular dynamics and density functional simulations are used to explain changes in Raman optical activity accompanying the formation of insulin fibrils.

Dependence of Adenine Raman Spectrum on Excitation Laser Wavelength: Comparison between Experiment and Theoretical Simulations

We acquired the Raman spectra of adenine in powder and aqueous phase using excitation lasers with 532, 633, and 785 nm wavelengths for the region between 300 and 1500 cm^-1. In comparison to the most distinct peak at 722 cm^-1, the peaks between 1200 and 1500 cm^-1 exhibited a characteristic increase in cross-section with decreasing excitation wavelength in both phases. This trend can be reproduced by different density functional theory (DFT) calculations for the adenine molecule in the gas phase as well as in the aqueous phase. Furthermore, from the calculation on the π-stacked dimer, hydrogen-bonded dimer, and trimer, we find that this trend toward excitation laser wavelength is not sensitive to the packing. When comparing the Raman spectra given by different excitation wavelength, one should take care in analyzing the cross-section, and present day DFT calculations are able to capture general trends in the excitation laser wavelength dependence of the Raman activity.

The computational prediction of Raman and ROA spectra of charged histidine tautomers in aqueous solution

Histidine is a key component of a number of enzymatic mechanisms, and undertakes many functionalities in biochemical systems.

Anharmonic Effects on Vibrational Spectra Intensities: Infrared, Raman, Vibrational Circular Dichroism, and Raman Optical Activity

The aim of this paper is 2-fold. First, we want to report the extension of our virtual multifrequency spectrometer (VMS) to anharmonic intensities for Raman optical activity (ROA) with the full inclusion of first- and second-order resonances for both frequencies and intensities in the framework of the generalized second-order vibrational perturbation theory (GVPT2) for all kinds of vibrational spectroscopies. Then, from a more general point of view, we want to present and validate the performance of VMS for the parallel analysis of different vibrational spectra for medium-sized molecules (IR, Raman, VCD, ROA) including both mechanical and electric/magnetic anharmonicity. For the well-known methyloxirane benchmark, careful selection of density functional, basis set, and resonance thresholds permitted us to reach qualitative and quantitative agreement between experimental and computed band positions and shapes. Next, the whole series of halogenated azetidinones is analyzed, showing that it is now possible to interpret different spectra in terms of mass, electronegativity, polarizability, and hindrance variation between closely related substituents, chiral spectroscopies being particular effective in this connection.

Conformational dynamics of carbohydrates: Raman optical activity of D-glucuronic acid and N-acetyl-D-glucosamine using a combined molecular dynamics and quantum chemical approach

As two biologically and medically relevant monosaccharides, the constituents of hyaluronic acid, d-glucuronic acid and N-acetyl-d-glucosamine, constitute perfect test cases for the development of carbohydrate-specific structural methods. These two molecules have been analysed by Raman optical activity (ROA), a spectroscopic technique exhibiting exquisite sensitivity to stereochemistry. We show that it is possible to support the experiment with a simulation approach combining density functional theory (DFT) and molecular dynamics (MD), both using explicit solvation. Thus, we have gained new insight into the crucial hydration effects that contribute to the conformational dynamics of carbohydrates and managed to characterize in detail the poorly understood vibrational nature of this class of biomolecules. Experimental and calculated ROA spectra of these two molecules are reported and excellent agreement has been found. More specifically, comparison has been made with the more commonly used gas phase and implicitly solvated calculation approaches, which offer poor or zero modelling of solvent interactions. The calculated spectra have been used to resolve the structural origins of the observed bands, a current challenge in the study of carbohydrates due to a lack of definitive vibrational assignments. We report and analyse major features in the fingerprint region of the ROA spectra, with recurrent structural and spectral features between the two monosaccharides observed.

The Raman optical activity of β-D-xylose: where experiment and computation meet

Besides its applications in bioenergy and biosynthesis, β-d-xylose is a very simple monosaccharide that exhibits relatively high rigidity. As such, it provides the best basis to study the impact of different solvation shell radii on the computation of its Raman optical activity (ROA) spectrum. Indeed, this chiroptical spectroscopic technique provides exquisite sensitivity to stereochemistry, and benefits much from theoretical support for interpretation. Our simulation approach combines density functional theory (DFT) and molecular dynamics (MD) in order to efficiently account for the crucial hydration effects in the simulation of carbohydrates and their spectroscopic response predictions. Excellent agreement between the simulated spectrum and the experiment was obtained with a solvation radius of 10 Å. Vibrational bands have been resolved from the computed ROA data, and compared with previous results on different monosaccharides in order to identify specific structure-spectrum relationships and to investigate the effect of the solvation environment on the conformational dynamics of small sugars. From the comparison with ROA analytical results, a shortcoming of the classical force field used for the MD simulations has been identified and overcome, again highlighting the complementary role of experiment and theory in the structural characterisation of complex biomolecules. Indeed, due to unphysical puckering, a spurious ring conformation initially led to erroneous conformer ratios, which are used as weights for the averaging of the spectral average, and only by removing this contribution was near perfect comparison between theory and experiment achieved.

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.

Recent advances in the use of vibrational chiroptical spectroscopic methods for stereochemical characterization of natural products

A comprehensive look into application of vibrational optical activity methods for conformational and configurational assignments in natural product molecules over the last 15 years is provided.

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.

Ramachandran Plot for Alanine Dipeptide as Determined from Raman Optical Activity

Accessible values of the φ and ψ torsional angles determining peptide main chain conformation are traditionally displayed in the form of Ramachandran plots. The number of experimental methods making it possible to determine such conformational distribution is limited. In the present study, Raman optical activity (ROA) spectra of Ac-Ala-NHMe were measured and fit by theoretical curves. This revealed the most favored conformers and a large part of the potential energy surface (PES) of this model dipeptide. Such experimental PES compares well to quantum chemical computations, whereas molecular dynamics (MD) modeling reproduces it less faithfully. The surface shape is consistent with the temperature dependence of the spectra, as observed experimentally and predicted by MD. Despite errors associated with spectral modeling and the measurement, the results are likely to facilitate future applications of ROA spectroscopy.

CH Stretching Region: Computational Modeling of Vibrational Optical Activity

Most organic compounds provide vibrational spectra within the CH stretching region, yet the signal is difficult to interpret because of multiple difficulties in experiment and modeling. To better understand various factors involved, the ability of several harmonic and anharmonic computational approaches to describe these vibrations was explored for α-pinene, fenchone, and camphor as test compounds. Raman, Raman optical activity (ROA), infrared absorption (IR), and vibrational circular dichroism (VCD) spectra were measured and compared to quantum chemical computations. Surprisingly, the harmonic vibrational approach reasonably well reproduced the measured spectral patterns, including the vibrational optical activity (VOA). The CH stretching, however, appeared to be more sensitive to the basis set and solvent variations than lower-frequency vibrations. For a higher accuracy in frequencies and spectral shapes, anharmonic corrections were necessary. Accurate harmonic and anharmonic force fields were obtained with the mPW2PLYP double-hybrid functional. A limited vibrational configuration interaction (LVCI) where the CH stretching motion was decoupled from other vibrations provided the best simulated spectra. A balanced harmonic oscillator basis set had to be used, containing also states indirectly interacting with fundamental vibrations. A simpler second-order perturbational approach (PT2) appeared less useful. The modeling provided unprecedented agreement with experimental vibrational frequencies; spectral shapes were reproduced less faithfully. The possibility of ab initio interpretation of the CH spectral region for relatively large molecules further broadens the application span of vibrational spectroscopy.

Determination of absolute configuration and conformation of a cyclic dipeptide by NMR and chiral spectroscopic methods

Increasing precision of contemporary computational methods makes spectroscopies such as vibrational (VCD) and electronic (ECD) circular dichroism attractive for determination of absolute configurations (AC) of organic compounds. This is, however, difficult for polar, flexible molecules with multiple chiral centers. Typically, a combination of several methods provides the best picture of molecular behavior. As a test case, all possible stereoisomers with known AC (RS, SR, SS, and RR) of the cyclic dipeptide cyclo(Arg-Trp) (CAT) were synthesized, and the performances of the ECD, infrared (IR), VCD, Raman, Raman optical activity (ROA), and nuclear magnetic resonance (NMR) techniques for AC determination were investigated. The spectra were interpreted with the aid of density functional theory (DFT) calculations. Folded geometries stabilized by van der Waals and electrostatic interactions between the diketopiperazine (DKP) ring and the indole group are predicted to be preferred for CAT, with more pronounced folding due to Arg-Trp stacking in the case of SS/RR-CAT. The RS/SR isomers prefer a twist-boat puckering of the DKP ring, which is relatively independent of the orientation of the side chains. Calculated conformer-averaged VCD and ECD spectra explain most of the experimentally observed bands and allow for AC determination of the tryptophan side-chain, whereas the stereochemical configuration of the arginine side-chain is visible only in VCD. NMR studies provide characteristic long-range (2)J(C,H) and (3)J(C,H) coupling constants, and nuclear Overhauser effect (NOE) correlations, which in combination with either ECD or VCD also allow for complete AC determination of CAT.

Exploring the active site structure of a photoreceptor protein by Raman optical activity

We have developed a near-infrared excited Raman optical activity (ROA) spectrometer and report the first measurement of near-infrared ROA spectra of a light-driven proton pump, bacteriorhodopsin. Our results demonstrate that a near-infrared excitation enables us to measure the ROA spectra of the chromophore within a protein environment. Furthermore, the ROA spectra of the all-trans, 15-anti and 13-cis, 15-syn isomers differ significantly, indicating a high structural sensitivity of the ROA spectra. We therefore expect that future applications of the near-infrared ROA will allow the experimental elucidation of the active site structures in other proteins as well as reaction intermediates.

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.