Lactic acid in solution: investigations of lactic acid self-aggregation and hydrogen bonding interactions with water and methanol using vibrational absorption and vibrational circular dichroism spectroscopies

The infrared vibrational absorption (VA) and vibrational circular dichroism (VCD) spectral features of L-(+)-lactic acid (LA) in CDCl3 solution are concentration dependent, showing evidence of oligomerization with increasing concentrations. To understand the observed spectra, geometry optimizations, vibrational frequencies, and VA and VCD intensities were evaluated for (LA)n with n=1-4 using density functional theory calculations at the B3LYP6-311++G(d,p), B3LYP/cc-pVTZ, and in some cases, B3LYP/aug-cc-pVTZ levels of theory. Comparisons with the experimental spectra indicate that the lowest energy LA dimer (AA), formed by two C Double Bond O...HO hydrogen bonds, is one of the dominating species in solution at room temperature. Possible contributions from the LA trimer and tetramer are also discussed. To model the VA and VCD spectra of LA in water and in methanol, both implicit polarizable continuum model and explicit hydrogen bonding considerations were used. For explicit hydrogen bonding, geometry optimizations of the AA-(water)n and AA-(methanol)n complexes, with n=2,4,6, were performed, and the corresponding VA and VCD spectra were simulated. Comparisons of the calculated and experimental VA and VCD spectra in the range of 1000-1800 cm(-1) show that AA-(water)n with n=6 best reproduces the experimental spectra in water. On the other hand, AA-(methanol)n with n=2 reproduces well the experimental results taken in methanol solution. In addition, we found evidence of chirality transfer, i.e., some vibrational bands of the achiral water subunits gain VCD strength upon complexation with the chiral LA solute. The study is the first to use VCD spectroscopy to probe the structures of LA aggregates and hydrogen bonding solvation clusters in the solution phase.     

Chirality transfer through hydrogen-bonding: experimental and ab initio analyses of vibrational circular dichroism spectra of methyl lactate in water

The infrared vibrational absorption (VA) and vibrational circular dichroism (VCD) spectra of methyl lactate were measured in the 1000-1800 cm(-1) region in the CCl(4) and H(2)O solvents, respectively. In particular, the chirality transfer effect, i.e. the H-O-H bending bands of the achiral water subunits that are hydrogen-bonded to the methyl lactate molecule exhibit substantial VCD strength, was detected experimentally. A series of density functional theory calculations using B3PW91 and B3LYP functionals with 6-311++G(d,p) and aug-cc-pVTZ basis sets were carried out to simulate the VA and VCD spectra of the methyl lactate monomer and the methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3. The population weighted VA and VCD spectra of the methyl lactate monomer are in good agreement with the experimental spectra in CCl(4). Implicit polarizable continuum model was found to be inadequate to account for the hydrogen-bonding effect in the observed VA and VCD spectra in H(2)O. The methyl lactate-(H(2)O)(n) complexes with n = 1, 2, 3 were used to model the explicit hydrogen-bonding. The population weighted VA and VCD spectra of the methyl lactate-H(2)O binary complex are shown to capture the main spectral features in the observed spectra in aqueous solution. The theoretical modeling shows that the extent of chirality transfer depends sensitively on the specific binding sites taken by the achiral water molecules. The observation of chirality transfer effect opens a new spectral window to detect and to model the hydrogen-bonding solvent effect on VCD spectra of chiral molecules.     

Matrix Isolation‐Vibrational Circular Dichroism Spectroscopy of 3‐Butyn‐2‐ol and its Binary Aggregates

Vibrational circular dichroism (VCD) spectroscopy has a unique specificity to chirality and is highly sensitive to the conformational equilibria of chiral molecules. On the other hand, the matrix‐isolation (MI) technique allows substantial control over sample compositions, such as the sample(s)/matrix ratio and the ratio among different samples, and yields spectra with very narrow bandwidths. We combined VCD spectroscopy with the MI technique to record MI‐VCD and MI‐vibrational absorption spectra of 3‐butyn‐2‐ol at different MI temperatures, which allowed us to investigate the conformational distributions of its monomeric and binary species. Good mirror‐imaged MI‐VCD spectra of opposite enantiomers were achieved. The related conformational searches were performed for the monomer and the binary aggregate and their vibrational absorption and VCD spectra were simulated. The well‐resolved experimental MI‐VCD bands provide the essential mean to assign the associated vibrational absorption spectral features correctly to a particular conformation in case of closely spaced bands. By varying the matrix temperature, we show that one can follow the self‐aggregation process of 3‐butyn‐2‐ol and confidently correlate the MI‐VCD spectral features with those obtained for a 0.1 M CCl₄ solution and as a neat liquid at room temperature. Comparison of the aforementioned experimental VCD spectra shows conclusively that there is a substantial contribution from the 3‐butyn‐2‐ol aggregate even at 0.1 M concentration. This spectroscopic combination will be powerful for studying self‐aggregation of chiral molecules, and chirality transfer from a chiral molecule to an interacting achiral molecule and in electron donor–acceptor chiral complexes.     

