Does the Residues Chirality Modify the Conformation of a Cyclo

Sep 5, 2017 - (46) The vibrational signature of the studied systems is obtained by means of infrared multiple photon dissociation (IRMPD) spectroscopy...
4 downloads 9 Views 1MB Size
Subscriber access provided by UNIV OF CAMBRIDGE

Article

Does the Residues Chirality Modify the Conformation of a Cyclo-Dipeptide? Vibrational Spectroscopy of Protonated Cyclo Diphenylalanine in the Gas Phase Ivan Alata, Ariel Francis Pérez-Mellor, Feriel Ben Nasr, Debora Scuderi, Vincent Steinmetz, Fabrice Gobert, Nejmeddine Jaïdane, and Anne Zehnacker-Rentien J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/acs.jpca.7b06159 • Publication Date (Web): 05 Sep 2017 Downloaded from http://pubs.acs.org on September 14, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry A is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Does the Residues Chirality Modify the Conformation of a Cyclo-Dipeptide? Vibrational Spectroscopy of Protonated Cyclo Diphenylalanine in the Gas Phase Ivan Alata1, Ariel Pérez-Mellor1, Feriel Ben Nasr1,2, Debora Scuderi3, Vincent Steinmetz3, Fabrice Gobert3, Nejm-Eddine Jaïdane2, Anne Zehnacker-Rentien1*. 1

Institut des Sciences Moléculaires d’Orsay (ISMO), CNRS, Univ. Paris-Sud, Université Paris-Saclay, F-91405 Orsay, France 2

Laboratoire de Spectroscopie Atomique Moléculaire et Applications (LSAMA) Université de Tunis El Manar, Tunis 1060, Tunisia

3

CNRS, Laboratoire de Chimie Physique (LCP), UMR8000, Orsay, F-91405 and Univ. Paris-Sud, Orsay, F-91405, France

Abstract

The structure of a protonated diketopiperazine dipeptide, cyclo diphenylalanine, is studied by means of Infra-Red Multiple Photon Dissociation (IRMPD) spectroscopy combined with quantum chemical calculations. Protonation exclusively occurs on the oxygen site and, in the most stable conformer, results to an intramolecular OH…π interaction, accompanied by a CH…π interaction. Higher-energy conformers with free OH and NH…π interactions are observed as well, due to kinetic trapping. Optimization of the intramolecular interactions involving the aromatic ring dictates the geometry of the benzyl substituents. Changing the chirality of one of the residue has consequences on the CH…π interaction, which is of CαH…π nature for LD while LL shows a CβH…π interaction. Higher-energy conformers also display some differences in the nature of the intramolecular interactions.

*e-mail : [email protected]

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

1. Introduction Chirality is one of the most striking characteristics of life and is ubiquitous in the chemistry of living organisms.1 The homochirality of life has been suggested to happen through initial symmetry breaking, followed by propagation and amplification,2 possibly by autocatalysis.3 It manifests itself by the incorporation of molecules of a specific chirality throughout life-related biological reactions. It is significant however that once the biological processes has stopped, homochirality tends to vanish. Amino acids racemization has indeed been observed in Pleistocene cave bear dentine and is used as a dating method.4 Chirality and stereochemistry play an important role in various biochemical reactions such as protein folding or enzymatic catalysis.1,5,6 Exogenous molecules like drugs also act differently depending on the chirality of their stereogenic centers.7 Numerous artificial drugs must therefore be of a specific chirality to be integrated into the complex biological system of living organisms. The different involvement of two enantiomers or diastereomers in biological reactions rests upon several factors, some of them being related to the environment or the solvent, others being intrinsic to the molecule itself. Gas phase studies enable studying the structural differences between diastereomers without the perturbations caused by the solvent.8-19 Studying each diastereomer in isolated conditions allows understanding the influence of chirality on intra and intermolecular interactions, which in turn govern the enantioselectivity of the reactive processes. Many biological molecules are protonated in physiological conditions, which influences their chemical and physical properties as well as the chemical processes they are involved in. Protonation can, for example, serve as a trigger for reversible structural interconversion from extended to helical forms of pyridine-derived oligoamides.20 The structure of protonated polypeptide often differs strongly from that of neutral ones. In particular, protonated polyalanines in the gas phase contrast with neutral ones and fail to adopt a helical conformation because the location of the charge at the N-terminus destabilizes the helix.21 Chirality-induced structural changes have been studied in small protonated peptides in isolated conditions. Limited differences have been observed between protonated or neutral LL and LD diphenylalanines or tetraphenylalanine.22-24 However, cyclic species seem to show increased sensitivity to stereochemical factors because of the constraints related to the ring. The structure of protonated model systems like 1-amino-2-indanol 25 or 4-hydroxyproline diastereomers15 are indeed sensitive to chirality, as also are neutral cyclic systems.13,16 This is why we have undertaken the study of cyclic dipeptides, both in the gas phase and in the condensed phase, with the aim of evidencing chirality-induced structural differences.26 Cyclic dipeptides result from the intramolecular peptide coupling of a linear α dipeptide.27-29 They all contain a Diketopyperazine (DKP) ring and are an important class of dipeptides with applications as active pharmaceutical ingredients.30,31 Cyclic dipeptides are often more stable than their linear counterpart in vivo and therefore promising in terms of therapeutic application.31 The DKP ring conformation depends on the nature of the residues and can be either planar, boat, or chair.32,33 34-36 A planar conformation is only observed for dipeptides molecules with small side groups like Glycine.37 The conformation remains planar for the slightly larger alanine (Ala) residue when the two residues are of opposite chirality, namely cyclo LAla-DAla. However, it becomes puckered for cyclo (LAla-LAla).33 Cyclic dipeptides containing a bulky residue like phenylalanine and tyrosine have a preferred conformation in which the aromatic ring in the side chain stacks over the DKP ring and the ring adopts a boat form.38 2 ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

