Modeling of protonation processes in acetohydroxamic acid

This work presents a theoretical study of acetohydroxamic acid and its protonation processes using ab initio methodology at the MP2(FC)/cc-pdVZ level. We find the amide form more stable than the imidic tautomer by less than 1.0 kcal mol(-)(1). For comparison with the experimental data, a three-dimensional conformational study is performed on the most stable tautomer (amide). From this study, the different barriers to rotation and inversion are determined and the intramolecular hydrogen bond between the OH group and the carbonyl oxygen is characterized. The electrostatic potential distribution shows three possible sites for electrophilic attack, but it is shown that only two of them, the carbonyl oxygen and the nitrogen atoms, are actual protonation sites. The protonation energy (proton affinity) is obtained from the results of the neutral and charged species. Proton affinities for the species charged on the carbonyl oxygen and the nitrogen atoms are estimated to be 203.4 and 194.5 kcal mol(-)(1), respectively. The development of a statistical model permits the quantification of DeltaG (gas-phase basicity) for the two protonation processes. In this way, the carbonyl oxygen protonated form is found to be more stable than that of the nitrogen atoms by 8.3 kcal mol(-)(1) at 1 atm and 298.15 K, due to the enthalpic contribution. As temperature increases, the proportion of the nitrogen protonated form increases slightly.     

An Investigation of Protonation Sites and Conformations of Protonated Amino Acids by IRMPD Spectroscopy

The protonation sites and structures of a series of protonated amino acids (Gly, Ala, Pro, Phe, Lys and Ser) are investigated by means of infrared multiple‐photon dissociation (IRMPD) spectroscopy and electronic‐structure calculations. The IRMPD spectra of the protonated species are recorded using the combination of a free‐electron laser (FEL) and an electrospray‐ion‐trap mass spectrometer. The structures of different possible isomers of these protonated species are optimized at the B3LYP/6‐311+G(d, p) level of theory and the IR spectra calculated using the same computational method. For every amino acid studied herein, the current results indicate that a proton is bound to the α‐amino nitrogen, except for lysine, in which the protonation site is the amino nitrogen in the side chain. According to the calculated and experimental IRMPD results, the structures of the protonated amino acids may be assigned unambiguously. For Gly, Ala, and Pro, in each of the most stable isomers the protonated amino group forms an intramolecular hydrogen bond with the adjacent carbonyl oxygen. In the case of Gly, the isomer containing a proton bound to the carbonyl oxygen is theoretically possible. However, it does not exist under the experimental conditions because it has a significantly higher energy (i.e. 26.6 kcal mol^?1) relative to the most stable isomer. For Ser and Phe, the protonated amino group forms two intramolecular hydrogen bonds with both the adjacent carbonyl oxygen and the side‐chain group in each of the most stable isomers. In protonated lysine, the protonated amino group in the side chain forms two hydrogen bonds with the α‐amino nitrogen and the carbonyl oxygen, which is a cyclic structure. Interestingly, for protonated lysine the zwitterionic structure is a local minimum energy isomer, but the experimental spectrum indicates that it does not exist under the experimental conditions. This is consistent with the fact that the zwitterionic isomer is 9.2 kcal mol^?1 higher in free energy at 298 K than the most stable isomer. The carbonyl stretching vibration in the range of 1760–1800 cm^?1 is especially sensitive to the structural change. In addition, IRMPD mechanisms for the protonated amino acids are also investigated.     

Corannulene as a Lewis base: computational modeling of protonation and lithium cation binding.

