Infrared Spectra of Protonated Neurotransmitters: Serotonin - The

Nov 22, 2010 - ... Monique van Scherpenzeel , Leo A. J. Kluijtmans , Ron A. Wevers ... Nitzan Mayorkas , Hanan Sachs , Markus Schütz , Shun-ichi Ishi...
0 downloads 0 Views 2MB Size
13268

J. Phys. Chem. A 2010, 114, 13268–13276

Infrared Spectra of Protonated Neurotransmitters: Serotonin Anita Lagutschenkov,† Judith Langer,† Giel Berden,‡ Jos Oomens,‡,§ and Otto Dopfer*,† Institut fu¨r Optik und Atomare Physik, Technische UniVersita¨t Berlin, Hardenbergstrasse 36, 10623 Berlin, Germany, FOM Institute for Plasma Physics Rijnhuizen, Edisonbaan 14, 3439 MN Nieuwegein, The Netherlands, and UniVersity of Amsterdam, Science Park 904, Amsterdam 1098XH, The Netherlands ReceiVed: September 29, 2010; ReVised Manuscript ReceiVed: NoVember 4, 2010

The gas-phase IR spectrum of the protonated neurotransmitter serotonin (5-hydroxytryptamine) was measured in the fingerprint range by means of IR multiple photon dissociation (IRMPD) spectroscopy. The IRMPD spectrum was recorded in a Fourier transform ion cyclotron resonance mass spectrometer coupled to an electrospray ionization source and an IR free electron laser. Quantum chemical calculations at the B3LYP and MP2 levels of theory using the cc-pVDZ basis set yield six low-energy isomers in the energy range up to 40 kJ/mol, all of which are protonated at the amino group. Protonation at the indole N atom or the hydroxyl group is substantially less favorable. The IRMPD spectrum is rich in structure and exhibits 22 distinguishable features in the spectral range investigated (530-1885 cm-1). The best agreement between the measured IRMPD spectrum and the calculated linear IR absorption spectra is observed for the conformer lowest in energy at both levels of theory, denoted g-1. In this structure, one of the three protons of the ammonium group points toward the indole subunit, thereby maximizing the intramolecular NH+-π interaction between the positive charge of the ammonium ion and the aromatic indole ring. This mainly electrostatic cation-π interaction is further stabilized by significant dispersion forces, as suggested by the substantial differences between the DFT and MP2 energies. The IRMPD bands are assigned to individual normal modes of the g-1 conformer, with frequency deviations of less than 29 cm-1 (average 90° or < -90°). The gauche isomers are further divided into those with positive and negative φ1 values, as indicated by g+ and g-, respectively. All gauche and trans isomers can further be classified by the orientation of the hydroxyl group described by φ3. The notation 1 and 2 differentiates isomers, in which the hydroxyl group is either in syn (1) or in anti (2) orientation with respect to the indole N atom. In general, there is good agreement between the relative energies (∆E) and free energies (∆G) of the serotoninH+ isomers calculated at the B3LYP and MP2 levels of theory (Figure 1). For example, for the four lowest-energy structures lying below 10 kJ/mol (all g isomers), the agreement is better than (1.5 kJ/mol. Both theoretical levels predict the g-1 isomer to be the global minimum on the potential energy surface. In this isomer, the ammonium group interacts preferentially with the phenol ring via NH+-π interaction. The g+1 isomer is slightly less stable than g-1, with an energy gap of 3.1 (4.6) kJ/mol at the B3LYP (MP2) level. In this isomer, the ammonium group interacts with the pyrrole ring via NH+-π interaction. The two isomers are separated by a substantial barrier of 4.8 kJ/mol (B3LYP) for g+1 f g-1 (i.e., 7.9 kJ/mol for g+1 r g-1), indicating that the internal rotation of the ammonium group above the indole plane is strongly hindered. Similar barrier heights have been reported previously, e.g. 5.3 kJ/mol for g+1 f g-1 at the B3LYP/6-31+G(d,p) level.14 As a general trend, the g(2 structures are less stable than the g(1 isomers by 5-10 kJ/mol. Apparently, the ion-dipole interaction between the positively charged ammonium group and the polar OH bond favors an orientation in which the hydroxyl group points away from the NH3+ unit. Interestingly, the energetic order predicted by B3LYP and MP2 differs for the g+2 and g-2 structures, which are, however, close in energy. B3LYP prefers g+2 slightly over g-2 by 0.2 kJ/mol, whereas MP2 yields an energetic preference for g-2 over g+2 by 1.7 kJ/mol. The two t isomers, t1 and t2, are considerably less stable than the corresponding g isomers by ∼20-30 kJ/mol, because they lack the cation-π interaction of the ammonium group with the indole ring. The barrier between the t1 isomer and g-1 is appreciable and amounts to 10.1 kJ/mol for t1 f g-1 at the B3LYP level (i.e., 33.7 kJ/mol for t1 r g-1). A similar barrier of 9.0 kJ/ mol was previously obtained at the B3LYP/6-31+G(d,p) level.14 Although there is in general good agreement between the relative energies of the various isomers calculated at the MP2 and B3LYP levels, there is a systematically larger energy difference between corresponding g and t isomers at the MP2 level. The relative energies of the t isomers are higher by 7 kJ/mol at the MP2 level than at the B3LYP level. This additional relative stabilization of the g isomers with respect to the t isomers at the MP2 level is attributed to dispersion interactions of the ammonium group with the aromatic indole ring, which are relevant only for the g isomers and are neglected in the B3LYP calculations.15 These additional dispersion forces of the intramolecular NH-π bond also lead to a shorter distance between the proton donor of the ammonium group and the aromatic ring at the MP2 level (by ∼0.09 Å for g-1). Similar effects have

Lagutschenkov et al. TABLE 2: Selected Bond Distances (in Å), Dihedral Angles (in Degrees), and Relative Energies and Free Energies (in kJ/mol) of the Protonated Serotonin Isomers g-1 and t1 Calculated at the B3LYP and MP2 Levels (Figures 1 and 2)a g-1

φ1 φ2 φ3 RNH · · · C3 RNH · · · C4 RNH · · · C9 ∆E ∆G

t1

B3LYP

MP2

B3LYP

MP2

-54.7 -94.2 175.4 2.430 2.468 2.312 0.0 0.0

-52.8 -91.6 173.8 2.322 2.485 2.220 0.0 0.0

-173.5 -107.4 178.3

-173.3 -107.6 177.3

23.6 21.3

30.7 28.5

a φ1 ) dihedral angle C3C10C11N2, φ2 ) dihedral angle C2C3C10C11, φ3 ) dihedral angle C4C5O1H.

