Hydration effects on the CN stretching mode frequencies of MeCN and MeSCN are investigated by carrying out ab initio calculations for a number of MeCN-water and MeSCN-water complexes with varying number of water molecules. It is found that the CN frequency shift induced by the hydrogen-bonding interactions with water molecules originate from two different ways to form hydrogen bonds with the nitrogen atom of the CN group. Considering the MeCN- and MeSCN-water cluster calculation results as databases, we first examined the validity of vibrational Stark effect relationship between the CN frequency and the electric field component parallel to the CN bond and found no strong correlation between the two. However, taking into account of additional electric field vector components is a simple way to generalize the vibrational Stark theory for the nitrile chromophore. Also, the electrostatic potential calculation method has been proposed and examined in detail. It turned out that the interactions of water molecules with nitrogen atom’s lone pair orbital and with nitrile π orbitals can be well described by the electrostatic potential calculation method. The present computational results will be of use to quantitatively simulate various linear and nonlinear vibrational spectra of nitrile compounds in solutions.

Nitrile- and thiocyanate-derivatized amino acids have been found to be useful IR probes for investigating their local electrostatic environments in proteins. To shed light on the CN stretch frequency shift and spectral lineshape change induced by interactions with hydrogen-bonding solvent molecules, we carried out both classical and quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulations for MeCN and MeSCN in water. These QM/MM and conventional force field MD simulation results were found to be inconsistent with the experimental results as well as with the high-level ab initio calculation results of MeCN-water and MeSCN-water potential energies. Thus, a new set of atomic partial charges of MeCN and MeSCN is obtained. By using the MD simulation trajectories and the electrostatic potential model recently developed, the CN and SCN stretching mode frequency trajectories were obtained and used to simulate the IR spectra. The CN frequency blueshifts of MeCN and MeSCN in water are estimated to be 9.0 and 1.9cm−1, respectively, in comparison with those of gas phase values. These values are found to be in reasonable agreement with the experimentally measured IR spectra of MeCN, MeSCN, β-cyano-L-alanine, and cyanylated cysteine in water and other polar solvents.

We have implemented the combined quantum mechanical (QM)/molecular mechanical (MM) molecular dynamics (MD) simulations of alanine dipeptide in water along with the polarizable and nonpolarizable classical MD simulations with different models of water. For the QM/MM MD simulation, the alanine dipeptide is treated with the AM1 or PM3 approximations and the fluctuating solute dipole moment is calculated by the Mulliken population analysis. For the classical MD simulations, the solute is treated with the polarizable or nonpolarizable AMBER and polarizable CHARMM force fields and water is treated with the TIP3P, TIP4P, or TIP5P model. It is found that the relative populations of right-handed α-helix and extended β and PII conformations in the simulation trajectory strongly depend on the simulation method. For the QM/MM MD simulations, the PM3/MM shows that the PII conformation is dominant, whereas the AM1/MM predicts that the dominant conformation is αR. Polarizable CHARMM force field gives almost exclusively PII conformation and other force fields predict that both α-helical and extended (β and PII) conformations are populated with varying extents. Solvation environment around the dipeptide is investigated by examining the radial distribution functions and numbers and lifetimes of hydrogen bonds. Comparing the simulated IR and vibrational circular dichroism spectra with experimental results, we concluded that the dipeptide adopts the PII conformation and PM3/MM, AMBER03 with TIP4P water, and AMBER polarizable force fields are acceptable for structure determination of the dipeptide considered in this paper.

A variety of spectroscopic probe molecules have been used to study the local electrostatic environment in proteins. Particularly, a few IR probes such as nitrile- and thiocyanate-derivatized amino acids were found to be quite useful not just because they are small but also because their IR absorption frequencies strongly depend on the strengths of hydrogen bonds with the surrounding protic solvent molecules. Recently, we experimentally demonstrated that azido-derivatized alanine is an excellent IR probe for studying structural change in protein in solution. The asymmetric stretching mode frequency of N3-group becomes blueshifted when it is dissolved in water. Such a blueshifting behavior upon hydrogen-bonding interaction with protic solvent molecules was commonly found in other IR probes containing a triple bond such as CN and SCN groups. In this paper, theoretical descriptions on the solvatochromic frequency shift and fluctuation of the azido stretch frequency are presented by carrying out ab initio calculations and both classical and quantum mechanical/molecular mechanical dynamics simulation studies for azidomethane and azidoalanine dipeptide dissolved in water. Two different methods developed here are based on the distributed multipole interaction models, and they are shown to be useful to describe site-specific hydrogen-bonding interaction-induced red- or blueshift of the azido stretch frequency. To test the validity of thus obtained interpolation formula, numerically simulated IR spectra of azidomethane and azidoalanine dipeptide in water are directly compared with experimental results. We anticipate that the present theoretical approaches will be of use in connecting experimentally measured azido stretch frequency to conformational change in protein containing this azido-derivatized alanine residue.