Conformational Space and Photochemistry of ... - ACS Publications

Subscriber access provided by HAIFA UNIV

Article

Conformational Space and Photochemistry of Tyramine Isolated in Argon and Xenon Cryomatrices Barbara Michela Giuliano, Sonia Melandri, Igor Reva, and Rui Fausto J. Phys. Chem. A, Just Accepted Manuscript • DOI: 10.1021/jp404707e • Publication Date (Web): 08 Sep 2013 Downloaded from http://pubs.acs.org on September 15, 2013

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 37

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

Conformational Space and Photochemistry of Tyramine Isolated in Argon and Xenon Cryomatrices

Barbara M. Giuliano,a Sonia Melandri,a,b* Igor Revaa and Rui Faustoa* a

Department of Chemistry, University of Coimbra, P-3004-535 Coimbra, Portugal b Department of Chemistry, University of Bologna, I-40126 Bologna, Italy

Abstract: The infrared spectra of tyramine monomers trapped in low temperature argon and xenon matrices were recorded. The presence of the flexible ethylamino side chain gives rise to a complex conformational surface that contains several minima of relatively low energies, some of them stabilized by a weak N−H…π hydrogen bond interaction between the amino group and the phenyl ring. The experimental infrared spectra confirm the presence of at least two stable conformers isolated in the matrices. Annealing experiments performed on the xenon matrix revealed a change in the relative population of the experimentally relevant conformers upon isolation in this polarizable matrix, compared to gas phase. The general interpretation of the spectra was based on harmonic and anharmonic quantum chemical calculations, undertaken at the DFT/B3LYP and MP2 levels of theory with the 6-311++G(d,p) basis set. The photochemical behaviour of the matrix-isolated compound upon narrowband UV irradiation was also investigated. Identification of ketene species in the spectra of the irradiated matrices suggests occurrence of a ring-opening reaction, which in the xenon matrix occurs concomitantly with the conformational isomerization of tyramine.

Keywords: Tyramine / Matrix isolation / IR spectroscopy / B3LYP and MP2 calculations / Conformational space / UV-Induced isomerization *

Corresponding authors e-mail: [email protected] (R. [email protected] (S. Melandri); +39 051 2099502 1 ACS Paragon Plus Environment

Fausto);

+351

239

854483;


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 2 of 37

1. Introduction Tyramine (4-hydroxyphenethylamine; TYR) is a biogenic amine derived from the decarboxylation of the amino acid tyrosine. It can be found in aged foods as a degradation product1,2 and can act as a releasing agent.3,4 The presence of the flexible ethylamino group as side chain of its phenolic aromatic ring is an essential structural feature of tyramine, leading to a complex conformational surface, which contains several minima of relatively low energy. This substance has already been investigated by microwave rotational spectroscopy and laser spectroscopy in the supersonic jet environment.5-9 Such experiments have evidenced the presence of a minimum of four, up to seven, different conformers in the gas phase, depending on the detection technique employed. In particular the dispersed fluorescence and IR spectra of tyramine are reported in ref. 7, where the vibrational spectra of seven different conformers were observed. All the gas phase experiments agree that the most stable conformers are stabilized by weak intramolecular hydrogen bond interactions and exhibit a folded geometry in which the amino hydrogen atoms interact with the π cloud of the aromatic ring. To the best of our knowledge, the infrared spectra of tyramine isolated in cryogenic matrices have never been reported before, nor its photochemistry has been investigated hitherto. In order to fill these gaps, we present in this article a study on the structure and infrared spectra of tyramine isolated in argon and xenon matrices, exploring its conformational landscape with the help of harmonic and anharmonic quantum chemical calculations, and analyzing its reactivity following narrowband UV irradiation.

2. Experimental and calculation methods

A sample of tyramine (99% purity) was purchased from Sigma-Aldrich and used without further purification. Matrices were prepared by co-deposition of tyramine vapors from a heated glass furnace containing a sample of the compound, and a large excess of the matrix gas (argon N60 from Air Liquide; scientific xenon 5.0 from Linde), onto the CsI optical substrate of the cryostat kept at 15 K and 30 K, respectively for argon and xenon matrices. An APD Cryogenics closed-cycle helium refrigeration system with a DE-202A expander was used. The temperature

2 ACS Paragon Plus Environment

Page 3 of 37

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


was measured directly at the sample holder by a silicon diode temperature sensor, connected to a digital controller (Scientific Instruments, model 9650-1), with an accuracy of 0.1 K. The IR spectra were recorded using a Nicolet 6700 Fourier transform infrared spectrometer, equipped with a deuterated triglycine sulphate (DTGS) detector and a Ge/KBr beam splitter, with 0.5 cm−1 resolution. In the photochemical experiments, narrowband tunable UV radiation was provided by the frequency doubled signal beam of the Quanta-Ray MOPO-SL optical parametric oscillator (fwhm ~0.2 cm−1, repetition rate 10 Hz, pulse energy ~3 mJ) pumped with a pulsed Nd:YAG laser. The quantum chemical calculations were performed with the Gaussian 03 software package.10 Equilibrium geometries for all studied species were fully optimized at the DFT level of theory, using the standard 6-311++G(d,p) basis set, and the three-parameter density functional B3LYP, which includes the Becke gradient exchange correction and the Lee, Yang and Parr correlation functional. The geometries and relative energies of the TYR minimum energy conformations were also computed using the MP2 method, with the 6-311++G(d,p) basis set. Both harmonic and anharmonic vibrational calculations were undertaken. The first were obtained at both the B3LYP and MP2 levels of theory. Anharmonic frequencies were obtained just at the B3LYP level of approximation. The theoretical normal modes were analyzed by carrying out potential energy distribution (PED) calculations. Transformations of the force constants with respect to the Cartesian coordinates to the force constants with respect to the molecule fixed internal coordinates allowed the PED analysis to be carried out as described by Schachtschneider and Mortimer.11 The internal coordinates used in this analysis were defined as recommended by Pulay et al.12 and are listed in Table S1 (Supporting Information).

3. Results and discussion

The molecule of tyramine (Figure 1) has four internal rotational degrees of freedom. Two of them are associated with the spatial orientation of the ethylamino side chain and originate from rotations around the CH2−CH2 and CH2−N bonds. The remaining two define the relative positions of the ethylamino moiety and hydroxyl group in relation to the aromatic ring, and correspond to internal rotations about the Caryl−CH2 and Caryl−O bonds, respectively.



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 4 of 37

