Model Simulations of the Thermal Dissociation of the TIK(H + ) 2 Tripeptide: Mechanisms and Kinetic Parameters.

Direct dynamics simulations, utilizing the RM1 semiempirical electronic structure theory, were performed to study the thermal dissociation of the doubly protonated tripeptide threonine-isoleucine-lysine ion, TIK(H ⁺ ) ₂ , for temperatures of 1250-2500 K, corresponding to classical energies of 1778-3556 kJ/mol. The number of different fragmentation pathways increases with increase in temperature. At 1250 K there are only three fragmentation pathways, with one contributing 85% of the fragmentation. In contrast, at 2500 K, there are 61 pathways, and not one dominates. The same ion is often formed via different pathways, and at 2500 K there are only 14 m/z values for the product ions. The backbone and side-chain fragmentations occur by concerted reactions, with simultaneous proton transfer and bond rupture, and also by homolytic bond ruptures without proton transfer. For each temperature the TIK(H ⁺ ) ₂ fragmentation probability versus time is exponential, in accord with the Rice-Ramsperger-Kassel-Marcus and transition state theories. Rate constants versus temperature were determined for two proton transfer and two bond rupture pathways. From Arrhenius plots activation energies E _a and A-factors were determined for these pathways. They are 62-78 kJ/mol and (2-3) × 10 ¹² s ^-1 for the proton transfer pathways and 153-168 kJ/mol and (2-4) × 10 ¹⁴ s ^-1 for the bond rupture pathways. For the bond rupture pathways, the product cation radicals undergo significant structural changes during the bond rupture as a result of hydrogen bonding, which lowers their entropies and also their E _a and A parameters relative to those for C-C bond rupture pathways in hydrocarbon molecules. The E _a values determined from the simulation Arrhenius plots are in very good agreement with the reaction barriers for the RM1 method used in the simulations. A preliminary simulation of TIK(H ⁺ ) ₂ collision-induced dissociation (CID), at a collision energy of 13 eV (1255 kJ/mol), was also performed to compare with the thermal dissociation simulations. Though the energy transferred to TIK(H ⁺ ) ₂ in the collisions is substantially less than the energy for the thermal excitations, there is substantial fragmentation as a result of the localized, nonrandom excitation by the collisions. CID results in different fragmentation pathways with a significant amount of short time nonstatistical fragmentation. Backbone fragmentation is less important, and side-chain fragmentation is more important for the CID simulations as compared to the thermal simulations. The thermal simulations provide information regarding the long-time statistical fragmentation.