Proton migration and tautomerism in protonated triglycine.

Proton migration in protonated glycylglycylglycine (GGG) has been investigated by using density functional theory at the B3LYP/6-31++G(d,p) level of theory. On the protonated GGG energy hypersurface 19 critical points have been characterized, 11 as minima and 8 as first-order saddle points. Transition state structures for interconversion between eight of these minima are reported, starting from a structure in which there is protonation at the amino nitrogen of the N-terminal glycyl residue following the migration of the proton until there is fragmentation into protonated 2-aminomethyl-5-oxazolone (the b(2) ion) and glycine. Individual free energy barriers are small, ranging from 4.3 to 18.1 kcal mol(-)(1). The most favorable site of protonation on GGG is the carbonyl oxygen of the N-terminal residue. This isomer is stabilized by a hydrogen bond of the type O-H.N with the N-terminal nitrogen atom, resulting in a compact five-membered ring. Another oxygen-protonated isomer with hydrogen bonding of the type O-H.O, resulting in a seven-membered ring, is only 0.1 kcal mol(-)(1) higher in free energy. Protonation on the N-terminal nitrogen atom produces an isomer that is about 1 kcal mol(-)(1) higher in free energy than isomers resulting from protonation on the carbonyl oxygen of the N-terminal residue. The calculated energy barrier to generate the b(2) ion from protonated GGG is 32.5 kcal mol(-)(1) via TS(6-->7). The calculated basicity and proton affinity of GGG from our results are 216.3 and 223.8 kcal mol(-)(1), respectively. These values are 3-4 kcal mol(-)(1) lower than those from previous calculations and are in excellent agreement with recently revised experimental values.     

Theoretical proton affinity and fluoride affinity of nerve agent VX

Proton affinity and fluoride affinity of nerve agent VX at all of its possible sites were calculated at the RI-MP2/cc-pVTZ//B3LYP/6-31G* and RI-MP2/aug-cc-pVTZ//B3LYP/6-31+G* levels, respectively. The protonation leads to various unique structures, with H(+) attached to oxygen, nitrogen, and sulfur atoms; among which the nitrogen site possesses the highest proton affinity of -ΔE ～ 251 kcal/mol, suggesting that this is likely to be the major product. In addition some H(2), CH(4) dissociation as well as destruction channels have been found, among which the CH(4) + [Et-O-P(═O)(Me)-S-(CH(2))(2)-N(+)(iPr)═CHMe] product and the destruction product forming Et-O-P(═O)(Me)-SMe + CH(2)═N(+)(iPr)(2) are only 9 kcal/mol less stable than the most stable N-protonated product. For fluoridization, the S-P destruction channel to give Et-O-P(═O)(Me)(F) + [S-(CH(2))(2)-N-(iPr)(2)](-) is energetically the most favorable, with a fluoride affinity of -ΔE ～ 44 kcal. Various F(-) ion-molecule complexes are also found, with the one having F(-) interacting with two hydrogen atoms in different alkyl groups to be only 9 kcal/mol higher than the above destruction product. These results suggest VX behaves quite differently from surrogate systems.     

On the unusual stability of valence anions of thymine based on very rare tautomers: A computational study

We characterized anionic states of thymine using various electronic structure methods, with the most accurate results obtained at the CCSD(T)/aug-cc-pVDZ level of theory followed by extrapolations to complete basis set limits. We found that the most stable anion in the gas phase is related to an imino-oxo tautomer, in which the N1H proton is transferred to the C5 atom. This valence anion, aT(c5)(nl), is characterized by an electron vertical detachment energy (VDE) of 1251 meV and it is adiabatically stable with respect to the canonical neutral nT(can) by 2.4 kcal/mol. It is also more stable than the dipole-bound (aT(dbs)(can)), and valence anion aT(val)(can) of the canonical tautomer. The VDE values for aT(dbs)(can)and T(val)(can) are 55 and 457 meV, respectively. Another, anionic, low-lying imino-oxo tautomer with a VDE of 2458 meV has a proton transferred from N3H to C5 aT(c5)(n3). It is less stable than aT(val)(can) by 3.3 kcal/mol. The mechanism of formation of anionic tautomers with the carbons C5 or C6 protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to C5. The six-member ring structure of the anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Within the PCM hydration model, the low-lying valence anions become adiabatically bound with respect to the canonical neutral; becomes the most stable, being followed by aT(c5)(nl), aT(c5)(n3), aT(can), and aT(c5)(nl).     