The systematic inspection of the PES of the molecule indicated that the optimized orientations of both the ethylamino and hydroxyl substituents remained essentially perpendicular and parallel to the ring plane, respectively, for all minimum energy conformations, independently of the conformation within the ethylamino moiety. Illustrative relaxed potential energy scans around the Caryl−CH2 and Caryl−OH bonds are shown in Figure 1. We choose, as example, the case where the conformations about the CH2−CH2 and CH2−N bonds of the ethylamino group are anti (Caryl−C−C−N dihedral angle is ~180°) and gauche (C−C−N−Lp is ~60°; Lp= lone electron pair of N). As it can be seen in this Figure, one full rotation around the Caryl−CH2 or Caryl−OH bond results in two minima, designated as Ag and Ag’ (first letter indicates the conformation about the CH2−CH2 bond; the second that about the CH2−N bond). If the conformation around the C−C−N−Lp dihedral angle was chosen as anti (~180°), the two minima would be mirror reflections of each other (Aa), corresponding either to reflection through the plane of the phenyl ring (flip of the ethylamine group), or in the plane bisecting the phenyl ring (flip of the OH group). According to the B3LYP calculations, the transformation between Ag and Ag’ via rotation about the Caryl−CH2 bond is associated with a barrier lower than 9 kJ mol−1, while that obtained via rotation of the OH group involves a larger barrier (about 15 kJ mol−1). It is important to note that transformations between the symmetry-equivalent minima can be accomplished either via internal rotation of the ethylamino fragment or, with the same effect, of the hydroxyl group. Then, to characterize a set of unique conformations of tyramine, one of these internal rotations does not need an explicit consideration. Since the OH group internal rotation is related with a higher energy barrier, in the further analysis we have oriented the OH group to have the cis value of the C2−C1−O−H dihedral angle. We kept this OH group orientation as a reference (or, equivalent to say, as an “asymmetry marker” of the phenyl ring) while analyzing all possible orientations of the ethylamino side chain. For a more detailed conformational characterization of the ethylamino fragment, we have performed a two-dimensional potential energy scan where the driving coordinates correspond to the dihedral angles related with the internal rotations around the CH2−CH2 and CH2−N bonds. These angles were incrementally fixed with a step of 20°, in all possible 324 combinations, while all other geometric parameters were fully optimized. The resulting potential energy surface is represented in Figure 2.


Page 5 of 37

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 complete potential energy surface contains 9 different minima. According to the adopted naming scheme, these minima are designated using two letters, appearing on the top and right sides of the surface representation in Figure 2. The first (capital) letter refers to the position of the amino group with respect to the phenyl ring, defined by the Caryl−C−C−N dihedral angle. The second (small) letter specifies the orientation of the amino group with respect to the ethylamino chain, defined by the C−C−N−Lp dihedral angle. This naming system is very similar to that implemented previously in the study of tryptamine,13 a molecule that also bears an ethylamino group connected to an asymmetric aromatic ring. The CH2−CH2 bond rotation gives rise to two types of conformers: the gauche conformers (denoted with the capital letters G and G’), adopting a “folded” position of the NH2 group with respect to the aromatic ring, and the anti forms (denoted with A) which assume an unfolded geometry. The internal rotation of the NH2 group produces three more orientations (denoted with g, a and g’) for each gauche or anti arrangement around the CH2−CH2 bond. Depending on the orientation of the amino group, some of the folded conformers may involve the presence of a weak N−H…π hydrogen bond interaction, while in the case of the unfolded conformations such interaction is impossible. The structures of the 9 different conformers of tyramine are presented in Figure 3. Their calculated energies are listed in Table 1. According to both the B3LYP and MP2 calculations, the lowest energy conformations are those stabilized by the intramolecular hydrogen bonding between one of the amino hydrogen atoms and the π electron cloud of the aromatic ring (Gg, G’g’, Ga and G’a; see also Figure 3). In the two folded conformations, where the nitrogen lone pair is pointing toward the phenyl ring (Gg’ and G’g), this interaction is not possible, and their energies relative to the most stable conformers are predicted to be much higher. The energies of the anti forms about the CH2−CH2 bond were predicted by the MP2 method to be only slightly lower than those of the Gg’ and G’g forms and considerably higher than those of the folded conformations bearing the stabilizing N−H…π interaction. At the B3LYP level of theory the relative energies of the anti forms about the CH2−CH2 bond were estimated to be smaller, with the Aa conformer being predicted to correspond to the second most stable form of tyramine. We can expect that the relative energies of the conformers are better predicted by the MP2 calculations, since this method is well-known to describe more consistently the intramolecular correlation energy than DFT based methods.5,14,15



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 6 of 37

From the MP2 calculated Gibbs energies at 298 K, the relative abundances of the different tyramine conformers in the gas phase equilibrium at that temperature were estimated (Table 1). The abundance of the conformers in a low-temperature matrix may, however, differ substantially from those, because of conformer interconversion either during deposition of the matrix or in the matrix itself. The probability of the occurrence of the conformational relaxation is related to the energy barriers between the conformers. In matrices, conformational cooling has been demonstrated to occur in compounds possessing barriers to intramolecular rotation of a few kJ mol−1.16-20 The B3LYP and MP2 predicted barriers for interconversions involving the nine conformers of tyramine can be seen in Figure 2 and in Scheme 1, respectively. As a general feature, the barriers for the NH2 group torsion are lower than those for the torsion of the ethylamino backbone. The conversion of Gg’ into the global minimum Gg shows the smallest calculated barrier, less than 12 kJ mol−1. This value is small enough to assume that the Gg’ form will be converted into the Gg form during matrix deposition. The same assumption can be made for the conversion of G’g into G’g’. Likewise, the barriers for the interconversions between the anti CH2−CH2 conformers (Ag, Ag’ and Aa), which are between 12 and 17 kJ mol-1, can also be considered low enough to lead us to expect that only one of these conformers (either the most stable Ag’ form or the Ag one which possesses the highest value of the dipole moment, according to the MP2 calculations, see Scheme 1) will be present in the low-temperature matrices after deposition. Similarly, the barriers for the amino group rotation converting the Ga form into the most stable Gg, and accordingly the G’a into the G’g’, have been calculated to be slightly higher (~12 kJ mol−1), which is in agreement with the presence of the weak N−H…π interaction in these conformers. Hence, the presence of Ga and G’a conformers in the matrices is not expected. Also, we will present in the following section, the experimental features for postulating the absence of these two conformers in the deposited matrices. Note that in alanine, like in tyramine, many conformers differing by the NH2 group orientation are possible. The corresponding barriers for the NH2 group internal rotation in alanine are of the order of ~10 kJ mol−1. The less stable forms of alanine readily collapse to the most stable conformer around this coordinate during deposition of the matrices and only two forms of alanine with different heavy atom backbone orientations could be trapped upon the matrix deposition.21


Page 7 of 37

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


