Driven molecular dynamics studies of the shared proton motion in the H5O2+·Ar cluster: the effect of argon tagging and deuteration on vibrational spectra

We report IR spectra of H5O2(+) and H5O2(+)·Ar and their deuterium isotopologues using ab initio molecular dynamics. The trajectories were propagated as microcanonical (NVE) ensembles at energies corresponding to temperatures 50 and 100 K. The potential energy surface is calculated on-the-fly at the MP2/aug-cc-pVDZ level of theory. The calculations show that adding an argon atom to H5O2(+) introduces symmetry breaking in the Zundel core ion, causes blueshift in the shared proton vibration by about 200 cm(-1), and leads to the splitting of the OH stretch vibrations into four bands. Driven molecular dynamics (DMD) method is used to assign the spectrum by coupling the dipole moment to an external electric field oscillating at frequency ω. The broad feature at 1100 cm(-1) in the H5O2(+)·Ar spectrum is ascribed to the large amplitude shared proton vibration coupled with torsion and wag modes. MD MP2 simulations predict the H/D redshift in the shared proton vibration and water bending vibration to be about 280 and 460 cm(-1), respectively, in good agreement with experimental observations.

Cation-π interactions in protonated phenylalkylamines

Phenylalkylamines of the general formula C6H5(CH2)nNH2 (n = 1-4) have been delivered to the gas phase as protonated species using electrospray ionization. The ions thus formed have been assayed by IRMPD spectroscopy in two different spectroscopic domains, namely, the 600-1800 and the 3000-3500 cm(-1) regions using either an IR free electron laser or a tabletop OPO/OPA laser source. The interpretation of the experimental spectra is aided by density functional theory calculations of candidate species and vibrational frequency analyses. Protonated benzylamine presents a relatively straightforward instance of a single stable conformer, providing a trial case for the adopted approach. Turning to the higher homologues, C6H5(CH2)nNH3(+) (n = 2-4), more conformations become accessible. For each C6H5(CH2)nNH3(+) ion (n = 2-4), the most stable geometry is characterized by cation-π interactions between the positively charged ammonium group and the aromatic π-electronic system, permitted by the folding of the polymethylene chain. The IRMPD spectra of the sampled ions confirm the presence of the folded structures by comparison with the calculated IR spectra of the various possible conformers. An inspection of the NH stretching region is helpful in this regard.

Vibrational spectroscopy and conformational structure of protonated polyalanine peptides isolated in the gas phase

The conformational structures of protonated polyalanine peptides, Ala(n)H(+), have been investigated in the gas phase for n = 3, 4, 5, and 7 using a combination of resonant-infrared multiphoton dissociation (R-IRMPD) spectroscopy in the NH and OH stretch regions and quantum chemical calculations. Agreement between theoretical IR and experimental R-IRMPD spectral features has enabled the assignment of specific hydrogen-bonded conformational motifs in the short protonated peptides and revealed their conformational evolution under elevated-temperature conditions, as a function of increasing chain length. The shortest peptide, Ala(3)H(+), adopts a mixture of extended and cyclic chain conformations, protonated, respectively, at a backbone carbonyl or the N-terminus. The longer peptides adopt folded, cyclic, and globular charge-solvated conformations protonated at the N-terminus, consistent with previous ion-mobility studies. The longest peptide, Ala(7)H(+), adopts a globular conformation with the N-terminus completely charge-solvated, demonstrating the emergence of "physiologically relevant" intramolecular interactions in the peptide backbone. The computed conformational relative free energies highlight the importance of entropic contributions in these peptides.

A THEORETICAL STUDY ON SECOND HARMONIC GENERATION HYPERPOLARIZABILITIES OF PHENYLALANINE POLYPEPTIDES

The second harmonic generation (SHG) hyperpolarizabilities of phenylalanine and homopolypeptides are investigated by configuration interaction among singly excited configurations (CIS) technique combined with the sum-over-states (SOS) method. The geometries of peptides containing phenylalanine ( Phe )n(n = 1–8) are optimized by B3LYP/6-31g(d) method, and they form the special structures like β-sheet (a common protein secondary structure). It is found that the energy gaps of various peptides are reduced and the hyperpolarizabilities are increased with the peptide chains lengthened. We discuss the origin of the second-order nonlinear optical response in phenylalanine homopolypeptides and confirm that the π → π* transitions in the aromatic residue of phenylalanine make the most important contributions to the second-order polarizability. Our results strongly suggest that the hyperpolarizabilities are dominated from the propagation direction of peptide chains.