previously been noted for dopamineH+,28 and to a lesser extent for neutral serotonin,12 demonstrating the substantial contribution of dispersion for these intramolecular cation-π interactions. Density functional calculations using functionals neglecting (B3LYP, B2LYP) and including dispersion (B3LYPD, B2LYP-D, M06-2X) confirm that the additional stabilization of 7 kJ/mol for the g isomers predicted at the MP2 level is indeed due to dispersion and not due to intramolecular basis set superposition error or possible overestimation of dispersion energy at the MP2 level (Table T1 in the Supporting Information). Further serotoninH+ structures with the excess proton attached to the N atom of the indole ring or at the O atom of the hydroxyl group are less stable by ∼100 and ∼170 kJ/mol than the g-1 isomer, respectively, at both levels of theory (Figure F1 in the Supporting Information). Protonation of primary amines at the amino group is clearly favored in the gas phase. Table 2 summarizes selected structural and energetic parameters for the most stable gauche and trans isomers of serotoninH+, g-1 and t1. The dihedral angles (φ1-φ3) describing the conformation of the alkyl side chain and the orientation of the hydroxyl group are listed along with selected bond lengths and relative energies. The corresponding data for all other calculated isomers are provided in Table T2 in the Supporting Information. All bond lengths of g-1 are given in Figure F2 in the Supporting Information. As mentioned above, the g-1 global minimum is significantly stabilized by the NH+-π cation-π interaction between the positively charged ammonium group and the aromatic indole ring, which involves a short intramolecular NH-C9 bond length of 2.312 (2.220) Å at the B3LYP (MP2) level. As a consequence of the substantial NH+-π interaction, the N-H bond of the proton donor of the ammonium group is substantially elongated (1.042 Å, B3LYP) as compared to the corresponding free N-H bond lengths (1.028 Å). The most stable isomer of neutral serotonin, denoted g+2(n) and shown in Figure F2 in the Supporting Information, has a structure in which one of the H atoms of the nearly neutral amino group interacts with the C2 atom of the pyrrole ring. It exhibits a much weaker NH-π interaction as compared to the g-1 isomer of serotoninH+, as evidenced by the relatively long NH-C3 bond lengths of 2.767 (2.629) Å. The IRMPD spectrum of serotoninH+ is compared in Figure 4 to the calculated linear IR absorption spectra of all considered isomers obtained at the B3LYP level. The spectra calculated at the MP2 level are similar in appearance and thus not discussed further here. This comparison reveals a single isomer to be mainly responsible for the observed IRMPD spectrum (Figure

IR Spectra of Protonated Serotonin

J. Phys. Chem. A, Vol. 114, No. 50, 2010 13273

Figure 5. Comparison of the IRMPD spectrum of serotoninH+ with the linear IR absorption spectra for the most stable gauche isomers of serotoninH+, g-1, and neutral serotonin, g+2(n), calculated at the B3LYP/cc-pVDZ level. Corresponding transitions are connected by dotted lines.

Figure 4. IRMPD spectrum of serotoninH+ and B3LYP/cc-pVDZ spectra of calculated isomers of serotoninH+ in energetic order (scaled by 0.98 with a fwhm of 30 cm-1). The calculated intensities are all on the same scale.

4), namely the most stable g-1 structure. Both the band positions and relative intensities of this isomer provide the closest match to the experimental spectrum. The other gauche isomers, namely g+1, g+2, and g-2, although close in energy and therefore potentially present in the ion cloud at room temperature, have a very intense transition predicted near 1400 cm-1, which is nearly absent in the IRMPD spectrum. This mode arises from a bending mode of the ammonium group (βNH3) slightly coupled to CH2 bending modes. Similarly, these three gauche isomers feature a relatively intense transition near 1050 cm-1, which is also not prominent in the IRMPD spectrum. Thus, significant contributions from gauche isomers other than g-1 can be excluded as major carriers of the IRMPD spectrum. This scenario is consistent with thermodynamic considerations. Assuming thermodynamic equilibrium at 300 K, the relative energies ∆E (∆G) calculated at the B3LYP level suggest a population ratio of 1:0.29:0.045:0.042 (1:0.41:0.058:0.045) for g-1, g+1, g+2, and g-2, which also yields a clear preference for the g-1 isomer (>65%). As the t2 isomer exhibits very weak activity in the 1300-1350 cm-1 range, in which strong resonances are detected in the IRMPD spectrum (bands P and Q), this structure can also be excluded as major carrier. The t1 isomer has in fact a spectrum similar to that of the g-1 isomer in the fingerprint range and can thus not be excluded by spectroscopic arguments. Nonetheless, on the basis of its high