A computational modeling of the protonation of corannulene at B3LYP/6-311G(d,p)//B3LYP/6-311G(d,p) and of the binding of lithium cations to corannulene at B3LYP/6-311G(d,p)//B3LYP/6-31G(d,p) has been performed. A proton attaches preferentially to one carbon atom, forming a sigma-complex. The isomer protonated at the innermost (hub) carbon has the best total energy. Protonation at the outermost (rim) carbon and at the intermediate (bridgehead rim) carbon is less favorable by ca. 2 and 14 kcal mol(-)(1), respectively. Hydrogen-bridged isomers are transition states between the sigma-complexes; the corresponding activation energies vary from 10 to 26 kcal mol(-)(1). With an empirical correction obtained from calculations on benzene, naphthalene, and azulene, the best estimate for the proton affinity of corannulene is 203 kcal mol(-)(1). The lithium cation positions itself preferentially over a ring. There is a small energetic preference for the 6-ring over the 5-ring binding (up to 2 kcal mol(-)(1)) and of the convex face over the concave face (3-5 kcal mol(-)(1)). The Li-bridged complexes are transition states between the pi-face complexes. Movement of the Li(+) cation over either face is facile, and the activation energy does not exceed 6 kcal mol(-)(1) on the convex face and 2.2 kcal mol(-)(1) on the concave face. In contrast, the transition of Li(+) around the corannulene edge involves a high activation barrier (24 kcal mol(-)(1) with respect to the lowest energy pi-face complex). An easier concave/convex transformation and vice versa is the bowl-to-bowl inversion with an activation energy of 7-12 kcal mol(-)(1). The computed binding energy of Li(+) to corannulene is 44 kcal mol(-)(1). Calculations of the (7)Li NMR chemical shifts and nuclear independent chemical shifts (NICS) have been performed to analyze the aromaticity of the corannulene rings and its changes upon protonation.     

Cyclization and rearrangement reactions of a(n) fragment ions of protonated peptides

a(n) ions are frequently formed in collision-induced dissociation (CID) of protonated peptides in tandem mass spectrometry (MS/MS) based sequencing experiments. These ions have generally been assumed to exist as immonium derivatives (-HN(+)═CHR). Using a quadrupole ion trap mass spectrometer, MS/MS experiments have been performed and the structure of a(n) ions formed from oligoglycines was probed by infrared spectroscopy. The structure and isomerization reactions of the same ions were studied using density functional theory. Overall, theory and infrared spectroscopy provide compelling evidence that a(n) ions undergo cyclization and/or rearrangement reactions, and the resulting structure(s) observed under our experimental conditions depends on the size (n). The a(2) ion (GG sequence) undergoes cyclization to form a 5-membered ring isomer. The a(3) ion (GGG sequence) undergoes cyclization initiated by nucleophilic attack of the carbonyl oxygen of the N-terminal glycine residue on the carbon center of the C-terminal immonium group forming a 7-membered ring isomer. The barrier to this reaction is comparatively low at 10.5 kcal mol(-1), and the resulting cyclic isomer (-5.4 kcal mol(-1)) is more energetically favorable than the linear form. The a(4) ion with the GGGG sequence undergoes head-to-tail cyclization via nucleophilic attack of the N-terminal amino group on the carbon center of the C-terminal immonium ion, forming an 11-membered macroring which contains a secondary amine and three trans amide bonds. Then an intermolecular proton transfer isomerizes the initially formed secondary amine moiety (-CH(2)-NH(2)(+)-CH(2)-NH-CO-) to form a new -CH(2)-NH-CH(2)-NH(2)(+)-CO- form. This structure is readily cleaved at the -CH(2)-NH(2)(+)- bond, leading to opening of the macrocycle and formation of a rearranged linear isomer with the H(2)C═NH(+)-CH(2)- moiety at the N terminus and the -CO-NH(2) amide bond at the C terminus. This rearranged linear structure is much more energetically favorable (-14.0 kcal mol(-1)) than the initially formed imine-protonated linear a(4) ion structure. Furthermore, the barriers to these cyclization and ring-opening reactions are low (8-11 kcal mol(-1)), allowing facile formation of the rearranged linear species in the mass spectrometer. This finding is not limited to 'simple' glycine-containing systems, as evidenced by the IRMPD spectrum of the a(4) ion generated from protonated AAAAA, which shows a stronger tendency toward formation of the energetically favorable (-12.3 kcal mol(-1)) rearranged linear structure with the MeHC═NH(+)-CHMe- moiety at the N terminus and the -CO-NH(2) amide bond at the C terminus. Our results indicate that one needs to consider a complex variety of cyclization and rearrangement reactions in order to decipher the structure and fragmentation pathways of peptide a(n) ions. The implications this potentially has for peptide sequencing are also discussed.     

