Gas phase basicities of polyfunctional molecules. Part 5: Non-aromatic sp2 nitrogen containing compounds.

This paper constitutes the fifth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435, Part 4: Carbonyl as basic site, Mass Spectrom Rev 2015, 34:493-534). This part is devoted to non-aromatic molecules characterized by a lone pair located on a sp2 nitrogen atom, it embraces functional groups such as imines, amidines, guanidines, diazenes, hydrazines, oximes, and phosphazenes. Specific examples are examined under five major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) imines, hydrazones, and oximes are considered. A second chapter describes the protonation energetic of aliphatic, conjugated, or cyclic amidines. Guanidines, polyguanides, and biomolecules containing guanidine were examined in the third chapter. A fourth chapter describes the particular case of the phosphazene molecules. Finally, diazenes and azides were considered in the last chapter. Experimental data were re-evaluated according to the presently adopted basicity scale, i.e., PA(NH3 ) = 853.6 kJ/mol, GB (NH3 ) = 819 kJ/mol. Structural and energetic information given by G4MP2 quantum chemistry computations on typical systems are presented. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 37:139-170, 2018.

Gas phase basicities of polyfunctional molecules. Part 6: Cyanides and isocyanides.

This paper gathers structural and thermochemical informations related to the gas-phase basicity of molecules containing cyanides (nitriles) and isocyanides (isonitriles) functional groups. It constitutes the sixth part of a general review devoted to gas-phase basicities of polyfunctional compounds. A large corpus of cyanides and isocyanides molecules is examined under seven major chapters. In the first one, a rapid overview of the definitions and methods leading to gas-phase basicity, GB, proton affinity, PA, and protonation entropy, Δp S°, is given. In the same chapter, several aspects of the gas phase chemistry of protonated cyanides and isocyanides are also presented. Chapters II-VI detail the protonation energetics of aliphatic, unsaturated, and heteroatom substituted (halogens, O, S, N, P) cyanides. A seventh chapter is devoted to isocyanides. Experimental data available in the literature (120 references) were reevaluated according to the presently adopted basicity scale that is the NIST database anchored to PA(NH3 ) = 853.6 kJ/mol and GB (NH3 ) = 819 kJ/mol. In this latter source, however, several erroneous values have been identified which were corrected in the present review. Structural and energetic information given by G4MP2 quantum chemistry computations on ca. 60 typical systems are presented. The present review includes the GB, PA, and Δp S° values of ca. 110 cyanides and isocyanides, and, for selected examples, is completed by a set of computed heats of formation (Δf H°) at 0 and 298 K.

Gas-phase lithium cation affinity of glycine.

The gas-phase lithium cation binding thermochemistry of glycine has been determined theoretically by quantum chemical calculations at the G4 level and experimentally by the extended kinetic method using electrospray ionization quadrupole time-of-flight tandem mass spectrometry. The lithium cation affinity of glycine, ∆(Li)H°(298)(GLY), i.e. the∆(Li)H°(298) of the reaction GlyLi(+)→ Gly + Li(+)) given by the G4 method is equal to 241.4 kJ.mol(-1) if only the most stable conformer of glycine is considered or to 242.3 kJ.mol(-1) if the 298K equilibrium mixture of neutral conformers is included in the calculation. The ∆(Li)H°(298)(GLY) deduced from the extended kinetic method is obviously dependent on the choice of the Li(+) affinity scale, thus∆(Li)H°(298)(GLY) is equal to 228.7±0.9(2.0) kJ.mol(- 1) if anchored to the recently re-evaluated lithium cation affinity scale but shifted to 235.4±1.0 kJ.mol(-1) if G4 computed lithium cation affinities of the reference molecules is used. This difference of 6.3 kJ.mol(-1) may originate from a compression of the experimental lithium affinity scale in the high ∆(Li)H°(298) region. The entropy change associated with the reaction GlyLi(+)→Gly + Li(+) reveals a gain of approximately 15 J.mol(-) 1.K(-1) with respect to monodentate Li(+) acceptors. The origin of this excess entropy is attributed to the bidentate interaction between the Li(+) cation and both the carbonyl oxygen and the nitrogen atoms of glycine. The computed G4 Gibbs free energy,∆(Li)G°(298)(GLY) is equal to 205.3 kJ.mol(-1), a similar result, 201.0±3.4 kJ.mol(-1), is obtained from the experiment if the∆(Li)G°(298) of the reference molecules is anchored on the G4 results.

