1
The dipeptide cyclic(glycyltryptophanyl) in the gas phase: A concerted action of density functional calculations, S0–S1 two-photon ionization, spectral UV/UV hole burning and laser photoelectron spectroscopy

2
The S0–S1 spectrum of the dipeptide cyclic(glycyltryptophanyl) (cGW) has been recorded by two-photon ionization (R2PI) spectroscopy of thermally evaporated, jet-cooled molecules.

3
The R2PI spectrum contains several vibronic transitions with relatively small spacing.

4
Density functional theory predicts more than 8 conformers.

5
Applying spectral UV/UV hole burning spectroscopy we found that the strongest S0–S1 transitions belong to a single conformer.

6
The transition at 35 058 cm−1 is attributed to the S1 origin and the other transitions to inter-ring modes.

7
To obtain further information on the conformer structure, resonant two-colour two-photon ionization photoelectron spectroscopy (R(1 + 1′)PI PES) has been performed via the two most intense S0–S1 transitions.

8
To our knowledge these are the first PE spectra of a dipeptide.

9
The spectra are broad with a smooth onset.

10
The lowest ionization onset lies at 7.709 eV, but is assumed to be not the adiabatic IE.

11
Theoretical calculations with the B3LYP/6-311++G** theory predict that in the most stable neutral conformer an N–H group of the dipeptide ring binds to the indole π system.

12
This structure has the highest IE of all conformers investigated and shows a strong geometry change upon ionization.

13
The fact that the IE of the experimentally observed cGW conformer is higher than that of 3-methylindole is taken as a signature of a repulsion between the two rings in the cation: This could be explained by the N–H group of the peptide ring now interacting with the positively charged indole π system.

14
The conformational assignment by R(1 + 1′)PI PES is tentative, because all of the lowest-energetic neutral conformers are expected to behave similarly in regard to ionization energies.

Introduction

15
The advantage of gas phase measurements is the applicability of spectroscopy methods such as resonant two-photon ionization (R2PI) spectroscopy combined with mass-selective detection, spectral UV/UV or IR/UV hole burning (UV/UV SHB or IR/UV SHB) spectroscopy and resonant two-colour two-photon ionization photoelectron spectroscopy (R(1 + 1′)PI PES).

16
Further advantages of gas phase investigations are usually sharp spectra which reveal detailed structural information.

17
Most importantly, however, gas phase data can be directly compared to theory.

18
For reviews on biological molecules in the gas phase see .ref. 1–4

19
Levy and co-workers were the first to perform spectroscopy of di- and tripeptides in the gas phase by using laser desorption with R2PI and mass-selective detection.5,6

20
For tryptophan-containing peptides they obtained sharp spectra with distinct lines despite the many vibrational degrees of freedom and the many possible conformers.

21
The spectroscopy of laser-desorbed small peptides was later improved by using double resonance methods.7,8

22
Neutral laser desorption and R2P ionization and laser fragmentation was also introduced as a new species selective inlet system into mass spectrometry.9–11

23
Thermal evaporation of natural unprotected peptides is impossible.12,13

24
Up to now only for some modified peptides in which the polar functional groups have been protected by methyl groups a thermal evaporation of intact molecules is reported.14,15

25
Cyclic dipeptides have no polar N or C terminal and hence, as we show here, can be vaporized intact by heating.

26
To our knowledge up to now no investigation of cyclic peptides in the gas phase has been performed.

27
UV/UV or IR/UV double resonance methods are widely used in spectroscopy of floppy molecules and clusters.2

28
The strategy to obtain structural information by using zero kinetic energy electron PES (ZEKE16), mass analyzed threshold ionization (MATI17) or high resolution conventional laser PES18 is not new.18–25

29
In these publications model systems, such as extended molecules with a chromophore and a flexible side chain, have been investigated.

30
Whereas for non-polar side groups sharp PE spectra are observed exclusively18–20 for polar groups sharp spectra are only expected for conformers where the interaction of the charge site and the polar group are minimized, i.e. unfolded, linear structures.18

31
In folded structures the polar group solvates the charge strongly and structural changes upon ionization become so large that it is impossible to reach the origin of the cation ground state and obtain vibrationally resolved spectra.22

32
Nevertheless, the IE shift together with theory may provide useful information to allow an identification of the neutral conformation.22

