Proton affinity of proline and modified prolines using the kinetic method: role of the conformation investigated by ab initio calculations

The proton affinities of proline, cis-3-methylproline and cis-3-ethylproline have been measured by the kinetic method using an ion trap instrument; the values obtained are 936, 940.5, and 943 kJ mol(-1), respectively. The experimental values are consistent with those obtained by high-level ab initio calculations (B3LYP/6-31+G*//B3LYP/6-31G* and B3P86/6-31+G*//B3LYP/6-31G*). Several conformations of neutral and protonated proline were considered, in particular the endo and exo ring structure and the position of the carboxyl group. These results show the importance of the position of the hydrogen atom of the carboxyl group in determining the most stable protonated proline structure.     

Conformational preference and cis-trans isomerization of 4(R)-substituted proline residues

We report here the conformational preference and prolyl cis-trans isomerization of 4(R)-substituted proline dipeptides, N-acetyl-N'-methylamides of 4(R)-hydroxy-L-proline and 4(R)-fluoro-L-proline (Ac-Hyp-NHMe and Ac-Flp-NHMe, respectively), studied at the HF/6-31+G(d), B3LYP/6-31+G(d), and B3LYP/6-311++G(d,p) levels of theory. The 4(R)-substitution by electron-withdrawing groups did not result in significant changes in backbone torsion angles as well as endocyclic torsion angles of the prolyl ring. However, the small changes in backbone torsion angles phi and psi and the decrease of bond lengths r(Cbeta-Cgamma) or r(Cgamma-Cdelta) appear to induce the increase of the relative stability of the trans up-puckered conformation and to alter the relative stabilities of transition states for prolyl cis-trans isomerization. Solvation free energies of local minima and transition states in chloroform and water were calculated using the conductor-like polarizable continuum model at the HF/6-31+G(d) level of theory. The population of trans up-puckered conformations increases in the order Ac-Pro-NHMe < Ac-Hyp-NHMe < Ac-Flp-NHMe in chloroform and water. The increase in population for trans up-puckered conformations in solution is attributed to the increase in population for the polyproline-II-like conformations with up puckering. The barriers DeltaGct++ to prolyl cis-to-trans isomerization for Ac-Hyp-NHMe and Ac-Flp-NHMe increase as the solvent polarity increases, as seen for Ac-Pro-NHMe. In particular, it was identified that the cis-trans isomerization proceeds through the clockwise rotation about the prolyl peptide bond for Ac-Hyp-NHMe and Ac-Flp-NHMe in chloroform and water, as seen for Ac-Pro-NHMe.     

Theoretical Study of Halogen Effect in Isomerization of 2-Halo-[9]-annulen Anion at the DFT Level

Cis-trans isomerization of [9]-annulenanion(1) and its 2-fluoro-,2-chloro-and 2-bromo-derivatives(2,3 and 4,respectively) were investigated at the HF/6-31G* and B3LYP/6-311++G** levels of theory.Cis,cis,cis,cis structures appear more stable than their corresponding cis,cis,cis,trans-isomers.The relative height of energy barriers for cis-trans isomerization is:2cis > 1cis > 3cis > 4cis.This trend for the reverse trans-cis isomerization follows the electronegativity of the substituent at C-2(2trans > 3trans > 4trans > 1trans).     

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.     

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.     

Theoretical studies on the butene double bond isomerization catalyzed by 5-H of 1-ethyl-3-methyl-imidazolium fluoride

The isomerization of 1-butene to trans-2-butene catalyzed by 5-H proton of 1-ethyl-3-methyl-imidazolium fluoride (EMImF) has been studied with density functional theory of quantum chemistry. The equilibrium states geometries and transition state geometry are optimized at the levels of B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p), respectively. The apparent activation barrier of isomerization is about 208 kJ/mol theoretically. It indicates that the 5-H proton on the imidazole ring of EMImF has certain catalytic activity to the butene double bond isomerization, which is similar to that of the 4-H proton. According to the data of intrinsic coordinate path, it can be determined that the isomerization is an elementary course and the hydrogen exchange of butene with EMImF is synergetic.     