Absolute configuration of C2-symmetric spiroselenurane: 3,3,3',3'-tetramethyl-1,1'-spirobi[3 H,2,1]benzoxaselenole

The enantiomers of 3,3,3',3'-tetramethyl-1,1'-spirobi[3 H,2,1]benzoxaselenole have been separated on a chiral preparative chromatographic column. The experimental vibrational circular dichroism (VCD) spectra have been obtained for both enantiomers in CH(2)Cl(2). The theoretical VCD spectra have been obtained by means of density functional theoretical calculations with the B3 LYP density functional. From a comparison of experimental and theoretical VCD spectra, the absolute configuration of an enantiomer with positive specific rotation in CH(2)Cl(2) at 589 nm is determined to be R. This conclusion has been verified by comparing results of experimental optical rotatory dispersion (ORD) and electronic circular dichroism (ECD) to predictions of the same properties using the B3 LYP functional for the title compound.     

Conformational Distributions of N‐Acetyl‐L‐cysteine in Aqueous Solutions: A Combined Implicit and Explicit Solvation Treatment of VA and VCD Spectra

The conformational distributions of N‐acetyl‐L ‐cysteine (NALC) in aqueous solutions at several representative pH values are investigated using vibrational absorption (VA), UV/Vis, and vibrational circular dichroism (VCD) spectroscopy, together with DFT and molecular dynamics (MD) simulations. The experimental VA and UV/Vis spectra of NALC in water are obtained under strongly acid, neutral, and strongly basic conditions, as well as the VCD spectrum at pH 7 in D₂O. Extensive searches are carried out to locate the most stable conformers of the protonated, neutral, deprotonated, and doubly deprotonated NALC species at the B3LYP/6‐311++G(d,p) level. The inclusion of the polarizable continuum model (PCM) modifies the geometries and the relative stabilities of the conformers noticeably. The simulated PCM VA spectra show significantly better agreement with the experimental data than the gas‐phase ones, thus allowing assignment of the conformational distributions and dominant species under each experimental condition. To further properly account for the discrepancies noted between the experimental and simulated VCD spectra, PCM and the explicit solvent model are utilized. MD simulations are used to aid the modelling of the NALC–(water)_N clusters. The geometry optimization, harmonic frequency calculations, and VA and VCD intensities are computed for the NALC–(water)_3,4 clusters at the B3LYP/6‐311++G(d,p) level without and with the PCM. The inclusion of both explicit and implicit solvation models at the same time provides a decisively better agreement between theory and experiment and therefore conclusive information about the conformational distributions of NALC in water and hydrogen‐bonding interactions between NALC and water molecules.     

Conformations of Serine in Aqueous Solutions as Revealed by Vibrational Circular Dichroism

Vibrational circular dichroism (VCD) spectroscopy is utilized to reveal the detailed conformational distributions of the dominant serine species in aqueous solutions under three representative pH conditions of 1.0, 5.7, and 13.0, together with vibrational absorption (VA) spectroscopy, density functional theory (DFT), and molecular dynamics simulation. The experimental VA and VCD spectra of serine in H₂O and D₂O in the fingerprint region under three pH values are obtained. DFT calculations at the B3LYP/6‐311++G(d,p) level are carried out for the protonated, zwitterionic, and deprotonated serine species. The lowest‐energy conformers of all three species are identified and their corresponding VA and VCD spectra simulated. A comparison between the gas‐phase simulations and the experimental VA and VCD spectra suggests that one or two of the most stable conformers of each species contribute predominantly to the observed data, although some discrepancies are noted. To account for the solvent effects, both the polarizable continuum model and the explicit solvation model are considered. Hydrogen‐bonded protonated, zwitterionic, and deprotonated serine–(water)₆ clusters are constructed based on radial distribution function analyses and molecular dynamics snapshots. Geometry optimization and VA and VCD simulations are performed for these clusters at the B3LYP/6‐311++G(d,p) level. Inclusion of the explicit water molecules is found to improve the agreement between theory and experiment noticeably in all three cases, thus enabling conclusive conformational distribution analyses of serine in aqueous solutions directly.     