Stabilization of very rare tautomers of 1-methylcytosine by an excess electron

We characterized valence anionic states of 1-methylcytosine using various electronic structure methods. We found that the most stable valence anion is related to neither the canonical amino-oxo nor a rare imino-oxo tautomer, in which a proton is transferred from the N4 to N3 atom. Instead, it is related to an imino-oxo tautomer, in which the C5 atom is protonated. This anion is characterized by an electron vertical detachment energy (VDE) of 2.12 eV and it is more stable than the anion based on the canonical tautomer by 1.0 kcal/mol. The latter is characterized by a VDE of 0.31 eV. Another unusual low-lying imino-oxo tautomer with a VDE of 3.60 eV has the C6 atom protonated and is 3.6 kcal/mol less stable than the anion of the canonical tautomer. All these anionic states are adiabatically unbound with respect to the canonical amino-oxo neutral, with the instability of 5.8 kcal/mol for the most stable valence anion. The mechanism of formation of anionic tautomers with carbon atoms protonated may involve intermolecular proton transfer or dissociative electron attachment to the canonical neutral tautomer followed by a barrier-free attachment of a hydrogen atom to the C5 or C6 atom. The six-member ring structure of anionic tautomers with carbon atoms protonated is unstable upon an excess electron detachment. Indeed the neutral systems collapse without a barrier to a linear or a bicyclo structure, which might be viewed as lesions to DNA or RNA. Within the PCM hydration model, the anions become adiabatically bound with respect to the corresponding neutrals, and the two most stable tautomers have a carbon atom protonated.     

Intriguing roles of reactive intermediates in dissociation chemistry of N-phenylcinnamides

In mass spectrometry of protonated N-phenylcinnamides, the carbonyl oxygen is the thermodynamically most favorable protonation site and the added proton is initially localized on it. Upon collisional activation, the proton transfers from the carbonyl oxygen to the dissociative protonation site at the amide nitrogen atom or the α-carbon atom, leading to the formation of important reactive intermediates. When the amide nitrogen atom is protonated, the amide bond is facile to rupture to form ion/neutral complex 1, [RC(6)H(4)CH[double bond, length as m-dash]CHCO(+)/aniline]. Besides the dissociation of the complex, proton transfer reaction from the α-carbon atom to the nitrogen atom within the complex takes place, leading to the formation of protonated aniline. The presence of electron-withdrawing groups favored the proton transfer reaction, whereas electron-donating groups strongly favored the dissociation (aniline loss). When the proton transfers from the carbonyl oxygen to the α-carbon atom, the cleavage of the C(α)-CONHPh bond results in another ion/neutral complex 2, [PhNHCO(+)/RC(6)H(4)CH[double bond, length as m-dash]CH(2)]. However, in this case, electron-donating groups expedited the proton transfer reaction from the charged to the neutral partner to eliminate phenyl isocyanate. Besides the cleavage of the C(α)-CONHPh bond, intramolecular nucleophilic substitution (a nucleophilic attack of the nitrogen atom at the β-carbon) and stepwise proton transfer reactions (two 1,2-H shifts) also take place when the α-carbon atom is protonated, resulting in the loss of ketene and RC(6)H(5), respectively. In addition, the H/D exchanges between the external deuterium and the amide hydrogen, vinyl hydrogens and the hydrogens of the phenyl rings were discovered by D-labeling experiments. Density functional theory-based (DFT) calculations were performed to shed light on the mechanisms for these reactions.     

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.     

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.     

Proton affinity of β-oxalylaminoalanine (BOAA): Incorporation of direct entropy correction into the single-reference kinetic method

A new version of the single-reference-extended kinetic method is presented in which direct entropy correction is incorporated. Results of calibration experiments with the monodentate base pyridine and the bidentate base ethylenediamine are presented for which the method provides proton affinities in excellent agreement with published values and reasonable predictions for the protonation entropies. The method is then used to determine the proton affinity and protonation entropy of the non-protein amino acid beta-oxalylaminoalanine (BOAA). The PA of BOAA is found to be 933.1 +/- 7.8 kJ/mol and a prediction for the protonation entropy of -39 J mol(-1) K(-1) is also obtained, indicating a significant degree of intramolecular hydrogen bonding in the protonated form. These results are supported by hybrid density functional theory calculations at the B3LYP/6-311++G**//B3LYP/6-31+G* level. They indicate that the preferred site of protonation is the alpha-nitrogen atom (PA = 935.0 kJ/mol) and that protonated BOAA has a strong hydrogen bond between the hydrogen on the alpha-amino group and one of the carbonyl oxygen atoms on the side chain.     

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).     

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.     