The rate enhancement for prolyl cis-to-trans isomerization of cyclic CPFC peptide is caused by an increase in the vibrational entropy of the transition state

The conformational preferences and prolyl cis-trans isomerization of oxidized and reduced Ac-Cys-Pro-Phe-Cys-NH2 (CPFC peptides) have been carried out using the ab initio HF/6-31+G(d) and hybrid density functional B3LYP/6-311++G(d,p) levels of theory. The most preferred conformations of oxidized and reduced CPFC peptides with the trans prolyl peptide bond have a type-I beta-turn for the Pro-Phe sequence in common. In particular, the transition states for both forms are stabilized by the intramolecular hydrogen bonds between the prolyl nitrogen and the N-H group of the Phe3 residue. The rotational barrier DeltaGct to the cis-to-trans isomerization for the oxidized CPFC peptide is calculated to be 19.37 kcal/mol at the B3LYP/6-311++G(d,p)//HF/6-31+G(d) level of theory, which is lower by 0.88 kcal/mol than that of the reduced CPFC peptide. This may indicate that the rate constant kc-->t of the prolyl cis-to-trans isomerization for the oxidized form is about 4 times larger than that of the reduced form, which is reasonably consistent with the value deduced from NMR experiments. In particular, the increase in vibrational entropy for the transition state of the oxidized form over that of the reduced form contributes to enhance the rate constant for the prolyl cis-to-trans isomerization of the oxidized form.     

Gas-phase reactions of protonated tryptophan

The gas phase reactions of protonated tryptophan have been examined in a quadrupole ion trap using a combination of collision induced dissociation, hydrogen-deuterium exchange, regiospecific deuterium labeling and molecular orbital calculations (at the B3LYP/6-31G* level of theory). The loss of ammonia from protonated tryptophan is observed as the primary fragmentation pathway, with concomitant formation of a [M + H - NH(3)](+) ion by nucleophilic attack from the C3 position of the indole side chain. Hydrogen-deuterium exchange and regiospecific deuterium labeling reveals that scrambling of protons in the C2 and C4 positions of the indole ring, via intramolecular proton transfer from the thermodynamically preferred site of protonation at the amino nitrogen, precedes ammonia loss. Molecular orbital calculations have been employed to demonstrate that the activation barriers to intramolecular proton transfer are lower than that for NH(3) loss.     

Protonation thermochemistry of selected hydroxy- and methoxycarbonyl molecules

The gas-phase basicities of a representative set of hydroxy- and methoxycarbonyl compounds (hydroxyacetone, 1, 3-hydroxybutanone, 2, 3-hydroxy-3-methylbutanone, 3, 1-hydroxy-2-butanone, 4, 4-hydroxy-2-butanone, 5, 5-hydroxy-2-pentanone, 6, methoxyacetone, 7, 3-methoxy-2-butanone, 8, 4-methoxy-2-butanone, 9, and 5-methoxy-2-pentanone, 10) were experimentally determined by the equilibrium method using Fourier transform ion cyclotron resonance and high-pressure mass spectrometry techniques. The latter method allows the measurement of proton transfer equilibrium constants at various temperatures and thus the estimate of both the proton affinities and the protonation entropies of the relevant species. Quantum chemical calculations at the G3 and the B3LYP/6-311+G(3df,2p)//6-31G(d) levels of theory were undertaken in order to find the most stable structures of the neutrals 1-10 and their protonated forms. Conformational and vibrational analyses have been done with the aim of obtaining a theoretical estimate of the protonation entropies.     

A novel mechanism of proton transfer in protonated peptides

The study presents quantum-chemical calculations on proton transfer in protonated N-acetylglycyl-N1-methylglycinamide (AGA) as a short oligopeptide model. All calculations employ the B3LYP functional and the 6-31++G** basis set. Two different mechanisms of proton transfer are discussed. The rate-determining step of the first mechanism exhibits an energy barrier of about 17.7 kcal mol-1, and it is represented by an isomerization of the proton around the double bond of the carbonyl group. The second mechanism is based on the large conformational flexibility of AGA, where all carbonyl oxygens cooperate. The rate-determining step of this mechanism exhibits an energy barrier of only 8.3 kcal mol-1.     