Ab initio study of the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2) and related C(4)H(10)(2+) isomers.

Ab initio calculations at the MP4(SDTQ)/6-311G//MP2/6-31G level were performed to study the structures and stabilities of the dimer of ethyl cation, (C(2)H(+)(5))(2), and related C(4)H(10)(2+) isomers. Two doubly hydrogen bridged diborane type trans 1 and cis 2 isomers were located as minima. The trans isomer was found to be more favorable than cis isomer by only 0.6 kcal/mol. Several other minima for C(4)H(10)(2+) were also located. However, the global energy minimum corresponds to C-H (C(4) position) protonated 2-butyl cation 10. Structure 10 was computed to be substantially more stable than 1 by 31.7 kcal/mol. The structure 10 was found to be lower in energy than 2-butyl cation 13 by 34.4 kcal/mol.     

The pentacyanocyclopentadienyl system: structures and energetics

In light of the recent silylation route to protonated pentacyanocyclopentadienyl (PCCp) anion reported by Richardson and Reed, PCCp anion, neutral radical, and protonated species have been investigated theoretically. The predicted adiabatic electron affinity for PCCp (ZPVE-corrected value in parentheses) is enormous, 5.56 (5.47) eV. As with the unsubstituted cyclopentadienyl radical, the PCCp radical exhibits a Jahn-Teller distortion from the D(5)(h) symmetry of the aromatic anion, yielding five equivalent (2)B(2) minima that should show uninhibited pseudorotation about a D(5)(h) conical intersection. Formation of the conjugate acid occurs via protonation of the anion at a nitrile nitrogen, which is favored over protonation at a ring carbon by 6.5 kcal/mol, with the preference explained by retention of aromaticity upon protonation at the nitrogen. Possible acid dimer structures have been investigated to evaluate the proposed polymeric acid structure of Richardson and Reed. Our predictions confirm their suggested polymeric structure, but we also present an alternative, self-contained dimer that should be competitive kinetically and thermodynamically. Vibrational frequencies and infrared intensities are predicted, to aid in the experimental identification of several of these species.     

Elimination of water from the carboxyl group of GlyGlyH+

The elimination of water from the carboxyl group of protonated diglycine has been investigated by density functional theory calculations. The resulting structure is identical to the b(2) ion formed in the mass spectrometric fragmentation of protonated peptides (therefore named "b2" in this study). The most stable geometry of the fragment ion ("b2") is an O-protonated diketopiperazine. However, its formation is kinetically disfavored as it requires a free energy of 58.2 kcal/mol. The experimentally observed N-protonated oxazolone is 3.0 kcal/mol less stable. The lowest energy pathway for the formation of the "b2" ion requires a free energy of 37.5 kcal/mol and involves the proton transfer from the amide oxygen of protonated diglycine to the hydroxyl oxygen. Fragmentation initiated by proton transfer from the terminal nitrogen has also a comparable free energy of activation (39.4 kcal/mol). Proton transfer initiating the fragmentation, from the highly basic terminal nitrogen or amide oxygen to the less basic hydroxyl oxygen is feasible at energies reached in usual mass spectrometric experiments. Amide N-protonated diglycine structures are precursors of mainly y(1) ions rather than "b2" ions. In the lowest energy fragmentation channels, proton transfer to the hydroxylic oxygen, bond breaking and formation of an oxazolone ring occur concertedly but asynchronously. Proton transfer to hydroxyl oxygen and cleavage of the corresponding C-O bond take place at the early stages of the fragmentation step, while ring closure to form an oxazolone geometry occurs at the later stages of the transition. The experimentally observed low kinetic energy release is expected to be due to the existence of a strongly hydrogen bonded protonated oxazolone-water complex in the exit channel. Whereas the threshold energy for "b2" ion formation (37.1 kcal/mol) is lower than for the y(1) ion (38.4 kcal/mol), the former requires a tight transition state with an activation entropy, DeltaS++ = -1.2 cal/mol.K and the latter has a loose transition state with DeltaS++ = +8.8 cal/mol.K. This leads to y(1) being the major fragment ion over a wide energy range.     