33
Only in one case of an aromatic molecule with a polar side chain, both folded and unfolded conformers showed sharp vibrationally resolved spectra.24,25

34
In this paper we report on spectroscopy of the cyclic(glycyltryptophanyl) (cGW) in the gas phase.

35
R2PI spectroscopy with mass selection has been used to record the S0–S1 spectrum.

36
UV/UV SHB was applied to analyze the number of conformers present in the molecular beam.

37
To obtain further information on the conformer structure, R(1 + 1′)PI PES was performed for the two most intense S0–S1 transitions.

38
Density functional theory (DFT) calculations on the B3LYP/6-311++G** level are used to interpret the number and structure of observed conformers, IEs and geometry shifts upon ionization.

39
A tentative conformer identification is given on the basis of a comparison between experiment and theory.

Experimental

40
The cyclic peptide cGW was purchased from Bachem.

41
It has been synthesized from pure L-tryptophan (W) and hence cGW has also L chirality.

42
The apparatus is described in detail elsewhere.18

43
Briefly it consists (Fig. 1) of a supersonic jet chamber (a), a field-free time-of-flight PE spectrometer (b) and a pulsed linear time-of-flight mass spectrometer (c).

44
A modified pulsed nozzle (General Valve, orifice diameter 300 μm) was heated to achieve a suitable sample vapor pressure.

45
The indole sample has been added to a sample container above the nozzle and heated to 325 K. The cGW sample, crystalline at room temperature, is placed directly in the nozzle head and heated to 525 K. Ar was used as carrier gas with a backpressure of 5 bar.

46
The jet passed the skimmer (1.5 mm orifice) and crossed the laser beam at right angles in the middle of the PE spectrometer.

47
The time-of-flight PE spectrometer is 900 mm long, μ-metal shielded and covered inside by graphite.

48
Only photoelectrons emitted into the small acceptance angle of the two-stage micro channel plate (MCP) detector were detected and their time measured with high precision in a start–stop logic.

49
Note that for photoelectron spectroscopy no ion extraction plates have been used within the PE spectrometer.

50
After ionization, the ions drift downwards with their supersonic jet velocity into a pulsed linear time-of-flight (linTOF) mass spectrometer (mass resolution mm × 100).

51
Hence R2PI mass spectra, mass-analyzed R2PI spectra and R(1 + 1′)PI PE spectra with simultaneous mass control can be recorded in the same apparatus.

52
Only for UV/UV SHB experiments three deflection plates were introduced at the ionization site in the PE spectrometer.

53
In this case the ions were formed in an electric field and deflected into the former electron detector.

54
The laser system consisted of two Nd:YAG laser-pumped and frequency-doubled tunable dye lasers (5–7 ns pulse width) and an additional excimer-pumped frequency-doubled dye laser which can be arbitrarily shifted in time with respect to the Nd:YAG-pumped dye lasers.

R2PI spectroscopy

55
For recording the R2PI spectrum ions were formed in the PE spectrometer and, after drifting down, mass-analyzed in the pulsed linear TOF mass spectrometer.

56
The S0–S1 spectrum was recorded by tuning the laser wavelength and recording the parent cations with a gated integrator.

57
For this purpose we used the Nd:YAG-pumped frequency doubled dye laser (resolution in the UV: < 1 cm−1, 100 μJ pulse−1, 1000 mm lens out of focus).

58
We found that the cooling efficiency of the supersonic jet very much depended on the expansion gas and improved with increasing mass of the carrier gas.

UV/UV SHB

59
For the SHB a three-plate ion optics (repeller plate, intermediate plate, grounded exit plate) was mounted in the PE spectrometer which deflected and accelerated the ions to the electron detector of the PE spectrometer.

60
By changing the MCP voltages suitable for ion detection mass and UV/UV SHB spectra have been measured.

61
The mass resolution mm of this linTOF is 1500.

62
This setup is needed to deflect the ions formed by the hole burning laser by using a pulsed electric field scheme.

63
The hole burning laser performed ground state depletion by R2PI 4 μs before the probe laser ionizes the remaining neutrals.

64
For ion discrimination coming from different lasers, ions stemming from the hole burning laser were for 4 μs accelerated backwards (+500 V at the second plate, repeller plate grounded) and neutralized on the repeller plate.