Conformational study and enantioselective, regiospecific syntheses of novel aminoxy trans-proline analogues derived from an acylnitroso Diels-Alder cycloaddition

The cis/trans isomerization of the proline amide bond has many implications in biological processes. The conformations of representative acylnitroso-derived proline analogues derived from cyclopentadiene were shown to exist exclusively as the E or trans conformation in CD2Cl2. The energetically favored conformations were determined using COSMO self-consistent reaction field calculations at the B3LYP/6-31G level of theory in addition to low temperature 1H NMR studies. The syntheses of the acylnitroso-derived peptides utilized two methods to selectively functionalize either of two chemically similar esters in the acylnitroso-derived amino acids. A novel transpeptidation of the amino acid that controlled the absolute stereochemistry in the acylnitroso Diels-Alder cycloaddition took advantage of an activated aminoxy amide linkage to control regiochemistry. Alternatively, an enantioselective and regiospecific enzymatic resolution of a racemic dimethyl ester provided a novel aminoxy acid.     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

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.     

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.     

Ab initio calculations and mass spectrometric determination of the gas-phase proton affinities of 4,4'-disubstituted 2,2'-bipyridines

The gas-phase proton affinities of 4,4'-di(R)-2,2'-bipyridines (R: H, Br, Cl, NO(2), Me) were determined by mass spectrometric measurements and by ab initio calculations at the HF/6-31G and MP2/6-31G levels of theory. The energy barriers for rotation about the central C-C bond were also studied computationally. Two minima were found for both unprotonated and protonated species, the global minima being at the trans planar and cis planar conformations, respectively. Local minima for the unprotonated compounds were at the cis nonplanar conformation and for the protonated compounds at the trans nonplanar. Two different proton affinity values were calculated for each compound by employing different conformations for the protonated species. The computational values were in good agreement with the experimental proton affinities. Substituents affect the proton affinity according to their ability to withdraw or to donate electrons, halogen and nitro-substituted bipyridines having a lower proton affinity and methyl-substituted bipyridine having a higher proton affinity than 2,2'-bipyridine itself.     

Coupling of intramolecular hydrogen bonding to the cis-to-trans isomerization of a proline imide bond of small model peptides

A relationship between intramolecular hydrogen bonding and the cis-trans isomerization of a proline imide bond for proline-containing short peptides were studied by proton NMR and infrared spectroscopy using DMSO-d6/CDCl3 mixed solvents. The percentage of the trans form increases with increasing fraction of CDCl3 in the mixed solvents except for compounds without possibility of intramolecular hydrogen bonding. Chemical shift variations of amide protons with solvent mixing ratios were found to be useful for judging whether the amide protons take part in the intramolecular hydrogen bonding to a considerable degree or not. These results and infrared spectra were used to specify intramolecularly hydrogen bonded structures of the peptides. Formation of the 10-membered or 13-membered hydrogen bonded ring which includes the carbonyl group precedent to the prolyl residue facilitates the cis-to-trans isomerization and these hydrogen bonded rings are strong enough to restrict the proline imide bond to the trans form in CDCl3 solution. On the other hand, a 7-membered hydrogen bonded ring is not so effective in restricting the proline imide bond.     

Aspartic acid side chain effect—Experimental and theoretical insight

Gas-phase H/D exchange and density functional theory study of the Asp and Glu side-chain carboxylic group intrinsic reactivity is reported. H/D exchange site specific treatment and some additional theoretical calculations showed that a side-chain carboxylic group may initiate proton transfer along with bond formation to one of its oxygens, i.e., possibility to initiate selective of cleavage peptide bond ("aspartic acid effect"). That finding is used to select aspartic acid cleavage mechanisms (side-chain proton transfer either to backbone carbonyl or to amide nitrogen) for further computational study. B3LYP/6-31G(d) and G3(MP2)//B3LYP potential energy profiles of both mechanisms on a model system CH3CO-Asp-NHCH3 were constructed. Although energy employed in low-energy collision induced dissociation suffices for both mechanisms thresholds, energy transferred to specific modes suggests a complex one-step mechanism of proton transfer (from the side-chain carboxylic group to the backbone amide group), bond formation (between deprotonated carboxylic group and carbon atom of the backbone carbonyl), and peptide bond cleavage as favorable.     