energy (20-30 kJ/mol above g-1) and the at most minor contributions of the lower-lying g+1, g+2, and g-2 isomers, it may be excluded for thermodynamic reasons. Comparison of the IRMPD spectrum with the calculated spectra may indeed suggest minor contributions from t2, g+1, and g+2, which would explain the enhanced intensities of bands S, P, and N, respectively. However, small discrepancies between experimental IRMPD intensities and linear IR absorption cross sections may also arise from deficiencies of the computational approach and/or deviations of IRMPD signals from a linear behavior. Finally, the isomers with protonation at the indole N atom and the O atom of the hydroxyl group can be eliminated from the list of major carriers both from spectroscopic and energetic points of view, as their predicted IR spectra are clearly different from the measured IRMPD spectrum, and their relative energies are more than 100 kJ/mol above that of the g-1 isomer. The comparison between experimental and theoretical spectra (Table 1 and Figure 5) is clearly in favor of attributing the IRMPD spectrum largely to the g-1 isomer. Several bands in the IRMPD spectrum correspond to single relatively intense calculated transitions of the g-1 isomer (e.g., bands A-D, F-I, L-O, T, U, and W in Figures 4 and 5), whereas other bands are due to more than one significant vibrational mode. The maximum deviation of the positions of the experimental band maxima from the calculated frequencies is 29 cm-1, with an average deviation of 13 cm-1, confirming the vibrational and isomer assignments given in Table 1. In addition, all modes with calculated oscillator strengths larger than ∼10 km/mol are visible in the experimental spectrum. The weak bands A and B at 595 and 640 cm-1 in the IRMPD spectrum of serotoninH+ are assigned to aromatic γCC vibrations of the g-1 isomer (see footnote d of Table 1 for the notation employed for vibrational modes). The weak feature C at 678 cm-1 corresponds to a ring mode, whereas band D is again attributed to a low-intensity aromatic γCC mode. The highintensity band E is composed of five closely overlapping aromatic γCH vibrations. The feature F at 840 cm-1 is assigned to an isolated aliphatic symmetric σCCN mode. However, the stretching modes within the ethylamine side chain strongly couple with each other and with twisting, wagging and torsional