65
A suitable pulsed voltage (750 V) was then applied to the repeller shortly after new ions had been formed by the probe laser.

66
The SHB spectrum was recorded by scanning the hole burning laser while the probe laser was set to the S0–S1 origin transition of cGW.

67
Each time when the neutral ground state of the cGW conformer was depopulated a dip appeared in the probe laser signal.

68
Note that in this mode the dip depth directly correlates with the transition moment of the transition.

69
Also a reversed scheme has been used to test the percentage of depletion for each transition of interest.

(R(1 + 1′)PI) PE spectroscopy

70
For R(1 + 1′)PI PE spectroscopy ionization was performed by two frequency-doubled dye lasers pumped by the same Nd:YAG laser.

71
The ionizing laser (laser 2) was time-delayed by 3 ns against laser 1 which performed the S0–S1 excitation.

72
The photon energy of laser 2 was not sufficient for a two-photon ionization.

73
The one-photon signal of laser 1 was lower than 5% of the total two-color ion signal.

74
Typical laser pulse energies were 10 μJ for laser 1 and 20 μJ for laser 2 (1000 mm lens, out of focus).

75
To optimize the electron energy resolution, the photon energy of laser 2 was chosen to keep the maximum photon excess energy below 200 meV.

76
For indole we achieved a PE energy resolution ΔE = 1 meV and an absolute accuracy of ±2 meV.

77
The accuracy for cGW was ±3 meV due to the smooth onset at the ionization onset.

Theoretical results

78
Our motivation to perform calculations on cGW was (i) to obtain the relative energies of the neutral and cation conformers, (ii) to calculate the vibrational frequencies for the lowest conformers, (iii) to calculate the vertical IEs and (iv) to explore the structural differences between the neutral and cation for the most stable conformers.

79
For conformer assignment we use the comparison of the experimental to calculated data, such as relative neutral conformer energies, vertical IEs and neutral to cation geometry shifts.

80
So far no ab initio or DFT calculations of a cyclic peptide are reported.

81
We used the density functional B3LYP26 because according to our experience18 it seems to better describe flexible molecules with polar groups than other methods.

82
It is well known that B3LYP systematically underestimates intermolecular electron correlation forces as present in folded molecules.

83
There could be several explanations.

84
One would be an accidental error compensation, another that electrostatic forces dominate the binding.

85
The greatest advantage of DFT methods is that they are fast enough to allow the use of big basis sets even for large molecules.

86
This is very important for folded molecules because it minimizes the basis set superposition error usually favoring folded molecules energetically.

87
DFT calculations have been performed with the Gaussian 98 program.27

88
The basis sets 3-21G, 6-31+G* and 6-311++G** were used.

89
The large basis set 6-311++G** was applied for neutral cGW to especially account for diffuse electron distributions.

90
By increasing the basis set from 3-21G to 6-31+G* we observe big changes in the relative conformer energies.

91
A saturation effect occurs when increasing the basis set from 6-31+G* to 6-311++G** (Table 1) which might be due to accidental error compensation for the smaller basis set.

92
We also calculated vibrational frequencies (B3LYP/6-31+G*) for the neutral conformers 1, 2 and 3 which are needed for calculating the Gibbs free energy and also used for a tentative vibrational assignment in section IV.2.

The conformers of neutral cGW

93
To find the local minima the conformational landscape of cGW has been explored using the program Tinker.28

94
It performs a general conformational search for the entire potential energy surface by using a hopping algorithm29 which is able to overcome barriers.

95
The program Tinker also minimizes each conformer after its generation.

96
The chosen MM3 force field30 is known to be quite successful in reproducing molecular IR spectra.

97
In order to make sure we do not miss a conformer, in addition also a potential energy surface scan was performed using the OPLS all atom force field.31

98
The results from the both scans were collected and the conformers processed with B3LYP by stepwise increasing the basis set size.

99
For discussions of the type of conformers of cGW in Fig. 2a the relevant rotational axes and the resultant structural possibilities are shown.

100
The conformers of cGW can be classified as follows: The peptide ring can be moved into three positions when rotating about the C1–C2 axis (Fig. 2a, right side, see also atomic numbering C0–C2).