Conformational isomerization about the CH2−CH2 bond involves rotation of heavy atoms. The associated barriers were calculated to be no less than 16 kJ mol−1 (Scheme 1), being the highest barriers for conformational isomerization in the studied molecule. This result can be used in favor of the presence of both gauche and anti-type conformers in the deposited matrices. Finally, the barrier for the rotation of the OH substituent was calculated to be around 15 kJ mol−1, independently of the conformation of the ethylamino substituent. Our recent studies on cytosine22 showed that two hydroxy conformers could convert to each other in the matrix even over the barrier of 30 kJ mol−1. Rotation of the OH group interconverts the Gg, Ga, Ag’ and Gg’ conformers into the G’g’, G’a, Ag and G’g forms, respectively. It is interesting to note that the conformers that differ only in the orientation of the OH group are almost isoenergetic (as can be seen from Table 1 and Scheme 1) but the conformer with the highest dipole moment should be favored in the matrix. Nevertheless since their calculated spectra are practically coincident, at our level of experimental resolution it is impossible to discriminate between these conformers. For this reason in the spectral analysis presented below, only one member of each of these classes of conformers will be considered and it will be the one with the highest dipole moment. We will then use G’g’ for the couple Gg/G’g’, Ag for Ag/Ag’, Ga for Ga/G’a and G’g for G’g/Gg’. After all these considerations, taking into account all the calculated values for the isomerizations barriers, and the dipole moments, the species that are expected to be predominant in the matrix after deposition are the Ag and G’g’ conformers. With the assumptions resulting from the analysis of the potential energy landscape of tyramine and taking into account the expected conformational conversions between conformers during matrix deposition, we can now roughly estimate the relative abundance (%) of the conformers in the as-deposited matrices, using the theoretically predicted (MP2) populations for the conformers in the gas phase before deposition presented in Table 1. The obtained value, G’g’:Ag = 71:29 (i.e., ~2.4), is, as it could be anticipated (see Figure 4 ), somewhat smaller than that found in ref. 6 for jet cooled tyramine using R2PI laser spectroscopy (in the spectrum shown in Figure 3 of ref. 6, the relative intensity of the bands assigned to the gauche conformers is ca. three times larger than the intensity of the bands assigned to the anti forms).



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 8 of 37

a. Band assignments Due to the characteristics of the molecular system under study (high flexibility of the side chain, relatively large molecule, similar spectra for the various conformers) the assignment of the main spectral features to individual species is difficult to undertake unequivocally. The most intense bands are split in various components that can be assigned to the presence of different conformers as well as to matrix site effects, which are particularly common for argon matrices. A further complication arises from the fact that the most conformation-sensitive vibrational modes have a highly anharmonic character. Therefore, anharmonic predictions of the spectra have been computed for the 7 most stable conformers of tyramine. In other cases,23,24 where the assignments are clearer, the B3LYP predicted frequencies have been found to reproduce better the experimental spectra than the MP2 ones, and therefore those calculations have been used to make the assignments also in the case of tyramine. In the case of ref. 7 the IR spectra of different conformers are obtained by IR dip spectroscopy or Laser Induced Fluorescence and the conformers are collected in groups based on the spectral similarities. The conformers belonging to the G or A groups show similarities in the CH out of plane bending in the benzene ring (around 550-1000 cm-1) while grouping a and g are based on the C-H methylene stretching region (2800-3000 cm-1). The main features of the observed spectra are well reproduced by B3LYP/aug-cc-pVDZ calculations. In our case the overlap of bands does not allow a straightforward assignment to a single conformation, nevertheless information on populations and populations changes can be extracted from the spectra. In this case the most informative region of the spectra is the fingerprint region. A detailed comparison of the B3LYP harmonic and anharmonic predictions will be provided in the next sections; shortly, the predicted harmonic frequencies, scaled by a common factor of 0.978, have been found to reproduce better the experimental spectra in the fingerprint region, and have been used to make the detailed assignment in this region of the spectra, while the anharmonic predictions have been very effective in reproducing the position of the OH stretching vibration for the various conformers.25 The experimental spectra of TYR isolated in argon and xenon matrices are shown in Figure 5. The simulated spectrum obtained by adding the frequency scaled calculated harmonic spectra of conformers G’g’ and Ag in a 2.4:1 ratio is also shown in this Figure. The sum of the predicted spectra of these two conformers is the minimum set of data necessary to satisfactorily reproduce 8 ACS Paragon Plus Environment

Page 9 of 37

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 recorded spectra. This result is compatible with the expected presence in the matrices of both G’g’ and Ag forms (vide supra). The small differences between the spectra obtained in Ar and Xe matrices are due to the different trapping sites which result in a slightly different splitting pattern. In addition, the different polarizability of the two matrices may induce a different stabilization order of the trapped conformers. The experimentally measured frequencies and their comparison with the B3LYP predictions, together with the PED analysis, for the G’g’ and Ag conformers are reported in Table 2. The calculated frequencies and PEDs for the remaining seven conformers are listed in Tables S2 and S3 (Supporting Information). A more detailed discussion about the vibrational assignments proposed for the different regions of the spectrum is given below.

(i) 3700–2800 cm−1 region (νOH, νNH and νCH stretching modes) The most intense band in this region is due to νOH, which appears as an intense sharp band at 3631 and 3617 cm−1 in Ar and Xe matrices, respectively (see figure 8). Since the OH group of tyramine is not involved in any intramolecular interaction, the νOH band appears at almost the same position as in matrix-isolated phenol.26 In argon, the band shows different components, possibly due to matrix site splitting effects (similarly as it is for phenol), in addition to the presence of different conformers of tyramine in the matrix. As can be seen in Table 2, the anharmonic predictions for the νOH mode are much closer to the observed frequencies and, in particular, the predicted vibrational shifts of the G’g’ and Ag conformers are in good agreement with the separation of the two observed components in the xenon matrix spectra. The intensity of the νNH2 bands is predicted to be very low and these bands are hardly discernible in the measured spectra. The calculations do not predict a strong effect on the νNH2 frequencies caused by the intramolecular N-H···π hydrogen bond. The shifts in the calculated frequencies for the conformers bearing this interaction in relation to the Gg’ and G’g conformers, in which this interaction is absent, are within a few wavenumbers. The proposed assignments for these modes are then tentative, due to their low intensity. The νCH modes can be roughly separated into two groups, i.e. the aromatic CH stretching, giving rise to bands above 3000 cm−1, and the aliphatic CH2 stretching vibrations, expected to absorb below 3000 cm−1. Although the theoretical predictions agree fairly well with the 9 ACS Paragon Plus Environment


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 10 of 37