Microhydration of Protonated Nα-Acetylhistidine: A Theoretical Approach.

Extensive exploration of the potential energy surfaces of protonated Nα-acetylhistidine hydrated by 0-3 molecules of water was performed. The methodology combined hierarchical and genealogical (Darwin family tree) approaches using polarizable AMOEBA force field and M06 functional. It is demonstrated that this mixed approach allows recovering a larger number of conformers than the number recovered by using any one of the two methods alone. Hydration enthalpies of protonated Nα-acetylhistidine and of model compounds have been computed using higher theoretical methods, up to the G4MP2 procedure. Excellent agreement with experiment is observed for successive hydration of methylamonium and imidazolium cations using MP2/6-311++G(2d,2p)//M06/6-311++G(d,p) and G4MP2 methods, thereby validating the theory levels used for hydrated protonated Nα-acetylhistidine. It is found that the first hydration enthalpy of protonated Nα-acetylhistidine is ca. 10 kJ mol(-1) lower than that of imidazolium, a result explained by the local environment of the positively charged imidazolium moiety.

Gas-phase basicities of polyfunctional molecules. Part 4: Carbonyl groups as basic sites.

This article constitutes the fourth part of a general review of the gas-phase protonation thermochemistry of polyfunctional molecules (Part 1: Theory and methods, Mass Spectrom Rev 2007, 26:775-835, Part 2: Saturated basic sites, Mass Spectrom Rev 2012, 31:353-390, Part 3: Amino acids, Mass Spectrom Rev 2012, 31:391-435). This fourth part is devoted to carbonyl containing polyfunctional molecules. After a short reminder of the methods of determination of gas-phase basicity and the underlying physicochemical concepts, specific examples are examined under two major chapters. In the first one, aliphatic and unsaturated (conjugated and cyclic) ketones, diketones, ketoalcohols, and ketoethers are considered. A second chapter describes the protonation energetic of gaseous acids and derivatives including diacids, diesters, diamides, anhydrides, imides, ureas, carbamates, amino acid derivatives, and peptides. Experimental data were re-evaluated according to the presently adopted basicity scale. Structural and energetic information given by G3 and G4 quantum chemistry computations on typical systems are presented.

From the mobile proton to wandering hydride ion: mechanistic aspects of gas-phase ion chemistry.

Structural characterization of molecular species by mass spectrometry supposes the knowledge of the type of ions generated and the mechanism by which they dissociate. In this context, a need for a rationalization of electrospray ionization(+)(-) mass spectra of small molecules has been recently expressed. Similarly, at the other end of the mass scale, efforts are currently made to interpret the major fragmentation processes of protonated and deprotonated peptides and their reduced forms produced in electron capture or electron transfer experiments. Most fragmentation processes of molecular and pseudo-molecular ions produced in the ion source of a mass spectrometer may be described by a combination of several key mechanistic steps: simple bond dissociation, formation of ion-neutral complex intermediates, hydrogen atom, hydride ion or proton migrations and nucleophilic attack. Selected crucial aspects of these elementary reactions, occurring inside positively charged ions, will be recalled and illustrated by examples taken in recent mass spectrometry literature. Emphasis will be given on the protonation process and its consequence in terms of structure and energetic.

Gas-phase structures and thermochemistry of neutral histidine and its conjugated acid and base.

Extensive exploration of the conformational space of neutral, protonated and deprotonated histidine has been conducted at the G4MP2 level. Theoretical protonation and deprotonation thermochemistry as well as heats of formation of gaseous histidine and its ionized forms have been calculated at the G4 level considering either the most stable conformers or an equilibrium population of conformers at 298 K. These theoretical results were compared to evaluated experimental determinations. Recommended proton affinity and protonation entropy deduced from these comparisons are PA(His) = 980 kJ mol(-1) and ΔpS(His) ∼ 0 J mol(-1) K(-1), thus leading to a gas-phase basicity value of GB(His) = 947.5 kJ mol(-1). Similarly, gas phase acidity parameters are ΔacidH(o)(His) = 1373 kJ mol(-1), ΔacidS(His) ∼ 10 J mol(-1) K(-1) and ΔacidG(o)(His) = 1343 kJ mol(-1). Computed G4 heats of formation values are equal to -290, 265 and -451 kJ mol(-1) for gaseous neutral histidine and its protonated and deprotonated forms, respectively. The present computational data correct, and complete, previous thermochemical parameter estimates proposed for gas-phase histidine and its acido-basic properties.