101
For each of these conformers the indole ring then can be switched into six positions by rotation around the C0–C1 axis (Fig. 2b, left side).

102
This gives theoretically 18 conformers.

103
The most stable conformers are formed by rotating the C1–C2 axis into the plane which is perpendicular to the indole ring (perpendicular conformer class).

104
The other tilted conformers are sterically unfavoured or unstable.

105
By our calculations 8 stable conformers have been found for neutral cGW.

106
The calculated relative energies with respect to the lowest conformer are presented in Table 1.

107
A full geometry optimization of the neutral conformers was performed by using the B3LYP/6-311++G** method.

108
The calculations show that the peptide ring (cyclic(glycylglycyl) or 2,5-piperazinedione) is not planar resulting in a double well potential (see Fig. 2b).

109
In the unperturbed cyclic(glycylglycyl) the minima are equal in energy and the barrier is expected to be low because of the high flexibility between the amide groups described in biochemical text books.

110
In cGW this double well is superimposed with a ring-to-ring interaction potential (schematic drawing in Fig. 2b).

111
Depending on the size and steepness of this ring-to-ring interaction potential, for each of the above-mentioned conformers either two sub-conformers a and b (this is the case for conformers 1 and 3), or only one conformer (a or b for conformers 2, 4, 5, 6, 7, 8) exist.

112
The peptide-ring sub-conformers are very similar in energy (energy difference 41 cm−1 for conformers 1a and 1b, and 18 cm−1 for conformers 3a and 3b, see Table 1).

113
The calculated isomerization barrier (B3LYP/6-31+G*) from conformer 1b to conformer 1a is 2 cm−1 and already the zero point vibration would be above the barrier.

114
The same is expected for conformers 3a and 3b for which the energetic difference is even smaller.

115
The accuracy of the reaction barrier calculation is not known.

116
However, even if the isomerization barrier between the a and the b forms would be larger than 2 cm−1 but lower than 300 cm−1, the more unstable form would be converted to the lower conformer in the supersonic expansion.32

117
In the following we assume that conformers 1b and 3b are experimentally not observable.

118
It is evident, however, that the resulting broad and anharmonic potential strongly reduces the vibrational frequencies in this coordinate, increases the vibrational density and by this favors these conformers on the Gibbs free energy landscape.

119
In Table 1 the results of the B3LYP density functional calculations of the 8 stable neutral conformers of the dipeptide cGW (see also a and b conformers) using the 3-21G, 6-31+G* and 6-311++G** basis sets are presented.

120
Relative conformer energies are given in cm−1.

121
The four lowest structures are shown in Fig. 3.

122
With the increase of the basis set from 3-21G to 6-31+G* the energetic order of some conformers changes.

123
When increasing the basis set from 6-31+G* to 6-311++G** the conformer order is conserved and the changes in the relative conformer energies become small, which might be due to an error compensation for the smaller basis set.

124
As discussed above, structures 1b and 3b are not stable by zero point vibration.

125
Conformer 2 which is characterized by a C–H⋯π interaction has a relative energy of 194 cm−1.

126
Conformer 3a which is similar to conformer 1a but with the indole ring switched by 180° lies 355 cm−1 above 1a.

127
In this conformer the N–H group seems to have better access to the indole π system.

128
The fact that it is less stable than conformer 1a is attributed to a positive partial charge on the carbon in position 4.

129
Only conformers with relative energies well below 365 cm−1 (∼ kT at 525 K) are expected to have measurable populations in the beam.

130
This also holds for conformer 2 but not for conformer 3.

131
To calculate the exact thermal population beside the energy difference of the electronic energies also the vibrational zero point energy and the vibrational density of states has to be considered.

132
Assuming harmonic potentials all this would be included in the Gibbs free energy.

133
For a temperature of 300 K it shows an increase of the energetic spacing between conformer 1 and 2 from 355 cm−1 to 675 cm−1.

134
This energy difference is expected to further increase when the strong anharmonicity of the inversion coordinate in conformer 1 would be taken into account.

135
Hence from theory we predict that a great entropy effect is forcing the dipeptide population into conformer 1.

136
The high energy of conformer 3 (355 cm−1) above conformer 1 should result in a marginal or vanishing population of conformer 3.

