Studying the Dissociation Behaviour of Ionized Non-covalent Complexes with a Cohesive Energetic and Structure Approach
Creators: name: Beneteau Renaud, Justin
Publisher: name: Université d'Ottawa / University of Ottawa
Year: 2014
Subjects: Mass Spectrometry
Non-covalent complex
Gas phase ions
RRKM unimolecular rate theory
Ion mobility spectrometry
Modelling dissociation energetics

description: This research explores the links between the structure and dissociation energetics of ionized non-covalent complexes. In chapter 3, a large series of similar non-covalent complexes were probed using electrospray tandem mass spectrometry (ESI-MS/MS) and RRKM modelling in order to identify any trends in the dissociation energetics based on charge state, overall size of the complex, or size of the substrate. Ion mobility spectrometry (IMS) in conjunction with molecular mechanics/molecular dynamics (MM/MD) was used to study the conformations of these non-covalent complexes in order to determine if the same trends identified in the energetics could be corroborated independently based on structure. The system of study consisted of varying lengths of the synthetic polymer, polymethylmethacrylate (PMMA) complexed with singly or doubly protonated diaminoalkanes (DAA) of varying length. The critical energies of dissociation (E0) increased as the length of the polymer increased and was not significantly affected by the length of the singly protonated DAA substrates. The E0 of dissociation of doubly protonated complexes was strongly influenced by the length of the DAA; longer DAA substrates had greater separation of charge which decreased coulombic repulsion within the complex resulting in higher E0 values. MM/MD low energy structures of all complexes were validated with experimental IMS measurements and showed that the arrangement between the polymer and DAA were similar for different singly protonated DAAs. When doubly protonated, the length of DAA was the most important factor in determining the overall structure of the complex. In chapter 4, a direct link is shown between the observed E0 dissociation energies and the molecular conformations for eight different peptide–saccharide complexes containing either a tri-saccharide (d-(+)-raffinose and d-panose) or tetra-saccharide (stachyose and maltotetraose) with a small peptide (FLEEL and FLEEV). The E0 values were highly related to the overall conformation adopted by the non-covalent complex in the gas phase. Complexes containing peptide FLEE(L/V) with the tri-saccharide raffinose or panose had similar E0 of dissociation (∼0.64 eV) and similar conformations based on MM/MD simulations and IMS drift times. Conversely, for complexes containing a FLEE(L/V) peptide with one of the isomeric tetra-saccharides; stachyose had a E0 ∼0.08 eV greater than maltotetraose. This difference of intermolecular interaction was also reflected by the IMS drift times; maltotetraose in complex with FLEEV or FLEEL had a 5.9% and 2.3% faster IMS drift time than stachyose respectively. This indicated that the molecular arrangement between maltotetraose and the peptides was more compact than the stachyose-peptide complexes. In chapter 5, RRKM modelling of breakdown diagrams is not possible when the reactant ion signal is overlapped by other isobaric species. Trimeric, non-covalent complexes that contained two PMMA molecules and a doubly protonated DAA, [(PMMAa)(DAA+2H)(PMMAb)]+2, have m/z signals that contain multiple different complexes having the same total number of polymer repeat units but differ in the length of the each polymer. In this situation, the applicability of using the simple kinetic method to gain insight into relative binding energies was explored. The major factors which determined the suitability of the kinetic method for this system were identified as the structural arrangement of the reactant ion complex, possible reverse activation barriers, and the evaluations of Δ(ΔS‡). MM/MD simulations coupled with IMS suggests that within the reactant ion, the DAA is almost equally shared between two PMMA oligomers and that the two PMMA oligomers interact predominately with the DAA, and not with each other. MS/MS of the trimeric reactant complexes proceeds by neutral loss of one polymer and is suggested to proceed with little or no reverse activation barrier based on the low coulombic repulsion factors. The IMS drift times of [(PMMAa)(DAA+2H)]+2 complexes that were generated directly by ESI-MS or by dissociation of a trimeric, [(PMMAa)(DAA+2H)(PMMAb)]+2 complex were found to be identical. This provides some evidence that Δ(ΔS‡) ≈ Δ(ΔS) and using a statistical mechanics approach, Δ(ΔS) ≈ 0. The effective temperature (Teff) variable in the kinetic method expression was found to decrease as a function of the size of the trimeric complex, suggesting that the population distribution of the dissociating ensemble of complexes narrows as size increases. Overall, when RRKM fitting is not possible, the simple kinetic method could provide relative energetic ranking of competing dissociations reactions however the Teff term contributed to the greatest uncertainty in obtaining absolute quantities. Fitting MS/MS breakdown diagrams of non-covalent complexes with multiple dissociation channels is difficult due to the number of total fitting variables. Building from the simple kinetic method, chapter 6 shows that the relationship between the natural logarithm of competing fragment ions and reciprocal collision energy yields a branching relationship that allows for the sign of Δ(ΔS‡) and Δ(E0) between the channels to be obtained. Furthermore, the relationships between the fitting variables of RRKM modelling are empirically related to the theoretical branching relationship characteristics. This allowed for the fitting variables of all dissociation channels to be expressed as a function of a single channel so that the theoretical branching relationship matches the experimental branching relationship. Using this method, RRKM fitting of a MS/MS breakdown diagram for APCI ionized anthracene determined the E0 and ∆S‡ was 4.69 ± 0.29 eV and -3 ± 17 J K-1; 4.21 ±0.29 eV and -19 ±15 J K-1; and 4.81 ± 0.29 eV and 36 ±22 J K-1 for hydrogen loss, acetylene loss and diacetylene loss respectively. With one exception, these values are within experimental error of the iPEPICO derived energetic values. In chapter 7, MS/MS of ammoniated triacylglycerides at multiple collision energies and computational analysis are used to explain the cause of uneven dissociation rates of the FAs from different positions on the glycerol backbone. The loss of sn-1 and sn-3 FAs are found to have lower activation energies than the loss of the sn-2 position FA, however the loss of the FA from the sn-2 position is more entropically favourable. Theoretical MS/MS breakdown curves were fit to experimental values using RRKM theory to estimate the E0 of dissociation of FAs from the three glycerol positions. The E0 values for cleavage from the sn-1 and sn-3 positions were found to be approximately 1.52 eV, while that for the sn-2 position was highly dependent on the identity of the FA at that position. Computational structures and energy analysis suggest that an important step in the dissociation of [TAG+NH4]+ is the loss of ammonia. In a model system, glyceryl tributyrate, the loss of NH3 produced two distinct [TAG+H]+ product structures sitting 148 kJ and 160 kJ in energy above the ammoniated structure. The [TAG+H]+ structure that leads to the loss of the sn-1(3) is 12 kJ lower than the [TAG+H]+ structure that leads to the loss of the sn-2 FA. From this, the loss of a neutral FA that follows sits only an additional 35–48 kJ above the [TAG+H]+ structures. In Chapter 8, singly deprotonated β-cyclodextrin monomers, [(β-CD-H+]-1, and doubly deprotonated dimers, [(β-CD)2-2H+]-2, are both present following ESI-MS and have the same monoisotopic m/z. Similar to chapter 5, this makes it difficult to generate an MS/MS breakdown diagrams that can be modelled with RRKM theory. IMS was used to mobility separate [(β-CD-H+]-1 and [(β-CD)2-2H+]-2 and was followed by MS/MS of the [(β-CycD)2-2H+]-2 ion. A second problem when generating a MS/MS breakdown diagram of non-covalent complexes that contain identical components is that the fragment ions could have an identical monoisotopic m/z as the reactant ion. MS/MS of [(β-CycD)2-2H+]-2 results in two [(β-CD-H+]-1 fragments. To overcome this, breakdown diagrams were then generated by monitoring the changes in the isotopic profile. The RRKM derived E0 for dissociation of [(β-CycD)2-H+]-1 and [(β-CycD)2-2H+]-2 were 1.85 ± 0.11eV and 1.79 ± 0.09eV, respectively, corresponding to a slight decrease in complex stability due to increased charge-charge repulsion in the dianion.