Calculation of Raman optical activity spectra of methyl-β-D-glucose incorporating a full molecular dynamics simulation of hydration effects

We report calculations of the Raman and Raman optical activity (ROA) spectra of methyl-β-D-glucose utilizing density functional theory combined with molecular dynamics (MD) simulations to provide an explicit hydration environment. This is the first report of such combination of MD simulations with ROA ab initio calculations. We achieve a significant improvement in accuracy over the more commonly used gas phase and polarizable continuum model (PCM) approaches, resulting in an excellent level of agreement with the experimental spectrum. Modeling the ROA spectra of carbohydrates has until now proven a notoriously difficult challenge due to their sensitivity to the effects of hydration on the molecular vibrations involving each of the chiral centers. The details of the ROA spectrum of methyl-β-D-glucose are found to be highly sensitive to solvation effects, and these are correctly predicted for the first time including those originating from the highly sensitive low frequency vibrational modes. This work shows that a thorough consideration of the role of water is pivotal for understanding the vibrational structure of carbohydrates and presents a new and powerful tool for characterizing carbohydrate structure and conformational dynamics in solution.     

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, that is, 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.     

Experimental and ab initio theoretical vibrational Raman optical activity of alanine

The vibrational Raman optical activity (ROA) spectra of l-alanine in water, 1 N NaOH and 1 N HCl between 720 and 1500 cm⁻¹ measured in backscattering are reported. Unlike the associated vibrational circular dichroism (VCD), the main ROA features are relatively insensitive to pH changes. Ab initio Raman and ROA intensities were evaluated using 6-31G and 6-31G* basis sets and found to agree remarkably well with the experimental parameters in the lower-frequency region.     

Infrared and vibrational CD spectra of partially solvated alpha-helices: DFT-based simulations with explicit solvent

Theoretical simulations are used to investigate the effects of aqueous solvent on the vibrational spectra of model alpha-helices, which are only partly exposed to solvent to mimic alpha-helices in proteins. Infrared absorption (IR) and vibrational circular dichroism (VCD) amide I' spectra for 15-amide alanine alpha-helices are simulated using density functional theory (DFT) calculations combined with the property transfer method. The solvent is modeled by explicit water molecules hydrogen bonded to the solvated amide groups. Simulated spectra for two partially solvated model alpha-helices, one corresponding to a more exposed and the other to a more buried structure, are compared to the fully solvated and unsolvated (gas phase) simulations. The dependence of the amide I spectra on the orientation of the partially solvated helix with respect to the solvent and effects of solvation on the amide I' of 13C isotopically substituted alpha-helices are also investigated. The partial exposure to solvent causes significant broadening of the amide I' bands due to differences in the vibrational frequencies of the explicitly solvated and unsolvated amide groups. The different degree of partial solvation is reflected primarily in the frequency shifts of the unsolvated (buried) amide group vibrations. Depending on which side of the alpha-helix is exposed to solvent, the simulated IR band-shapes exhibit significant changes, from broad and relatively featureless to distinctly split into two maxima. The simulated amide I' VCD band-shapes for the partially solvated alpha-helices parallel the broadening of the IR and exhibit more sign variation, but generally preserve the sign pattern characteristic of the alpha-helical structures and are much less dependent on the alpha-helix orientation with respect to the solvent. The simulated amide I' IR spectra for the model peptides with explicitly hydrogen-bonded water are consistent with the experimental data for small alpha-helical proteins at very low temperatures, but overestimate the effects of solvent on the protein spectra at ambient temperatures, where the peptide-water hydrogen bonds are weakened by thermal motion.     

Two Spectroscopies in One: Interference of Circular Dichroism and Raman Optical Activity

Previously, we and other laboratories have reported an unusual and strong Raman optical activity (ROA) induced in solvents by chiral dyes. Various theories of the phenomenon appeared, but they were not capable of explaining fully the observed ROA band signs and intensities. In this work, an analysis based both on the light scattering theory and dedicated experiments provides a more complete understanding. For example, double-cell magnetic circular dichroism and magnetic ROA experiments with copper-porphyrin complex show that the induced chirality is observed without any contact of the solvents with the complex. The results thus indicate that a combination of electronic circular dichroism (ECD) with the polarized Raman scattering is responsible for the effect. The degree of circularity of solvent vibrational bands is a principal molecular property participating in the event. The insight and the possibility to predict the chirality transfer promise future applications in spectroscopy, chemical analysis and polarized imaging.     