137
Comparing conformer 1 and 3 the entropy argument is not helpful because also conformer 3 has a broad double minimum in the peptide ring inversion coordinate.

138
It is well known that barriers lower than 300 cm−1 can be overcome during the supersonic jet expansion.32

139
To compare theoretical data to the experimental data it would be important to calculate the barrier heights between conformers 2 and 1 as well as 3 and 1.

140
Unfortunately, because of the complexity of the coordinates and the size of the system, it is impossible to calculate such isomerization barriers.

The cGW cation

141
We calculated vertical (6-311++G** basis set) and adiabatic (3-21G and 6-31+G* basis sets) IEs for comparison with the PE results in the hope to be able to assign the structure of the observed conformer.

142
The cation geometry optimization was performed with the 3-21G basis set for all 8 conformers, and with the 6-31+G* basis set for seven conformers.

143
The B3LYP/6-31+G* calculation for conformer 7 did not converge presumably because the potential surface is very flat.

144
We checked the dependence of the cation structural changes on the basis set, but found that the 3-21G structures are very similar to the structures found with the 6-31+G* basis set.

145
The calculated data are collected in Table 2 and the comparison of neutral and cation structures is shown in Fig. 4.

146
We also calculated the vertical IEs of the b conformers of structures 1 and 3, although it is clear that they are not populated due to their low isomerization barriers (see above).

147
The differences between the vertical IEs of conformers a and b nevertheless give hints as to how steep the ion ground state potential at the vertical geometry might be.

148
When optimizing the geometry in the cation ground state the sub-conformers a and b formed by the peptide ring inversion converge for conformers 1a and 1b as well as 3a and 3b to a single minimum.

149
This is explained by a large repulsive peptide ring-to-indole ring interaction in the cation which forces the system into one minimum.

150
A calculation of the geometrical average of the changes of all atomic coordinates between neutral and cation (termed RMS value in Table 2) shows that, except one conformer, all conformers undergo large structural changes upon ionization.

151
Fig. 4 shows the neutral-to-cation geometry changes of conformers 1, 2 and 6.

152
For this the neutral and the cation structure are overlaid with the indole ring and conformational changes indicated by arrows and lines.

153
In all cases the structural changes are so large that it might be difficult to reach the adiabatic IE.

154
For conformers 1, 2 (Fig. 4a,b) and 3 (not shown) the large geometry change is explained by the difference between the attractive (neutral) versus the repulsive (cation) N–H⋯π (conformers 1 and 3) or C–H⋯π (conformer 2) interaction.

155
In Fig. 4c the neutral and the cationic structures of conformer 6, the most stable cation conformer, are compared.

156
In this conformer the CO group (note the negative partial charge on O) has access to the indole chromophore.

157
Whereas conformer 6 is rather unstable in the neutral due to the repulsive CO⋯π interaction it becomes the most stable conformer when the π system is positively charged.

Experimental results and discussion

158
In the following we present R2PI mass, R2PI S0–S1, UV/UV SHB and R(1 + 1′)PI PE spectra of the dipeptide cGW.

The R2PI mass spectrum

159
For recording the R2PI mass spectrum of cGW (Fig. 5) the second harmonic of a Nd:YAG laser-pumped dye laser was tuned to the S0–S1 origin transition.

160
The mass spectrum was recorded by using the PE spectrometer as flight tube for the linTOF.

161
Mass 234 Da corresponds to the parent radical cation.

162
Two fragments with masses 131 Da and 117 Da are observed.

163
These masses correspond to the neutral fragments 3-methylindole (131 Da) and indole (117 Da) which are then ionized by the laser pulse.

164
In the case of a cationic fragmentation of tryptophan instead of mass 117 Da mass 130 Da would be expected as the most intense fragment.33

165
Measurements with different laser intensities and wavelengths confirm that both masses are neutral fragments containing the indole chromophore.

166
The ionization of the neutral fragments is a minor contribution in the mass spectrum at low photon energies but increases when coming closer to the indole S0–S1 resonance.

167
Because of the mass selection there should be no influence of the fragment ionization on the R2PI spectrum of cGW.

168
Nevertheless we find a small peak in the cGW spectrum at exactly the energy of the S0–S1 origin transition of indole which does not show the usual a1 and b1 progressions of cGW transitions.