Protonated carbonyl sulfide: prospects for the spectroscopic observation of the elusive HSCO+ isomer

Spurred by the apparent conflict between ab initio predictions and infrared spectroscopic evidence regarding the relative stability of isomers of protonated carbonyl sulfide, key stationary points on the isomerization surface of HOCS(+) have been examined via systematic extrapolations of ab initio energies. Electron correlation has been accounted for using second-order M?ller-Plesset perturbation theory and coupled cluster theory through triple excitations [CCSD, CCSD(T), and CCSDT] in conjunction with the correlation consistent hierarchy of basis sets, cc-pVXZ (X=D,T,Q,5,6). HSCO(+) is predicted to lie lower in energy than HOCS(+) by 4.86 kcal mol(-1), computed using the focal point extrapolation scheme of Allen and co-workers [J. Chem. Phys. 99, 4638 (1993)] with corrections for anharmonic zero-point vibrational energy, core correlation, non-Born-Oppenheimer, and scalar relativistic effects. A transition state has been located, constituting the barrier to isomerization of HSCO(+) to HOCS(+), lying 68.9 kcal mol(-1) higher in energy than HSCO(+). This is well above predicted exothermicity [DeltaH(r) (o)(0 K)=48.1 kcal mol(-1), cc-pVQZ CCSD(T)] for the reaction considered in the experiments (HSCO(+)+H(2)-->OCS+H(3) (+)). Though proton tunneling will lead to a lower effective barrier, this prediction is consistent with the lack of HSCO(+) in electrical discharges in H(2)OCS, since the relative populations of HOCS(+) and HSCO(+) will depend on the experimental details of the protonation route rather than the relative thermodynamic stability of the isomers. Anharmonic vibrational frequencies and vibrationally corrected rotational constants from cc-pVTZ CCSD(T) cubic and quartic force constants are provided, to aid in the spectroscopic observation of the energetically favorable but apparently elusive HSCO(+) isomer.     

The gas-phase basicity and proton affinity of 1,3,5-cycloheptatriene--energetics,structure and interconversion of dihydrotropylium ions

The hitherto unknown gas-phase basicity and proton affinity of 1,3,5-cycloheptatriene (CHT) have been determined by Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry. Several independent techniques were used in order to exclude ambiguities due to proton-induced isomerisation of the conjugate cyclic C(7)H(9)(+) ions, [CHT + H](+). The gas-phase basicity obtained by the thermokinetic method, GB(CHT) = 799 +/- 4 kJ mol(-1), was found to be identical, within the limits of experimental error, with the values measured by the equilibrium method starting with protonated reference bases, and with the values resulting from the measurements of the individual forward and reverse rate constants, when corrections were made for the isomerised fraction of the C(7)H(9)(+) population. The experimentally determined gas-phase basicity leads to the proton affinity of cycloheptatriene, PA(CHT) = 833 +/- 4 kJ mol(-1), and the heat of formation of the cyclo-C(7)H(9)(+) ion, deltaH(f)(0)([CHT + H](+)) = 884 +/- 4 kJ mol(-1). Ab initio calculations are in agreement with these experimental values if the 1,2-dihydrotropylium tautomer, [CHT + H((1))](+), generated by protonation of CHT at C-1, is assumed to be the conjugate acid, resulting in PA(CHT) = 825 +/- 2 kJ mol(-1) and deltaH(f)(0)(300)([CHT + H((1))](+)) = 892 +/- 2 kJ mol(-1). However, the calculations indicate that protonation of cycloheptatriene at C-2 gives rise to transannular C-C bond formation, generating protonated norcaradiene [NCD + H](+), a valence tautomer being 19 kJ mol(-1) more stable than [CHT + H((1))](+). The 1,4-dihydrotropylium ion, [CHT + H((3))](+), generated by protonation of CHT at C-3, is 17 kJ mol(-1) less stable than [CHT + H((1))](+). The bicyclic isomer [NCD + H](+) is separated by relatively high barriers, 70 and 66 kJ mol(-1) from the monocyclic isomers, [CHT + H((1))](+) and [CHT + H((3))](+), respectively. Therefore, the initially formed 1,2-dihydrotropylium ion [CHT + H((1))](+) does not rearrange to the bicyclic isomer [NCD + H](+) under mild protonation conditions.     