This work aims at determining the structural and spectroscopic differences due to chirality for the two diastereomers of a protonated DKP dipeptide built on phenylalanine (Phe). The studied systems, namely, protonated cyclo (LPhe-LPhe), referred to as c-LLH+, and cyclo (LPhe-DPhe), referred to as c-LDH+, are shown in Figure 1. The gas-phase structure of the neutral peptides has been studied already, as well as that of the cyclo (LPhe-LPhe) crystal.26,36 The gas-phase neutral molecules display a boat conformation, with one Phe extended towards the amide nitrogen while the other one is stacked over the DKP ring. Due to the aromatic nature of phenylalanine, dispersion forces are important in the stabilization of the peptide and may favor OH… π,39-41 NH…π,42 or CH…π43,44 non-covalent interactions, which have been suggested to play a role in chiral recognition.19,45 The two diastereomer dipeptides are studied by combining Fourier transform ion cyclotron resonance (FTICR) mass spectrometry with laser spectroscopy.46 The vibrational signature of the studied systems is obtained by means of infrared multiple photon dissociation (IRMPD) spectroscopy and compared to quantum chemical calculations. Both fingerprint (1000-1800 cm-1) and OH or NH stretch (3100-3700 cm-1) regions are explored. IRMPD spectroscopy has been shown indeed to be a powerful tool for studying the structure of peptides, for example protonated peptides47-50, their complexes with alkali metals47 or divalent and trivalent metals,51,52 their cation or anions,53 as well as their reaction products54 or their derivatives with post-translational modifications.55-57 Studies of DKP peptides are limited to the protonated model systems, cyclo GlyGly58,59, cyclo AlaGly,60 and linear PhePhe in interaction with Ca2+.61

2. Methods 2.1. Experimental section The experiment setup has been described in details already and is only briefly summarized here.46 It involved a 7T Fourier FT-ICR hybrid mass spectrometer (Bruker, Apex Qe) coupled to a Free Electron Laser (FEL) at the Centre Laser Infrarouge d’Orsay (CLIO).62 The protonated species were generated by electrospraying a 100 μM solution obtained by diluting a 1 mM stock solution into a slightly acidified (1% formic acid) mixture of methanol and water (50:50). The cyclic dipeptides were obtained from GeneCust- Luxembourg (98% purity) and used without further purification. Infrared spectra were obtained by monitoring the fragmentation efficiency φ= -ln(P/(F + P)) as a function of the IR wavelength, with F being the sum of the abundances of the fragment ions produced by IRMPD and P that of the parent. The 900-2000 cm-1 region was covered resorting to the CLIO FEL.62 The CLIO FEL beam was mildly focused by a 2000 mm Ag-protected spherical concave mirror, with a typical beam waist at the centre of the FT-ICR of the order of 1 mm diameter. The CLIO FEL operated at 25 Hz with 8 µs long bunch pulses containing 0.5-3 ps micropulses separated by 16 ns. The spectral bandwidth (full width at half maximum FWHM) was about 7 cm-1 with pulse energy of 1600 to 900 mW from 900 to 2000 cm-1. The wavelength calibration was ensured by simultaneously recording the spectrum of polystyrene. The irradiation time was 1s. Additional broadband CO2 laser synchronized with CLIO (Universal Laser system 10 W at cw operation centered at λ = 10.6 μm) was employed to sustain the fragmentation. The CO2 laser pulse length was adjusted to 3.8 ms to avoid photo-dissociation of the molecule by the CO2 laser alone.

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The vibrational spectrum in the OH and NH stretch region was recorded resorting to an IR Optical Parametric Oscillator laser (LaserVision OPO, 25 Hz, 5-10 mJ/pulse). The irradiation time was 3 s and the CO2 laser pulse length was adjusted to 30 ms.

2.2. Computational details Exploration of the potential energy surface (PES) was performed using the OPLS-2005 and MMFFs force fields with the advanced conformational search implemented in the MacroModel suite, a part of the Schrödinger package.63 All the peptides structures found thereby with energy below 20 kJ/mol were optimized within the frame of the DFT theory at the b97-d3 /TZVPP level of theory, using the resolution of the identity approximation for the evaluation of the electron-repulsion integrals in order to decrease the calculation time.64 The conformers of the most stable family, namely, those protonated on oxygen (vide infra), were re-optimized using the dispersion-corrected functional B3LYP-D3 associated to the 6-311G++(d, p) basis set. The B3LYP-D3 functional was chosen here, as it was in previous studies,65,66 because it combines the good frequency description of the hybrid functional B3LYP67 and inclusion of empirical dispersion corrections, which might be important for peptides containing aromatic residues.68,69 Indeed, calculation of the interaction energy between neutral or protonated alanine and benzene indicate that methods taking dispersion into account are necessary to reproduce the OH…π, NH…π, and CH…π, interactions.70 To make sure that the whole PES was fully explored, structures were also built from the most stable neutral conformers found previously.26 A proton was added to each of the two oxygens, with the two possible orientations of the OH, and the resulting structure was optimized as described above. The vibrational spectra were simulated by convoluting the harmonic frequencies obtained at the same level of calculation by a Lorentzian line shape (FWHM 10 cm-1). A scaling factor of 0.970 in the 900-1700 cm-1 region and 0.955 in the 3000-3800 cm-1 region was applied to correct the frequencies for anharmonicity and basis set incompleteness. All calculations were performed with the Gaussian 09 packages71 or with the Turbomole 6.6 software.72

3. Results and discussion 3.1. CID fragments MS2 Collision-Induced Dissociation (CID) spectra of c-LLH+ and c-LDH+, obtained in the hexapole collision cell with a collision voltage of -6.5 V, are displayed in Figure 2. Both systems show similar fragmentation pattern and efficiency. For an activation time of 0.2 s, loss of CO (m/z 267) and formation of the iminium ion of Phe (m/z 120) are observed as the major fragmentation channels. For longer activation time, sequential fragmentation is observed. Subsequent loss of NH3 and CO from m/z 267 results to the m/z 250 and m/z 222 fragments. The fragmentation pattern described above is typical of DKP dipeptides.73-75 A mechanism has been proposed for model protonated DKP based on Glycine and Alanine, which involves either fragmentation from the most stable Oprotonated tautomer, in the case of cyclo AlaAla, or isomerization to the N-protonated tautomer with subsequent fragmentation, in the case of cyclo GlyGly.76 Minor differences in the fragments relative intensity are observed in c-LDH+ relative to cLLH+. While the primary fragments (m/z 267, m/z 120) are slightly more intense in c-LLH+, the

4 ACS Paragon Plus Environment

Page 4 of 22

Page 5 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

secondary fragments are slightly more intense in c-LDH+ .This observation indicates that c-LDH+ has slightly larger fragmentation efficiency.

3.2. IRMPD Spectra