169
We believe that this peak is due to an intermolecular charge transfer (CT) effect where after ionization the charge hops from indole to the cGW which has the lower IE.

Mass selective R2PI spectroscopy: the S0–S1 transition

170
To record the R2PI spectrum of cGW the second harmonic of a Nd:YAG laser-pumped dye laser was tuned and the intensity of the parent ion monitored as a function of the wavelength.

171
In the R2PI spectrum shown in Fig. 6 an intense first transition is followed by several transitions which have decreasing intensities with increasing energy.

172
The spacing is not regular which may be a hint for the presence of several conformers in the beam.

173
UV/UV SHB spectroscopy proves (see section III.3) that this is not the case: All intense transitions correspond to one single conformer.

174
In the following we make use of this information and assign transitions a1, b1, c1, d1 and presumably Y to vibrations of this conformer.

175
A tentative conformer identification has been made in section IV.4.

176
The R2PI spectrum of cGW shown in Fig. 6 covers 450 cm−1.

177
No transitions were observed below 35 000 cm−1.

178
Above 35 300 cm−1 the indole vibrations and their many combination bands with vibrations a, b, c and d are expected.

179
The justification to confine our investigation to a range of 450 cm−1 is, that the four most stable calculated conformers have an attractive peptide ring-to-indole ring interaction and therefore should have similar S0–S1 origin energies.

180
The S0–S1 origin transition of cGW is situated at 35 058 cm−1.

181
It lies 180 cm−1 below the S0–S1 transition of indole34 and 178 cm−1 above the S0–S1 transition of 3-methylindole.35

182
The S1 energy in cGW is expected to be influenced by three parameters: (i) The electronic effect caused by the increase of the side chain of indole from 3-methylindole to cGW, (ii) the solvation of indole by the dipeptide ring and (iii) vibrational zero point energy effects.

183
Concerning (i) we expect in analogy to the alkylbenzene series (relative shift with respect to methylbenzene: +110 cm−1 for ethylbenzene, +105 cm−1 for propylbenzene, +99 cm−1 for butylbenzene36) a blue shift in comparison to 3-methylindole when increasing the chain size.

184
Concerning (ii): In a structure where the N–H group can interact with the indole π system (see most stable conformer 1 in Fig. 2a), a red shift of the S0–S1 transition is expected when the π system is excited, i.e. enlarged,18 because the N–H⋯π interaction is stabilized further.

185
However, (iii) this shift may be partially compensated by a change in the zero point vibrational energy if some vibrations increase strongly with the excitation to S1.

186
Interestingly the origin of S1 in cGW lies 353 cm−1 above the S1 of the lowest and 163 cm−1 above the second conformer of linear NH2–G–W–COOH.5,6,8

187
The fact that our red-shift is not as strong as in linear peptides may be interpreted as follows: The peptide ring formation forces the functional groups CO and N–H of the peptide backbone into positions which do not allow ideal solvation of the indole ring.

188
This lengthy discussion shows how difficult it is to interpret the tryptophan S0–S1 absorption in peptides with regard to structural changes of the environment.

189
The relatively intense transitions to one quantum of the low-frequency modes a1 (30 cm−1) and b1 (52 cm−1) indicate some S0–S1 geometry change in the corresponding coordinates.

190
Interestingly the transition a2 is not observed.

191
Only the combination band a1b1 and tentatively a small transition b2 is observed.

192
Because the intense indole vibrations lie above 300 cm−1,34 the vibrations a, b, c (108 cm−1) and d (174 cm−1) are attributed to vibrations involving a relative motion between the two rings.

193
For conformer 1a the frequency calculation of the neutral ground state gives the following vibrations below 90 cm−1: 21 cm−1 (global ring–ring bending), 39 cm−1 (ring–ring torsion along C1–C2), 51 cm−1 (boat to inverse boat motion in the peptide ring) and 59 cm−1 (higher order bending of the peptide ring, the angle at C0–C1–C2 and the indole ring).

194
The frequency of 51 cm−1 for the ring inversion mode is an artefact of the harmonic potential approximation.

195
Because of the small energy difference of 41 cm−1 between conformers 1a and 1b, at the energy of the first vibrational quantum the real potential width is much broader as assumed by the harmonic approximation.

196
This should lead to a very low experimental frequency.