13274

J. Phys. Chem. A, Vol. 114, No. 50, 2010

modes of the ethylamine side chain. Band G at 936 cm-1 is attributed to a further isolated ring mode, whereas band H at 961 cm-1 again is an aliphatic σCCN mode of the ethylamine side chain. The intensity of band H in the IRMPD spectrum appears somewhat enhanced compared to the calculated value. Band I at 1027 cm-1 arises from an aliphatic CC stretch mode between C3 of the indole ring and C10 of the side chain. The feature K at 1066 cm-1 is composed of the two close lying CH2 twisting vibrations located at C10 and C11 of the ethylamine side chain with nearly equal weight. Band L at 1094 cm-1 is described as a coupled aromatic βNH/CH mode of the indole ring. The intense band M is assigned to the isolated δOH mode of the hydroxyl group. Band N at 1220 cm-1 is dominated by the τCH2 torsional mode located at C10 of the ethylamine side chain, with very weak contributions of an overlapping aromatic βCH vibration. Band O at 1263 cm-1 is assigned to a single aromatic βCH mode. In contrast, band P at 1304 cm-1 consists of two overlapping modes, the more intense being the σCO vibration and the less intense being the τCH2 torsional mode at C11 of the ethylamine side chain. Four different vibrational modes produce band Q at 1342 cm-1. Two of them are the βCH2 wagging modes at C10 and C11, and the third one is attributed to a further aromatic βCH vibration. However, the by far strongest contribution comes from an aromatic CC stretch (corresponding to ν14 in Wilson notation)33 predicted at 1371 cm-1 with high intensity. The somewhat larger red shift of about -30 cm-1 observed for this intense mode between the IRMPD and calculated frequency is ascribed to the IRMPD process, which shifts intense modes to lower frequency if other vibrations with significant intensity occur in the spectrum at slightly lower frequency.32 The weak shoulder R at 1417 cm-1 is probably the most prominent indicator for the assignment of g-1 to the IRMPD spectrum. As mentioned above, this band should be the most intense feature in this spectral range, in the case of significant contributions from other gauche isomers (in disagreement with the experimental observation). On the basis of the g-1 isomer, band R is assigned to two overlapping vibrations, namely a coupled σCC and σCN mode of the indole ring and the βCH2 scissoring mode at C10 of the side chain. The broad band S at 1454 cm-1 is a prominent feature in the IRMPD spectrum and also in the spectrum calculated for g-1. It consists of three almost equally intense fundamental transitions, namely, the inphase and out-of-phase βNH3 symmetric bending modes (umbrella mode) coupled with the βCH2 scissoring mode at C11, and an aromatic σCC mode located at the phenol ring (ν19a). The shoulder T at 1499 cm-1 arises from a similar aromatic σCC mode (ν19b) also located at the phenol ring. In contrast, band U at 1542 cm-1 is due to an aromatic σCC mode of the pyrrole ring. The intense band V at 1583 cm-1 is assigned to three overlapping transitions: the two asymmetric βNH3 bending modes of the ammonium group and one aromatic σCC mode located at the phenol ring (ν8b). The highest-frequency band in the measured spectral range is W at 1624 cm-1 and is attributed to a further aromatic σCC mode of the phenol ring (ν8a). Comparison of the properties of neutral serotonin with those of serotoninH+ establishes the effects of protonation on its geometric and electronic structure. The calculations predict an energetic preference for gauche isomers for both neutral and protonated serotonin, because they are stabilized through the intramolecular NH(+)-π interaction with the aromatic indole ring. The major difference is that the neutral species prefers NH-π bonding to the pyrrole ring, whereas the protonated species interacts more effectively via NH+-π bonding with the phenol ring (Figure F2 in the Supporting Information). In the