Ab initio studies of solvent effect on molecular conformation and vibrational spectra of diacetamide

The conformational equilibria and vibrational spectra of diacetamide have been investigated by ab initio molecular orbital studies using the basis sets 6-31g(d,p) and 6-31++g(d,p) at Hartree-Fock and MP2 levels. The vibrational spectra of diacetamide have been satisfactorily interpreted taking into consideration the agreement between the calculated harmonic vibrational frequencies, infrared and Raman band intensities and shifts in deuterated molecules with those observed. The solvent effects were investigated by the self-consistent reaction field (SCRF) theory. The effect of solvent on the conformational equilibria and vibrational spectra is discussed. The calculated changes in the geometry and vibrational spectra on going from the gas phase to the solvent medium are in accord with the increasing weight of the dipolar resonance structure of the amide group in more polar solvents.     

Understanding the Signatures of Secondary‐Structure Elements in Proteins with Raman Optical Activity Spectroscopy

A prerequisite for the understanding of functional molecules like proteins is the elucidation of their structure under reaction conditions. Chiral vibrational spectroscopy is one option for this purpose, but provides only indirect access to this structural information. By first‐principles calculations, we investigate how Raman optical activity (ROA) signals in proteins are generated and how signatures of specific secondary‐structure elements arise. As a first target we focus on helical motifs and consider polypeptides consisting of twenty alanine residues to represent α‐helical and 3₁₀‐helical secondary‐structure elements. Although ROA calculations on such large molecules have not been carried out before, our main goal is the stepwise reconstruction of the ROA signals. By analyzing the calculated ROA spectra in terms of rigorously defined localized vibrations, we investigate in detail how total band intensities and band shapes emerge. We find that the total band intensities can be understood in terms of the reconstructed localized vibrations on individual amino acid residues. Two different basic mechanisms determining the total band intensities can be established, and it is explained how structural changes affect the total band intensities. The band shapes can be rationalized in terms of the coupling between the localized vibrations on different residues, and we show how different band shapes arise as a consequence of different coupling patterns. As a result, it is demonstrated for the chiral variant of Raman spectroscopy how collective vibrations in proteins can be understood in terms of well‐defined localized vibrations. Based on our calculations, we extract characteristic ROA signatures of α helices and of 3₁₀‐helices, which our analysis directly relates to differences in secondary structure.     

Theoretical study of the Raman optical activity spectra of 3(10)-helical polypeptides

Raman optical activity (ROA) spectroscopy is a promising analytical method for studying the structure and conformation of polypeptides and proteins in solution. However, the structural information obtained from such vibrational spectra is only indirect and theoretical studies are often necessary to identify how the structure determines the observed spectra. One particular target is the identification and discrimination of different helical secondary structure elements. Herein, a theoretical investigation of the ROA spectra of a series of 3(10)-helical polypeptides is presented. In particular, the effect of chain length, C(α)-substitution pattern, the introduction of larger aliphatic side chains, and the variation of their conformation on the ROA spectra is studied. To extract general principles from these calculations, the positions, intensities, and shapes of the ROA bands are analyzed in terms of localized modes, which makes it possible to identify possible ROA signatures of 3(10) -helical structures, but also provides fundamental insight into the generation of ROA signals in complex polypeptides. Finally, the calculated spectra can be compared to the previously reported ROA spectrum of a specifically designed 3(10) -helical heptapeptide. This allows most of the features in the experimental spectrum to be assigned.     

A combined theoretical and experimental study of the structure and vibrational absorption, vibrational circular dichroism, Raman and Raman optical activity spectra of the L-histidine zwitterion