197
Excitation of this motion would, however, only be expected if the peptide ring would have a strong coupling to the indole ring, which is not expected.

198
For a tentative assignment of the vibration in S1 we therefore only take the global ring–ring vibrations into account.

199
The ring–ring interaction in the lowest-energetic conformer 1 is attractive.

200
One can expect that this attraction increases in the S1 state due to the increase of the size of the π system upon excitation.

201
As a consequence the frequencies of the ring–ring motions would increase.

202
Hence the S0 vibrations with the energy 21 cm−1 and 39 cm−1 might correspond to the S1 vibrations a1 and b1.

203
The type of these vibrational motions also agrees with the concept that the N–H⋯π interaction changes when exciting from S0 to S1.

204
The S1 vibrations observed in the S0–S1 spectra of linear NH2–G–W–COOH are 24 cm−1 and 38 cm−1 and not too far off our values although the structures are very different.5,6,8

205
This indicates that the first two or three low-frequency vibrations of all similar-sized chromophore–tail molecules lie in the range of 20–50 cm−1.

206
For conformer 2 for example the first three calculated S0 vibrational frequencies are with 26 cm−1, 41 cm−1 and 56 cm−1 very similar to those of conformer 1.

207
Hence the S0 frequencies provide no arguments for conformer assignment.

208
The S0–S1 R2PI spectrum was recorded with 5 bar Ar nozzle backpressure.

209
Experiments with He and Ne as carrier gas resulted in less perfect vibrationally cooling.

Spectral UV/UV hole burning spectroscopy (UV/UV SHB)

210
To clarify the number of conformers in our molecular beam, we used R2PI UV/UV SHB.

211
As mentioned in the experimental section, for the SHB a three-plate ion optics was mounted and the PE spectrometer used as a linTOF mass spectrometer.

212
This setup was needed in order to remove and neutralize the prompt ions formed by the hole burning laser by applying a pulsed electric field.

213
The spectrum was recorded by scanning the hole burning laser while the probe laser was set to the S0–S1 origin transition.

214
Each time when the neutral ground state of the cGW conformer was depopulated by the hole burning laser a dip appeared in the probe laser parent ion signal.

215
Note that with this scheme the transition moments of the individual peaks determine the depth of the SHB dip.

216
In Fig. 7a the UV/UV SHB spectrum of cGW and in (b) for comparison the R2PI spectrum are shown.

217
We confined our investigation on the low-energetic section of the R2PI spectrum of cGW because (i) it contains the most intense transitions and (ii) the lowest-energetic four conformers are expected to lie in this range because of their similar peptide ring-to-indole ring interaction.

218
In Fig. 7a all transitions except E show a strong hole burning effect and hence correspond to one conformer.

219
To investigate the hole burning for transition E in detail, but also to probe accidental overlap of conformer transitions, we carefully checked the degree of hole burning for each transition.

220
Fig. 8 shows the degree of depletion of the probe signal by the hole burning laser for the first six transitions of the R2PI spectrum.

221
The conditions have been the same as for Fig. 7a, except that now the hole burning laser is kept fix at transition A and the probe laser is scanned to the individual transitions A-F.

222
With this scheme the depletion percentage should be the same for all transitions.

223
The transition intensities recorded with the probe laser have been measured with the hole burning laser “off” (left trace) and “on” (right trace).

224
For each “off” and “on” sequence of each transition 15 measurements (each averaged over 20 laser shots) are displayed.

225
The hole burning reduces the probe signal by 50% for all transitions.

226
This shows that no origin transition of a minor conformer component is accidentally covered by a vibronic transition of the most stable conformer.

227
Note that theoretically more than 50% of the neutral ground state levels can be depopulated37.

Photoelectron spectroscopy

228
Photoelectron spectroscopy is in principle a double-resonance spectroscopy which is able to provide conformer- and state-selective radical cation spectra.

229
The initial hope was that we could sort out and assign the different conformers of cGW by PE spectroscopy alone.

230
To achieve a resolution of 1 meV it is necessary to keep the electron kinetic energy below 200 meV.38

231
To fulfil this condition a second frequency-doubled dye laser which is pumped by the second harmonic of the same Nd:YAG laser was used for ionization of excited S1 cGW.