The IRMPD spectra of LLCH+ and LDCH+ are presented in Figure 3. Their aspect is similar on the overall spectral range; they show however minor differences. The fingerprint region is dominated by an intense and congested feature in the region of the β(OH) and aliphatic β(CH) bending modes, which looks identical for the two diastereomers, with two maxima at at 1270-1310 cm-1 in c-LLH+ and 1280-1312 cm-1 in c-LLH+. The β(NH) range displays two bands of weak intensity at 1440 cm-1 and +

1493 in c-LDH+, whose counterpart seems to overlap in c-LLH , at 1452-1492 cm-1. Last, the ν(CO) stretch region is characterized by two intense bands located at 1674-1734 cm-1 in the c-LDH+ +

spectrum, and at 1681-1750 cm-1 in that of c-LLH . These bands are typical of the amide I vibration mode of protonated DKP rings.59 In the fingerprint region, the main differences between the two diastereomers is the blue shift of several bands of c-LLH+ relative to c-LDH+. The high frequencies region is characterized by two features identical for the two diastereomers: First, a broad and intense absorption that peaks at ~ 3250 cm-1 and extends from 3050 to 3350 cm-1. Second, an intense and narrow peak appears at 3415 cm-1. The main differences appear on the one hand in the ν(OH) stretch region where the spectrum of c-LDH+ shows two peaks at 3592 and 3564 cm-1, while that of c-LLH+ only shows a broadened peak at 3567 cm-1. Last, a peak appears at 3380 +

cm-1 in the spectrum of c-LDH+ while it is absent in that of c-LLH .

3.3. Theoretical results and assignment The nomenclature used for the molecule is shown in Figure 1. The peptide bond atoms will be referred to with a subscript “p” when they are on the protonated part. The carbon atoms of the residue will be denoted by their position relative to the peptide bond, i.e. Cα or Cβ. Last, the ipso carbon of the benzene rings will be denoted by the subscript “i”. Two protonation sites can be considered, namely, protonation of the amide oxygen and protonation of the amide nitrogen. For both c-LLH+ and c-LDH+ the most stable structures correspond to protonation of the oxygen atom. The most stable conformer protonated on the amide nitrogen is higher in energy, by about 13 kcal/mol, than the most stable conformer protonated on the amide oxygen. It is therefore unlikely that it is populated in our experimental condition at room temperature. The most stable N-protonated structures are shown in Figure S1 of the Supplementary Information. Moreover, the calculated spectra of the N-protonated species do not show satisfactory agreement with the experimental results, as shown in Figures S2 and S3 of the Supplementary Information. These conclusions are in line with proton affinity measurements resting on the Cooks kinetic method77 that indicate a larger stability of the O-protonated form of cyclic peptides due to the formation of a resonance stabilized cation.78 This results also parallel those deduced from the IRMPD spectra of the most simple DKP dipeptide, cyclo GlyGly.58

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The most stable conformers of c-LDH+ and c-LLH+ are presented in Figure 4 and the relative Gibbs energies given in Table 1. The energetic characteristics obtained at the B97-D3/TZVPP level are identical to those at the B3LYP-D3/6-311++g(d,p) level and are summarized in Table S1 of the Supplementary Information. In what follows, we will describe the B3LYP-D3 only. The most stable conformation of c-LDH+, CLD0, is stabilized by two interactions involving the aromatic ring. An ionic OpH…π interaction is accompanied by a neutral CαpH…π interaction. These interactions result to short distances of 2.51 (2.53) Å between OpH (CpH) and the center of the aromatic ring. The CLD1 conformation (1.6 kcal/mol) also shows an ionic OpH…π interaction, but the CαpH…π interaction is missing and is replaced by a NpH…π interaction. CLD2 and CLD3 both involve intramolecular NpH…π and CαH…π interactions but differ by the orientation of the OH group. Despite the resonance stabilization of the cation, the NpH…π interaction is less stabilizing than the OpH…π, which explains that CLD2 and CLD3 are higher in energy than the conformers that involve an OpH…π interaction, by a few kcal/mol. The NpH…π interaction is characterized by a distance between NpH and the center of the ring of more than 3 Å (3.06 Å for CLD2 for example). The two peptide bonds of the molecule are now different because of the dissymmetry introduced by protonation of one of the amide oxygens. This manifests itself by different NH bond lengths: 1.013 Å for the NH on the neutral side, vs. 1.17 to 1.025 Å on the protonated side depending on whether NpH is involved in an hydrogen bond or not. However, the amide bond remains planar or quasi planar (dihedral angle < 5°) in all cases. In the most stable conformer CLD0, the non-protonated amide bond is very similar to what is observed or calculated for the neutral molecule. The DKP ring structure is slightly puckered. It shows a flattened boat structure with a puckering angles of ~-10°, similar to that of the neutral molecule (-11°).36 One of the aromatic rings is in a flagpole position, which corresponds to an NCαpCβpCip dihedral angle of ~60°, very close to what is calculated for the neutral molecule (65°). The ~60° angle allows CH…π interaction to take place. In contrast, the NCαCβCi dihedral angle calculated for the Phe in extended position is of the order of 180°, very different from the −58° calculated in the neutral (unpublished results). This is due to a rotation of the phenyl that allows optimization of the OpH…π interaction. The other conformers differ from the most stable by the DKP structure, much closer to planarity, and the residue conformation. The NCαCβCi dihedral angle of ~180° or ~60° correspond to OpH…π (LCD1) or CαH…π (CLD2 and CLD3) interactions, respectively, while the NCαpCβpCip dihedral angle of ~55° corresponds to an NpH…π interaction. +

The results obtained for c-LLH parallel those described for c-LDH+ with however minor differences. CLL3 differs from CLD3 and stands out by the nature of the CH…π and NH…π interactions. While the NH…π interaction involves NpH on the protonated side in CLD3, it involves NH on the neutral side in CLL3 and is therefore weaker (distance between NH and the center of the ring of 3.3 Å). Conversely, CLD3 displays a CαH…π interaction while CLL3 shows a CβH…π interaction. The calculated Gibbs energy of the most stable structures, namely CLD0 and CLL0, indicates that CLL0 is more stable by 0.17 kcal/mol. at the B3LYP-D3/6-311++g(d,p) level. The calculated difference is within the error and one cannot draw definitive conclusions from these numbers. The experimental IRMPD spectrum of c-LDH+ and the spectra simulated for CLD0, CLD1, CLD2 and CLD3 at the B3LYP-D3/6-311++g(d,p) level of theory are presented in Figure 5. A tentative assignment of the main IRMPD features of c-LDH+ is proposed in Table 2 by comparing the experimental vibrational frequencies with those calculated for the IR active vibrational modes. Comparison with the B97-D3/TZVPP results is given in Figures S4 and S5 of the Supplementary 6 ACS Paragon Plus Environment