Protonated tert-pentyl dication (C5H122+, isopentane dication)

Structures of the tert-pentyl cation (C(5)H(11)(+)) and its protonated dication (C(5)H(12)(2+), isopentane dication) were studied using ab initio methods at the MP2/cc-pVTZ level. Both C-C and C-H hyperconjugatively stabilized structures 1 and 2 , respectively, were found to be minima on the potential energy surface (PES) of the tert-pentyl cation. Structure 1 was computed to be about as stable as structure 2 (slightly more stable by 0.5 kcal mol(-1)). Inter-conversion between 1 and 2 through transition state 3 has a kinetic barrier of only 1.5 kcal mol(-1). The C-H protonated form (H(3)C)(2)C(+)CH(2)CH(4)(+)4 was found to be the global minimum for the protonated tert-pentyl dication. Charges and (13)C NMR chemical shifts of the dication 4 were calculated and compared to those of monocation 1 to study the effect of the additional charge in the dication.     

Polyglycine conformational analysis: calculated vs experimental gas-phase basicities and proton affinities

Structures of neutral and protonated polyglycines (Gly(n) and Gly(n)H(+) with n = 1-6) in the vicinity of global energy minima were calculated using the density functional theory at the B3LYP/6-311++G** (A) and B3LYP/6-31+G** (B) levels. Ninety-three structures were chosen for conformation and protonation studies. Geometries of the peptides are found to vary from open chains to multiple rings. Intramolecular hydrogen bonding is deduced to be the driving force for conformational stability. The preferred protonation sites are shown to be the terminal nitrogen atom and its adjacent amide oxygen atom. Structural series are developed according to geometrical form, hydrogen bonding, and protonation site. Physical factors that influence the relative electronic and thermodynamic stabilities of different structural series are examined. To obtain ab initio values of highest quality for gas-phase basicity (GB) and proton affinity (PA), electronic energies for n = 1-6 and thermal corrections to Gibbs free energy and enthalpy for n = 1-3 were calculated at level A, supplemented by thermal corrections for n = 4-6 at level B. Calculated GB and PA values are compared with mass spectral results obtained by the kinetic method (KM) and reaction bracketing (RB). The KM results and the ab initio values derived from structurally compatible pairs of lowest free energies are generally in good agreement, but the RB results for GB are lower by 2-8 kcal/mol for n = 2-6. Several reaction pathways are proposed to elucidate the experimental results. On the basis of theoretical structures consistent with the measurements, it is concluded that KM mostly samples the neutral and protonated structures of highest populations at thermal equilibrium, whereas RB targets those with sterically most accessible sites for protonation and deprotonation.     

The maximum number of carbonyl groups around an Ru6C polyhedral cluster: hexanuclear ruthenium carbonyl carbides