A combined theoretical and experimental study of the vibrational absorption (VA)/IR, vibrational circular dichroism (VCD), Raman and Raman optical activity (ROA) spectra of l-histidine in aqueous solution has been undertaken to answer the questions (i) what are the species present and (ii) which conformers of the species are present under various experimental conditions. The VA spectra of l-histidine have been measured in aqueous solution and the spectral bands which can be used to identify both species (cation, zwitterion, anion) and conformer of the species have been identified and subsequently used to identify the species (zwitterion) and conformer (gauche minus minus, gauche minus plus for the side chain dihedral angles) present in solution at pH 7.6. The VCD spectral intensities have been used subsequently in combination with further theoretical studies to confirm the conclusions that have been arrived at by only analyzing the VA/IR spectra. Finally a comparison of measured Raman and ROA spectra of l-histidine with Raman and ROA spectral simulations for the conformers and species derived from the combined VA/IR and VCD experimental and theoretical work is presented as a validation of the conclusions arrived at from VA/IR and VCD spectroscopy. The combination of VA/IR and VCD with Raman and ROA is clearly superior and both sets of experiments should be performed.     

Raman and theoretical study of the solvent effects on the sizable intramolecular charge transfer in the push-pull 5-(dimethylamino)-5'-nitro-2,2'-bithiophene

In this paper, we analyze the degree of intramolecular charge transfer in a push-pull pi-conjugated system, 5-(dimethylamino)-5'-nitro-2,2'-bithiophene, from changes in frequencies and relative intensities of its strongest Raman scatterings in a bunch of solvents with different polarities. Density functional theory (DFT) was used as a support of the experimental study. Solvent effects on the molecular and electronic structures and on the vibrational properties were estimated by performing B3LYP/6-31G calculations within the framework of the polarized continuum model (PCM) developed by Tomasi. Calculations reveal that the molecule is highly polarized in the ground state and behaves as a very efficient photoinduced push-pull system. The polarization of the molecule strongly increases with solvent polarity and determines that the profile of the Raman spectra greatly changes from one solvent to another and in going to the solid. The strongest Raman scattering associated with the nu(sym)(NO(2)) stretching undergoes a downshift of 48 cm(-1) in passing from CCl(4) to the solid. DFT calculations provide a comprehensive interpretation of the evolution of the Raman spectra with solvent polarity.     

Effects of hydrogen bond and solvent polarity on the C=O stretching of bis(2-thienyl)ketone in solution

The optimized structural parameters, the absorption and the resonance Raman spectra have been investigated for the bis(2-thienyl)ketone in gas phase, in cyclohexane, methanol, and acetonitrile solvents by means of time dependent density functional theory calculations, the solvent electronic polarization effect on the solvation shift is examined and in well accordance with the calculation. The effect of increasing the polarity of the solvent is well represented by the polarizable continuum model, both for the absorption spectra and resonance Raman intensities. The Raman spectra of the C=O stretching mode, which is sensitive to the intermolecular interaction for bis(2-thienyl)ketone dissolved in solvents, were systematically studied. It was found that the hydrogen bond effect plays an important role in reducing the carbonyl stretching wavenumbers. The results of Raman shifts were interpreted through the dilution effect, solvation effects, and hydrogen bond-forming effects. Furthermore, the excitation profiles of several important Raman bands of bis(2-thienyl)ketone molecule in different solvents have been critically analyzed. The solvent effects on structural and symmetry properties of the molecule in S2 electronic state as well as the short-time photo relaxation dynamics have been discussed.     

Theoretical studies of UO2(H2O)n(2+), NpO2(H2O)n(+), and PuO2(H2O)n(2+) complexes (n=4-6) in aqueous solution and gas phase