Lagutschenkov et al. condensed phase, the preferential configuration of (protonated) serotonin changes drastically due to the effects of the environment. The interaction with solvent molecules and counterions is stronger than the intramolecular NH-π interaction and leads for both species to a preferential stabilization of trans conformers via intermolecular NH hydrogen bonding, although the gauche isomers are clearly calculated to be the global minima on the potential energy surface of the isolated species. Interestingly, the present calculations at the B3LYP and MP2 level using the cc-pVDZ basis set yield a global minimum structure for neutral serotonin, denoted g-2(n), which slightly differs from the global minimum derived previously at lower level calculations.10,12 Calculations at the B3LYP/6-31G* and MP2/6-31G* yield a global minimum, denoted Gpy(out), which is 0.7 or 2.5 kJ/mol more stable than the Gpy(up) conformation, which corresponds to g-2(n). At the B3LYP/cc-pVDZ and MP2/cc-pVDZ level, the g-2(n) isomer is lower in energy than the Gpy(out) structure by 2.7 and 0.5 kJ/mol, respectively. These subtle differences illustrate the difficulties in predicting the energetic order of close lying isomers of these flexible molecules. Experimental evidence for the relative stability of the neutral serotonin isomers has been obtained from fluorescence and ionization yields,12 which may, however, slightly depend on the conformation and ionization efficiency of the molecule and thus do not provide quantitative information about their stability. In the following, we compare the properties of the most stable gauche structures of isolated serotonin(H+), g-1 and g+2(n), as obtained by the present quantum chemical calculations. Detailed structural parameters for both species are given in Figure F2 in the Supporting Information. The energy difference between both structures corresponds to the proton affinity of 967.6 (972.85) kJ/mol at the B3LYP (MP2) level. No experimental value appears to be available for comparison. As already mentioned, protonation drastically enhances the strength of the intramolecular NH-π interaction leading to substantially shorter contacts between the NH proton and the aromatic carbon atoms. Protonation at the N-terminus also leads to a small average elongation of the N-H bonds and simultaneously to an elongation of the neighboring N-C bond. All other bond length changes are less significant. The structural changes induced by protonation of serotonin translate directly into the vibrational properties and the corresponding IR spectra, which are compared in Figure 5 in the fingerprint range. Clearly, protonation has a profound effect on both the positions and the IR intensities of the vibrational modes. The vibrational frequencies and IR intensities of g+2(n) and g-1 obtained at the B3LYP level are compared in Table T3 in the Supporting Information, along with the assignment of the normal modes. As expected, the N-H bending modes experience the largest impact upon protonation in the frequency range investigated. All three N-H bend fundamentals have much larger IR intensities for the protonated species (band V) due to the large positive partial charge of the NH3+ group. In particular, the intense symmetric N-H umbrella modes (here split into two modes due to coupling with CH2 scissoring modes) at 1439 and 1445 cm-1 are characteristic for the charged g-1 species (band S), whereas the IR spectrum of the neutral g+2(n) molecule has no intense absorption in this spectral range. On the other hand, the large intensity of the aromatic C-C stretch mode of g+2(n) at 1499 cm-1 (ν19b, 103 km/mol) is largely reduced upon protonation (1508 cm-1, 19 km/mol, band T). As there is no IR spectrum of isolated serotonin available in the literature in the fingerprint range, comparison of the IR spectrum

IR Spectra of Protonated Serotonin calculated for g+2(n) with experiment is not possible at the present stage. The NBO population analysis for g+2(n) and g-1, is detailed in Figure F3 in the Supporting Information. As expected, the ethylamine side chain carries only little charge (0.03 e, B3LYP) in neutral serotonin and nearly the total positive charge (+0.97 e) in serotoninH+. The latter one is mainly localized on the ammonium group (+0.62 e) and to lesser extent in the adjacent CH2 units (+0.30 and +0.05 e). The charge on the aromatic ring is not affected upon protonation, consistent with the lack of hyperconjugation in aromatic alkane derivatives. The large positive partial charge on the NH3+ group gives rise to the substantial charge-enhancement of the NH+-π interaction in the g-1 isomer of serotoninH+ as compared to the neutral molecule. Finally, the fragmentation process of serotoninH+ is considered. The observed IRMPD fragment of serotoninH+ is m ) 160 u, corresponding to the loss of NH3. In contrast to dopamineH+,28 for which fragmentation upon IRMPD occurs in a sequential fashion into various fragment channels, IRMPD of serotoninH+ yields only one single fragment. Predominant loss of NH3 upon collisional activation of serotoninH+ has previously been observed by other groups.20 Interestingly, Chang and Yeung report the elimination of the methylamine group (CH2NH2) to be the dominant fragmentation channel (m ) 146 u).34 However, it is likely that, due to their insufficient mass resolution, they were indeed observing the fragmentation of the serotonin+ radical cation rather than serotoninH+, which was shown to eliminate CH2NHn with n ) 1 and 2.35 A variety of possible structures for the m ) 160 u fragment ion observed upon IRMPD of g-1 is shown in Figure F4 of the Supporting Information, along with their relative energies and dissociation energies for NH3 elimination. However, it is difficult at the present stage to identify the actual structure of the fragment ion observed in the experiment. 4. Concluding Remarks The conformation and intramolecular cation-π interactions of isolated serotoninH+ were investigated by IRMPD spectroscopy and quantum chemical calculations. Comparison of the linear IR absorption spectra calculated for various gauche and trans conformers of serotoninH+ and the experimental IRMPD spectrum in the fingerprint region yields the best agreement for the gauche conformer g-1, which is calculated to be the lowestenergy isomer at the B3LYP and MP2 level. In the g-1 isomer, protonation occurs at the N-terminus of the ethylamine side chain, allowing for an efficient intramolecular NH+-π cation-π interaction of the ammonium group and the phenol moiety of the aromatic indole ring. Other isomers are higher in energy and have predicted IR spectra, which differ from the observed IRMPD spectrum, suggesting that g-1 is the major carrier of the experimental spectrum. Although the cation-π interaction in the gauche isomers of serotoninH+ are largely governed by electrostatic and inductive contributions arising from the substantial positive charge localized mainly at the NH3+ group (∼20-30 kJ/mol), dispersion contributions to the NH+-π interaction (∼7 kJ/mol) are inferred from the different stabilization energies obtained at the B3LYP and MP2 energies. Interestingly, this difference is only ∼5 kJ/mol for dopamineH+, suggesting that the dispersion forces are larger in serotoninH+ due to the interaction with the more extended aromatic π-electron system. The major differences between serotonin and its protonated form are that the excess charge strongly enhances the NH(+)-π interaction and that protonation changes the