Page 6 of 22

Page 7 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Information. The agreement with the experimental data is slightly better at the B3LYP-D3/6311++g(d,p) level and only the latter results will be discussed. The highest-energy peaks observed in the experimental spectrum at 3564 and 3592 cm-1 are typical of free ν(OH) stretches. They appear in the same frequency range as the ν(OH) stretch observed by Wang et al. in the model protonated DKP, cyclo (GlyGly)H+, at 3585 cm-1.58 The 3564 and 3592 cm-1 bands can therefore be assigned to the free ν(OpH) stretch of CLD3 and CLD2, respectively. Their relatively weak intensity results from the low population of CLD2 and CLD3. It must be noted that the band assigned to the ν(OH) stretch dominates the spectrum of cyclo (GlyGly)H+, the other bands being one order of magnitude less intense. The other bands expected for CLD2 and CLD3 are therefore not expected to be detectable and the other features observed in the spectrum must be due to CLD0 and CLD1.79 It must be noted that a Maxwell-Boltzmann distribution at room temperature predicts negligible population of CLD3 and CLD2 compared to CLD0 (see Table 1), which should be the only conformer substantially populated. Kinetic trapping explains the presence of thermodynamicallydisfavored species, as already observed in mass spectrometry experiments involving an electrospray source.80-82 Another possibility would be incorrect description of the energy ordering, possibly due to imperfect description of dispersion. The doublet at 3380-3415 cm-1 is characteristic of the ν(NH) stretches of the protonated DKP ring. In cyclo (GlyGly)H+, the amide ν(NH) stretches appear at 3370 and 3430 cm-1 for the protonated and neutral amide, respectively. The peaks at 3380-3415 cm-1 are attributed to their counterpart in CLD0, possibly superimposed with the ν(NH) stretch of CLD1, which is located at the same position a that CLD0. The broad and intense band between 3350 and 3050 cm-1 is assigned to the ν(OpH) stretch of CLD0, possibly superimposed with that of CLD1. The broadening and the red shift of the band are due to the interaction between OpH and the π electron of the aromatic ring. The fingerprint region is not very specific: the superposition of the spectra simulated for the two most stable conformers yield satisfactory agreement with the experiment. Contribution of LD2 and LD3 is also compatible with the observed spectrum. The CO stretch appears at 1734 cm-1 while the other intense band at 1674 cm-1 is due to the ν(CpNp) stretch, which gains intensity due to the charge. These bands are typical of the amide I in protonated DKP rings.59 Indeed, this spectral region does not look the same for protonation on N (see Supplementary Information). The two partly resolved bands around 1480 cm-1 contains the contribution of the β(NH) and β(CH) bend, which gains intensity via the OH…π interaction. The intense feature at ~1280 cm-1 contains the contribution of multiple β(CH) bends, the most intense being the β(CαH). It should be noted that the bands due to the CO stretches are not as strong in the experimental spectrum as theoretically predicted, as often observed in peptides an DNA basis.24,83 The results for cLLH +, shown in Figure 6, parallel those for c-LDH+ and we will only stress the differences between them. The fingerprint region only shows minor differences that are difficult to interpret. In the high-energy region, the bands assigned to the ν(OH) stretches of CLL2 and CLL3 are closer to each other, which explains that the feature appearing as a distinct doublet in c-LDH+ is now a single band with a shoulder. The NH stretch region displays a single band, which contrasts with the doublet observed for c-LDH +. This is due to the fact that the OH…π interaction is weaker in CLL0 than in CLD0. As a result, the ν(OH) stretch is less intense and less shifted and overlaps with the ν(NH)

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

stretching mode located on the NH far from the protonation site, which is therefore hidden and does not appear in the spectrum.

4. Conclusion The diastereomers of a protonated cyclic dipeptide built on phenylalanine have been structurally characterized by IRMPD spectroscopy in two spectral regions, 1000-1800 cm-1 and 3000-3700 cm-1, combined with DFT calculations. Experimental findings, especially in the fingerprint region, and theoretical results unambiguously show that protonation happens on the amide oxygen for both c-LDH+ and c-LLH+. In particular, the amide I stretch region is typical of the DKP ring protonated on the oxygen. The experimental spectra can be satisfactorily reproduced only by taking into account the presence of the four most stable conformers, kinetically trapped in our experimental conditions. The DKP ring structure is slightly puckered and takes a dissymmetrical geometry due to protonation. The most stable form shows a DKP ring geometry close to that of the neutral form, namely a flattened boat geometry, for both c-LDH+ and c-LLH+. The chirality of the residues does not have dramatic influence on the molecular structure. However, c-LDH+ and c-LLH+ differ by the strength of the intramolecular OpH…π hydrogen bonds, which results in differences in the ν(NH)/ν(OH) stretch region. The higher energy conformers display a wide range of intramolecular hydrogen bonds, namely, OpH…π (CLD1 and CLL1) or CαH…π (CLD2, CLD3 and CLL2) accompanied by an NpH…π interaction. CLL3 stands out and differ from its CLD3 counterpart by showing CβH…π and NH…π interactions. The fragmentation efficiency is slightly larger for c-LDH+ than for c-LLH+, which is compatible with slightly smaller formation energy for CLD0 than for CLL0. However, the difference is within the error and one cannot rule out kinetic effects. 84-86 Experiments in cryogenic ion traps are in progress for better differentiation of the diastereomer cyclic dipeptides.87,88

5. Supporting Information In this section are shown the geometry and energetics of higher-energy conformers of c-LDH+ and c-LLH+, protonated on N, as well as their simulated IR spectra. This section also gives the spectra and energetics of the conformers discussed above calculated at the b97-D3/TZVPP level of theory.