Extensive ab initio calculations both in gas phase and solution have been carried out to study the equilibrium structure, vibrational frequencies, and bonding characteristics of various actinyl (UO2(2+), NpO2(+), and PuO2(2+)) and their hydrated forms, AnO2(H2O)n(z+) (n=4, 5, and 6). Bulk solvent effects were studied using a continuum method. The geometries were fully optimized at the coupled-cluster singles + doubles (CCSD), density-functional theory (DFT), and M?ller-Plesset (MP2) level of theories. In addition vibrational frequencies have been obtained at the CCSD as well as MP2/DFT levels. The results show that both the short-range and long-range solvent effects are important. The combined discrete-continuum model, in which the ionic solute and the solvent molecules in the first and second solvation shells are treated quantum mechanically while the solvent is simulated by a continuum model, can predict accurately the bonding characteristics. Moreover, our values of solvation free energies suggest that five- and six-coordinations are equally preferred for UO2(2+), and five-coordinated species are preferred for NpO2(+) and PuO2(2+). On the basis of combined quantum-chemical and continuum treatments of the hydrated complexes, we are able to determine the optimal cavity radii for the solvation models. The coupled-cluster computations with large basis sets were employed for the vibrational spectra and equilibrium geometries both of which compare quite favorably with experiment. Our most accurate computations reveal that both five- and six-coordination complexes are important for these species.     

Preferential hydration of lysozyme in water/glycerol mixtures: a small-angle neutron scattering study

In solution small-angle neutron scattering has been used to study the solvation properties of lysozyme dissolved in water/glycerol mixtures. To detect the characteristics of the protein-solvent interface, 35 different experimental conditions (i.e., protein concentration, water/glycerol fraction in the solvent, content of deuterated compounds) have been considered and a suitable software has been developed to fit simultaneously the whole set of scattering data. The average composition of the solvent in the close vicinity of the protein surface at each experimental condition has been derived. In all the investigated conditions, glycerol resulted especially excluded from the protein surface, confirming that lysozyme is preferentially hydrated. By considering a thermodynamic hydration model based on an equilibrium exchange between water and glycerol from the solvation layer to the bulk, the preferential binding coefficient and the excess solvation number have been estimated. Results were compared with data previously derived for ribonuclease A in the same mixed solvent: even if the investigated solvent compositions were very different, the agreement between data is noticeable, suggesting that a unique mechanism presides over the preferential hydration process. Moreover, the curve describing the excess solvation number as a function of the solvent composition shows the occurrence of a region of maximal hydration, which probably accounts for the changes in protein stability detected in the presence of cosolvents.     

A frequency-weighted sum rule for vibrational Raman optical activity

The sums of frequency-weighted vibrational Raman optical activity (ROA) intensities of isotopic molecules are shown to be related through the elements of the inverse force constant matrix. These relations could be used to determine force constants from experimental ROA intensities.     

Demonstration of the ring conformation in polyproline by the Raman optical activity

Raman and Raman optical activity (ROA) spectra of poly-L-proline were recorded in a wide frequency range and analyzed with respect to the proline side chain conformation. The analysis was based on comparison to ab initio simulations of spectral band positions and intensities. The presence of two conformer states of the five-member ring was found, approximately equally populated in the polypeptide. Additionally, Raman and ROA spectral shapes indicated that the peptide adopts the polyproline II helical conformation, in both aqueous and TFE solutions. The helix, however, is perturbed by fluctuations, which affects the vibrational coupling among amino acid residues and broadens the ROA bands. Contributions of the side and main peptide chains to the polyproline ROA intensities have comparable magnitudes. Thus understanding of the origins of both signals is important for determination of the peptide structure by ROA.     

Role of hydration in determining the structure and vibrational spectra of L-alanine and N-acetyl L-alanine N′-methylamide in aqueous solution: a combined theoretical and experimental approach