232
The R2PI PE spectra of cGW ionized via the vibrationless origin transition of S1 and the transition S1 a1 are shown in Fig. 9a and b.

233
Both spectra show a broad ionization onset with no vibrational structure.

234
According to the calculations of the neutral and cation structures which show a great geometry change upon ionization (Fig. 4a) we believe that we do not observe the origin transition to the cation ground state, but excite to a dense higher excited vibrational manifold.

235
Due to the smooth onset the IE can only be determined with an accuracy of ±3 meV.

236
The fact that the ionization onset via the S1a1 state is 6 cm−1 lower than when ionized via the S1 000 state (7.715 eV ±3 meV) is attributed to a better Franck–Condon factor to lower vibrational quantum numbers in the cation ground state.

237
The PE spectrum of indole (Fig. 9c) taken in the same machine shows that we could resolve vibrational structures if their spacing were larger than 10 cm−1.

238
The vibrations of the indole cation ground state have been recently investigated in detail by ZEKE spectroscopy39 and are not discussed here.

239
The IE value of indole determined by ZEKE spectroscopy allows us to calibrate our indole and cGW PE spectra with an accuracy of ±2 meV.

240
Unfortunately there is a background in the peptide PE spectra which stems from ionization of neutral fragments, preferentially indole, which is wavelength-dependent.

241
It is small for the wavelength range of the first two S1 transitions of cGW, but increases as the closer the wavelength comes to the indole S0–S1 origin.

242
Without coincidence techniques it is in principle impossible to distinguish between fragment and peptide photoelectrons.

243
We therefore confined our PE measurements to the first two S0–S1 transitions.

244
In Fig. 9 the IE of 3-methylindole (7.514 eV40) is also marked.

245
The IEs of substituted indoles should monotonically decrease with chain length (alkylbenzene41).

246
The fact that with respect to 3-methylindole the IE of cGW is blue-shifted and the smooth onset indicates a relatively strong repulsive interaction between the two rings in the cation, as predicted by theory.

247
Unfortunately the calculation of absolute IEs is relatively inaccurate and the measured IE corresponds neither directly to the calculated adiabatic IE nor the vertical IE.

248
The fact that conformer 1 has the highest IE of all conformers agrees well with the observed blue shift with respect to 3-methylindole.

Conclusion and tentative conformer assignment

249
DFT predicts that one conformer (conformer 1) is predominantly stable and further favoured by entropy due to a low-energy peptide ring inversion mode.

250
This is qualitatively in agreement with our R2PI and R2PI SHB spectra in which we attribute the most intense transition in the spectra to one conformer.

251
In conformer 1 the N–H group of the dipeptide ring interacts with the π system.

252
This interaction is attractive and is expected to increase with S1 excitation.

253
This is qualitatively confirmed by comparison to the calculated low-frequency S0 vibrations.

254
The fact that the experimentally observed conformer has an IE which is 195 meV higher than 3-methylindole shows that the solvation of the charge in the indole by the peptide ring is repulsive.

255
The most stable neutral conformer has the highest vertical IE of all 8 conformers and shows a very strong geometry change upon ionization.

256
Hence conformer 1 fulfils all experimental findings.

257
We admit that also conformer 2 and 3 follow the pattern of a attractive ring–ring interaction in the neutral cGW and a repulsive ring–ring interaction in the cationic cGW and that our arguments for the conformer assignment are very weak.

258
For a confident conformer assignment it would be helpful to have several conformers in the molecular beam and be able to compare their IEs.

259
Because the relative energies of calculated IEs are more accurate this would allow a more confident assignment.

260
In conclusion, we propose that conformer 1 is the cGW conformer we observe in our molecular beam.

261
Our work shows that R2PI and the R(1 + 1′)PI PE spectra predominantly provide information about the structural changes in the nearest environment of the chromophore.

262
In cGW, however, the structural changes upon ionization are so great that with conventional high resolution PES a resolution of vibrational structures in the PE spectrum is impossible.

263
In general a similar neutral-to-cation Franck–Condon problem might exist for peptides or amino acids which have polar groups near the chromophore.

264
The IR/UV double resonance technique may provide additional clues for conformer identification and especially allow us to distinguish between conformers 1 (N–H⋯π interaction) and 2 (C–H⋯π interaction).