Octahedral, trigonal prismatic, and capped square pyramidal structures have been optimized for the Ru(6)C(CO)(n) clusters (15 ≤ n ≤ 20) using density functional theory. The experimentally known very stable Ru(6)C(CO)(17) is predicted to have an octahedral structure in accord with experiment as well as the Wade-Mingos rules. The stability of Ru(6)C(CO)(17) is indicated by its high carbonyl dissociation energy of ~37 kcal mol(-1) and the high energy of ~33 kcal mol(-1) required for disproportionation into Ru(6)C(CO)(18) + Ru(6)C(CO)(16). Theoretical calculations predict a doubly carbonyl bridged octahedral Ru(6)C(CO)(17) structure to be ~0.7 kcal mol(-1) more stable than the experimentally observed singly bridged structure. A trigonal prismatic Ru(6)C(CO)(19) cluster isoelectronic with the known Co(6)C(CO)(15)(2-) dianion does not appear to be viable as indicated by a low carbonyl dissociation energy of 8.8 kcal mol(-1) and a required energy of only 4.9 kcal mol(-1) for disproportionation into Ru(6)C(CO)(20) + Ru(6)C(CO)(18). The predicted instability of Ru(6)C(CO)(n) (n ≥ 18) derivatives suggests a maximum of 17 external carbonyl groups around a stable polyhedral Ru(6)C structure.     

Is there a generalized reverse anomeric Effect? substituent and solvent effects on the configurational equilibria of neutral and protonated N-arylglucopyranosylamines and N-aryl-5-thioglucopyranosylamines

The effects of substitution and solvent on the configurational equilibria of neutral and protonated N-(4-Y-substituted-phenyl) peracetylated 5-thioglucopyranosylamines (Y = OMe, H, CF(3), NO(2)) 1-4 and N-(4-Y-substituted-phenyl) peracetylated glucopyranosylamines (Y = OMe, H, NO(2)) 9-11 are described. The configurational equilibria were determined by direct integration of the resonances of the individual isomers in the (1)H NMR spectra after equilibration of both alpha- and beta-isomers. The equilibrations of the neutral compounds 1-4 in CD(3)OD, CD(3)NO(2), and (CD(3))(2)CO were achieved by HgCl(2) catalysis and those of the neutral compounds 9-11 in CD(2)Cl(2) and CD(3)OD by triflic acid catalysis. The equilibrations of the protonated compounds in both the sulfur series (solvents, CD(3)OD, CD(3)NO(2), (CD(3))(2)CO, CDCl(3), and CD(2)Cl(2)) and oxygen series (solvents, CD(2)Cl(2) and CD(3)OD) were achieved with triflic acid. The substituent and solvent effects on the equilibria are discussed in terms of steric and electrostatic effects and orbital interactions associated with the endo-anomeric effect. A generalized reverse anomeric effect does not exist in neutral or protonated N-aryl-5-thioglucopyranosylamines and N-arylglucopyranosylamines. The anomeric effect ranges from 0.85 kcal mol(-)(1) in 2 to 1.54 kcal mol(-)(1) in 10. The compounds 1-4 and 9-11 show an enhanced endo-anomeric effect upon protonation, ranging from 1.73 kcal mol(-)(1) in 2 to 2.57 kcal mol(-)(1) in 10. We estimate the increase in the anomeric effect upon protonation of 10 to be approximately 1.0 kcal mol(-)(1). However, this effect is offset by steric effects due to the associated counterion which we estimate to be approximately 1.2 kcal mol(-)(1). The values of K(eq)(axial-equatorial) in protonated 1-4 increase in the order OMe < H < CF(3) < NO(2), in agreement with the dominance of steric effects (due to the counterion) over the endo-anomeric effect. The values of K(eq)(axial-equatorial) in protonated 9-11 show the trend OMe > H < NO(2) that is explained by the balance of the endo-anomeric effect and steric effects in the individual compounds. The trends in the values of the C(1)-H(1) coupling constants for 1-4 and the corresponding deacetylated compounds 5-8 as a function of substituent and alpha- or beta-configuration are discussed in terms of the Perlin effect and the interplay of the endo- and exo-anomeric effects.     

N-acylamidine palladium(II) complexes. Synthesis,structures, and catalytic activity for Suzuki-Miyaura cross-coupling reactions: experiments and DFT calculations