A 'meta effect' in the fragmentation reactions of ionised alkyl phenols and alkyl anisoles.

The competition between benzylic cleavage (simple bond fission [SBF]) and retro-ene rearrangement (RER) from ionised ortho, meta and para RC(6) H(4) OH and RC(6) H(4) OCH(3) (R = n-C(3) H(7) , n-C(4) H(9) , n-C(5) H(11) , n-C(7) H(15) , n-C(9) H(19) , n-C(15) H(31) ) is examined. It is observed that the SBF/RER ratio is significantly influenced by the position of the substituent on the aromatic ring. As a rule, phenols and anisoles substituted by an alkyl group in meta position lead to more abundant methylene-2,4-cyclohexadiene cations (RER fragmentation) than their ortho and para homologues. This 'meta effect' is explained on the basis of energetic and kinetic of the two reaction channels. Quantum chemistry computations have been used to provide estimate of the thermochemistry associated with these two fragmentation routes. G3B3 calculation shows that a hydroxy or a methoxy group in the meta position destabilises the SBF and stabilises the RER product ions. Modelling of the SBF/RER intensities ratio has been performed assuming two single reaction rates for both fragmentation processes and computing them within the statistical RRKM formalism in the case of ortho, meta and para butyl phenols. It is clearly demonstrated that, combining thermochemistry and kinetics, the inequality (SBF/RER)(meta) < (SBF/RER)(ortho) < (SBF/RER)(para) holds for the butyl phenols series. It is expected that the 'meta effect' described in this study enables unequivocal identification of meta isomers from ortho and para isomers not only of alkyl phenols and alkyl anisoles but also in other alkyl benzene series.

Gas phase basicities of polyfunctional molecules. Part 3: Amino acids.

The present article is the third part of a general overview of the gas-phase protonation thermochemistry of polyfunctional molecules (first part: Mass Spectrom. Rev., 2007, 26:775-835, second part: Mass Spectrom. Rev., 2011, in press). This review is devoted to the 20 proteinogenic amino acids and is divided in two parts. In the first one, the experimental data obtained during the last 30 years using the equilibrium, thermokinetic and kinetic methods are presented. A general re-assignment of the values originating from these various experiments has been done on the basis of the commonly accepted Hunter & Lias 1998 gas-phase basicity scale in order to provide an homogeneous set of data. In the second part, theoretical investigations on gaseous neutral and protonated amino acids are reviewed. Conformational landscapes of both types of species were examined in order to provide theoretical protonation thermochemistry based on the truly identified most stable conformers. Proton affinities computed at the presently highest levels of theory (i.e. composite methods such as Gn procedures) are presented. Estimates of thermochemical parameters calculated using a Boltzmann distribution of conformers at 298K are also included. Finally, comparison between experiment and theory is discussed and a set of evaluated proton affinities, gas-phase basicities and protonation entropies is proposed.

Use of the HS-PTR-MS for online measurements of pyrethroids during indoor insecticide treatments.

A high-sensitivity proton transfer reaction mass spectrometer (HS-PTR-MS) has been used to study the temporal evolution of pesticide concentrations in indoor environments. Because of the high time variability of the indoor air concentrations during household pesticide applications, the use of this online high time resolution instrument is found relevant. Four pyrethroid pesticides of the latest generation that are commonly found in electric vaporizer refills, namely, transfluthrin, empenthrin, tetramethrin, and prallethrin, were considered. A controlled pesticide generation system was settled and coupled to a HS-PTR-MS analyzer, and a calibration procedure based on the fragmentation patterns of the protonated molecules was performed. To illustrate the functionality of the method, measurements of the concentration-time profiles of transfluthrin contained in an electric vaporizer were carried out in a full-scale environmental room under air exchange rate-controlled conditions. This study demonstrates that the HS-PTR-MS technique can provide online and high time-resolved measurements of semi-volatile organic compounds such as pyrethroid insecticides.

Isomerization and dissociation of n-butylbenzene radical cation.