6. Acknowledgments The research described here has been supported by the “Investissements d’Avenir” LabEx PALM contract (ANR-10-LABX-0039-PALM) and by the ANR contract (ESBODYR ANR-14-CE06-0019-03). Dr. Philippe Maître and the mass spectrometry platform SMAS of the Laboratoire de Chimie Physique (Univ. Paris-Sud) are gratefully acknowledged. We thank Dr J. M. Ortega and the CLIO team for technical assistance. We acknowledge the TGE FT-ICR ((FR 3624) for financial support and the use of the computing center MésoLUM of the LUMAT research federation (FR LUMAT 2764).

7. References

8 ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(1) Bonner, W. A. Chirality and Life. Origins Life Evol. Biosphere 1995, 25, 175-190. (2) Podlech, J. Origin of Organic Molecules and Biomolecular Homochirality. Cellular and Molecular Life Sciences 2001, 58, 44-60. (3) Soai, K.; Sato, I. Asymmetric Autocatalysis and its Application to Chiral Discrimination. Chirality 2002, 14, 548-554. (4) De Torres, T.; Ortiz, J. E.; Garcia, M. J.; Llamas, J. F.; Canoira, L.; De la Morena, M. A. G.; Julia, R. Geochemical Evolution of Amino Acids in Dentine of Pleistocene Bears. Chirality 2001, 13, 517-521. (5) Kwiecinska, J. I.; Cieplak, M. Chirality and Protein Folding. J. Phys.: Condens. Matter 2005, 17, S1565-S1580. (6) Nakamura, M.; Maki, S.; Amano, Y.; Ohkita, Y.; Niwa, K.; Hirano, T.; Ohmiya, Y.; Niwa, H. Firefly Luciferase Exhibits Bimodal Action Depending on the Luciferin Chirality. Biochem. Biophys. Res. Commun. 2005, 331, 471-475. (7) Nau, C.; Strichartz, G. R. Drug Chirality in Anesthesia. Anesthesiology 2002, 97, 497502. (8) Zehnacker, A. Chirality Effects in Gas-Phase Spectroscopy and Photophysics of Molecular and Ionic Complexes: Contribution of Low and Room Temperature Studies. Int. Rev. Phys. Chem. 2014, 33, 151-207. (9) Zehnacker, A.; Suhm, M. A. Chirality Recognition Between Neutral Molecules in the Gas Phase. Angew. Chem. Int. Ed. 2008, 47, 6970-6992. (10) Dean, J. C.; Buchanan, E. G.; James, W. H.; Gutberlet, A.; Biswas, B.; Ramachandran, P. V.; Zwier, T. S. Conformation-Specific Spectroscopy and Populations of Diastereomers of a Model Monolignol Derivative: Chiral Effects in a Triol Chain. J. Phys. Chem. A 2011, 115, 8464-8478. (11) Abo-Riziq, A. G.; Bushnell, J. E.; Crews, B.; Callahan, M. P.; Grace, L.; de Vries, M. S. Discrimination Between Diastereoisomeric Dipeptides by IR-UV Double Resonance Spectroscopy and ab initio Calculations. Int. J. Quantum Chem 2005, 105, 437-445. (12) Gloaguen, E.; Pagliarulo, F.; Brenner, V.; Chin, W.; Piuzzi, F.; Tardivel, B.; Mons, M. Intramolecular Recognition in a Jet-cooled Short Peptide Chain: gamma-turn Helicity Probed by a Neighbouring Residue. Phys. Chem. Chem. Phys. 2007, 9, 4491-4497. (13) Alauddin, M.; Gloaguen, E.; Brenner, V.; Tardivel, B.; Mons, M.; Zehnacker-Rentien, A.; Declerck, V.; Aitken, D. J. Intrinsic Folding Proclivities in Cyclic β-Peptide Building Blocks: Configuration and Heteroatom Effects Analyzed by Conformer-Selective Spectroscopy and Quantum Chemistry. Chem. Eur. J. 2015, 21. (14) James, W. H., III; Baquero, E. E.; Shubert, V. A.; Choi, S. H.; Gellman, S. H.; Zwier, T. S. Single-Conformation and Diastereomer Specific Ultraviolet and Infrared Spectroscopy of Model Synthetic Foldamers: alpha/beta-Peptides. J. Am. Chem. Soc. 2009, 131, 6574-6590. (15) Crestoni, M. E.; Chiavarino, B.; Scuderi, D.; Di Marzio, A.; Fornarini, S. Discrimination of 4-Hydroxyproline Diastereomers by Vibrational Spectroscopy of the Gaseous Protonated Species. J. Phys. Chem. B 2012, 116, 8771-8779. (16) Mahjoub, A.; Chakraborty, A.; Lepere, V.; Le Barbu-Debus, K.; Guchhait, N.; Zehnacker, A. Chirality-Dependent Hydrogen Bond Direction in Jet-cooled (S)-1,2,3,4-tetrahydro-3isoquinoline methanol (THIQM): IR-ion Dip Vibrational Spectroscopy of the Neutral and the Ion. Phys. Chem. Chem. Phys. 2009, 11, 5160-5169. (17) Sen, A.; Bouchet, A.; Lepere, V.; Le Barbu-Debus, K.; Scuderi, D.; Piuzzi, F.; ZehnackerRentien, A. Conformational Analysis of Quinine and its Pseudo Enantiomer Quinidine: a Combined Jet-cooled Spectroscopy and Vibrational Circular Dichroism Study. J. Phys. Chem. A2012, 116, 83348344. (18) Sen, A.; Lepere, V.; Le Barbu-Debus, K.; Zehnacker, A. How do Pseudoenantiomers Structurally Differ in the Gas Phase? An IR/UV Spectroscopy Study of Jet-Cooled Hydroquinine and Hydroquinidine. Chem. Phys. Chem : a European Journal of Chemical Physics and Physical Chemistry 2013, 14, 3559-3568. 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(19) Scuderi, D.; Le Barbu-Debus, K.; Zehnacker, A. The Role of Weak Hydrogen Bonds in Chiral Recognition. Phys. Chem. Chem. Phys. 2011, 13, 17916-17929. (20) Kolomiets, E.; Berl, V.; Lehn, J. M. Chirality Induction and Protonation-Induced Molecular Motions in Helical Molecular Strands. Chemistry-a European Journal 2007, 13, 5466-5479. (21) Kaleta, D. T.; Jarrold, M. F. Disrupting Helix Formation in Unsolvated Peptides. J. Phys. Chem. B 2001, 105, 4436-4440. (22) Dunbar, R. C.; Steill, J. D.; Oomens, J. Conformations and Vibrational Spectroscopy of Metal-Ion/Polylalanine Complexes. Int. J. Mass spectrom. 2010, 297, 107-115. (23) Dunbar, R. C.; Steill, J. D.; Oomens, J. Chirality-Induced Conformational Preferences in Peptide-Metal Ion Binding Revealed by IR Spectroscopy. J. Am. Chem. Soc. 2011, 133, 1212-1215. (24) Lepere, V.; Le Barbu-Debus, K.; Clavaguera, C.; Scuderi, D.; Piani, G.; Simon, A.-L.; Chirot, F.; MacAleese, L.; Dugourd, P.; Zehnacker, A. Chirality-Dependent Structuration of Protonated or Sodiated Polyphenylalanines: IRMPD and Ion Mobility Studies. Phys. Chem. Chem. Phys. 2016, 18, 1807-1817. (25) Bouchet, A.; Klyne, J.; Piani, G.; Dopfer, O.; Zehnacker, A. Diastereo-Specific Conformational Properties of Neutral, Protonated and Radical Cation Forms of (1R,2S)-cis and (1R,2R)-trans Amino-Indanol by Gas Phase Spectroscopy. Phys. Chem. Chem. Phys. 2015, 17, 2580925821. (26) Perez-Mellor, A.; Zehnacker, A. Vibrational Circular Dichroism of a 2,5Diketopiperazine (DKP) Peptide: Evidence for Dimer Formation in Cyclo LL or LD Diphenylalanine in the Solid State. Chirality 2017, 29, 89-96. (27) Borthwick, A. D. 2,5-Diketopiperazines: Synthesis, Reactions, Medicinal Chemistry, and Bioactive Natural Products. Chem. Rev. 2012, 112, 3641-3716. (28) Capasso, S.; Vergara, A.; Mazzarella, L. Mechanism of 2,5-Dioxopiperazine Formation. J. Am. Chem. Soc. 1998, 120, 1990-1995. (29) Basiuk, V. A.; Gromovoy, T. Y. The Gas Solid-Phase 2,5-Dioxopiperazine Synthesis Cyclization of Vaporous Dipeptides on Silica Surface. Collect. Czech. Chem. Commun. 1994, 59, 461466. (30) Prasad, C. Bioactive Cyclic Dipeptides. Peptides 1995, 16, 151-164. (31) Bellezza, I.; Peirce, M. J.; Minelli, A. Cyclic Dipeptides: from Bugs to Brain. Trends in Molecular Medicine 2014, 20, 551-558. (32) Karle, I. L. In The Peptides Analysis, Synthesis, Biology, Modern Techniques of Conformational Structural, and Configurational Analysis; Gross, E., Ed.; Elsevier: 1981; Vol. 4, p 1-54. (33) Sletten, E. Conformation of Cyclic Dipeptides. Crystal and Molecular Structures of Cyclo-D-Alanyl-L-Alanyl and Cyclo-L-Alanyl-L-Alanyl (3,6-Dimethylpiperazine-2,5-dione). J. Am. Chem. Soc. 1970, 92, 172-177. (34) Deslauriers, R.; Grzonka, Z.; Schaumburg, K.; Shiba, T.; Walter, R. C-13 Nuclear Magnetic-Resonance Studies of Conformations of Cyclic Dipeptides. J. Am. Chem. Soc. 1975, 97, 5093-5100. (35) Lin, C.-F.; Webb, L. E. Crystal Structures and Conformations of the Cyclic Dipeptides cyclo-(Glycyl-L-Tyrosyl) and cyclo-(L-Seryl-L-Tyrosyl) Monohydrate. J. Am. Chem. Soc. 1973, 95, 68036811. (36) Gdaniec, M.; Liberek, B. Structure of cyclo(-L-Phenylalanyl-L-phenylalanyl). Acta Crystallogr., Sect. C: Cryst. Struct. Commun. 1986, 42, 1343-1345. (37) Corey, R. B. The Crystal Structure of Diketopiperazine. J. Am. Chem. Soc. 1938, 60, 1598−1604. (38) Webb, T. H.; Wilcox, C. S. Enantioselective and Diastereoselective Molecular Recognition of Neutral Molecules. Chem. Soc. Rev. 1993, 22, 383-395. (39) Pribble, R. N.; Zwier, T. S. Size-Specific Infrared-Spectra of Benzene-(H2O)(n) Clusters (n=1 through 7) - Evidence for Noncyclic (H2O)(n) Structures. Science 1994, 265, 75-79.

10 ACS Paragon Plus Environment

Page 10 of 22

Page 11 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(40) Fredericks, S. Y.; Jordan, K. D.; Zwier, T. S. Theoretical Characterization of the Structures and Vibrational Spectra of Benzene-(H2O)n (n=1-3) Clusters. J. Phys. Chem. 1996, 100, 7810-7821. (41) Seurre, N.; Sepiol, J.; Lahmani, F.; Zehnacker-Rentien, A.; Le Barbu-Debus, K. Vibrational Study of the S0 and S1 States of 2-Naphthyl-1-ethanol/(water)(2) and 2-Naphthyl-1ethanol/(Methanol)(2) Complexes by IR/UV Double Resonance Spectroscopy. Phys. Chem. Chem. Phys. 2004, 6, 4658-4664. (42) Gloaguen, E.; Valdes, H.; Pagliarulo, F.; Pollet, R.; Tardivel, B.; Hobza, P.; Piuzzi, F.; Mons, M. Experimental and Theoretical Investigation of the Aromatic-Aromatic Interaction in Isolated Capped Dipeptides. J. Phys. Chem. A 2009, 114, 2973-2982. (43) Tsuzuki, S.; Fujii, A. Nature and Physical Origin of CH...π Interaction: Significant Difference From Conventional Hydrogen Bonds. Phys. Chem. Chem. Phys. 2008, 10, 2584-2594. (44) Tsuzuki, S.; Honda, K.; Fujii, A.; Uchimaru, T.; Mikami, M. CH...π Interactions in Methane Clusters With Polycyclic Aromatic Hydrocarbons. Phys. Chem. Chem. Phys. 2008, 10, 28602865. (45) Le Barbu-Debus, K.; Broquier, M.; Mahjoub, A.; Zehnacker-Rentien, A. Chiral Recognition in Jet-Cooled Complexes of (1R,2S)-(+)-cis-1-Amino-2-indanol and Methyl Lactate: on the Importance of the CH...π Interaction. Phys. Chem. Chem. Phys. 2009, 11, 7589-7598. (46) Bakker, J. M.; Besson, T.; Lemaire, J.; Scuderi, D.; Maitre, P. Gas-Phase Structure of a π-allyl-palladium Complex: Efficient Infrared Spectroscopy in a 7 T Fourier Transform Mass Spectrometer. J. Phys. Chem. A 2007, 111, 13415-13424. (47) Semrouni, D.; Balaj, O. P.; Calvo, F.; Correia, C. F.; Clavaguera, C.; Ohanessian, G. Structure of Sodiated Octa-Glycine: IRMPD Spectroscopy and Molecular Modeling. J. Am. Soc. Mass. Spectrom. 2010, 21, 728-738. (48) Martens, J. K.; Compagnon, I.; Nicol, E.; McMahon, T. B.; Clavaguera, C.; Ohanessian, G. Globule to Helix Transition in Sodiated Polyalanines. J. Phys. Chem. Letters 2012, 3, 3320-3324. (49) Gaigeot, M. P. Theoretical Spectroscopy of Floppy Peptides at Room Temperature. A DFTMD Perspective: Gas and Aqueous Phase. Phys. Chem. Chem. Phys. 2010, 12, 3336-3359. (50) Sediki, A.; Snoek, L. C.; Gaigeot, M. P. N-H+ Vibrational Anharmonicities Directly Revealed from DFT-Based Molecular Dynamics Simulations on the Ala(7)H(+) Protonated Peptide. Int. J. Mass spectrom. 2011, 308, 281-288. (51) Dunbar, R. C.; Steill, J. D.; Polfer, N. C.; Berden, G.; Oomens, J. Peptide Bond Tautomerization Induced by Divalent Metal Ions: Characterization of the Iminol Configuration. Angew. Chem. Int. Ed. 2012, 51, 4591-4593. (52) Prell, J. S.; Flick, T. G.; Oomens, J.; Berden, G.; Williams, E. R. Coordination of Trivalent Metal Cations to Peptides: Results from IRMPD Spectroscopy and Theory. J. Phys. Chem. A 2010, 114, 854-860. (53) Osburn, S.; Berden, G.; Oomens, J.; O'Hair, R. A. J.; Ryzhov, V. Structure and Reactivity of the N-Acetyl-Cysteine Radical Cation and Anion: Does Radical Migration Occur? J. Am. Soc. Mass. Spectrom. 2011, 22, 1794-1803. (54) Bythell, B. J.; Maitre, P.; Paizs, B. Cyclization and Rearrangement Reactions of a(n) Fragment Ions of Protonated Peptides. J. Am. Chem. Soc. 2010, 132, 14766-14779. (55) Correia, C. F.; Balaj, P. O.; Scuderi, D.; Maitre, P.; Ohanessian, G. Vibrational Signatures of Protonated, Phosphorylated Amino Acids in the Gas Phase. J. Am. Chem. Soc. 2008, 130, 3359-3370. (56) Lanucara, F.; Chiavarino, B.; Crestoni, M. E.; Scuderi, D.; Sinha, R. K.; Maitre, P.; Fornarini, S. S-nitrosation of Cysteine as Evidenced by IRMPD spectroscopy. Int. J. Mass spectrom. 2012, 330, 160-167. (57) Gregori, B.; Guidoni, L.; Chiavarino, B.; Scuderi, D.; Nicol, E.; Frison, G.; Fornarini, S.; Crestoni, M. E. Vibrational Signatures of S-Nitrosoglutathione as Gaseous, Protonated Species. J. Phys. Chem. B 2014, 118, 12371-12382. 11 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(58) Wang, D.; Gulyuz, K.; Stedwell, C. N.; Polfer, N. C. Diagnostic NH and OH Vibrations for Oxazolone and Diketopiperazine Structures: b(2) from Protonated Triglycine. J. Am. Soc. Mass. Spectrom. 2011, 22, 1197-1203. (59) Zou, S.; Oomens, J.; Polfer, N. C. Competition Between Diketopiperazine and Oxazolone Formation in Water Loss Products from Protonated ArgGly and GlyArg. Int. J. Mass spectrom. 2012, 316, 12-17. (60) Perkins, B. R.; Chamot-Rooke, J.; Yoon, S. H.; Gucinski, A. C.; Somogyi, A.; Wysocki, V. H. Evidence of Diketopiperazine and Oxazolone Structures for HA b(2)(+) Ion. J. Am. Chem. Soc. 2009, 131, 17528-17529. (61) Dunbar, R. C.; Steill, J. D.; Oomens, J. Encapsulation of Metal Cations by the PhePhe Ligand: A Cation-π Ion Cage. J. Am. Chem. Soc. 2011, 133, 9376-9386. (62) Prazeres, R.; Glotin, F.; Insa, C.; Jaroszynski, D. A.; Ortega, J. M. Two-Colour Operation of a Free-Electron Laser and Applications in the Mid-Infrared. Eur. Phys. J. D 1998, 3, 87-93. (63) MacroModel Version 9.8; ed. Schrödinger, LLC: New York, NY, 2010. MacroModel version 9.8; ed. Schrödinger, LLC: New York, NY, 2010 2010. (64) Weigend, F.; Haser, M. RI-MP2: First Derivatives and Global Consistency. Theor. Chem. Acc. 1997, 97, 331-340. (65) Altnoeder, J.; Bouchet, A.; Lee, J. J.; Otto, K. E.; Suhm, M. A.; Zehnacker-Rentien, A. Chirality-dependent balance between hydrogen bonding and London dispersion in isolated (+-)-1indanol clusters. Phys. Chem. Chem. Phys. 2013, 15, 10167-10180. (66) Semrouni, D.; Clavaguera, C.; Dognon, J. P.; Ohanessian, G. Assessment of density functionals for predicting the infrared spectrum of sodiated octa-glycine. Int. J. Mass spectrom. 2010, 297, 152-161. (67) Halls, M. D.; Velkovski, J.; Schlegel, H. B. Harmonic frequency scaling factors for Hartree-Fock, S-VWN, B-LYP, B3-LYP, B3-PW91 and MP2 with the Sadlej pVTZ electric property basis set. Theor. Chem. Acc. 2001, 105, 413. (68) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132. (69) Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456-1465. (70) Mohan, N.; Vijayalakshmi, K. P.; Koga, N.; Suresh, C. H. Comparison of Aromatic NH ...π, OH ...π, and CH ...π Interactions of Alanine Using MP2, CCSD, and DFT Methods. J. Comput. Chem. 2010, 31, 2874-2882. (71) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.et al.; Gaussian Inc.: Wallingford CT, 2009, p Gaussian 09, Revision A.02. (72) Ahlrichs, R.; Bar, M.; Haser, M.; Horn, H.; Kolmel, C. Electronic-Structure Calculations on Workstation Computers - the Program System Turbomole. Chem. Phys. Lett. 1989, 162, 165-169. (73) Bratakos, S. M.; Sinanoglou, V. J.; Matsoukas, M. T.; Siapi, E.; Papahatjis, D. P.; Riganakos, K.; Zoumpoulakis, P. Fragmentation Patterns of Aromatic 2,5-diketopiperazines Using Liquid chromatography/Mass Spectrometry. Curr. Anal. Chem. 2016, 12, 439-449. (74) Chen, Y. H.; Liou, S. E.; Chen, C. C. Two-Step Mass Spectrometric Approach for the Identification of Diketopiperazines in Chicken Essence. Eur. Food Res. and Tech. 2004, 218, 589-597. (75) Vandergreef, J.; Tas, A. C.; Nijssen, L. M.; Jetten, J.; Hohn, M. Identification and Quantitation of Diketopiperazines by Liquid-Chromatography Mass-Spectrometry, Using a Moving Belt Interface. J. Chromatogr. 1987, 394, 77-88. (76) Shek, P. Y. I.; Lau, J. K.-C.; Zhao, J.; Grzetic, J.; Verkerk, U. H.; Oomens, J.; Hopkinson, A. C.; Siu, K. W. M. Fragmentations of Protonated cyclic-GlycylGlycine and cyclic-AlanylAlanine. Int. J. Mass spectrom. 2012, 316, 199-205.

12 ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

(77) Brodbeltlustig, J. S.; Cooks, R. G. Determination of Relative Gas-Phase Basicities by the Proton-Transfer Equilibrium Technique and the Kinetic Method in a Quadrupole Ion-Trap. Talanta 1989, 36, 255-260. (78) Nold, M. J.; Cerda, B. A.; Wesdemiotis, C. Proton Affinities of the N- and C-Terminal Segments Arising upon the Dissociation of the Amide Bond in Protonated Peptides. J. Am. Soc. Mass. Spectrom. 1999, 10, 1-8. (79) Bodo, E.; Ciavardini, A.; Giardini, A.; Paladini, A.; Piccirillo, S.; Rondino, F.; Scuderi, D. Infrared Multiple Photon Dissociation Spectroscopy of Ciprofloxacin: Investigation of the Protonation Site. Chem. Phys. 2012, 398, 124-128. (80) Voronina, L.; Rizzo, T. R. Spectroscopic Studies of Kinetically Trapped Conformations in the Gas Phase: the Case of Triply Protonated Bradykinin. Phys. Chem. Chem. Phys. 2015, 17, 2582825836. (81) Xia, H. X.; Attygalle, A. B. Effect of Electrospray Ionization Source Conditions on the Tautomer Distribution of Deprotonated p-Hydroxybenzoic Acid in the Gas Phase. Anal. Chem. 2016, 88, 6035-6043. (82) Corinti, D.; De Petris, A.; Coletti, C.; Re, N.; Chiavarino, B.; Crestoni, M. E.; Fornarini, S. Cisplatin Primary Complex with L-Histidine Target Revealed by IR Multiple Photon Dissociation (IRMPD) Spectroscopy. Chem. Phys. Chem 2017, 18, 318-325. (83) Salpin, J. Y.; Scuderi, D. Structure of Protonated Thymidine Characterized by Infrared Multiple Photon Dissociation And Quantum Calculations. Rapid Communications in Mass Spectrometry 2015, 29, 1-7. (84) Horeau, A.; Nouaille, A. Micromethod for Determining Configuration of Secondary Alcohols by Kinetic Reduction - Use of Mass-Spectrography. Tetrahedron Lett. 1990, 31, 2707-2710. (85) Vedejs, E.; Jure, M. Efficiency in Nonenzymatic Kinetic Resolution. Angew. Chem. Int. Ed. 2005, 44, 3974-4001. (86) Liang, Y. J.; Bradshaw, J. S.; Dearden, D. V. The Thermodynamic Basis for Enantiodiscrimination: Gas-Phase Measurement of the Enthalpy and Entropy of Chiral Amine Recognition by Dimethyldiketopyridino-18-crown-6. J. Phys. Chem. A 2002, 106, 9665-9671. (87) Soorkia, S.; Broquier, M.; Gregoire, G. Conformer- and Mode-Specific Excited State Lifetimes of Cold Protonated Tyrosine Ions. J. Phys. Chem. Letters 2014, 5, 4349-4355. (88) Feraud, G.; Broquier, M.; Dedonder, C.; Jouvet, C.; Gregoire, G.; Soorkia, S. Excited State Dynamics of Protonated Phenylalanine and Tyrosine: Photo-Induced Reactions Following Electronic Excitation. J. Phys. Chem. A 2015, 119, 5914-5924.

13 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 22

Figure and tables Table 1: Relative electronic energies (ΔE) and Gibbs energies (ΔG) at room temperature (298 K ) in kcal/mol calculated at the B3LYP-D3/6-311++g(d,p) level of theory. Calculated (NH) and (OH) distances (Å) and relevant angles (°).

ΔE (kcal/mol) B3LYP-D3/6311++g(d,p)

ΔG (kcal/mol) Boltzmann B3LYP-D3/6factor 311++g(d,p) (% )

CαpCpNpCα dihedral angle (°)

CαCNCαp dihedral angle (°)

NCαCβCi dihedral angle (°)

NCαpCβpCip dihedral angle (°)

d(Np…H) d(OpH) (A°) (A°)

CLD0

0.0

0.0

93

-16

-11

+174

61

1.017

0.988

CLD1

3.8

1.6

6

9

2

-170

56

1.022

0.985

CLD2

4.1

2.8

1

-16

-2

61

54

1.024

0.968

CLD3

4.8

3.8