A single chiroptical spectroscopic method may not be able to establish the absolute configurations of diastereomers: dimethylesters of hibiscus and garcinia acids

Electronic circular dichroism (ECD), optical rotatory dispersion (ORD), and vibrational circular dichroism (VCD) spectra of hibiscus acid dimethyl ester have been measured and analyzed in combination with quantum chemical calculations of corresponding spectra. These results, along with those reported previously for garcinia acid dimethyl ester, reveal that none of these three (ECD, ORD, or VCD) spectroscopic methods, in isolation, can unequivocally establish the absolute configurations of diastereomers. This deficiency is eliminated when a combined spectral analysis of either ECD and VCD or ORD and VCD methods is used. It is also found that the ambiguities in the assignment of absolute configurations of diastereomers may also be overcome when unpolarized vibrational absorption is included in the spectral analysis.     

Enantiopure, monodisperse alleno-acetylenic cyclooligomers: effect of symmetry and conformational flexibility on the chiroptical properties of carbon-rich compounds

A series of enantiopure, monodisperse alleno-acetylenic cyclooligomers were synthesized. The single-crystal X-ray structures of the cyclic trimer and hexamer were resolved, providing insights into the symmetry of these molecules. Electronic circular dichroism (ECD), optical rotatory dispersion (ORD), Raman spectroscopy, and vibrational circular dichroism (VCD) data were analyzed with the aid of theoretical calculations. This multidimensional approach ultimately provided general guidelines that are useful for designing carbon-rich compounds with intense chiroptical properties.     

Determination of the absolute configurations of natural products via density functional theory calculations of vibrational circular dichroism, electronic circular dichroism, and optical rotation: the iridoids plumericin and isoplumericin

The absolute configurations (ACs) of the iridoid natural products, plumericin (1) and isoplumericin (2), have been re-investigated using vibrational circular dichroism (VCD) spectroscopy, electronic circular dichroism (ECD) spectroscopy, and optical rotatory dispersion (ORD). Comparison of DFT calculations of the VCD spectra of 1 and 2 to the experimental VCD spectra of the natural products, (+)-1 and (+)-2, leads unambiguously to the AC (1R,5S,8S,9S,10S)-(+) for both 1 and 2. In contrast, comparison of time-dependent DFT (TDDFT) calculations of the ECD spectra of 1 and 2 to the experimental spectra of (+)-1 and (+)-2 does not permit definitive assignment of their ACs. On the other hand, TDDFT calculations of the ORD of (1R,5S,8S,9S,10S)-1 and -2 over the range of 365-589 nm are in excellent agreement with the experimental data of (+)-1 and (+)-2, confirming the ACs derived from the VCD spectra. Thus, the ACs initially proposed by Albers-Sch?nberg and Schmid are shown to be correct, and the opposite ACs recently derived from the ECD spectra of 1 and 2 by Els?sser et al. are shown to be incorrect. As a result, the ACs of other iridoid natural products obtained by chemical correlation with 1 and 2 are not in need of revision.     

A configurational and conformational study of aframodial and its diasteriomers via experimental and theoretical VA and VCD spectroscopies

In this work we present the experimental and theoretical vibrational absorption (VA) and the theoretical vibrational circular dichroism (VCD) spectra for aframodial. In addition, we present the theoretical VA and VCD spectra for the diasteriomers of aframodial. Aframodial has four chiral centers and hence has 2⁴ = 16 diasteriomers, which occur in eight pairs of enantiomers. In addition to the four chiral centers, there is an additional chirality due to the helicity of the entire molecule, which we show by presenting 12 configurations of the 5S,8S,9R,10S enantiomer of aframodial. The VCD spectra for the diasteriomers and the 12 configurations of one enantiomer are shown to be very sensitive not only to the local stereochemistry at each chiral center, but in addition, to the helicity of the entire molecule. Here one must be careful in analyzing the signs of the VCD bands due to the ‘non-chiral’ chromophores in the molecule, since one has two contributions; one due to the inherent chirality at the four chiral centers, and one due to the chirality of the side chain groups in specific conformers, that is, its helicity. Theoretical simulations for various levels of theory are compared to the experimental VA recorded to date. The VCD spectra simulations are presented, but no experimental VCD and Raman spectra have been reported to date, though some preliminary VCD measurements have been made in Stephens’ lab in Los Angeles. The flexible side chain is proposed to be responsible for the small size of the VCD spectra of this molecule, even though the chiral part of the molecule is very rigid and has four chiral centers. In addition to VCD and Raman measurements, Raman optical activity (ROA) measurements would be very enlightening, as in many cases bands which are weak in both the VA and VCD, may be large in the Raman and/or ROA spectra. The feasibility of using vibrational spectroscopy to monitor biological structure, function and activity is a worthy goal, but this work shows that a careful theoretical analysis is also required, if one is to fully utilize and understand the experimental results. The reliability, reproduceability and uniqueness of the vibrational spectroscopic experiments and the information which can be gained from them is discussed, as well as the details of the computation of VA, VCD and Raman (and ROA) spectroscopy for molecules of the complexity of aframodial, which have multiple chiral centers and flexible side chains. Festschrift in Honor of Philip J. Stephens’ 65th Birthday.     

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.     