Gas-phase structure of protonated histidine and histidine methyl ester: combined experimental mass spectrometry and theoretical ab initio study

Gas-phase H/D exchange experiments with CD3OD and D2O and quantum chemical ab initio G3(MP2) calculations were carried out on protonated histidine and protonated histidine methyl ester in order to elucidate their bonding and structure. The H/D exchange experiments show that both ions have three equivalent fast hydrogens and one appreciably slower exchangeable hydrogen assigned to the protonated amino group participating in a strong intramolecular hydrogen bond (IHB) with the nearest N(sp2) nitrogen of the imidazole fragment and to the distal ring NH-group, respectively. It is taken for granted that the proton exchange in the IHB is much faster than the H/D exchange. Unlike in other protonated amino acids (glycine, proline, phenylalanine, tyrosine, and tryptophan) studied earlier, the exchange rate of the carboxyl group in protonated histidine is slower than that of the amino group. The most stable conformers and the enthalpies of neutral and protonated histidine and its methyl ester are calculated at the G3(MP2) level of theory. It is shown that strong intramolecular hydrogen bonding between the amino group and the imidazole ring nitrogen sites is responsible for the stability and specific properties of the protonated histidine. It is found that the proton fluctuates between the amino and imidazole groups in the protonated form across an almost vanishing barrier. Proton affinity (PA) of histidine calculated by the G3(MP2) method is 233.2 and 232.4 kcal mol(-1) for protonation at the imidazole ring and at the amino group nitrogens, respectively, which is about 3-5 kcal mol(-1) lower than the reported experimental value.     

Acidity of adenine and adenine derivatives and biological implications. A computational and experimental gas-phase study

The gas-phase acidities of adenine, 9-ethyladenine, and 3-methyladenine have been investigated for the first time, using computational and experimental methods to provide an understanding of the intrinsic reactivity of adenine. Adenine is found to have two acidic sites, with the N9 site being 19 kcal mol(-1) more acidic than the N10 site; the bracketed acidities are 333 +/- 2 and 352 +/- 4 kcal mol(-1), respectively. Because measurement of the less acidic site can be problematic, we benchmarked the adenine N10 measurement by bracketing the acidity of 9-ethyladenine, which has the N9 site blocked and allows for exclusive measurement of the N10 site. The acidity of 9-ethyladenine brackets to 352 +/- 4 kcal mol(-1), comparable to that of the N10 site of the parent adenine. Calculations and experiments with 3-methyladenine, a harmful mutagenic nucleobase, uncovered the surprising result that the most commonly written tautomer of 3-methyladenine is not the most stable in the gas phase. We have found that the most stable tautomer is the "N10 tautomer" 10, as opposed to the imine tautomer 3. The bracketed acidity of 10 is 347 +/- 4 kcal mol(-1). Since 10 is not a viable species in DNA, 3 is a likely tautomer; calculations indicate that this form has an extremely high acidity (320-323 kcal mol(-1)). The biological implications of these results, particularly with respect to enzymes that cleave alkylated bases from DNA, are discussed.     

A theoretical study of potassium cation binding to glycylglycine (GG) and alanylalanine (AA) dipeptides

By combining Monte Carlo conformational search technique with high-level density functional calculations, the geometry and energetics of K(+) interaction with glycylglycine (GG) and alanylalanine (AA) were obtained for the first time. The most stable K(+)-GG and K(+)-AA complexes are in the charge-solvated (CS) form with K(+) bound to the carbonyl oxygens of the peptide backbone, and the estimated 0 K binding affinities (DeltaH(0)) are 152 and 157 kJ mol(-1), respectively. The K(+) ion is in close alignment with the molecular dipole moment vector of the bound ligand, that is, electrostatic ion-dipole interaction is the key stabilizing factor in these complexes. Furthermore, the strong ion-dipole interaction between K(+) and the amide carbonyl oxygen atom of the peptide bond is important in determining the relative stabilities of different CS binding modes. The most stable zwitterionic (ZW) complex involves protonation at the amide carbonyl oxygen atom and is approximately 48 kJ mol(-1) less stable than the most stable CS form. The usefulness of proton affinity (PA) as a criterion for estimating the relative stability of ZW versus CS binding modes is examined. The effect of chain length and the nature of metal cations on cation-dipeptide interactions are discussed. Based on results of this study, the interaction of K(+) with longer peptides consisting of aliphatic amino acids are rationalized.     

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.     

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.     