experimental findings, these modes cannot distinguish between the different conformations because this region of the spectrum is very congested. (ii) 1850–1000 cm−1 region The main vibrational modes observed in this region, though giving rise to very intense bands, are not specific for distinguishing the various conformers. Nevertheless, the general agreement between the calculated (particularly in the case of the scaled harmonic calculations) and experimental data is rather good. As for νOH, these bands appear with several components (especially in Ar matrix), which can be ascribed to matrix site splitting and presence of different conformers in the matrix. The band exhibiting the largest peak intensity in the argon matrix spectrum, observed at 1518 cm−1 (1514 cm−1 in Xe), originates from the all-in-phase in-plane bending CH vibration of the aromatic ring and is predicted at very close values for all conformers. The νCO mode is observed at ca. 1260 (Ar) and 1254 (Xe) cm−1, also in good agreement with the theoretical predictions. The intense band at 1170 cm−1 (1169 cm−1 in Xe) is assigned to the δOH in-plane bending mode. Compared to these three bands, all the remaining features observed in this spectral region are, in agreement with the theoretical predictions, very much less intense (see Table 2 for detailed assignments). (iii) 1000–400 cm−1 region Unlike the other regions of the spectrum, the differences of the calculated frequencies for the various conformers in this spectral region are adequate to be resolved and allow tentative assignment of bands to different conformers. Figure 6 shows a closer view of the 900-750 cm−1 region of the experimental spectra, and compares it with the predicted spectra of the G’g’ and Ag conformers. To help in the visual comparison of the calculated and measured spectra, the relative intensity of the calculated bands has been scaled accordingly with the theoretical predictions. The most intense absorption in this region is due to the NH2 wagging mode, predicted at ca. 800 cm−1. The band at 767 cm−1 is clearly distinguishable from the other components and shall be ascribed to the Ag’/Ag conformers exclusively. In turn, the bands at ca. 825 cm−1 are also clearly separated from the others, and are distinctive of the Gg/G’g’ conformers. The bands between 775 and ca. 815 cm−1, on the other hand, are more difficult to doubtlessly assign based only on the comparison of the experimental spectra under discussion and the calculated data. 10 ACS Paragon Plus Environment

Page 11 of 37

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


However, as will be shown in the next sections, annealing and UV-irradiation experiments in xenon matrix allowed us also to make additional safe detailed band assignments to individual conformers in this range (see Table 2). It is interesting to compare tyramine amino group vibrations with those of the model compound n-propylamine. The latter compound, like tyramine, may also form two types of conformers around the CCCN backbone, gauche and anti. Like in tyramine, the wagging vibrations of the gauche forms of matrix-isolated n-propylamine were found at higher wavenumbers than those due to anti-forms.27 Unlike tyramine, n-propylamine cannot establish intramolecular H-bond interactions. Interestingly, the gauche forms of n-propylamine have their NH2 wagging modes at lower frequencies than in gauche-tyramine. The increase of frequency of NH2 wagging modes usually occurs when the corresponding groups are involved in H-bond interactions.16 The observed increase of NH2 wagging frequency in tyramine, comparing to n-propylamine, is consistent with the presence of a weak stabilizing interaction in gauchetyramine. CH out-of-plane (oop) vibrations in phenol produce a very strong IR-band around 750 −1

cm . Due to proximity with NH2 wagging, CH oop vibrations could also contribute around 800 cm−1 in tyramine, along with NH2 wag. This may contribute to the complexity of observed spectrum around 800 cm−1. We refer also to spectra of a similar molecule, propylbenzene, isolated in argon matrix. Instead of NH2 group, propylbenzene has a CH3 group at the same position as in tyramine. It shows a strong infrared band around 748 cm−1, and nothing particularly strong around 800 cm−1. This can be a proof that in our case the bands around 800 cm−1 are due to the amino group wagging modes. Another interesting analogy with tyramine, is that both gauche and anti conformers of propylbenzene were trapped in argon matrices at 12 K, without any conformational cooling. This should reinforce our interpretation, confirming that it should be possible to trap also gauche (Gx) and anti (Ax) conformers of tyramine. b. Annealing of the matrices Annealing experiments were performed in both argon and xenon matrices to help in the interpretation of the spectra, in particular to confirm the presence of different conformers in the matrices. Due to its physical characteristics, annealing of the argon matrix (deposited at 15 K) was 11 ACS Paragon Plus Environment


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 12 of 37

carried out only up to 30 K, to avoid destruction of the matrix upon further heating. The xenon matrix was deposited at 30 K and annealed up to 60 K without remarkable alteration of the baseline, which is an indication of the physical state of the matrix, and without signals of tyramine aggregation in the recorded spectrum. In the argon matrix, annealing did not produce any band intensity changes, while in the experiments performed with xenon, changes in the relative intensity of the components of the most intense bands started to occur between 40 and 42 K (Figures 7-9). A similar annealing behavior was observed for rearrangements of conformers of trimethylphosphate (TMP).28 The rotations of methoxy groups in TMP (similar in size to the ethylamino fragment in tyramine) have calculated barriers of ca. 10 kJ mol−1 and internal conversion around these barriers in TMP in xenon matrices was observed at temperatures above 40 K. Figure 7 shows the observed changes in the region of the NH2 wagging vibration, compared with the frequency scaled harmonic predicted spectra for the G’g’ and Ag conformers. Relative intensities are scaled according to the calculated relative population of the conformers. Other spectral regions, including the OH stretching region, are shown in Figures 8 and 9. In these cases comparison with anharmonic calculations is also presented. Unfortunately, the anharmonic prediction in the NH2 wagging region is still very shifted from the experimental values and does not aid the interpretation. Indeed, the amino group wagging vibrations are known to be problematic for anharmonic simulation.29,30 The large deviation from the experimental value for this mode is in agreement with anharmonic calculations performed on glycine.31 The components that increase upon annealing fit well the calculated spectrum of the anti (Ax) conformers. These conformers are predicted by the calculations to be less stable than the gauche forms. Hence, based on the calculated relative energies, the reverse conversion could be expected, since annealing of the matrix must lead the conformational system to approach the low temperature equilibrium, i.e., as a general rule, to conversion of the less stable conformers into the more stable ones. These results suggest that, at least in the highly polarizable xenon matrix, the anti conformers become more stable than the gauche forms. Nevertheless, the putative conflict between the experimental and calculated data may also be a result of insufficient accuracy of calculations to predict either relative conformational energies or vibrational spectra. Conclusions about the relative stability of the two conformers in the argon matrix could not be extracted due to the lack of observation of any spectral changes in the annealing experiments 12 ACS Paragon Plus Environment

Page 13 of 37

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


