Modeling the infrared and circular dichroism spectroscopy of a bridged cyclic diamide

Near-Ultraviolet Circular Dichroism and Two-Dimensional Spectroscopy of Polypeptides

A fully quantitative theory of the relationship between protein conformation and optical spectroscopy would facilitate deeper insights into biophysical and simulation studies of protein dynamics and folding. In contrast to intense bands in the far-ultraviolet, near-UV bands are much weaker and have been challenging to compute theoretically. We report some advances in the accuracy of calculations in the near-UV, which were realised through the consideration of the vibrational structure of the electronic transitions of aromatic side chains.

Fundamental frequency from classical molecular dynamics

We give a theoretical validation for calculating fundamental frequencies of a molecule from classical molecular dynamics (MD) when its anharmonicity is small enough to be treated by perturbation theory. We specifically give concrete answers to the following questions: (1) What is the appropriate initial condition of classical MD to calculate the fundamental frequency? (2) From that condition, how accurately can we extract fundamental frequencies of a molecule? (3) What is the benefit of using ab initio MD for frequency calculations? Our analytical approaches to those questions are classical and quantum normal form theories. As numerical examples we perform two types of MD to calculate fundamental frequencies of H2O with MP2/aug-cc-pVTZ: one is based on the quartic force field and the other one is direct ab initio MD, where the potential energies and the gradients are calculated on the fly. From those calculations, we show comparisons of the frequencies from MD with the post vibrational self-consistent field calculations, second- and fourth-order perturbation theories, and experiments. We also apply direct ab initio MD to frequency calculations of C-H vibrational modes of tetracene and naphthalene. We conclude that MD can give the same accuracy in fundamental frequency calculation as second-order perturbation theory but the computational cost is lower for large molecules.

Vibronic structure in the far-UV electronic circular dichroism spectra of proteins

The Franck–Condon effect is considered and the vibrational structure of the π_nbπ* transition of the peptide backbone is incorporated into matrix method calculations of the electronic circular dichroism (CD) spectra of proteins in the far-ultraviolet. We employ the state-averaged CASPT2 method to calculate the ground and π_nbπ* excited state geometries and frequencies of N-methylacetamide (NMA), which represents the peptide chromophore. The results of these calculations are used to incorporate vibronic levels of the excited states into the matrix method calculation. The CD spectra of a set of 49 proteins, comprising a range of structural types, are calculated to assess the influence of the vibrational structure. The calculated spectra of α-helical proteins are better resolved using the vibronic parameters and correlation between the experimental and the calculated intensity of less regular β structure proteins improves over most wavelengths in the far-UV. No obvious improvement is observed in the calculated spectra of regular β-sheet proteins. Our high-level ab initio calculations of the vibronic structure of the π_nbπ* transition in NMA have provided some further insight into the physical origins of the nature of protein CD spectra in the far-UV.

Electronic circular dichroism of the chiral rigid tricyclic dilactam with nonplanar tertiary amide groups

Electronic circular dichroism (ECD) of the spirocyclic dilactam 5,8-diazatricyclo[6,3,0,0(1,5)]undecane-4,9-dione has been measured in the extended wavelength range (170-260 nm) utilizing far-UV CD instrumentation including synchrotron radiation light source. The data of this model of two nonplanar tertiary amide groups interacting within the rigid chiral environment provided new information particularly about the shorter wavelength π-π* transition region below 190 nm. The interpretation using TDDFT calculations confirmed that effects of amide nonplanarity follow our previous observations on monolactams as far as amide n-π* transitions are concerned. ECD band in the n-π* transition region of the nonplanar diamide exhibits an identical bathochromic shift and its sign remains tied to the sense of nonplanar deformation in the same way. As far as n-π* transitions are concerned amide nonplanarity acts as a local phenomenon independently reflecting sum properties of single amide groups. On the other hand, CD bands associated with π-π* transitions (found between ∼170 to 210 nm) form an exciton-like couplet with the sign pattern determined by mutual orientation of the associated electric transition moments. This sign pattern follows predictions pertaining to a coupled oscillator. The influence of amide nonplanarity on π-π* transitions is only minor and concentrates into the shorter wavelength lobe of the π-π* couplet. The detailed analysis of experimental ECD with the aid of TDDFT calculations shows that there is only little interaction between effects of inherent chirality caused by nonplanarity of amide groups and amide-amide coupling. Consequently these two effects can be studied nearly independently using ECD. In addition, the calculations indicate that participation of other type of transitions (n-σ*, π-σ* or Rydberg type transitions) is only minor and is concentrated below 180 nm.

Fully anharmonic IR and Raman spectra of medium-size molecular systems: accuracy and interpretation

Computation of full infrared (IR) and Raman spectra (including absolute intensities and transition energies) for medium- and large-sized molecular systems beyond the harmonic approximation is one of the most interesting challenges of contemporary computational chemistry. Contrary to common beliefs, low-order perturbation theory is able to deliver results of high accuracy (actually often better than those issuing from current direct dynamics approaches) provided that anharmonic resonances are properly managed. This perspective sketches the recent developments in our research group toward the development of a robust and user-friendly virtual spectrometer rooted in second-order vibrational perturbation theory (VPT2) and usable also by non-specialists essentially as a black-box procedure. Several examples are explicitly worked out in order to illustrate the features of our computational tool together with the most important ongoing developments.

Anharmonic vibrations of N-H in Cl(-)(N-methylacetamide)1(H2O)(0-2)Ar2 cluster ions. Combined IRPD experiments and BOMD simulations

Infrared Predissociation (IRPD) spectra of Cl(-)(NMA)1(H2O)0-2Ar2 combined with Born-Oppenheimer Molecular Dynamics (BOMD) IR spectra have been acquired, providing the structure and dynamics of these systems. We show that the chloride ion is bound to the hydrogen of the amide N-H group, forming a strong ionic hydrogen bond, weakening the N-H stretch, and shifting it to lower frequency. The presence of water molecules enhances the ionic hydrogen bond by binding to the amide carbonyl oxygen of NMA and shifts the N-H stretch further to lower frequency. The BOMD IR spectra can recapture all, but about 100 cm(-1), of the 600 to 700 cm(-1) shifts due to the strong N-H stretch anharmonicities observed in experiments. This residual error was found to be due to the lack of zero point energy in the classical treatment of motion in the BOMD method.

Computational study of the structure and electronic circular dichroism spectroscopy of blue copper proteins

The calculation of the electronic circular dichroism (CD) spectra of the oxidized form of the blue copper proteins plastocyanin and cucumber basic protein and the relationship between the observed spectral features and the structure of the active site of the protein is investigated. Excitation energies and transition strengths are computed using multireference configuration interaction, and it is shown that computed spectra based on coordinates from the crystal structure or a single structure optimized in quantum mechanics/molecular mechanics (QM/MM) or ligand field molecular mechanics (LFMM) are qualitatively incorrect. In particular, the rotational strength of the ligand to metal charge transfer band is predicted to be too small or have the incorrect sign. By considering calculations on active site models with modified structures, it is shown that the intensity of this band is sensitive to the nonplanarity of the histidine and cysteine ligands coordinated to copper. Calculation of the ultraviolet absorption and CD spectra based upon averaging over many structures drawn from a LFMM molecular dynamics simulation are in good agreement with experiment, and superior to analogous calculations based upon structures from a classical molecular dynamics simulation. This provides evidence that the LFMM force field provides an accurate description of the molecular dynamics of these proteins.