J. Phys. Chem. A, Vol. 114, No. 50, 2010 13275 preferred binding motif from NH-π bonding to the pyrrole ring to NH+-π bonding to the phenol ring. Clearly, further experimental information about the strength of the NH+-π interaction could come from IR spectra of serotoninH+ recorded in the NH stretch range using optical parametric oscillator lasers.36 Acknowledgment. This work was supported by Technische Universita¨t Berlin and Deutsche Forschungsgemeinschaft (DO 729/3). The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under Grant Agreement No. 226716. This work is part of the research program of FOM, which is financially supported by the Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO). We gratefully acknowledge the excellent assistance by the FELIX staff (B. Redlich et al.). Supporting Information Available: Geometries, energies, charge distributions, fragment ion structures, and vibrational properties of various serotoninH+ isomers and neutral serotonin evaluated at various levels of theory. This material is available free of charge via the Internet at http://pubs.acs.org. References and Notes (1) Mazak, K.; Doczy, V.; Kokosi, J.; Noszal, B. Chem. BiodiVersity 2009, 6 (4), 578. Frazer, A.; Hensler, J. G. In Basic Neurochemistry: Molecular, Cellular and Medical Aspects; Siegel, G. J., Agranoff, B. W., Albers, R. W., Fisher, S. K., Uhler, M. D., Eds.; Lippincott-Raven: Philadelphia, PA, 1999. (b) Hoyer, D.; Hannon, J. P.; Martin, G. R. Pharmacol., Biochem. BehaV. 2002, 71 (4), 533. (2) Decloedt, E. H.; Stein, D. J. Neuropsychiatr. Dis. Treat. 2010, 6, 233. (3) Soloff, P. H.; Price, J. C.; Meltzer, C. C.; Fabio, A.; Frank, G. K.; Kaye, W. H. Biol. Psychiatry 2007, 62 (6), 580. (4) Paterson, D. S.; Trachtenberg, F. L.; Thompson, E. G.; Belliveau, R. A.; Beggs, A. H.; Darnall, R.; Chadwick, A. E.; Krous, H. F.; Kinney, H. C. J. Am. Med. Assoc. 2006, 296 (17), 2124. (5) Pratuangdejkul, J.; Schneider, B.; Launay, J. M.; Kellermann, O.; Manivet, P. Curr. Med. Chem. 2008, 15 (30), 3214. (6) Chattopadhyay, A.; Rukmini, R.; Mukherjee, S. Biophys. J. 1996, 71 (4), 1952. (7) Beene, D. L.; Brandt, G. S.; Zhong, W. G.; Zacharias, N. M.; Lester, H. A.; Dougherty, D. A. Biochemistry 2002, 41 (32), 10262. (8) Pisterzi, L. F.; Almeida, D. R. P.; Chass, G. A.; Torday, L. L.; Papp, J. G.; Varro, A.; Csizmadia, I. G. Chem. Phys. Lett. 2002, 365 (56), 542. (9) Mazurek, A. P.; Weinstein, H.; Osman, R.; Topiol, S.; Ebersole, B. J. Int. J. Quantum Chem. 1984, 183. (10) van Mourik, T.; Emson, L. E. V. Phys. Chem. Chem. Phys. 2002, 4 (23), 5863. (11) Bayari, S.; Saglam, S.; Ustundag, H. F. J. Mol. Struct.-THEOCHEM 2005, 726 (1-3), 225. (12) LeGreve, T. A.; Baquero, E. E.; Zwier, T. S. J. Am. Chem. Soc. 2007, 129 (13), 4028. (13) Alagona, G.; Ghio, C.; Nagy, P. I. J. Chem. Theory Comput. 2005, 1 (5), 801. (14) Pratuangdejkul, J.; Jaudon, P.; Ducrocq, C.; Nosoongnoen, W.; Guerin, G. A.; Conti, M.; Loric, S.; Launay, J. M.; Manivet, P. J. Chem. Theory Comput. 2006, 2 (3), 746. (15) Alagona, G.; Ghio, C. J. Mol. Struct.-THEOCHEM 2006, 769 (13), 123. (16) LeGreve, T. A.; Clarkson, J. R.; Zwier, T. S. J. Phys. Chem. A 2008, 112 (17), 3911. LeGreve, T. A.; James, W. H.; Zwier, T. S. J. Phys. Chem. A 2009, 113 (2), 399. (17) Ison, R. R.; Roberts, G. C. K.; Partington, P. J. Pharm. Pharmacol. 1972, 24 (1), 82. (18) Bugg, C. E.; Thewalt, U. Science 1970, 170 (3960), 852. Amit, A.; Mester, L.; Klewe, B.; Furberg, S. Acta Chem. Scand., Ser. A: Phys. Inorg. Chem. 1978, 32 (3), 267. Karle, L.; Dragonet, Ks.; Brenner, S. A. Acta Crystallogr. 1965, 19, 713. (19) Thewalt, U.; Bugg, C. E. Acta Crystallogr., Sect. B: Struct. Crystallogr. Cryst. Chem. 1972, 28, 82. (20) Thomson, B. A.; Iribarne, J. V.; Dziedzic, P. J. Anal. Chem. 1982, 54 (13), 2219. Steiner, W. E.; Clowers, B. H.; English, W. A.; Hill, H. H.