performed in this matrix. However, we can assume that the barrier is large enough to prevent the conformational interconversion, for those forms that are trapped in the matrix, in the temperature range accessible for experiments with argon. c. Narrow band UV irradiation experiments Freshly deposited argon and xenon matrices of tyramine were irradiated using narrowband tunable UV laser light, through a quartz window. In gas phase, tyramine exhibits an absorption band in the UV spectrum with origin at about 282 nm.7 In the performed experiments, irradiation of the deposited samples was then undertaken starting at 290 nm and ending at 270 nm, with steps of 1 nm. The irradiation time for each step was 2 minutes. Taking into account previous results in related compounds,32-35 for matrix-isolated tyramine the expected photochemical processes should be mainly photoisomerizations (Scheme 2). Depending on the reaction site of the molecule, isomerization processes in tyramine may lead to three classes of products: isomerization of the ethylamine side chain can deplete one of the gauche or anti conformers in favor of the other one (path A in Scheme 2); isomerization of the phenol ring can lead either to production of a Dewar isomer (path B) or an open ring conjugated ketene (path C). The formation of a Dewar isomer could be easily excluded, since no new bands were observed in the OH stretching region of the spectra of the irradiated matrices (the OH stretching vibrations for the Dewar forms are calculated at least 16 cm-1 away from the ones of the other conformers of TYR). On the other hand, new bands in the region 2100-2150 cm-1 started to appear at the early stages of irradiation at 290 nm, in the spectra of both argon and xenon irradiated tyramine matrices. The presence of these bands is an unequivocal indication of the formation of a ketene species as can be deduced from the infrared spectra of ketenes reported in the literature.36,37 For ketenes, the shapes of absorption bands of the same species are different in different environment, and even in the same environment; the discussion of this effect is not trivial, and requires a dedicated study which is beyond the scope of this paper. Ketenes have been found to be main photoproducts resulting from the opening of aromatic 6-membered rings bearing oxygen containing substituents.32-35 The ketene marker bands in the region 2100-2150 cm-1 continuously increased along the performed irradiation experiments. Figure 10 shows this region of the spectra 13 ACS Paragon Plus Environment


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 37

obtained during a series of irradiations of the argon and xenon matrices. As can be seen from the total intensity of the new bands, the rate of the ring-opening reaction is quite low. Irradiation of the xenon matrix at 287 nm has shown, in addition to the bands associated to the ketene species, results similar to those obtained during the annealing experiments. Some components of the most intense bands, assigned to the Ag conformer, increased, while bands ascribed to the G’g’ form, decreased. It can then be concluded that in the xenon matrix the photoisomerization of gauche form into anti form occurs simultaneously with the ring-opening reaction. On the other hand, irradiation of the argon matrix did not lead to a clear observation of the gauche→anti isomerization process, with all the components of the recorded bands being observed to decrease simultaneously. In argon matrix, the single reliably observed photoprocess was then the formation of the isomeric open ring conjugated ketene.

4. Conclusions

The infrared spectra of tyramine isolated in low temperature argon and xenon matrices were described for the first time. The spectra of the as-deposited matrices are compatible with the presence of at least two types of conformers: one has been assigned to the most stable gauche forms about the CH2−CH2 bond (G’g’+Gg) and the other one to the most stable anti type conformers (Ag’+Ag). The most stable conformer in gas phase (Gg) has a folded structure stabilized by the presence of a weak N−H…π hydrogen bond interaction. Observation of both gauche and anti forms in the as-deposited matrices is consistent with the calculated relative values for the different conformational isomerization barriers, with that associated with the rotation around the CH2−CH2 bond, involving rotation of heavy atoms, being predicted to be the highest among all interconversion barriers. Annealing experiments performed in the xenon matrix, compared with harmonic and anharmonic predictions for the most stable conformers, seem to confirm the presence in the deposited matrix of the gauche and anti conformers and the occurrence of the gauche→anti isomerization reaction upon temperature increase. This last result indicates a change in the relative energies of the two experimentally relevant conformers upon isolation in the polarizable xenon matrix, compared to gas phase, with the more polar anti form being stabilized in the matrix environment.


Page 15 of 37

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


New bands formed after UV irradiation of the matrix-isolated tyramine (both in argon and xenon matrices) are consistent with the formation of a conjugated ketene species, generated via ring-opening reaction, which has already been showed to occur for similar 6-membered ring compounds bearing an oxygen-containing substituent.32-35 Concomitant occurrence of the gauche→anti isomerization reaction was observed in the xenon matrix. For both matrices, no evidence of photoinduced formation of Dewar isomers of tyramine has been found.

Acknowledgments. The research leading to these results has received funding from the European Community’s Seventh Frame work Programme under Grant Agreement No. 228334, and the Portuguese “Fundação para a Ciência e a Tecnologia” (FCT) Projects PTDC/QUIQUI/111879/2009

and

PTDC/QUI-QUI/118078/2010,

FCOMP-01-0124-FEDER-021082,

cofunded by QREN-COMPETE-UE. B.M.G. acknowledges FCT for the postdoctoral grant No. SFRH/BPD/44689/2008. Supporting Information: Internal coordinates used in the normal modes analyses for tyramine (Table S1) as well as calculated harmonic and anharmonic wavenumbers, absolute intensities and potential energy distribution for seven tyramine conformers at the B3LYP/6-311++G(d,p) level of theory (Tables S2 and S3) and complete reference 10 are given as Supporting Information. This material is available free of charge via the Internet at http://pubs.acs.org.



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 16 of 37

References 1. 2. 3. 4. 5.

6. 7.

8.

9. 10.

11.

12.

13.

14.

15.

Karovičová, J.; Kohajdová, Z. Biogenic Amines in Food. Chem. Pap. 2005, 59, 70-79. Silla Santos, M. H. Biogenic Amines: Their Importance in Foods. Int. J. Food Microbiol. 1996, 29, 213-231. Snyder, S. H. Drugs and the Brain; Scientific American Library: New York, 1999. Mycek, M. J.; Harvey, R. A.; Champe, P. C. Farmacologia; Zanichelli: Bologna, 2000. Melandri, S.; Maris, A. Intramolecular Hydrogen Bonds and Conformational Properties of Biogenic Amines: A Free-jet Microwave Study of Tyramine. Phys. Chem. Chem. Phys. 2004, 6, 2863-2866. Yoon, I.; Seo, K.; Lee, S.; Lee, Y.; Kim, B. Conformational Study of Tyramine and its Water Clusters by Laser Spectroscopy. J. Phys. Chem. A 2007, 111, 1800-1807. Makara, K.; Misawa, K.; Miyazaki, M.; Mitsuda, H.; Ishiuchi, S.; Fujii, M. Vibrational Signature of the Conformers in Tyramine Studied by IR Dip and Dispersed Fluorescence Spectroscopies. J. Phys. Chem. A 2008, 112, 13463-13469. Martinez, S. J.; Alfano, J. C.; Levy, D. H. The Electronic Spectroscopy of Tyrosine and Phenylalanine Analogs in a Supersonic Jet: Basic Analogs. J. Mol. Spectrosc. 1993, 158, 82-92. Robertson, E. G.; Simons, J. P.; Mons, M. Structural and Vibrational Assignment of pmethoxyphenyl -ethylamine Conformers. J. Phys. Chem. A 2001, 105, 9990-9992. 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, et al. Gaussian 03, revision C.02; Gaussian, Inc.: Wallingford, CT, 2004. Complete reference is given as Supporting Information Schachtschneider, J. H.; Mortimer, F. S. Vibrational Analysis of Polyatomic Molecules. VI. FORTRAN IV Programs for Solving the Vibrational Secular Equation and for the LeastSquares Refinement of Force Constants. Project No. 31450. Structural Interpretation of Spectra, Emeryville, CA, 1969. Pulay, P.; Fogarasi, G.; Pang, F.; Boggs, J. E. Systematic Ab Initio Gradient Calculation of Molecular Geometries, Force Constants and Dipole-Moment Derivatives. J. Am. Chem. Soc. 1979, 101, 2550-2560. Carney, J. R.; Zwier, T. S. The Infrared and Ultraviolet Spectra of Individual Conformational Isomers of Biomolecules: Tryptamine. J. Phys. Chem. A 2000, 104, 86778688. Riley, K. E.; Hobza, P. Assessment of the MP2 Method, along with Several Basis Sets, for the Computation of Interaction Energies of Biologically Relevant Hydrogen Bonded and Dispersion Bound Complexes. J. Phys. Chem. A 2007, 111, 8257-8263. Frija, L. M. T.; Reva, I.; Ismael, A.; Coelho, D. V.; Fausto, R.; Cristiano, M. L. S. Sigmatropic Rearrangements in 5-allyloxytetrazoles. Org. Biomol. Chem. 2011, 9, 60406054.


