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

• Published in: The journal of physical chemistry. A, Molecules, spectroscopy, kinetics, environment, & general theory Molecules, spectroscopy, kinetics, environment, & general theory, 2016-10, Vol.120 (42), p.8211-8227
• Main Authors: Homayoon, Zahra, Pratihar, Subha, Dratz, Edward, Snider, Ross, Spezia, Riccardo, Barnes, George L, Macaluso, Veronica, Martin Somer, Ana, Hase, William L
• Format: Article
• Language: English
• Subjects: Chemical Physics; Physics
• ISSN: 1089-5639; 1520-5215
• DOI: 10.1021/acs.jpca.6b05884

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+)2, 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+)2 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) × 1012 s–1 for the proton transfer pathways and 153–168 kJ/mol and (2–4) × 1014 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+)2 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+)2 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 simulation