Fragmentation mechanisms of ionized butylbenzene to give m/z 91 and m/z 92 fragment ions have been examined at the G3B3 and G3MP2B3 levels of theory. It is shown that the energetically favored pathways lead to tropylium, Tr(+), and methylene-2,4-cyclohexadiene, MCD(•+), ions. Formation of m/z 91 benzyl ions, Bz(+), by a simple bond fission (SBF) process, needs about 30 kJ/mol more energy than Tr(+). Possible formation of C(7)H(8)(•+) ions of structures different from the retro-ene rearrangement (RER) product, MCD(•+), has been also considered. Comparison with experimental data of this "thermometer" system is done through a kinetic modeling using Rice-Ramsperger-Kassel-Marcus (RRKM) and orbiting transition state (OTS) rate constant calculations on the G3MP2B3 0 K energy surface. The results agree with previous experimental observation if (i) the competitive formation of Tr(+) and Bz(+) is taken into account in the m/z 91 pathway, and (ii) the stepwise character of the RER fragmentation is introduced in the m/z 92 fragmentation route.

Gas-phase basicities of polyfunctional molecules. Part 2: Saturated basic sites.

The present article is the second part of a general overview of the gas-phase protonation thermochemistry of polyfunctional molecules. The first part of the review (Mass Spectrom. Rev., 2007, 26:775-835) was devoted to the description of the physico-chemical concepts and of the methods of determination, both experimental and theoretical, of gas-phase basicity. Several clues concerning the structural and energetic aspects of the protonation of isolated species have been emphasized. In the present article, specific examples are examined. The field of investigation is limited to molecules containing a "saturated" basic site, that is, nitrogen or oxygen atoms engaged in simple σ bonds with their neighboring. Aliphatic, cyclic and aromatic poly-amines, amino alcohols, alcohols, ethers, and hydroxyl-ethers, are successively presented.

Fragmentation reactions of phenoxide anions: deprotonated Dinoterb and related structures.

Dinoterb (6-t-butyl-2,4-dinitrophenol), 1, Dinoseb (6-secbutyl-2,4-dinitrophenol), 2, TBP (2-t-butylphenol), 3, and DNP (2,4-dinitrophenol), 4, have been analyzed by electrospray ionization in the negative mode (ESI-N) - tandem mass spectrometry. Nominal laboratory collision energy was varied from zero to 60 eV during the experiments. Apparent fragmentation energies were estimated from a parametric fitting of the collision efficiency curves. In parallel, fragmentation mechanisms of the deprotonated molecules [M-H](-) were explored using quantum chemistry modeling at the B3LYP/6-31 + G(d,p) level. A major fragmentation of the [M-H](-) ions of Dinoterb and Dinoseb is elimination of an alcohol molecule. This reaction is shown to involve one oxygen atom originating from a nitro group rather than the phenoxide moiety. Eliminations of NO, C(4) and CH(2) = C(CH(3))(2), i.e. reactions involving significant rearrangements, constitute the major part of the other fragmentation pathways observed from [3-H](-) and [4-H](-) ions.

Acid-base thermochemistry of gaseous oxygen and sulfur substituted amino acids (Ser, Thr, Cys, Met).

Acid-base thermochemistry of isolated amino acids containing oxygen or sulfur in their side chain (serine, threonine, cysteine and methionine) have been examined by quantum chemical computations. Density functional theory (DFT) was used, with B3LYP, B97-D and M06-2X functionals using the 6-31+G(d,p) basis set for geometry optimizations and the larger 6-311++G(3df,2p) basis set for energy computations. Composite methods CBS-QB3, G3B3, G4MP2 and G4 were applied to large sets of neutral, protonated and deprotonated conformers. Conformational analysis of these species, based on chemical approach and AMOEBA force field calculations, has been used to identify the lowest energy conformers and to estimate the population of conformers expected to be present at thermal equilibrium at 298 K. It is observed that G4, G4MP2, G3B3, CBS-QB3 composite methods and M06-2X DFT lead to similar conformer energies. Thermochemical parameters have been computed using either the most stable conformers or equilibrium populations of conformers. Comparison of experimental and theoretical proton affinities and Δ(acid)H shows that the G4 method provides the better agreement with deviations of less than 1.5 kJ mol(-1). From this point of view, a set of evaluated thermochemical quantities for serine, threonine, cysteine and methionine may be proposed: PA = 912, 919, 903, 938; GB = 878, 886, 870, 899; Δ(acid)H = 1393, 1391, 1396, 1411; Δ(acid)G = 1363, 1362, 1367, 1382 kJ mol(-1). This study also confirms that a non-negligible ΔpS° is associated with protonation of methionine and that the most acidic hydrogen of cysteine in the gas phase is that of the SH group. In several instances new conformers were identified thus suggesting a re-examination of several IRMPD spectra.

The shape of gaseous n-butylbenzene: assessment of computational methods and comparison with experiments.

Conformational landscape of neutral and ionized n-butylbenzene has been examined. Geometries have been optimized at the B3LYP/6-31G(d), B3LYP/6-31+G(d,p), B3LYP-D/6-31+G(d,p), B2PLYP/6-31+G(d,p), B2PLYP-D/6-31+G(d,p), B97-D/6-31+G(d,p), and M06-2X/6-31+G(d,p) levels. This study is complemented by energy computations using 6-311++G(3df,2p) basis set and CBS-QB3 and G3MP2B3 composite methods to obtain accurate relative enthalpies. Five distinguishable conformers have been identified for both the neutral and ionized systems. Comparison with experimentally determined rotational constants shows that the best geometrical parameters are provided by B3LYP-D and M06-2X functionals, which include an explicit treatment of dispersion effects. Composite G3MP2B3 and CBS-QB3 methods, and B2PLYP-D, B3LYP-D, B97-D, and M06-2X functionals, provide comparable relative energies for the two sets of neutral and ionized conformers of butyl benzene. An exception is noted however for conformer V(+) the stability of which being overestimated by the B3LYP-D and B97-D functionals. The better stability of neutral conformers I, III, and IV, and of cation I(+) , demonstrated by our computations, is in perfect agreement with conclusions based on micro wave, fluorescence, and multiphoton ionization experiments.

Acid-base thermochemistry of gaseous aliphatic α-aminoacids.

Acid-base thermochemistry of isolated aliphatic amino acids (denoted AAA): glycine, alanine, valine, leucine, isoleucine and proline has been examined theoretically by quantum chemical computations at the G3MP2B3 level. Conformational analysis on neutral, protonated and deprotonated species has been used to identify the lowest energy conformers and to estimate the population of conformers expected to be present at thermal equilibrium at 298 K. Comparison of the G3MP2B3 theoretical proton affinities, PA, and ΔH(acid) with experimental results is shown to be correct if experimental thermochemistry is re-evaluated and adapted to the most recent acidity-basicity scales. From this point of view, a set of evaluated proton affinities of 887, 902, 915, 916, 919 and 941 kJ mol(-1), and a set of evaluated ΔH(acid) of 1433, 1430, 1423, 1423, 1422 and 1426 kJ mol(-1), is proposed for glycine, alanine, valine, leucine, isoleucine and proline, respectively. Correlations with structural parameters (Taft's σ(α) polarizability parameter and molecular size) suggest that polarizability of the side chain is the major origin of the increase in PA and decrease in ΔH(acid) along the homologous series glycine, alanine, valine and leucine/isoleucine. Heats of formation of gaseous species AAA, AAAH(+) and [AAA-H](-) were computed at the G3MP2B3 level. The present study provides previously unavailable Δ(f)H°(298) for the ionized species AAAH(+) and [AAA-H](-). Comparison with Benson's estimate, and correlation with molecular size, show that several experimental Δ(f)H°(298) values of neutral or gaseous AAA might be erroneous.

Aromatic substitution reactions between ionized benzene derivatives and neutral methyl isocyanide.

Aromatic substitution reactions of selected ionized benzenes derivatives (chlorobenzene, nitrobenzene, dimethylphtalates) and phenoxy cation using neutral methyl isocyanide as a nucleophile are shown to efficiently occur in the gas phase. Nitrilium ions are produced in high abundance during these processes. These reactions have been performed in the hexapole collision cell of a large-scale hybrid tandem mass spectrometer. Computed 298 K enthalpy diagrams at the B3LYP/6-31+G(d,p) level of theory confirm the exothermic formation of the N-methylbenzonitrilium ions starting with ionized chloro- and nitrobenzene molecular ions. In this last case, two other exothermic processes are also detected: (i) an oxygen atom transfer yielding ionized nitrosobenzene and neutral methyl isocyanate and (ii) a loss of carbon monoxide from the ion/molecule reaction product generated when metastably generated phenoxy cations (produced in the hexapole collision cell) react with methyl isocyanide. Using extended theoretical calculations, several reaction pathways have been derived. The behavior of the three isomeric dimethyl phthalates has been investigated in the same way.

Reactions of ionized methyl benzoate with methyl isocyanide in the gas phase: nucleophilic aromatic substitutions vs hydrogen migrations.

The chemistry leading to the competitive eliminations of H, CH(3), and OCOCH(3) from adducts of ionized methyl benzoate and neutral methyl isocyanide has been explored using density functional theory molecular orbital calculations. The energies of the various reactants and transition structures were estimated at the B3LYP/6-31+G(d,p) level of theory. Nucleophilic aromatic substitution is proposed to account for the H and OCOCH(3) eliminations. The corresponding sigma-complex intermediates, B(1ipso) and B(1ortho), are stable species lying in deep energy wells situated 70 and 120 kJ/mol, respectively, below the reactants, ionized methyl benzoate and methyl isocyanide. The latter complex, B(1ortho), may be also at the origin of a multistep rearrangement involving hydrogen migrations and methyl elimination from the original methoxy group of the benzoate moiety.

Gas-phase protonation thermochemistry of glutamic acid.

Proton affinity, PA(Glu), and protonation entropy (i.e., the difference Delta(p)S(o)(Glu) = S(o)(GluH(+)) - S(o)(Glu)) of glutamic acid have been experimentally determined by the extended kinetic method using electrospray ionization triple quadrupole-time-of-flight (ESI-Q-TOF) tandem mass spectrometry. The values deduced from these experiments are PA(Glu) = 945.3 +/- 2.8(5.8) kJ x mol(-1) and Delta(p)S(o)(Glu) = -28 +/- 4(9) J x mol(-1) x K(-1) thus leading to a gas-phase basicity, GB(Glu), of 904.4 +/- 3.0(6.4) kJ x mol(-1) (uncertainties are standard deviation and, in parentheses, 95% confidence limit). Theoretical calculations performed at the G3MP2B3 level provide information on the structures, conformations, and energetics of the neutral and protonated species. Thermochemical data are calculated at this level and include a correction to the computation of the entropy associated with hindered rotation. When the lowest energy conformers of protonated and neutral glutamic acid are considered the following values are calculated: PA(Glu) = 948.1 kJ x mol(-1) and Delta(p)S(o)(Glu) = -31.3 J x mol(-1) x K(-1). Using G3MP2B3 data to estimate the gas-phase distribution of conformers at 298 K, the averaged molar quantities becomes PA(Glu) = 949.8 kJ x mol(-1) and Delta(p)S(o)(Glu) = -36.0 J x mol(-1) x K(-1). Both computations give comparable GB(Glu) = 906.4-906.7 kJ x mol(-1).

Gas phase protonation thermochemistry of phenylalanine and tyrosine.

Gas phase basicities of phenylalanine and tyrosine, GB(Phe) = 892.0 +/- 1.3(2.6) kJ.mol(-1) and GB(Tyr) = 894.9 +/- 2.8(5.9) kJ.mol(-1) (uncertainties are standard deviation and, in parentheses, 95% confidence limit), have been experimentally determined by the extended kinetic method using ESI-TQ tandem mass spectrometry. Proton affinities deduced from these experiments, PA(Phe) = 931.3 +/- 1.1(2.3) kJ.mol(-1) and PA(Tyr) = 934.8 +/- 2.5(5.2) kJ.mol(-1), are perfectly reproduced by theoretical calculations performed at the B3LYP/6-311++G(3df,2p)//B3LYP/6-31+G(d,p) level. An entropy loss of approximately -25 J.mol(-1).K(-1) occurs upon protonation of both Phe and Tyr. The origin of this entropy change is attributed (i) to the change in strength of the interaction between the amino group and the aromatic moiety in the neutral and protonated forms and (ii) to the larger entropy of mixing associated with the population of neutral conformers with respect to their protonated counterparts. Previous neglect of the protonation entropy term has led to underestimated tabulated PA values; the evaluated values proposed in the present study are PA(Phe) = 932 +/- kJ.mol(-1) and PA(Tyr) = 935 +/- kJ.mol(-1).