N-Acylamidines 2 are easily prepared by acylation of amidines 1. Upon treatment with PdCl(2)(PhCN)(2), they form 2:1 PdCl(2) complexes 3 with trans configuration, acting as monodentate ligands via the nitrogen atom remote from the carbonyl group. The structures of the complexes 3a-c in the solid state were obtained by X-ray crystallography. As studied by DFT calculations on 2:1 model complexes, many isomeric structures 6 with different conformations and configurations are associated with local minima on the energy hyperface within an energy range of 4 kcal/mol. The complexes 3 show very high catalytic activity in Suzuki-Miyaura cross-coupling reactions, either as isolated crystals or prepared in situ without isolation.     

Proton grease: an acid accelerated molecular rotor

A molecular rotor was designed that rotates 7 orders of magnitude faster upon protonation. The quinoline rotor is based on a rigid N-arylimide framework that displays restricted rotation due to steric interaction between the quinoline nitrogen and imide carbonyls. At rt (23 °C), the rotor rotates slowly (t(1/2) = 26 min, ΔG(?) = 22.2 kcal/mol). However, upon addition of 3.5 equiv of acid the rotor rotates rapidly (t(1/2) = 2.0 × 10(-4) s, ΔG(?) = 12.9 kcal/mol). Mechanistic studies show that this dramatic acid catalyzed change is due to stabilization of the planar transition state by the formation of an intramolecular hydrogen bond between the protonated quinoline nitrogen (N(+)-H) and an imide carbonyl (O═C). The acid catalyzed acceleration is reversible and can be stopped by addition of base.     

Effect of the N-terminal basic residue on facile Cα-C bond cleavages of aromatic-containing peptide radical cations

Fragmentation of radical cationic peptides [R(G)(n-2)X(G)(7-n)]˙(+) and [R(G)(m-2)XG]˙(+) (X = Phe or Tyr; m = 2-5; n = 2-7) leads selectively to a(n)(+) product ions through in situ C(α)-C peptide backbone cleavage at the aromatic amino acid residues. In contrast, substituting the arginine residue with a less-basic lysine residue, forming [K(G)(n-2)X(G)(7-n)]˙(+) (X = Phe or Tyr; n = 2-7) analogs, generates abundant b-y product ions; no site-selective C(α)-C peptide bond cleavage was observed. Studying the prototypical radical cationic tripeptides [RFG]˙(+) and [KFG]˙(+) using low-energy collision-induced dissociation and density functional theory, we have examined the influence of the basicity of the N-terminal amino acid residue on the competition between the isomerization and dissociation channels, particularly the selective C(α)-C bond cleavage viaβ-hydrogen atom migration. The dissociation barriers for the formation of a(2)(+) ions from [RFG]˙(+) and [KFG]˙(+)via their β-radical isomers are comparable (33.1 and 35.0 kcal mol(-1), respectively); the dissociation barrier for the charge-induced formation of the [b(2)- H]˙(+) radical cation from [RFG]˙(+)via its α-radical isomer (39.8 kcal mol(-1)) was considerably higher than that from [KFG]˙(+) (27.2 kcal mol(-1)). Thus, the basic arginine residue sequesters the mobile proton to promote the charge-remote selective C(α)-C bond cleavage by energetically hindering the competing charge-induced pathways.     

Doubly charged protonated a ions derived from small peptides

Protonated a(2) and a(3) (therefore doubly charged) ions in which both charges lie on the peptide backbone are formed in collision-induced dissociations of [La(III)(peptide)(CH(3)CN)(m)](3+) complexes. Abundant (a(3)+H)(2+) ions are formed from triproline (PPP) and peptides with a proline residue at the N-terminus; these peptides are the most effective in producing ions of the type (a(2)+H)(2+) and (a(3)+H)(2+). A systematic study of the effect of the location of the proline residue and other residues of aliphatic amino acids on the generation of protonated a ions is reported. Density functional theory calculations at B3LYP/6-311++G(d,p) gave the proton affinity of the a(3) ion derived from PPP to be 167.6 kcal mol(-1), 2.6 kcal mol(-1) higher than that of water. The protonated a(2) ions of diglycine and diproline and a(3) ions of triglycine have lower proton affinities and are only observed in lower abundances, possibly due to proton transfer to water in ion-molecule reactions.     

