Dataset Information

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules.

ABSTRACT: Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural representation of the potential energy surface. In this paper, a nuclear quantum vibrational perturbation (QVP) method is described for efficiently determining the instantaneous vibrational frequency of a chromophore in molecular dynamics simulations. Computational efficiency is achieved through the use of (a) discrete variable representation of the vibrational wave functions, (b) a perturbation theory to evaluate the vibrational energy shifts due to solvent dynamic fluctuations, and (c) a combined QM/MM potential for the systems. It was found that first-order perturbation is sufficiently accurate, enabling time-dependent vibrational frequencies to be obtained on the fly in molecular dynamics. The QVP method is illustrated in the mode-specific linear and 2D-IR spectra of the H-Cl stretching frequency in the HCl-water clusters and the carbonyl stretching vibration of acetone in aqueous solution. To further reduce computational cost, a hybrid strategy was proposed, and it was found that the computed vibrational spectral peak position and line shape are in agreement with experimental results. In addition, it was found that anharmonicity is significant in the H-Cl stretching mode, and hydrogen-bonding interactions further enhance anharmonic effects. The present QVP method complements other computational approaches, including path integral-based molecular dynamics, and represents a major improvement over the electrostatics-based spectroscopic mapping procedures.

SUBMITTER: Xue RJ

PROVIDER: S-EPMC5793873 | biostudies-literature | 2017 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules.

Xue Rui-Jie RJ Grofe Adam A Yin He H Qu Zexing Z Gao Jiali J Li Hui H

Journal of chemical theory and computation 20161207 1

Linear and two-dimensional infrared (IR) spectroscopy of site-specific probe molecules provides an opportunity to gain a molecular-level understanding of the local hydrogen-bonding network, conformational dynamics, and long-range electrostatic interactions in condensed-phase and biological systems. A challenge in computation is to determine the time-dependent vibrational frequencies that incorporate explicitly both nuclear quantum effects of vibrational motions and an electronic structural repre ...[more]

PMID: 28068771

Similar Datasets

Project description:The vibrational spectroscopy of hydration water in dilauroylphosphatidylcholine lipid multi-bilayers is investigated using molecular dynamics simulations and a mixed quantum/classical model for the OD stretch spectroscopy of dilute HDO in H(2)O. FTIR absorption spectra, and isotropic and anisotropic pump-probe decay curves have been measured experimentally as a function of the hydration level of the lipid multi-bilayer, and our goal is to make connection with these experiments. To this end, we use third-order response functions, which allow us to include non-Gaussian frequency fluctuations, non-Condon effects, molecular rotations, and a fluctuating vibrational lifetime, all of which we believe are important for this system. We calculate the response functions using existing transition frequency and dipole maps. From the experiments it appears that there are two distinct vibrational lifetimes corresponding to HDO molecules in different molecular environments. In order to obtain these lifetimes, we consider a simple two-population model for hydration water hydrogen bonds. Assuming a different lifetime for each population, we then calculate the isotropic pump-probe decay, fitting to experiment to obtain the two lifetimes for each hydration level. With these lifetimes in hand, we then calculate FTIR spectra and pump-probe anisotropy decay as a function of hydration. This approach, therefore, permits a consistent calculation of all observables within a unified computational scheme. Our theoretical results are all in qualitative agreement with experiment. The vibrational lifetime of lipid-associated OD groups is found to be systematically shorter than that of the water-associated population, and the lifetimes of each population increase with decreasing hydration, in agreement with previous analysis. Our theoretical FTIR absorption spectra successfully reproduce the experimentally observed red-shift with decreasing lipid hydration, and we confirm a previous interpretation that this shift results from the hydrogen bonding of water to the lipid phosphate group. From the pump-probe anisotropy decay, we confirm that the reorientational motions of water molecules slow significantly as hydration decreases, with water bound in the lipid carbonyl region undergoing the slowest rotations.