ab initio Calculation of the stereoelectronic structure of 1-(1-trichlorogermylethyl)-pyrrolidin-2-one*

RHF/6-31G(d) and MP2/6-31G(d) calculations were carried out to study the stereoelectronic structure of 1-(1-trichlorogermylethyl)pyrrolidin-2-one with a pentacoordinated germanium atom. These results were compared with the X-ray diffraction structural analysis data. Upon formation of the Ge ← O coordination bond in this molecule, the electron density of all the atoms of the coordination polyhedron of the germanium atom, including the oxygen atom, increases, especially the axial chlorine atom, while the electron density of the germanium, nitrogen, and carbonyl group carbon atoms decreases. Different polarization of all three valence p-orbitals of each Cl atom of this molecule was established. ³⁵Cl nuclear quadrupole resonance spectrum parameters were evaluated. The molecule also has stable form, in which the germanium atom is tetracoordinated. The total energy of this form is 2.7 kcal/mol higher than for the structure with a pentacoordinated germanium atom.     

Proton and sodium cation affinities of harpagide: a computational study

The aim of this work was to estimate the proton and sodium cation affinities of harpagide (Har), an iridoid glycoside responsible for the antiinflammatory properties of the medicinal plant Harpagophytum. Monte Carlo conformational searches were performed at the semiempirical AM1 level to determine the most stable conformers for harpagide and its protonated and Na+-cationized forms. The 10 oxygen atoms of the molecule were considered as possible protonation and cationization sites. Geometry optimizations were then refined at the DFT B3LYP/6-31G level from the geometries of the most stable conformers found. Final energetics were obtained at the B3LYP/6-311+G(2d,2p)//B3LYP/6-31G level. The proton and sodium ion affinities of harpagide have been estimated at 223.5 and 66.0 kcal/mol, respectively. Since harpagide mainly provides HarNa+ ions in electrospray experiments, the DeltarG298 associated with the reaction of proton/sodium exchange between Har and methanol, MeOHNa+ + HarH+ --> MeOH2+ + HarNa+ (1), has been calculated; it has been estimated to be 1.9 kcal/mol. Complexing a methanol molecule to each reagent and product of reaction 1 makes the reaction become exothermic by 1.7 kcal/mol. These values are in the limit of the accuracy of the method and do not allow us to conclude definitely whether the reaction is endo- or exothermic, but, according to these very small values, the cation exchange reaction is expected to proceed easily in the final stages of the ion desolvation process.     

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.     

Formation and stability of enolates of acetamide and acetate anion: an Eigen plot for proton transfer at alpha-carbonyl carbon

Second-order rate constants were determined in D(2)O for deprotonation of acetamide, N,N-dimethylacetamide, and acetate anion by deuterioxide ion and for deprotonation of acetamide by quinuclidine. The values of k(B) = 4.8 x 10(-8) M(-1) s(-1) for deprotonation of acetamide by quinuclidine (pK(BH) = 11.5) and k(BH) = 2-5 x 10(9) M(-1) s(-1) for the encounter-limited reverse protonation of the enolate by protonated quinuclidine give pK(a)(C) = 28.4 for ionization of acetamide as a carbon acid. The limiting value of k(HOH) = 1 x 10(11) s(-1) for protonation of the enolate of acetate anion by solvent water and k(HO) = 3.5 x 10(-9) M(-1) s(-1) for deprotonation of acetate anion by HO(-) give pK(a)(C) approximately 33.5 for acetate anion. The change in the rate-limiting step from chemical proton transfer to solvent reorganization results in a downward break in the slope of the plot of log k(HO) against carbon acid pK(a) for deprotonation of a wide range of neutral alpha-carbonyl carbon acids by hydroxide ion, from -0.40 to -1.0. Good estimates are reported for the stabilization of the carbonyl group relative to the enol tautomer by electron donation from alpha-SEt, alpha-OMe, alpha-NH(2), and alpha-O(-) substituents. The alpha-NH(2) and alpha-OMe groups show similar stabilizing interactions with the carbonyl group, while the interaction of alpha-O(-) is only 3.4 kcal/mol more stabilizing than for alpha-OH. We propose that destabilization of the enolate intermediates of enzymatic reactions results in an increasing recruitment of metal ions by the enzyme to provide electrophilic catalysis of enolate formation.