Determination of the Absolute Configuration of Pentacoordinate Chiral Phosphorus Compounds in Solution by Using Vibrational Circular Dichroism Spectroscopy and Density Functional Theory

Vibrational circular dichroism (VCD) spectroscopic measurements and density functional theory (DFT) calculations have been used to obtain the absolute structural information about four sets of diastereomers of pentacoordinate spirophosphoranes derived separately from l‐ (or d‐ ) valine and l‐ (or d‐ ) leucine for the first time. Each compound contains three stereogenic centers: one at the phosphorus center and two at the amino acid ligands. Extensive conformational searches for the compounds have been carried out and their vibrational absorption (VA) and VCD spectra have been simulated at the B3LYP/6‐311++G** level. Although both VA and VCD spectra are highly sensitive to the structural variation of the apical axis, that is, the O? P? O or N? P? O arrangement, the rotamers generated by the aliphatic amino side chains show little effect on both. The dominant experimental VCD features in the 1100–1500 cm^?1 region were found to be controlled by the chirality at the phosphorus center, whereas those at the C?O stretching region are determined by the chirality of the amino acid ligands. The good agreement between the experimental VA and VCD spectra in CDCl₃ solution and the simulated ones allows us to assign the absolute configurations of these pentacoordinate phosphorus compounds with high confidence. This study shows that the VCD spectroscopy complemented with DFT calculations is a powerful and reliable method for determining the absolute configurations and dominating conformers of synthetic phosphorus coordination complexes in solution.     

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.     

Vibrational Circular Dichroism Spectroscopy of Three Multidentate Nitrogen Donor Ligands: Conformational Flexibility and Solvent Effects

A series of multidentate nitrogen donor ligands have been synthesized and characterized and their conformational distributions in solution have been investigated. Vibrational absorption (VA) and vibrational circular dichroism (VCD) spectroscopy, complemented with DFT calculations, have been used to probe the conformations of these important ligands in solution directly. These three ligands demonstrate very different conformational flexibility; the pyridine subunits and amine groups may adopt a number of different conformations. Experimental VA and VCD data measured in CDCl₃ have been compared to the theoretical spectra of all possible most stable conformers. Solvent effects have been taken into account by using the implicit polarizable continuum model and explicit solvation model. The explicit hydrogen‐bonding solvation model is important for explaining the VCD sign‐reverse phenomenon in the amide I region. Good agreement has been achieved between experimental and predicted spectra for all three ligands; thus allowing detailed examination of the related conformational structures and distributions in solution.     

Spectroscopic investigation of the structures of dialkyl tartrates and their cyclodextrin complexes