13276

J. Phys. Chem. A, Vol. 114, No. 50, 2010

Rapid Commun. Mass Spectrom. 2004, 18 (8), 882. Bourcier, S.; Hoppilliard, Y. Rapid Commun. Mass Spectrom. 2009, 23 (1), 93. (21) Bourcier, S.; Benoist, J. F.; Clerc, F.; Rigal, O.; Taghi, M.; Hoppilliard, Y. Rapid Commun. Mass Spectrom. 2006, 20 (9), 1405. (22) Macleod, N. A.; Simons, J. P. Mol. Phys. 2006, 104 (20-21), 3317. Macleod, N. A.; Simons, J. P. Phys. Chem. Chem. Phys. 2004, 6, 2821. Vaden, T. D.; de Boer, T.; MacLeod, N. A.; Marzluff, E. M.; Simons, J. P.; Snoek, L. C. Phys. Chem. Chem. Phys. 2007, 9 (20), 2549. (23) Eyler, J. R. Mass Spectrom. ReV. 2009, 28 (3), 448. Polfer, N. C.; Oomens, J. Mass Spectrom. ReV. 2009, 28 (3), 468. Knorke, H.; Langer, J.; Oomens, J.; Dopfer, O. Astrophys. J. Lett. 2009, 706, L66. Seydou, M.; Gregoire, G.; Liquier, J.; Lemaire, J.; Schermann, J. P.; Desfrancois, C. J. Am. Chem. Soc. 2008, 130 (12), 4187. Fridgen, T. D. Mass Spectrom. ReV. 2009, 28 (4), 586. (24) Zhao, D. W.; Langer, J.; Oomens, J.; Dopfer, O. J. Chem. Phys. 2009, (18), 131. (25) Rizzo, T. R.; Stearns, J. A.; Boyarkin, O. V. Int. ReV. Phys. Chem. 2009, 28 (3), 481. Nagornova, N. S.; Rizzo, T. R.; Boyarkin, O. V. J. Am. Chem. Soc. 2010, 132 (12), 4040. Fujihara, A.; Matsumoto, H.; Shibata, Y.; Ishikawa, H.; Fuke, K. J. Phys. Chem. A 2008, 112 (7), 1457. Guidi, M.; Lorenz, U. J.; Papadopoulos, G.; Boyarkin, O. V.; Rizzo, T. R. J. Phys. Chem. A 2009, 113 (5), 797. Stearns, J. A.; Guidi, M.; Boyarkin, O. V.; Rizzo, T. R. J. Chem. Phys. 2007, 127 (15), 154322. (26) MacAleese, L.; Maitre, P. Mass Spectrom. ReV. 2007, 26 (4), 583. Dopfer, O. J. Phys. Org. Chem. 2006, 19, 540. Lagutschenkov, A.; Springer, A.; Lorenz, U. J.; Maitre, P.; Dopfer, O. J. Phys. Chem. A 2010, 114 (5), 2073. Lagutschenkov, A.; Sinha, R. K.; Maitre, P.; Dopfer, O. J. Phys. Chem. A 2010, 114 (42), 11053. (27) O. Dopfer, Report for Project FELIX-003 (2009). (28) Lagutschenkov, A.; Langer, J.; Berden, G.; Oomens, J.; Dopfer, O. Phys. Chem. Chem. Phys., in press. (29) Valle, J. J.; Eyler, J. R.; Oomens, J.; Moore, D. T.; van der Meer, A. F. G.; von Helden, G.; Meijer, G.; Hendrickson, C. L.; Marshall, A. G.;