In this work we have utilized recent density functional theory Born-Oppenheimer molecular dynamics simulations to determine the first principles locations of the water molecules in the first solvation shell which are responsible for stabilizing the zwitterionic structure of L-alanine. Previous works have used chemical intuition or classical molecular dynamics simulations to position the water molecules. In addition, a complete shell of water molecules was not previously used, only the water molecules which were thought to be strongly interacting (H-bonded) with the zwitterionic species. In a previous work by Tajkhorshid et al. (J Phys Chem B 102:5899) the L-alanine zwitterion was stabilized by 4 water molecules, and a subsequent work by Frimand et al. (Chem Phys 255:165) the number was increased to 9 water molecules. Here we found that 20 water molecules are necessary to fully encapsulate the zwitterionic species when the molecule is embedded within a droplet of water, while 11 water molecules are necessary to encapsulate the polar region with the methyl group exposed to the surface, where it migrates during the MD simulation. Here we present our vibrational absorption, vibrational circular dichroism and Raman and Raman optical activity simulations, which we compare to the previous simulations and experimental results. In addition, we report new VA, VCD, Raman and ROA measurements for L-alanine in aqueous solution with the latest commercially available FTIR VA/VCD instrument (Biotools, Jupiter, FL, USA) and Raman/ROA instrument (Biotools). The signal to noise of the spectra of L-alanine measured with these new instruments is significantly better than the previously reported spectra. Finally we reinvestigate the causes for the stability of the P_π structure of the alanine dipeptide, also called N-acetyl-L-alanine N′-methylamide, in aqueous solution. Previously we utilized the B3LYP/6-31G* + Onsager continuum level of theory to investigate the stability of the NALANMA4WC Han et al. (J Phys Chem B 102:2587) Here we use the B3PW91 and B3LYP hybrid exchange correlation functionals, the aug-cc-pVDZ basis set and the PCM and CPCM (COSMO) continuum solvent models, in addition to the Onsager and no continuum solvent model. Here by the comparison of the VA, VCD, Raman and ROA spectra we can confirm the stability of the NALANMA4WC due to the strong hydrogen bonding between the four water molecules and the peptide polar groups. Hence we advocate the use of explicit water molecules and continuum solvent treatment for all future spectral simulations of amino acids, peptides and proteins in aqueous solution, as even the structure (conformer) present cannot always be found without this level of theory. Festschift in Honor of Philip J. Stephens’ 65th Birthday. During the proof stage of this article a very relevant article has been published by M. Losada and Y. Xu titled “Chirality transfer through hydrogen-bonding: Experimental and ab initio analyses of vibrational circular dichroism spectra of methyl lactate in water” in Phys Chem Chem Phys 2007, 9: 3127–3135. In that work they confirm that the effects of water are seen in the VCD spectra and hence it is fundamental to include explicit water molecules in modeling studies of the vibrational spectra of biomolecules in aqueous solution.     

Transfer and Amplification of Chirality Within the “Ring of Fire” Observed in Resonance Raman Optical Activity Experiments

We report extremely strong chirality transfer from a chiral nickel complex to solvent molecules detected as Raman optical activity (ROA). Electronic energies of the complex were in resonance with the excitation‐laser light. The phenomenon was observed for a wide range of achiral and chiral solvents. For chiral 2‐butanol, the induced ROA was even stronger than the natural one. The observations were related to so‐called quantum (molecular) plasmons that enable a strong chiral Rayleigh scattering of the resonating complex. According to a model presented here, the maximal induced ROA intensity occurs at a certain distance from the solute, in a three‐dimensional “ring of fire”, even after rotational averaging. Most experimental ROA signs and relative intensities could be reproduced. The effect might significantly increase the potential of ROA spectroscopy in bioimaging and sensitive detection of chiral molecules.     

Comparison of the Electronic and Vibrational Optical Activity of a Europium(III) Complex

The geometry and the electronic structure of chiral lanthanide(III) complexes are traditionally probed by electronic methods, such as circularly polarised luminescence (CPL) and electronic circular dichroism (ECD) spectroscopy. The vibrational phenomena are much weaker. In the present study, however, significant enhancements of vibrational circular dichroism (VCD) and Raman optical activity (ROA) spectral intensities were observed during the formation of a chiral bipyridine–Eu^III complex. The ten‐fold enhancement of the vibrational absorption and VCD intensities was explained by a charge‐transfer process and the dominant effect of the nitrate ion on the spectra. A much larger enhancement of the ROA and Raman intensities and a hundred‐fold increase of the circular intensity difference (CID) ratio were explained by the resonance of the λ=532 nm laser light with the ⁷F₀ → ⁵D₀ transitions. This phenomenon is combined with a chirality transfer, and mixing of the Raman and luminescence effects involving low‐energy ⁷F states of europium. The results thus indicate that the vibrational optical activity (VOA) may be a very sensitive tool for chirality detection and probing of the electronic structure of Eu^III and other coordination compounds.