Page 17 of 37

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


16.

17.

18.

19.

20.

21.

22.

23. 24.

25. 26.

27.

28. 29.

Reva, I. D.; Plokhotnichenko, A. M.; Stepanian, S. G.; Ivanov, A. Yu.; Radchenko, E. D.; Sheina, G. G.; Blagoi, Yu. P. The Rotamerization of Conformers of Glycine Isolated in Inert Gas Matrices. An Infrared Spectroscopic Study. Chem. Phys. Lett. 1995, 232, 141-148. Erratum. The Rotamerization of Conformers of Glycine Isolated in Inert Gas Matrices. An Infrared Spectroscopic Study (Chem. Phys. Letters 232 (1995) 141). Chem. Phys. Lett. 1995, 235, 617. Gómez-Zavaglia, A.; Fausto, R. Conformational Study of Sarcosine as Probed by MatrixIsolation FT-IR Spectroscopy and Molecular Orbital Calculations. Vib. Spectrosc. 2003, 33, 105-126. Reva, I. D.; Ilieva S. V.; Fausto, R. Conformational Isomerism in Methyl Cyanoacetate: A Combined Matrix-Isolation Infrared Spectroscopy and Molecular Orbital Study. Phys. Chem. Chem. Phys. 2001, 3, 4235-4241. Reva, I. D.; Stepanian, S. G.; Adamowicz, L.; Fausto, R. Missing Conformers. Comparative Study of Conformational Cooling in Cyanoacetic Acid and Methyl Cyanoacetate Isolated in Low Temperature Inert Gas Matrixes. Chem. Phys. Lett. 2003, 374, 631-638. Nunes, C. M., Lapinski, L., Fausto, R. and Reva, I. Near-IR Laser Generation of a HighEnergy Conformer of L-alanine and the Mechanism of its Decay in a Low-Temperature Nitrogen Matrix J. Chem. Phys., 2013, 138, 125101. Reva, I., Nowak, M. J., Lapinski, L., Fausto, R. Spontaneous Tunneling and Near-InfraredInduced Interconversion Between the Amino-Hydroxy Conformers of Cytosine J. Chem. Phys. 2012, 136, 064511-8. Reva, I., Lapinski, L., Fausto, R. Infrared Spectra of Methyl Isocyanate Isolated in Ar, Xe and N2 Matrices. J. Mol. Struct. 2010, 976, 333-341. Breda, S., Reva, I., Fausto, R. Molecular Structure and Vibrational Spectra of 2(5H)-Furanone and 2(5H)-Thiophenone Isolated in Low Temperature Inert Matrix. J. Mol. Struct. 2008, 887, 75-86. Sharma, A., Reva, I., Fausto, R. Conformational Switching Induced by Near-Infrared Laser Irradiation. J. Am. Chem. Soc. 2009, 131, 8752-8753. Giuliano, B. M.; Reva, I.; Lapinski, L.; Fausto, R. Infrared Spectra and Ultraviolet-Tunable Laser Induced Photochemistry of Matrix-Isolated Phenol and Phenol-d5. J. Chem. Phys. 2012, 136, 024505. Sato, N., Hamada, Y., Tsuboi, M. Vibrational and Conformational Analysis of nPropylamine by Means of i.r. Spectroscopy and Ab Initio MO Calculations. Spectrochim. Acta A, 1987, 43, 943-954. Reva, I., Simão, A., Fausto, R. Conformational Properties of Trimethyl Phosphate Monomer. Chem. Phys. Lett. 2005, 406, 126-136. Shin-ya, K., Takahashi, O., Katsumoto, Y., Ohno, K. Intramolecular CH⋯π and CH⋯O Interactions in the Conformational Stability of Benzyl Methyl Ether Studied by MatrixIsolation infrared Spectroscopy and Theoretical Calculations. J. Mol. Struct. 2007, 827, 155-164.



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

30.

31.

32.

33.

34. 35.

36.

37.

Page 18 of 37

Wojciechowski, P., Helios, K., Michalska, D. Theoretical Anharmonic Raman and Infrared Spectra with Vibrational Assignments for Monoﬂuoroaniline Isomers. Vib. Spectrosc. 2011, 57, 126-134. Bludský, O.; Chocholoušová, J.; Vacek, J.; Huisken, F.; Hobza, P. Anharmonic Treatment of the Lowest-Energy Conformers of Glycine: A Theoretical Study. J. Chem. Phys. 2000, 113, 4629-7. Breda, S.; Lapinski, L.; Reva, I.; Fausto, R. 4,6-Dimethyl-α-pyrone: a Matrix Isolation Study of the Photochemical Generation of Conjugated Ketene, Dewar Valence Isomer and 1,3-Dimethyl-Cyclobutadiene. J. Photochem. Photobiol. A-Chem. 2004, 162, 139-151. Kuş, N.; Breda, S.; Reva, I. D.; Tasal, E.; Ogretir, C.; Fausto, R. FTIR Spectroscopic and Theoretical Study of the Photochemistry of Matrix-Isolated Coumarin. Photochem. Photobiol. 2007, 83, 1237-1253. Erratum: FTIR Spectroscopic and Theoretical Study of the Photochemistry of Matrixisolated Coumarin Photochem. Photobiol. 2007, 83, 1541-1542. Kuş, N.; Reva, I.; Bayarı, S.; Fausto, R. In situ Photoproduction of Dichlorodibenzo-pDioxin From Non-Ionic Triclosan Isolated in Solid Argon. J. Mol. Struct. 2012, 1007, 8894. Kappe C. O., Wong M. W., Wentrup C. Acetylketene: Conformational Isomerism and Photochemistry. Matrix Isolation Infrared and Ab Initio Studies J. Org. Chem. 1996, 60, 1686-1695. Zuhse R. H., M. W., Wentrup C. Photochemistry of Deuterated Acetylketenes: Matrix Isolation Infrared Spectroscopic and ab Initio Studies J. Phys. Chem. 1996, 100, 39173922.