Chirality effects on the IRMPD spectra of basket resorcinarene/nucleoside complexes

The IRMPD spectra of the ESI-formed proton-bound complexes of the R,R,R,R- and S,S,S,S-enantiomers of a bis(diamido)-bridged basket resorcin[4]arene (R and S) with cytosine (1), cytidine (2), and cytarabine (3) were measured in the region 2800-3600 cm(-1). Comparison of the IRMPD spectra with the corresponding ONIOM (B3LYP/6-31(d):UFF)-calculated absorption frequencies allowed the assessment of the vibrational modes that are responsible for the observed spectroscopic features. All of the complexes investigated, apart from [R⋅H⋅3](+), showed similar IRMPD spectra, which points to similar structural and conformational landscapes. Their IRMPD spectra agree with the formation of several isomeric structures in the ESI source, wherein the N(3)-protonated guest establishes noncovalent interactions with the host amidocarbonyl groups that are either oriented inside the host cavity or outside it between one of the bridged side-chains and the upper aromatic nucleus. The IRMPD spectrum of the [R⋅H⋅3](+) complex was clearly different from the others. This difference is attributed to the effect of intramolecular hydrogen-bonding interactions between the C(2')-OH group and the aglycone oxygen atom of the nucleosidic guest upon repulsive interactions between the same oxygen atom and the aromatic rings of the host.

Characterizing the intramolecular H-bond and secondary structure in methylated GlyGlyH+ with H2 predissociation spectroscopy

We report vibrational predissociation spectra of the four protonated dipeptides derived from glycine and sarcosine, GlyGlyH(+)•(H(2))(1,2), GlySarH(+)•(D(2))(2), SarGlyH(+)•(H(2))(2), and SarSarH(+)•(D(2))(2), generated in a cryogenic ion trap. Sharp bands were recovered by monitoring photoevaporation of the weakly bound H(2) (D(2)) molecules in a linear action regime throughout the 700-4200 cm(-1) range using a table-top laser system. The spectral patterns were analyzed in the context of the low energy structures obtained from electronic structure calculations. These results indicate that all four species are protonated on the N-terminus, and feature an intramolecular H-bond involving the amino group. The large blue-shift in the H-bonded N-H fundamental upon incorporation of a methyl group at the N-terminus indicates that this modification significantly lowers the strength of the intramolecular H-bond. Methylation at the amide nitrogen, on the other hand, induces a significant rotation (~110°) about the peptide backbone.

Infrared spectroscopy of the alanine dipeptide analog in liquid water with DFT-MD. Direct evidence for P(II)/beta conformations

Following our previous work [J. Phys. Chem. B. Lett., 2009, 113, 10059], DFT-based molecular dynamics (DFTMD) simulations of 2-Ala peptide (i.e. Ac-Ala-NHMe dialanine peptide analog with methyl group caps at the extremities) immersed in liquid water at room temperature are reported. Our goal here is the theoretical calculation of the infrared spectrum of aqueous 2-Ala, in order to provide a definitive understanding of the average conformation adopted by this peptide in the liquid phase, taking into account solute and solvent at the same theoretical level of representation. We find that the experimental Amide I-II band predominantly results from a mixture of partially unfolded P(II) and unfolded beta conformational equilibrium of aqueous 2-Ala at room temperature.

IRPD spectroscopy and ensemble measurements: effects of different data acquisition and analysis methods

Three different commonly used infrared photodissociation (IRPD) spectroscopy acquisition and analysis methods are described, and results from these methods are compared using the same dataset for an extensively hydrated metal cation, La(3+)(H(2)O)(36). Using the first-order laser-induced photodissociation rate constant as an IRPD intensity has several advantages over photodissociation yield and depletion/appearance methods in that intensities can be more directly compared with calculated infrared absorption spectra, and the intensities can be readily corrected for changes in laser power or irradiation times used for optimum data acquisition at each frequency. Extending IRPD spectroscopy to large clusters can be complicated when blackbody infrared radiative dissociation competes strongly with laser-induced photodissociation. A new method to obtain IRPD spectra of single precursor ions or ensembles of precursor ions that is nearly equivalent to the photodissociation rate constant method for single precursor ions is demonstrated. The ensemble IRPD spectra represent the "average" structure of clusters of a given size range, and this method has the advantage that spectra with improved signal-to-noise ratios can be obtained with no increase in data acquisition time. Results using this new method for a precursor ensemble consisting of La(3+)(H(2)O)(35-37) are compared with results for La(3+)(H(2)O)(36).

Theoretical spectroscopy of floppy peptides at room temperature. A DFTMD perspective: gas and aqueous phase

Theoretical spectroscopy is mandatory for a precise understanding and assignment of experimental spectra recorded at finite temperature. We review here room temperature DFT-based molecular dynamics simulations for the purpose of interpreting finite temperature infrared spectra of peptides of increasing size and complexity, in terms of temperature-dependent conformational dynamics and flexibility, and vibrational anharmonicities (potential energy surface anharmonicities, vibrational mode couplings and dipole anharmonicities). We take examples from our research projects in order to illustrate the main key-points and strengths of dynamical spectra modeling in that context. The calculations are presented in relation to room temperature gas phase IR-MPD experiments and room temperature liquid phase IR absorption experiments. These illustrations of floppy polypeptides have been chosen in order to convey the following ideas: temperature-dependent spectra modeling is pivotal for a precise understanding of gas phase spectra recorded at room temperature, including conformational dynamics and vibrational anharmonicities; harmonic spectroscopy (as commonly performed in the literature) can be misleading and even erroneous for a proper interpretation of spectra recorded at finite temperature; taking into account vibrational anharmonicities is pivotal for a proper interplay between theory and experiments; amide I-III bands are not necessarily the most relevant fingerprints for unraveling the local structures of peptides and more complex systems; liquid phase simulations have unraveled relationships between the zwitterionic properties of the peptide bonds and infrared signatures. The review presents a state-of-the-art account of the domain and offers perspectives and new developments for future still more challenging applications.

Infrared consequence spectroscopy of gaseous protonated and metal ion cationized complexes

In this article, the new and exciting techniques of infrared consequence spectroscopy (sometimes called action spectroscopy) of gaseous ions are reviewed. These techniques include vibrational predissociation spectroscopy and infrared multiple photon dissociation spectroscopy and they typically complement one another in the systems studied and the information gained. In recent years infrared consequence spectroscopy has provided long-awaited direct evidence into the structures of gaseous ions from organometallic species to strong ionic hydrogen bonded structures to large biomolecules. Much is being learned with respect to the structures of ions without their stabilizing solvent which can be used to better understand the effect of solvent on their structures. This review mainly covers the topics with which the author has been directly involved in research: structures of proton-bound dimers, protonated amino acids and DNA bases, amino acid and DNA bases bound to metal ions and, more recently, solvated ionic complexes. It is hoped that this review reveals the impact that infrared consequence spectroscopy has had on the field of gaseous ion chemistry.

Vibrational spectroscopy of bare and solvated ionic complexes of biological relevance

The low density of ions in mass spectrometers generally precludes direct infrared (IR) absorption measurements. The IR spectrum of an ion can nonetheless be obtained by inducing photodissociation of the ion using a high-intensity tunable laser. The emergence of free electron lasers (FELs) and recent breakthroughs in bench-top lasers based on nonlinear optics have now made it possible to routinely record IR spectra of gas-phase ions. As the energy of one IR photon is insufficient to cause dissociation of molecules and strongly bound complexes, two main experimental strategies have been developed to effect photodissociation. In infrared multiple-photon dissociation (IR-MPD) many photons are absorbed resonantly and their energy is stored in the bath of vibrational modes, leading to dissociation. In the "messenger" technique a weakly bound van der Waals atom is detached upon absorption of a single photon. Fundamental, historical, and practical aspects of these methods will be presented. Both of these approaches make use of very different methods of ion preparation and manipulation. While in IR-MPD ions are irradiated in trapping mass spectrometers, the "messenger" technique is generally carried out in molecular beam instruments. The main focus of this review is the application of IR spectroscopy to biologically relevant molecular systems (amino acids, peptides, proteins). Particular issues that will be addressed here include gas-phase zwitterions, the (chemical) structures of peptides and their collision-induced dissociation (CID) products, IR spectra of gas-phase proteins, and the chelation of metal-ligand complexes. Another growing area of research is IR spectroscopy on solvated clusters, which offer a bridge between the gas-phase and solution environments. The development of state-of-the-art computational approaches has gone hand-in-hand with advances in experimental techniques. The main advantage of gas-phase cluster research, as opposed to condensed-phase experiments, is that the systems of interest can be understood in detail and structural effects can be studied in isolation. It will be shown that IR spectroscopy of mass-selected (bio)molecular systems is now well-placed to address specific questions on the individual effect of charge carriers (protons and metal ions), as well as solvent molecules on the overall structure.

The ab initio assigning of the vibrational probing modes of tryptophan: linear shifting of approximate anharmonic frequencies vs. multiplicative scaling of harmonic frequencies

To gain insight into the prospects for a few-dimensional ab initio quantum-mechanical description of the vibrational motions of conformationally flexible molecular systems, the NH-, NH(2)-, CO- and OH-stretching and COH-bending vibrations of the most stable tryptophan conformations have been probed using simple one- and two-dimensional anharmonic Hamiltonians and potential energy functions evaluated by means of the standard RI-MP2, CCSD(T) and DFT-D quantum chemical procedures. Although strongly dependent on the procedure used, the calculated vibrational spectral patterns have been found to be in a robust one-to-one harmony with their experimental counterparts, thus proving the adequacy of the theory used for the reliable assignment of the experimental data. Therefore, the approach appears to be a suitable tool for assigning the vibrational probing modes even of systems which are too large to be tractable by the standard normal-coordinate analysis.

Alanine polypeptide structural fingerprints at room temperature: what can be gained from non-harmonic Car-Parrinello molecular dynamics simulations

Structural infrared fingerprints of neutral gas phase alanine peptides of increasing size and complexity (dipeptide, octapeptide, and beta-strand peptide) are characterized through DFT-based Car-Parrinello molecular dynamics simulations. Harmonic and nonharmonic vibrational signatures are calculated from the time correlation of the dipole moment of the gas phase peptide in a direct way (without any approximation) respectively from low temperature (20 K) and room temperature (300 K) molecular dynamics. Our main purpose is to answer the two following questions: (i) Is the direct inclusion of temperature for the calculation of infrared spectra mandatory for the comprehension of the vibrational signatures experimentally recorded at room temperature? (ii) To what extent is the amide I, II, and III domain sensitive enough to the local structure of the peptides, to provide vibrational signatures that can be definitely used to assess the peptide conformation at 300 K?

Conformational preferences of an amyloidogenic peptide: IR spectroscopy of Ac-VQIVYK-NHMe

The (306)VQIVYK(311) sequence in the tau peptide is essential for the formation of intracellular amyloid fibrils related to Alzheimer's disease, where it forms interdigitating cross-beta-structures. The inherent conformational preferences of the capped hexapeptide Ac-VQIVYK-NHMe were characterized in the gas phase. IR/UV double-resonance spectroscopy of the peptide isolated in a cold molecular beam was used to probe the conformation of the neutral peptide. The influence of protonation at the lysine side chain was investigated at 298 K by characterizing the protonated peptide ion, Ac-VQIVYK(H(+))-NHMe, with IRMPD spectroscopy in the fingerprint and amide I/II band region in an FTICR mass spectrometer. The conformations for both neutral and protonated peptides were predicted by an extensive conformational search procedure followed by cluster analysis and then DFT calculations. Comparison of the experimental and computed IR spectra, with consideration of the relative energies, was used to assign the dominant conformations observed. The neutral peptide adopts a beta-hairpin-like conformation with two loosely extended peptide chains, demonstrating the preference of the sequence for extended beta-strand-like structures. In the protonated peptide, the lysine NH3(+) disrupts this extended conformation by binding to the backbone and induces a transition to a random-coil-like structure.

Quantum Chemical Benchmark Energy and Geometry Database for Molecular Clusters and Complex Molecular Systems (www.begdb.com): A Users Manual and Examples

Our previous benchmark CCSD(T)/ complete basis set limit calculations were collected into a database named begdb - Benchmark Energy and Geometry DataBase. Web-based interface to this database was prepared and is available at www.begdb.com. Users can browse, search and plot the data online or download structures and energy tables.