Structures of three dialkyl tartrates, namely, dimethyl tartrate, diethyl tartrate, and diisopropyl tartrate, in CCl4, dimethyl sulfoxide (DMSO)/DMSO-d6, and H2O/D2O solvents have been investigated using vibrational absorption (VA), vibrational circular dichroism (VCD), and optical rotatory dispersion (ORD). VA, VCD, and ORD spectra are found to be dependent on the solvent used. Density functional theory (DFT) calculations are used to interpret the experimental data in CCl4 and DMSO. The trans-COOR conformer with hydrogen bonding between the OH group and the C=O group attached to the same chiral carbon is dominant for dialkyl tartrates both in vacuum and in CCl4. The experimental VA, VCD, and ORD data of dialkyl-D-tartrates in CCl4 correlated well with those predicted for dimethyl-(S,S)-tartrate molecule as both isolated and solvated in CCl4. In DMSO solvent, dialkyl tartrate molecules favor formation of intermolecular hydrogen bonding with DMSO molecules. Clusters of dimethyl-(S,S)-tartrate, with one molecule of dimethyl-(S,S)-tartrate hydrogen bonded to two DMSO molecules, are used for the DFT calculations. A trans-COOR cluster and a trans-H cluster are needed to obtain a reasonable agreement between the predicted and experimental data of dimethyl tartrate in DMSO solvent. VA, VCD, and optical rotations are also measured for dialkyl tartrate-cyclodextrin complexes. It is noted that these properties are barely affected by complexation of dialkyl tartrates with cyclodextrins, indicating weak interaction between tartrates and cyclodextrin. Binding constants of alpha-CD and beta-CD with diethyl L-tartrate in both H2O and DMSO have been determined using isothermal titration calorimetry technique. The smaller binding constants (less than 100) confirmed the weak interaction between tartrates and cyclodextrin in the solution state.     

光学活性和手性光谱的溯源和发展

概述了19世纪以来光学活性和手性光谱的发现和发展,着重对旋光色散(ORD)、电子圆二色(ECD)和振动圆二色(VCD)光谱的发展背景和基本原理作出介绍。其中特别提及华人科学家徐光宪和徐云洁在手性光谱发展历程中的杰出贡献。     

Vibrational circular dichroism spectroscopy of two chiral binaphthyl diphosphine ligands and their palladium complexes in solution

BINAP (2,2'-diphenylphosphino-1,1'-binaphthyl) is a unique binaphthyl diphosphine ligand with axial chirality. The palladium complexes of BINAP and of its derivative TOLBINAP have found extensive applications in the field of asymmetric syntheses. The conformational changes in the BINAP and TOLBINAP ligands before and after coordination with palladium have been investigated using density functional theory, vibrational absorbance (VA) and vibrational circular dichroism (VCD) spectroscopy. VA and VCD spectra of these two chiral ligands and their corresponding palladium complexes have been recorded in CDCl(3) solution. Extensive conformational searches have been carried out for both the ligands and the associated palladium complexes. Coordination with palladium has been found to introduce structural rigidity to the ligands. The calculated VA and VCD spectra of the ligands and complexes in the gas phase show substantial differences to the experimental data. Incorporation of the implicit polarisable continuum solvation model has provided much better agreement between theory and experiment, especially for the complexes, allowing clear identification of the species and conformations. This and the high specificity of VCD spectral signatures to chirality and to conformations suggest the potential applications of VCD spectroscopy for following these important catalytic species in solution reactions directly.     

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.     

Solvent effects on IR and VCD spectra of helical peptides: DFT-based static spectral simulations with explicit water

Simulations of IR and VCD spectra are carried out for model alpha-helical, 3(10)-helical, and 3(1)-helical (polyProII-like) oligopeptides, with up to 21 amide groups, and including explicit consideration of effects of directly hydrogen-bonded solvent (water). Parameters used were obtained from ab initio density functional theory (DFT) computations of force field, atomic polar and axial tensors for oligopeptides of 5 to 7 amides, whose structures were constrained in (phi,psi) to target the secondary structure type but otherwise fully optimized. By comparison with experimental data as well as with calculations for identical but isolated (gas phase) peptides, the computed effects of an inner shell of aqueous solvent on the vibrational spectra of helical oligopeptides are illustrated. The interaction with solvent causes significant frequency shifts of the amide bands, but only minor changes in the characteristic IR intensity distributions and splittings and the VCD band shapes. Better agreement with experimental band shapes is achieved for the alpha-helical amide I' (N-deuterated) VCD by inclusion of explicit solvent in the calculations. Some improvements are also observed in theoretical VCD predictions for 13C labeled alpha-helical peptides when solvated models are used. However, the qualitative isotopic splitting patterns are preserved and just shifted in frequency due to consistent, solvent independent interamide coupling constants. The critical match of experiment and theory for relative positions of transitions in peptides with specifically separated 13C=O labels, including and neglecting solvent, confirms the stability of the coupling interactions. Despite these solvation effects, the calculated VCD band shape of the amide I mode is shown to be a reliable conformational probe, since it remains basically insensitive to frequency shifts caused by environment. Thus theoretical VCD simulations, even vacuum calculations, are shown to provide useful spectral predictions for solution-phase peptides.