Arachno, nido, and closo aromatic isomers of the Li6B6H6 molecule

We analyzed chemical bonding in low-lying isomers of the recently computationally predicted B(6)H(6)Li(6) molecule. According to our calculations the benzene-like B(6)H(6)Li(6) (D(2h), (1)A(1g)) arachno structure with the planar aromatic B(6)H(6)(6-) anion is the most stable one. A nido isomer with two aromatic B(6)H(6)(4-) (pentagonal pyramid) and Li(3)(+) (triangular) moieties, which can be considered as derived from the global minimum structure through a two-electron intramolecular transfer from B(6)H(6)(6-) to three Li(+) cations, was found to be 10.7 kcal/mol higher in energy. A closo isomer with three aromatic moieties (octahedral B(6)H(6)(2-) and two Li(3)(+)) was found to be 31.3 kcal/mol higher in energy than the global minimum. Another isomer with three aromatic moieties (two B(3)H(3)(2-) and Li(3)(+)) was found to be substantially higher in energy (74.4 kcal/mol). Thus, the intramolecular electron transfers from the highly charged B(6)H(6)(6-) anion to cations are not favorable for the B(6)H(6)Li(6) molecule, even when a formation of three-dimensional aromatic B(6)H(6)(2-) anion and two sigma-aromatic Li(3)(+) cations occurs in the closo isomer.     

Benchmark ab initio proton affinity and gas-phase basicity of α-alanine based on coupled-cluster theory and statistical mechanics

We determine the proton affinity (PA) and gas-phase basicity (GB) of amino acid α-alanine at a chemically accurate level by performing explicitly-correlated CCSD(T)-F12b/aug-cc-pVDZ geometry optimizations and normal mode vibrational frequency calculations as well as CCSD(T)-F12b/aug-cc-pVTZ energy computations at the possible neutral and protonated geometries. Temperature effects at 298.15 K considering translational, rotational, and vibrational enthalpy and entropy corrections are obtained via standard statistical mechanics utilizing the molecular geometries and the harmonic vibrational energy levels. Both the amino nitrogen (N) and the carbonyl oxygen (O) atoms are proven to be potential protonation sites and a systematic conformational search reveals 3 N- and 9 O-protonated conformers in the 0.00–7.88 and 25.43–30.43 kcal/mol energy ranges at 0 K, respectively. The final computed PA and GB values at (0)298.15 K in case of N-protonation are (214.47)216.80 and 207.07 kcal/mol, respectively, whereas the corresponding values for O-protonation are (189.04)190.63 and 182.31 kcal/mol. The results of the benchmark high-level coupled-cluster computations are utilized to assess the accuracy of several lower-level cost-effective methods such as MP2 and density functional theory with various functionals (SOGGA11-X, M06-2X, PBE0, B3LYP, M06, TPSS).     

Structure of the [M + H - H2O]+ ion from tetraglycine: a revisit by means of density functional theory and isotope labeling

Collision-induced dissociations of protonated (18)O-labeled tetraglycines labeled separately at either the first or the second amide bond established that water loss from the backbone occurs from the N-terminal residue. Density functional theory at B3LYP/6-311++G(d,p) predicted that the low-energy [G(4) + H - H(2)O](+) product ion is an N(1)-protonated 3,5-dihydro-4H-imidazol-4-one. The ion at the lowest energy, III, is 24.8 kcal mol(-1) lower than the protonated oxazole structure, II, proposed by Bythell et al. (J. Phys. Chem A2010, 114, 5076-5082). In addition, structure III has a predicted IR spectrum that provides a better match with the published experimental IRMPD spectrum than that of structure II.