Page 19 of 37

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


Table 1. Zero-point corrected relative energies (∆EZPE/ kJ mol−1), dipole moments (µ/ Debye), relative Gibbs energies at 298 K (∆G298/ kJ mol−1) and relative populations (P298, %; at T= 298 K) of the nine different conformers of TYR, calculated at the B3LYP and MP2 levels of approximation with the 6-311++G(d,p) basis set. Conformer a

B3LYP Gg G’g’ Ga G’a Aa Ag’ Ag Gg’ G’g

c

0.00 0.26 0.52 0.42 0.17 0.63 0.73 5.39 5.60

∆G298

P298b

MP2

MP2

MP2

0.67 2.65 2.23 1.67 2.02 0.90 2.57 2.70 2.72

e

23.9 20.8 11.2 11.3 9.1 12.9 6.7 1.9 2.2

µ

∆EZPE MP2

B3LYP

d

0.00 0.53 1.81 1.95 5.25 4.77 5.16 7.60 7.73

0.65 2.47 2.09 1.57 1.86 0.87 2.46 2.37 2.37

a

Conformers structures are shown in Figure 3. Estimated using the calculated ∆G298 values. c Absolute energy: –441.397289 Eh. d Absolute energy: –440.109036 Eh. e Absolute Gibbs energy: –440.144844 Eh. b


0.00 0.34 1.88 1.85 2.40 1.52 3.14 6.26 5.87


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 20 of 37

Page 21 of 37

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 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 22 of 37

Page 23 of 37

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 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 24 of 37

Captions to the Figures and Scheme Figure 1. Relaxed B3LYP/6-311++G(d,p) potential energy scans around the C8−C7−C4−C3 (left) and C2−C1−O−H21 (right) dihedral angles in tyramine. Along the Caryl−CH2 scan (left), the optimized C2−C1−O−H, N−C8−C7−C4 and Lp-N−C8−C7 dihedral angles stayed around cis, anti and gauche orientations, respectively. Along the Caryl−OH scan (right), the C8−C7−C4−C3, N−C8−C7−C4 and Lp-N−C8−C7 dihedral angles stayed around perpendicular, anti and gauche orientations, respectively. “Lp” stands for lone pair of electrons of the nitrogen atom, and its orientation is defined to be in the plane bisecting the C8−N9−H19 and C8−N9−H20 planes. Center: optimized structure of Gg conformer with atom numbering. Figure 2. Relaxed potential energy surface of tyramine for rotation around CH2−CH2 and CH2−NH2 bonds calculated at the B3LYP/6-311++G(d,p) level. The names of conformers correspond to the combinations of letters (G/A/G’) and (g/a/g’) shown at the edges of the surface. The global energy minimum is Gg and its energy was chosen as the relative zero level. The isoenergy levels are presented by continuous thin lines and are spaced by 3 kJ mol−1. Two additional isoenergy levels are shown at 2 (bold line) and 16 (dotted line) kJ mol−1. “Lp” stands for “lone pair” of the nitrogen atom, and its orientation was chosen to be in the plane bisecting the two CCNH planes. Figure 3. The nine conformers of tyramine and their optimized geometries. Note that positions of the conformers in this figure correspond to their positions on the two-dimensional potential energy surface shown in Figure 2. The dashed lines specify the amino group hydrogen atoms involved in the stabilizing interaction with the aromatic ring. Figure 4. Calculated Boltzmann relative populations of G- and A-type conformers as a function of temperature, according to the MP2 calculated Gibbs energies. Figure 5. Experimental infrared spectra of TYR isolated in freshly deposited Ar (a) and Xe (b) matrices, at 15 and 30 K, respectively; (c) Simulated theoretical spectrum which includes population-weighted contributions of the G’g’ (71%) and Ag (29%) conformers. Before simulation, all frequencies were calculated at the B3LYP/6-311++G(d,p) level of theory and scaled by a factor of 0.978. Then, the absorption bands were simulated using Lorentzian functions centered at the calculated (scaled) frequencies and with the full width at half maximum (fwhm) equal to 2 cm−1. Figure 6. Part of the experimental infrared spectra of TYR isolated in freshly deposited Ar (a) and Xe (b) matrices, at 15 and 30 K, respectively; (c) Simulated infrared spectra of the G’g’ (black line) and Ag (red line) conformers. For the details of simulation see caption of Figure 5. Figure 7.


Page 25 of 37

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


Experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) simulated spectra for the G’g’ (black line) and Ag (red line) conformers. For the details of simulation see caption of Figure 5. Figure 8. (a) OH stretching region of the experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) anharmonic prediction of the νOH for the G’g’ (black line) and Ag (red line) conformers. Absorption bands were simulated using Lorentzian profiles centered at the B3LYP calculated (non-scaled) anharmonic frequency and having fwhm = 2 cm−1. Infrared intensities were taken from the harmonic B3LYP calculation and then scaled by 0.71 and 0.29, respectively according to the calculated relative population of the conformers. Figure 9. (a) Experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) Simulated anharmonic spectra for the G’g’ (black line) and Ag (red line) conformers. For details of the simulation see caption of Figure 8. Figure 10. Region of the C=C=O antisymmetric stretching vibration in the infrared spectra of TYR isolated in an argon (a) and xenon (b) matrices, in the freshly deposited matrix (black line) and after UV irradiations. Scheme 1. Relative energies (∆E/ kJ mol−1) of transition states between the nine conformers of TYR calculated at the MP2/6-311++G(d,p) level of theory. Scheme 2. Possible reaction pathways after UV irradiation of matrix isolated TYR.



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

Relative Energy /

Page 26 of 37

0

90

15 12 9 6 3 0 180 3 6 9 12 15

Relative Energy / kJ mol −1

kJ mol −1

Ag'

Ag

270

0

Ag

Ag'

15 12 9 6 3 90 0 3 6 9 12 15

270

180

C8-C7-C4-C3 dihedral angle / degree

C2-C1-O-H21 dihedral angle / degree

Figure 1. Relaxed B3LYP/6-311++G(d,p) potential energy scans around the C8−C7−C4−C3 (left) and C2−C1−O−H21 (right) dihedral angles in tyramine. Along the Caryl−CH2 scan (left), the optimized C2−C1−O−H, N−C8−C7−C4 and Lp-N−C8−C7 dihedral angles stayed around cis, anti and gauche orientations, respectively. Along the Caryl−OH scan (right), the C8−C7−C4−C3, N−C8−C7−C4 and Lp-N−C8−C7 dihedral angles stayed around perpendicular, anti and gauche orientations, respectively. “Lp” stands for lone pair of electrons of the nitrogen atom, and its orientation is defined to be in the plane bisecting the C8−N9−H19 and C8−N9−H20 planes. Center: optimized structure of Gg conformer with atom numbering.