4-Methylpseudoproline derived from α-methylserine - synthesis and conformational studies

This paper presents the synthesis and solution conformational studies of the tripeptides Fmoc-Ala-(R)-(αMe)Ser(Ψ(H,H)Pro)-Ala-OBu(t) (6a) and Fmoc-Ala-(S)-(αMe)Ser(Ψ(H,H)Pro)-Ala-OBu(t) (6b). Additionally, the X-ray structure of 6a is given. NMR analysis corroborated by theoretical calculations (XPLOR) shows that in both peptides the amide bond between pseudoproline and the preceding amino acid is in the trans conformation. The same amide bond geometry was observed in the crystal state of 6a. The latter is additionally influenced by the presence of two symmetrically independent molecules in an asymmetric unit. Both molecules adopt a conformation which resembles β-turn type II, stabilized by hydrogen bonding. The conformational preferences and prolyl cis-trans isomerization of Ac-(αMe)Ser(Ψ(H,H)Pro)-NHMe (7) were explored at the IEFPCM/B3LYP/6-31+G(d) level of theory in vacuum, water and chloroform. It has been shown that the trans isomer predominates in water solutions and the cis isomer is preferred in chloroform. The conformation of 7 is down-puckered independently of the geometry of the amide bonds, with lower puckering in the transition state of the cis-trans isomerization.     

Proton affinity of canavanine and canaline, oxyanalogues of arginine and ornithine, from the extended kinetic method

The absolute proton affinities of the nonprotein amino acids canavanine and canaline have been determined using the extended kinetic method in an electrospray ionization quadrupole ion trap instrument. Canavanine results from the substitution of an oxygen atom for the delta-CH2 group in the side chain of the protein amino acid arginine, whereas canaline results from a similar substitution at the delta-CH2 group in the side chain of ornithine. Absolute proton affinities of 1001+/-9 and 950+/-7 kJ/mol are obtained for canavanine and canaline, respectively. For canaline, this proton affinity is in excellent agreement with theoretical predictions obtained using the hybrid density functional theory method B3LYP/6-311++G**//B3LYP/6-31+G*. For canavanine, theory predicts a somewhat larger proton affinity of 1015 kJ/mol. Oxygen atom substitution in these nonprotein amino acids results in a decrease in their proton affinities of 40-50 kJ/mol compared to arginine and ornithine.     

Gas-phase reactivity of 2,7-dimethyl-[1,2,4]-triazepine thio derivatives toward Cu+ cation: a DFT study

The gas-phase interactions of 2,7-dimethyl-[1,2,4]-triazepine and its thio derivatives with Cu+ were studied through the use of high-level density functional theory (DFT) calculations. The structure of all possible tautomers and their conformers was optimized at the B3LYP/6-31G(d) level of theory. Final energies were obtained at the B3LYP/6-311+G(2df,2p) level. It has been found that the direct association of Cu+ occurs at the oxygen atom attached to position 3 in the case of the dioxo derivative and at the sulfur atom in all other cases. For the dithio derivatives, the global minimum of the PES corresponds to the structure in which the metal ion bridges between the heteroatom at position 3 and the nitrogen atom at position 4 of the corresponding enolic tautomer, forming a four-membered ring structure; for the dioxo derivative, this conformer competes with the ketone tautomer. Moreover, the isomerization processes leading from the most stable adduct to the other stable conformers were investigated. Among all the considered compounds, the 3,5-dithiotriazepines-Cu+ is found to be the one that associates Cu+ more tightly in the gas phase. The calculated Cu+ binding energies show a good correlation with the experimental proton affinities.