Project description:In this paper, we describe an automated integration-free path-integral (AIF-PI) method, based on Kleinert's variational perturbation (KP) theory, to treat internuclear quantum-statistical effects in molecular systems. We have developed an analytical method to obtain the centroid potential as a function of the variational parameter in the KP theory, which avoids numerical difficulties in path-integral Monte Carlo or molecular dynamics simulations, especially at the limit of zero-temperature. Consequently, the variational calculations using the KP theory can be efficiently carried out beyond the first order, i.e., the Giachetti-Tognetti-Feynman-Kleinert variational approach, for realistic chemical applications. By making use of the approximation of independent instantaneous normal modes (INM), the AIF-PI method can readily be applied to many-body systems. Previously, we have shown that in the INM approximation, the AIF-PI method is accurate for computing the quantum partition function of a water molecule (3 degrees of freedom) and the quantum correction factor for the collinear H(3) reaction rate (2 degrees of freedom). In this work, the accuracy and properties of the KP theory are further investigated by using the first three order perturbations on an asymmetric double-well potential, the bond vibrations of H(2), HF, and HCl represented by the Morse potential, and a proton-transfer barrier modeled by the Eckart potential. The zero-point energy, quantum partition function, and tunneling factor for these systems have been determined and are found to be in excellent agreement with the exact quantum results. Using our new analytical results at the zero-temperature limit, we show that the minimum value of the computed centroid potential in the KP theory is in excellent agreement with the ground state energy (zero-point energy) and the position of the centroid potential minimum is the expectation value of particle position in wave mechanics. The fast convergent property of the KP theory is further examined in comparison with results from the traditional Rayleigh-Ritz variational approach and Rayleigh-Schrödinger perturbation theory in wave mechanics. The present method can be used for thermodynamic and quantum dynamic calculations, including to systematically determine the exact value of zero-point energy and to study kinetic isotope effects for chemical reactions in solution and in enzymes.

Project description:I.r. spectroscopy has been applied to the study of hydrogen-bonding of the unique ester carbonyl group of acylchymotrypsins in the oxyanion hole of the enzyme. This catalytic device provides electrophilic stabilization of negative charge in the transition states and tetrahedral intermediates along the reaction pathway. The use of 13C isotope substitution of the ester carbonyl group reinforces the previous observation [White & Wharton (1990) Biochem. J. 270, 627-637] that the ester carbonyl group is significantly polarized in the ground state by hydrogen bonding in the oxyanion hole. I.r. difference spectra of [carbonyl-12C]-minus [carbonyl-13C]-cinnamoyl-chymotrypsin as well as each of these acylenzymes minus free enzyme are reported. These spectra show that the contribution of protein perturbation (i.e. spectral features that arise from the enzyme which is distorted on acylation) in [carbonyl-12C]cinnamoyl-chymotrypsin minus free enzyme spectra is significant. The contribution of the perturbation components of the spectra is pH-dependent and can represent up to 50% of the total absorbance in the spectral region from 1690 to 1740 cm-1. Use of the isotopic difference method has allowed problems associated with protein perturbation to be eliminated. Similar difference spectra are presented for dihydrocinnamoyl-chymotrypsin. In this case the effect of perturbation is very marked and leads to the cancellation of the band assigned to the non-bonded conformation of the acyl group which has previously only been observed at higher pH. The isotopic difference method again proves reliable and shows that the frequency difference previously used to calculate the ground-state electronic strain induced by the oxyanion-hole catalytic device is not affected by the perturbation, although the amplitudes of the spectral features are different. A study of the deacylation of cinnamoyl-chymotrypsin in water and deuterium oxide using both u.v. and i.r. spectroscopies has confirmed that the use of deuterium oxide as solvent has no serious effect on the deacylation behaviour of the enzyme. I.r. bands assigned to nonproductive and productive conformers decline identically during deacylation, which shows that the conformers are in dynamic exchange on the reaction time-scale.

Dataset Information

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules.

Publications

Perturbation Approach for Computing Infrared Spectra of the Local Mode of Probe Molecules.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets