Binding Energies of Proton-Bound Dimers of Imidazole and n-Acetylalanine Methyl Ester Obtained by Blackbody Infrared Radiative Dissociation

The dissociation kinetics of protonated n-acetyl-L-alanine methyl ester dimer (AcAlaME(d)), imidazole dimer, and their cross dimer were measured using blackbody infrared radiative dissociation (BIRD). Master equation modeling of these data was used to extract threshold dissociation energies (E(o)) for the dimers. Values of 1.18 +/- 0.06, 1.11 +/- 0.04, and 1.12 +/- 0.08 eV were obtained for AcAlaME(d), imidazole dimer, and the cross dimer, respectively. Assuming that the reverse activation barrier for dissociation of the ion-molecule complex is negligible, the value of E(o) can be compared to the dissociation enthalpy (DeltaH(d) degrees ) from HPMS data. The E(o) values obtained for the imidazole dimer and the cross dimer are in agreement with HPMS values; the value for AcAlaME(d) is somewhat lower. Radiative rate constants used in the master equation modeling were determined using transition dipole moments calculated at the semiempirical (AM1) level for all dimers and compared to ab initio (RHF/3-21G*) calculations where possible. To reproduce the experimentally measured dissociation rates using master equation modeling, it was necessary to multiply semiempirical transition dipole moments by a factor between 2 and 3. Values for transition dipole moments from the ab initio calculations could be used for two of the dimers but appear to be too low for AcAlaME(d). These results demonstrate that BIRD, in combination with master equation modeling, can be used to determine threshold dissociation energies for intermediate size ions that are in neither the truncated Boltzmann nor the rapid energy exchange limit.     

Binding energies of protonated betaine complexes: a probe of zwitterion structure in the gas phase

The dissociation kinetics of proton-bound dimers of betaine with molecules of comparable gas-phase basicity were investigated using blackbody infrared radiative dissociation (BIRD). Threshold dissociation energies were obtained from these data using master equation modeling. For bases that have comparable or higher gas-phase basicity, the binding energy of the protonated base.betaine complex is ~1.4 eV. For molecules that are ~2 kcal/mol or more less basic, the dissociation energy of the complexes is ~1.2 eV. The higher binding energy of the former is attributed to an ion-zwitterion structure which has a much larger ion-dipole interaction. The lower binding energy for molecules that are ~2 kcal/mol or more less basic indicates that an ion-molecule structure is more favored. Semiempirical calculations at both the AM1 and PM3 levels indicate the most stable ion-molecule structure is one in which the base interacts with the charged quaternary ammonium end of betaine. These results indicate that the measurement of binding energies of neutral molecules to biological ions could provide a useful probe for the presence of zwitterions and salt bridges in the gas phase. From the BIRD data, the gas-phase basicity of betaine obtained from the kinetic method is found to be 239.2 +/- 1.0 kcal/mol. This value is in excellent agreement with the value of 239.3 kcal/mol (298 K) from ab initio calculations at the MP2/6-31+g** level. The measured value is slightly higher than those reported previously. This difference is attributed to entropy effects. The lower ion internal energy and longer time frame of BIRD experiments should provide values closer to those at standard temperature.     

Dissociation energetics and mechanisms of leucine enkephalin (M+H)+ and (2M+X)+ ions (X=H, Li, Na, K, and Rb) measured by blackbody infrared radiative dissociation

The dissociation kinetics of protonated leucine enkephalin and its proton and alkali metal bound dimers were investigated by blackbody infrared radiative dissociation in a Fourier-transform mass spectrometer. From the temperature dependence of the unimolecular dissociation rate constants, Arrhenius activation parameters in the zero-pressure limit are obtained. Protonated leucine enkephalin dissociates to form b(4) and (M-H(2)O)(+) ions with an average activation energy (E(a)) of 1.1 eV and an A factor of 10(10.5) s(-1). The value of the A factor indicates that these dissociation processes are rearrangements. The b(4) ions subsequently dissociate to form a(4) ions via a process with a relatively high activation energy (1.3 eV), but one that is entropically favored. For the cationized dimers, the thermal stability decreases with increasing cation size, consistent with a simple electrostatic interaction in these noncovalent ion-molecule complexes. The E(a) and A factors are indistinguishable within experimental error with values of ~1.5 eV and 10(17) s(-1), respectively. Although not conclusive, results from master equation modeling indicate that all these BIRD processes, except for b(4) --> a(4), are in the rapid energy exchange limit. In this limit, the internal energy of the precursor ion population is given by a Boltzmann distribution and information about the energetics and dynamics of the reaction are obtained directly from the measured Arrhenius parameters.     

Dissociation of heme-globin complexes by blackbody infrared radiative dissociation: Molecular specificity in the gas phase?

The temperature dependence of the unimolecular kinetics for dissociation of the heme group from holo-myoglobin (Mb) and holo-hemoglobin alpha-chain (Hb-alpha) was investigated with blackbody infrared radiative dissociation (BIRD). The rate constant for dissociation of the 9 + charge state of Mb formed by electrospray ionization from a "pseudo-native" solution is 60% lower than that of Hb-alpha at each of the temperatures investigated. In solutions of pH 5.5-8.0, the thermal dissociation rate for Mb is also lower than that of HB-alpha (Hargrove, M. S. et al. J. Biol. Chem.1994, 269, 4207-4214). Thus, Mb is thermally more stable with respect to heme loss than Hb-alpha both in the gas phase and in solution. The Arrhenius activation parameters for both dissociation processes are indistinguishable within the current experimental error (activation energy 0.9 eV and pre-exponential factor of 10(8-10) s(-1)). The 9+ to 12+ charge states of Mb have similar Arrhenius parameters when these ions are formed from pseudo-native solutions. In contrast, the activation energies and pre-exponential factors decrease from 0.8 to 0.3 eV and 10(7) to 10(2) s(-1), respectively, for the 9 + to 12 + charge states formed from acidified solutions in which at least 50% of the secondary structure is lost. These results demonstrate that gas-phase Mb ions retain clear memory of the composition of the solution from which they are formed and that these differences can be probed by BIRD.     

Activation of Peptide ions by blackbody radiation: factors that lead to dissociation kinetics in the rapid energy exchange limit

Unimolecular rate constants for blackbody infrared radiative dissociation (BIRD) were calculated for the model protonated peptide (AlaGly)(n) (n = 2-32) using a variety of dissociation parameters. Combinations of dissociation threshold energies ranging from 0.8 to 1.7 eV and transition entropies corresponding to Arrhenius preexponential factors ranging from very "tight" (A(infinity) = 10(9.9) s(-1)) to "loose" (A(infinity) = 10(16.8) s(-1)) were selected to represent dissociation parameters within the experimental temperature range (300-520 K) and kinetic window (k(uni) = 0.001-0.20 s(-1)) typically used in the BIRD experiment. Arrhenius parameters were determined from the temperature dependence of these values and compared to those in the rapid energy exchange (REX) limit. In this limit, the internal energy of a population of ions is given by a Boltzmann distribution, and kinetics are the same as those in the traditional high-pressure limit. For a dissociation process to be in this limit, the rate of photon exchange between an ion and the vacuum chamber walls must be significantly greater than the dissociation rate. Kinetics rapidly approach the REX limit either as the molecular size or threshold dissociation energy increases or as the transition-state entropy or experimental temperature decreases. Under typical experimental conditions, peptide ions larger than 1.6 kDa should be in the REX limit. Smaller ions may also be in the REX limit depending on the value of the threshold dissociation energy and transition-state entropy. Either modeling or information about the dissociation mechanism must be known in order to confirm REX limit kinetics for these smaller ions. Three principal factors that lead to the size dependence of REX limit kinetics are identified. With increasing molecular size, rates of radiative absorption and emission increase, internal energy distributions become relatively narrower, and the microcanonical dissociation rate constants increase more slowly over the energy distribution of ions. Guidelines established here should make BIRD an even more reliable method to obtain information about dissociation energetics and mechanisms for intermediate size molecules.     

Direct determination of the ionization energies of PtC, PtO, and PtO2 with VUV radiation

Photoionization efficiency curves were measured for gas-phase PtC, PtO, and PtO2 using tunable vacuum ultraviolet (VUV) radiation at the Advanced Light Source. The molecules were prepared by laser ablation of a platinum tube, followed by reaction with CH4 or N2O and supersonic expansion. These measurements provide the first directly measured ionization energy for PtC, IE(PtC) = 9.45 +/- 0.05 eV. The direct measurement also gives greatly improved ionization energies for the platinum oxides, IE(PtO) = 10.0 +/- 0.1 eV and IE(PtO2) = 11.35 +/- 0.05 eV. The ionization energy connects the dissociation energies of the neutral and cation, leading to greatly improved 0 K bond dissociation energies for the neutrals: D0(Pt-C) = 5.95 +/- 0.07 eV, D0(Pt-O) = 4.30 +/- 0.12 eV, and D0(OPt-O) = 4.41 +/- 0.13 eV, as well as enthalpies of formation for the gas-phase molecules DeltaH(0)(f,0)(PtC(g)) = 701 +/- 7 kJ/mol, DeltaH(0)(f,0)(PtO(g)) = 396 +/- 12 kJ/mol, and DeltaH(0)(f,0)(PtO2(g)) = 218 +/- 11 kJ/mol. Much of the error in previous Knudsen cell measurements of platinum oxide bond dissociation energies is due to the use of thermodynamic second law extrapolations. Third law values calculated using statistical mechanical thermodynamic functions are in much better agreement with values obtained from ionization energies and ion energetics. These experiments demonstrate that laser ablation production with direct VUV ionization measurements is a versatile tool to measure ionization energies and bond dissociation energies for catalytically interesting species such as metal oxides and carbides.     

Protonation processes and electronic spectra of histidine and related ions

A full structural assignment of the neutral, protonated, and deprotonated histidine conformers in the gas phase is presented. A total of 3024 unique trial structures were generated by all combinations of internal single-bond rotamers of these species and optimized at the B3LYP/6-311G* level and further optimized at the B3LYP/6-311++G** level. A set of unique conformers is found, and their relative energies, free energies, dipole moments, rotational constants, electron affinities, ionization energies, and harmonic frequencies are determined. The population ratio of histidine and its tautomer is 1:0.16 at 298 K. Massive conformational changes are observed due to protonation and deprotonation, and the intramolecular H-bonds are characterized with the atoms in molecules theory. The calculated proton dissociation energy, gas-phase acidity, proton affinity, and gas-phase basicity are in excellent agreement with the experiments. The deprotonation and protonation of gaseous histidine both occur on the imidazole ring, explaining the versatile biofunctions of histidine in large biomolecules. The UV spectra of neutral and singly and doubly protonated histidine are investigated with the TDDFT/B3LYP/6-311+G(2df,p) calculations. The S0-S1, S0-S2, and S0-S3 excitations of histidine are mixed pipi*/npi* transitions at 5.37, 5.44, and 5.69 eV, respectively. The three excitation energies for histidine tautomer are 4.85, 5.47, and 5.52 eV, respectively. The three excitations for protonated histidine are mainly npi* transitions at 5.45, 5.67, and 5.82 eV, respectively. The S0-S1 excitation of protonated histidine produces ImH-CbetaH2-CalphaH(COOH)-NH2+, while the S0-S2 and S0-S3 transitions produce ImH-CbetaH2-CalphaH(NH2)-(COOH)+. These data may help to understand the mechanisms of the UV fragmentation of biomolecules.     

Blackbody infrared radiative dissociation of bradykinin and its analogues: energetics, dynamics, and evidence for salt-bridge structures in the gas phase

Blackbody infrared radiative dissociation (BIRD) spectra of singly and doubly protonated bradykinin and its analogues are measured in a Fourier-transform mass spectrometer. Rate constants for dissociation are measured as a function of temperature with reaction delays up to 600 s. From these data, Arrhenius activation parameters in the zero-pressure limit are obtained. The activation parameters and dissociation products for the singly protonated ions are highly sensitive to small changes in ion structure. The Arrhenius activation energy (E(a)) and pre-exponential (or frequency factor, A) of the singly protonated ions investigated here range from 0.6 to 1.4 eV and 10(5) to 10(12) s(-1), respectively. For bradykinin and its analogues differing by modification of the residues between the two arginine groups on either end of the molecule, the singly and doubly protonated ions have average activation energies of 1.2 and 0.8 eV, respectively, and average A values of 10(8) and 10(12) s(-1), respectively, i.e., the presence of a second charge reduces the activation energy by 0.4 eV and decreases the A value by a factor of 10(4). This demonstrates that the presence of a second charge can dramatically influence the dissociation dynamics of these ions. The doubly protonated methyl ester of bradykinin has an E(a) of 0.82 eV, comparable to the value of 0.84 eV for bradykinin itself. However, this value is 0.21 +/- 0.08 eV greater than that of singly protonated methyl ester of bradykinin, indicating that the Coulomb repulsion is not the most significant factor in the activation energy of this ion. Both singly and doubly protonated Lys-bradykinin ions have higher activation energies than the corresponding bradykinin ions indicating that the addition of a basic residue stabilizes these ions with respect to dissociation. Methylation of the carboxylic acid group of the C-terminus reduces the E(a) of bradykinin from 1.3 to 0.6 eV and the A factor from 1012 to 105 s(-1). This modification also dramatically changes the dissociation products. Similar results are observed for [Ala(6)]-bradykinin and its methyl ester. These results, in combination with others presented here, provide experimental evidence that the most stable form of singly protonated bradykinin is a salt-bridge structure.     

Packing Structures of Thiophene Dimers and Their Effects on Excitation Energies of Thiophene Dimers

IntroductionPi(π) conjugated polymers have received exten-sive experimental and theoretical attention because oftheir fascinating electronic and optical properties[1—3].Recently, many new experiment results[4—6]have re-vealed that the inter-chain inter…     

Dissociative electron attachment to gas phase valine: a combined experimental and theoretical study

Using a crossed electron/molecule beam technique the dissociative electron attachment (DEA) to gas phase L-valine, (CH(3))(2)CHCH(NH(2))COOH, is studied by means of mass spectrometric detection of the product anions. Additionally, ab initio calculations of the structures and energies of the anions and neutral fragments have been carried out at G2MP2 and B3LYP levels. Valine and the previously studied aliphatic amino acids glycine and alanine exhibit several common features due to the fact that at low electron energies the formation of the precursor ion can be characterized by occupation of the pi* orbital of the carboxyl group. The dominant negative ion (M-H)(-) (m/Z=116) is observed at electron energies of 1.12 eV. This ion is the dominant reaction product at electron energies below 5 eV. Additional fragment ions with m/Z=100, 72, 56, 45, 26, and 17 are observed both through the low lying pi* and through higher lying resonances at about 5.5 and 8.0-9.0 eV, which are characterized as core excited resonances. According to the threshold energies calculated here, rearrangements play a significant role in the formation of DEA fragments observed from valine at subexcitation energies.     

Revisiting the photophysical properties and excited singlet-state dipole moments of several coumarin derivatives

The solvent effects on the electronic absorption and fluorescence emission spectra of several coumarins derivatives, containing amino, N,N-dimethyl-amino, N,N-diethyl-amino, hydroxyl, methyl, carboxyl, or halogen substituents at the positions 7, 4, or 3, were investigated in eight solvents with various polarities. The first excited singlet-state dipole moments of these coumarins were determined by various solvatochromic methods, using the theoretical ground-state dipole moments which were calculated by the AM1 method. The first excited singlet-state dipole moment values were obtained by the Bakhshiev, Kawski-Chamma-Viallet, Lippert-Mataga, and Reichardt-Dimroth equations, and were compared to the ground-state dipole moments. In all cases, the dipole moments were found to be higher in the excited singlet-state than in the ground state because of the different electron densities in both states. The red-shifts of the absorption and fluorescence emission bands, observed for most compounds upon increasing the solvent polarity, indicated that the electronic transitions were of π-π* nature.     

Weak metal-metal bonding in small manganese cluster ions, Mn(N)+ (N < or = 7)

The binding energies of manganese cluster ions Mn(N)+ (N = 5-7) were determined by the photodissociation experiments in the near-infrared and visible-photon-energy ranges. The bond dissociation energies of Mn(N)+, D0(Mn(N-1)+...Mn), were obtained to be 1.70+/-0.08, 1.04+/-0.10, and 1.46+/-0.11 eV, respectively, for N = 5, 6, and 7 from the threshold energies for the two-atom loss processes and the bond dissociation energies of Mn3(+) and Mn4(+) reported previously [A. Terasaki et al., J. Chem. Phys. 117, 7520 (2002)]. Correspondingly, binding energies per atom are obtained to be 0.99+/-0.03, 1.00+/-0.03, and 1.06+/-0.03 eV/at. for N = 5, 6, and 7, respectively. A gradual increase in the binding energy from N = 2 to N = 7 shows an increasing contribution of nonbonding 3d orbitals to the bonding via weak hybridization with valence 4s orbitals as the cluster size increases. These binding energies per atom are still much smaller than the bulk cohesive energy of manganese (2.92 eV/at.), and this finding indicates exceptionally weak metal-metal bonds in this size range.     

Gas-phase oxidation of Cm+ and Cm2+ --thermodynamics of neutral and ionized CmO

Fourier transform ion cyclotron resonance mass spectrometry was employed to study the products and kinetics of gas-phase reactions of Cm (+) and Cm (2+); parallel studies were carried out with La (+/2+), Gd (+/2+) and Lu (+/2+). Reactions with oxygen-donor molecules provided estimates for the bond dissociation energies, D[M (+)-O] (M = Cm, Gd, Lu). The first ionization energy, IE[CmO], was obtained from the reactivity of CmO (+) with dienes, and the second ionization energies, IE[MO (+)] (M = Cm, La, Gd, Lu), from the rates of electron-transfer reactions from neutrals to the MO (2+) ions. The following thermodynamic quantities for curium oxide molecules were obtained: IE[CmO] = 6.4 +/- 0.2 eV; IE[CmO (+)] = 15.8 +/- 0.4 eV; D[Cm-O] = 710 +/- 45 kJ mol (-1); D[Cm (+)-O] = 670 +/- 40 kJ mol (-1); and D[Cm (2+)-O] = 342 +/- 55 kJ mol (-1). Estimates for the M (2+)-O bond energies for M = Cm, La, Gd, and Lu are all intermediate between D[N 2-O] and D[OC-O] - that is, 167 kJ mol (-1) < D[M (2+)-O] < 532 kJ mol (-1) - such that the four MO (2+) ions fulfill the thermodynamic requirement for catalytic oxygen-atom transport from N2O to CO. It was demonstrated that the kinetics are also favorable and that the CmO (2+), LaO (2+), GdO (2+), and LuO (2+) dipositive ions each catalyze the gas-phase oxidation of CO to CO2 by N2O. The CmO 2 (+) ion appeared during the reaction of Cm (+) with O 2 when the intermediate, CmO (+), was not collisionally cooled - although its formation is kinetically and/or thermodynamically unfavorable, CmO 2 (+) is a stable species.     

Computational comparison of the stacking interactions between the aromatic amino acids and the natural or (cationic) methylated nucleobases

The strongest gas-phase MP2/6-31G*(0.25) stacking energies between the aromatic amino acids and the natural or methylated nucleobases were considered. The potential energy surfaces of dimers were searched as a function of the vertical separation, angle of rotation and horizontal displacement between monomers stacked according to their centers of mass. Our calculations reveal that the stacking interactions of adducts for a given nucleobase are dependent on the methylation site (by up to 20 kJ mol(-1)), where the relative magnitudes of the interactions are determined by the dipole moments of the adducts and the proton affinities of nucleobase methylation sites. Nevertheless, the differences in the (gas-phase) stacking of methylated adducts are small compared with the differences between the stacking of the corresponding natural and methylated nucleobases. Indeed, methylation increases the stacking energy by up to 40 kJ mol(-1) (or 135%). Although immersing the dimers in different solvents decreases the gas-phase stacking energies with an increase in the polarity of the environment, base methylation still has a significant effect on the nucleobase stacking ability in solvents with large dipole moments, and, perhaps more importantly, environments that mimic enzyme active sites. Our results shed light on the workings of DNA repairs enzymes that selectively remove a wide variety of alkylated nucleobases over the natural bases.     

Reactivity and infrared spectroscopy of gaseous hydrated trivalent metal ions

Hydrated trivalent rare earth metal ions containing yttrium and all naturally abundant lanthanide metals are formed using electrospray ionization, and the structures and reactivities of these ions containing 17-21 water molecules are probed using blackbody infrared radiative dissociation (BIRD) and infrared action spectroscopy. With the low-energy activation conditions of BIRD, there is an abrupt transition in the dissociation pathway from the exclusive loss of a single neutral water molecule to the exclusive loss of a small protonated water cluster via a charge-separation process. This transition occurs over a narrow range of cluster sizes that differs by only a few water molecules for each metal ion. The effective turnover size at which these two dissociation rates become equal depends on metal ion identity and is poorly correlated with the third ionization energies of the isolated metals but is well correlated with the hydrolysis constants of the trivalent metal ions in bulk aqueous solution. Infrared action spectra of these ions at cluster sizes near the turnover size are largely independent of the specific identity of the trivalent metal ion, suggesting that any differences in the structures of the ions present in our experiment are subtle.     

Theoretical study of the electronic gas-phase spectrum of glycine, alanine, and related amines and carboxylic acids

A theoretical study on the origin of the common electronic excitations in amino acids is presented, focusing on the excited states of glycine, alanine and the related substructures formic acid, acetic acid, propionic acid, ammonia, methylamine, and ethylamine. Special attention is given to the valence excitation from the nonbonding lone-pair on the carboxylic oxygen atom to the antibonding pi-orbital (n(O) --> pi*(CO)) and the first Rydberg excitation from the nonbonding lone-pair on the nitrogen atom (n(N) --> 3s). From extensive calculations on formic acid and methylamine, different basis sets and electron correlation treatments are benchmarked using a hierarchy of coupled cluster (CC) methods, consisting of CCS, CC2, CCSD, CCSDR(3), and CC3, in combination with augmented correlation consistent basis sets. The dependence of the excitation energies on the size of the backbone structure in the two groups of molecules is investigated, and 0-0 transition energies for the n(O) --> pi*(CO) and n(N) --> 3s transitions are calculated for the smallest molecules. Excellent agreement with experimental values is found where secure experimental assignments are available. A few outstanding problems in the experimental assignments found in the literature are described for both the carboxylic acids and the amines. Final predictions for vertical excitation energies are given for all molecules, including glycine and alanine where no gas-phase experimental results are available. Finally, calculations on protonated amino acids are presented showing an isolation of the n(O) --> pi*(CO) from higher lying states by as much as 1.9 eV for alanine.     

Total dissociative electron attachment cross sections of selected amino acids

Total dissociative electron attachment cross sections are presented for the amino acids, glycine, alanine, proline, phenylalanine, and tryptophan, at energies below the first ionization energy. Cross section magnitudes were determined by observation of positive ion production and normalization to ionization cross sections calculated using the binary-encounter-Bethe method. The prominent 1.2 eV feature in the cross sections of the amino acids and the closely related HCOOH molecule is widely attributed to the attachment into the -COOH pi* orbital. The authors discuss evidence that direct attachment to the lowest sigma* orbital may instead be responsible. A close correlation between the energies of the core-excited anion states of glycine, alanine, and proline and the ionization energies of the neutral molecules is found. A prominent feature in the total dissociative electron attachment cross section of these compounds is absent in previous studies using mass analysis, suggesting that the missing fragment is energetic H-.     

Binding energies of water to lithiated valine: Formation of solution-phase structure in vacuo

Dissociation kinetics for loss of a water molecule from hydrated ions of lithiated valine, alanine ethyl ester and betaine are determined using blackbody infrared radiative dissociation at temperatures between -60 and 110 degrees C. From master equation modeling of these data, values of the threshold dissociation energy are obtained for clusters containing one through three water molecules. By comparing the values for valine with its two isomers, one a model for the nonzwitterion structure, the other a model for the zwitterion structure, information about the structure of valine in these hydrated clusters is inferred. Structures, relative energies, and water binding energies for these ions are also calculated at the B3LYP/6-31++G** level of theory. With one water molecule, both experiment and theory indicate that valine is not a zwitterion and that the lithium ion coordinates with the amino nitrogen and the carbonyl oxygen (NO coordinated) and the water molecule interacts directly with the lithium ion. With two water molecules, the zwitterion and nonzwitterion structures are nearly isoenergetic, but the experiment clearly indicates a NO-coordinated nonzwitterion structure. With three water molecules, both the experimental data and theory indicate that the lithium ion binds to the carboxylate group of valine, i.e., valine is zwitterionic with three water molecules. The agreement between the experimentally determined and calculated binding energies is good for all the clusters, with deviations of <== 0.12 eV.     

Structures of cationized proline analogues: evidence for the zwitterionic form

The structures of lithiated and sodiated alpha-methyl-proline (alpha-Me-Pro) and structural isomers, both with and without a water molecule, are investigated using blackbody infrared radiative dissociation (BIRD) and density functional theory. From the BIRD kinetic data measured as a function of temperature, combined with master equation modeling of these data, threshold dissociation energies for the loss of a water molecule from these clusters are obtained. These energies are 77.5 +/- 0.5 and 53 +/- 1 kJ/mol for lithiated and sodiated alpha-Me-Pro, respectively. For the nonzwitterionic isomer, proline methyl ester, these values are 3.0-4.5 kJ/mol higher. These results provide compelling experimental evidence that alpha-Me-Pro is zwitterionic in these clusters. Theory at the temperature corrected B3LYP/6-311++G**//B3LYP/6-31++G** level indicates that the salt-bridge or zwitterionic forms of lithiated and sodiated alpha-Me-Pro are between 17 and 23 kJ/mol lower in energy than the nonzwitterionic or charge-solvated forms and that attachment of a single water molecule does not significantly change the structure or the relative energies of these clusters. The proton affinity of proline is 8 kJ/mol higher than that of alpha-Me-Pro, indicating that lithiated and sodiated singly hydrated proline should also be zwitterionic.     

DFT studies using supercells and projector-augmented waves for structure, energetics, and dynamics of glycine, alanine, and cysteine

A large variety of gas phase conformations of the amino acids glycine, alanine, and cysteine is studied by numerically efficient semi-local gradient-corrected density functional theory calculations using a projector-augmented wave scheme and periodic boundary conditions. Equilibrium geometries, conformational energies, dipole moments, vibrational modes, and IR optical spectra are calculated from first principles. A comparison of our results with values obtained from quantum-chemistry methods with localized basis sets and nonlocal exchange-correlation functionals as well as with experimental data is made. For conformations containing strong intramolecular hydrogen bonds deviations in their energetic ordering occur, which are traced back to different treatments of spatial nonlocality in the exchange-correlation functional. However, even for these structures, the comparison of calculated and measured vibrational frequencies shows satisfying agreement.