Page 27 of 37

G

360

A

30

16

24

3

12

12

9

21 27

16

24 g' 27

6

6

120

3

16

27

3

180

30

9

21

16

240

9 18

18

12 27

G' 6

18 9

300 CCN(Lp) dihedral angle / degree

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


12

27

12

24 18

16 6

21

a

3

3 18

9

9

24

27

24 27

18

16

9 21

60

3

24 30

16 18

6

3

9

16 6

12 21

18

g

30

0 0

60

120 180 240 300 CCCN dihedral angle / degree

360

Figure 2. Relaxed potential energy surface of tyramine for rotation around CH2−CH2 and CH2−NH2 bonds calculated at the B3LYP/6-311++G(d,p) level. The names of conformers correspond to the combinations of letters (G/A/G’) and (g/a/g’) shown at the edges of the surface. The global energy minimum is Gg and its energy was chosen as the relative zero level. The isoenergy levels are presented by continuous thin lines and are spaced by 3 kJ mol−1. Two additional isoenergy levels are shown at 2 (bold line) and 16 (dotted line) kJ mol−1. “Lp” stands for “lone pair” of the nitrogen atom, and its orientation was chosen to be in the plane bisecting the two CCNH planes.



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 28 of 37

Gg’

Ag’

G’g’

Ga

Aa

G’a

Gg

Ag

G’g

Figure 3. The nine conformers of tyramine and their optimized geometries. Note that positions of the conformers in this figure correspond to their positions on the two-dimensional potential energy surface shown in Figure 2. The dashed lines specify the amino group hydrogen atoms involved in the stabilizing interaction with the aromatic ring.


Page 29 of 37

G-type conformers A-type conformers

100

80 Relative Population/ %

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


60

40

20

0 0

50

100

150

200

250

300

Temperature/ K

Figure 4 Calculated Boltzmann relative populations of G- and A-type conformers as a function of temperature, according to the MP2 calculated Gibbs energies.



(a)

Absorbance

0,35

Relative Intensity

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 30 of 37

0,00

(b) 0,35

0,00

(c) 35

0 3600

3200

2800 1600

1200

800

Wavenumber / cm−1

Figure 5. Experimental infrared spectra of TYR isolated in freshly deposited Ar (a) and Xe (b) matrices, at 15 and 30 K, respectively; (c) Simulated theoretical spectrum which includes population-weighted contributions of the G’g’ (71%) and Ag (29%) conformers. Before simulation, all frequencies were calculated at the B3LYP/6-311++G(d,p) level of theory and scaled by a factor of 0.978. Then, the absorption bands were simulated using Lorentzian functions centered at the calculated (scaled) frequencies and with the full width at half maximum (fwhm) equal to 2 cm−1.


Page 31 of 37

(a)

Absorbance

0,2

0,0

(b) 0,2

Relative Intensity

0,0

(c)

G'g' 20

Ag 0 875

850

825

800

775

750

Wavenumber / cm−1

Absorbance

Figure 6. Part of the experimental infrared spectra of TYR isolated in freshly deposited Ar (a) and Xe (b) matrices, at 15 and 30 K, respectively; (c) Simulated infrared spectra of the G’g’ (black line) and Ag (red line) conformers. For the details of simulation see caption of Figure 5.

(a)

0,2

0,1

0,0

Relative Intensity

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


G'g' Ag

20

10

G'g' G'g'

(b)

G'g'

Ag

Ag

0 875

850

825

800

775

750

Wavenumber / cm−1

Figure 7. Experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) simulated spectra for the G’g’ (black line) and Ag (red line) conformers. For the details of simulation see caption of Figure 5.



Absorbance

(a) 0,2

0,1

0,0

Relative Intensity

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 32 of 37

G'g'

(b)

20

10

Ag

0 3650

3625

3600

Wavenumber / cm

−1

Figure 8. (a) OH stretching region of the experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) anharmonic prediction of the νOH for the G’g’ (black line) and Ag (red line) conformers. Absorption bands were simulated using Lorentzian profiles centered at the B3LYP calculated (non-scaled) anharmonic frequency and having fwhm = 2 cm−1. Infrared intensities were taken from the harmonic B3LYP calculation and then scaled by 0.71 and 0.29, respectively according to the calculated relative population of the conformers.


Page 33 of 37

Absorbance

(a) 0,2 0,1 0,0

Relative Intensity

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


G'g'

G'g'

20

(b)

Ag Ag

10

G'g'

0 1250

1200

1175

1150

1125

−1

Wavenumber / cm

Figure 9. (a) Experimental infrared spectra of TYR in a freshly deposited Xe matrix at 30 K (black line); and after annealing at 42 K (red line); (b) Simulated anharmonic spectra for the G’g’ (black line) and Ag (red line) conformers. For details of the simulation see caption of Figure 8.



0.02 283 nm 284 nm 285 nm

(a)

0.01

Absorbance

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 34 of 37

0.00

0.02

284 nm 286 nm 287 nm

(b)

0.01

0.00 2200

2150

2100

Wavenumber / cm

2050 −1

Figure 10. Region of the C=C=O antisymmetric stretching vibration in the infrared spectra of TYR isolated in an argon (a) and xenon (b) matrices, in the freshly deposited matrix (black line) and after UV irradiations.


Page 35 of 37

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


Scheme 1. Relative energies (∆E/ kJ mol−1) of transition states between the nine conformers of TYR calculated at the MP2/6-311++G(d,p) level of theory.a

26.93

11.13

12.34

11.72

Gg’ 8.61 [2.70]

Ag’ 5.86 [0.90]

G’g’ 0.49 [2.65]

21.06

15.94

22.34

26.37

Ga 1.67 [2.23]

22.14

15.99

23.29

Aa 6.23 [2.02]

14.44

23.29

G’a 1.74 [1.67]

14.12

16.16

15.86

Gg 0.00b [0.67]

Ag 5.95 [2.57]

G’g 9.16 [2.72]

11.13

22.08

26.93

21.27

12.34

a

22.34

26.37

11.72

Data enclosed in squares represent names, relative energies (kJ mol−1), and net dipole moments (Debye, in brackets) of all the minima, positioned in the same way as on the map shown in Figure 2. Data between squares represent the energies of the respective transition states. Values on the perimeter show the transition states energies between the minima on the opposite edges and, to emphasize it, they are repeated on right/left and top/bottom. b Absolute energy: –440.287276 Eh.



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

NH2

NH2

C H O

OH

A

B NH2

H2N

OH OH

Scheme 2. Possible reaction pathways after UV irradiation of matrix isolated TYR.


Page 36 of 37

Page 37 of 37

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


Table of Contents Image


Conformational Space and Photochemistry of ... - ACS Publications

Recommend Documents