Lagutschenkov et al. Blakney, G. T. ReV. Sci. Instrum. 2005, 76 (2), 7. Oepts, D.; Vandermeer, A. F. G.; Vanamersfoort, P. W. Infrared Phys. Technol. 1995, 36 (1), 297. (30) Dunning, T. H. J. Chem. Phys. 1989, 90 (2), 1007. (31) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Montgomery, J. A., Jr.; Vreven, T.; Kudin, K. N.; Burant, J. C.; Millam, J. M.; Iyengar, S. S.; Tomasi, J.; Barone, V.; Mennucci, B.; Cossi, M.; Scalmani, G.; Rega, N.; Petersson, G. A.; Nakatsuji, H.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Klene, M.; Li, X.; Knox, J. E.; Hratchian, H. P.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Ayala, P. Y.; Morokuma, K.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Zakrzewski, V. G.; Dapprich, S.; Daniels, A. D.; Strain, M. C.; Farkas, O.; Malick, D. K.; Rabuck, A. D.; Raghavachari, K.; Foresman, J. B.; Ortiz, J. V.; Cui, Q.; Baboul, A. G.; Clifford, S.; Cioslowski, J.; Stefanov, B. B.; Liu, G.; Liashenko, A.; Piskorz, P.; Komaromi, I.; Martin, R. L.; Fox, D. J.; Keith, T.; Al.Laham, M. A.; Peng, C. Y.; Nanayakkara, A.; Challacombe, M.; Gill, P. M. W.; Johnson, B.; Chen, W.; Wong, M. W.; Gonzalez, C.; Pople, J. A. Gaussian 03, revision C.02; Gaussian, Inc.: Wallingford, CT, 2004. (32) Oomens, J.; Sartakov, B. G.; Meijer, G.; Von Helden, G. Int. J. Mass Spectrom. 2006, 254 (1-2), 1. (33) Wilson, E. B. Phys. ReV. 1934, 45 (10), 706. (34) Chang, S. Y.; Yeung, E. S. Anal. Chem. 1997, 69 (13), 2251. (35) Cardoso, A. M.; Correla, A. J. F. Eur. J. Mass Spectrom. 1999, 5 (1), 11. (36) Pasker, F.; Solca, N.; Dopfer, O. J. Phys. Chem. A 2006, 110, 12793. Andrei, H.-S.; Solca, N.; Dopfer, O. ChemPhysChem 2006, 7, 107. Andrei, H.-S.; Solca, N.; Dopfer, O. J. Phys. Chem. A 2005, 109, 3598. Dopfer, O. Z. Phys. Chem. 2005, 219, 125. Solca, N.; Dopfer, O. Phys. Chem. Chem. Phys. 2004, 6, 2732. Solca, N.; Dopfer, O. Eur. Phys. J. D 2002, 20, 469. Bieske, E. J.; Dopfer, O. Chem. ReV. 2000, 100, 3963.

JP109337A