Dynamics of Protonated Peptide Ion Collisions with Organic Surfaces: Consonance of Simulation and Experiment

Landing Proteins on Graphene Trampoline Preserves Their Gas-Phase Folding on the Surface

Molecule-surface collisions are known to initiate dynamics that lead to products inaccessible by thermal chemistry. These collision dynamics, however, have mostly been examined on bulk surfaces, leaving vast opportunities unexplored for molecular collisions on nanostructures, especially on those that exhibit mechanical properties radically different from those of their bulk counterparts. Probing energy-dependent dynamics on nanostructures, particularly for large molecules, has been challenging due to their fast time scales and high structural complexity. Here, by examining the dynamics of a protein impinging on a freestanding, single-atom-thick membrane, we discover molecule-on-trampoline dynamics that disperse the collision impact away from the incident protein within a few picoseconds. As a result, our experiments and ab initio calculations show that cytochrome c retains its gas-phase folded structure when it collides onto freestanding single-layer graphene at low energies (∼20 meV/atom). The molecule-on-trampoline dynamics, expected to be operative on many freestanding atomic membranes, enable reliable means to transfer gas-phase macromolecular structures onto freestanding surfaces for their single-molecule imaging, complementing many bioanalytical techniques.

Surface-induced Dissociation Mass Spectrometry as a Structural Biology Tool

Native mass spectrometry (nMS) is evolving into a workhorse for structural biology. The plethora of online and offline preparation, separation, and purification methods as well as numerous ionization techniques combined with powerful new hybrid ion mobility and mass spectrometry systems has illustrated the great potential of nMS for structural biology. Fundamental to the progression of nMS has been the development of novel activation methods for dissociating proteins and protein complexes to deduce primary, secondary, tertiary, and quaternary structure through the combined use of multiple MS/MS technologies. This review highlights the key features and advantages of surface collisions (surface-induced dissociation, SID) for probing the connectivity of subunits within protein and nucleoprotein complexes and, in particular, for solving protein structure in conjunction with complementary techniques such as cryo-EM and computational modeling. Several case studies highlight the significant role SID, and more generally nMS, will play in structural elucidation of biological assemblies in the future as the technology becomes more widely adopted. Cases are presented where SID agrees with solved crystal or cryoEM structures or provides connectivity maps that are otherwise inaccessible by "gold standard" structural biology techniques.

Exploring the Effects of Methylation on the CID of Protonated Lysine: A Combined Experimental and Computational Approach

We report the results of experiments, simulations, and DFT calculations that focus on describing the reaction dynamics observed within the collision-induced dissociation of l-lysine-H⁺ and its side-chain methylated analogues, N_ε-methyl-l-lysine-H⁺ (Me₁-lysine-H⁺), N_ε,N_ε-dimethyl-l-lysine-H⁺ (Me₂-lysine-H⁺), and N_ε,N_ε,N_ε-trimethyl-l-lysine-H⁺ (Me₃-lysine-H⁺). The major pathways observed in the experimental measurements were m/z 130 and 84, with the former dominant at low collision energies and the latter at intermediate to high collision energies. The m/z 130 peak corresponds to loss of N(CH₃)_nH_3-n, while m/z 84 has the additional loss of H₂CO₂ likely in the form of H₂O + CO. Within the time frame of the direct dynamics simulations, m/z 130 and 101 were the most populous peaks, with the latter identified as an intermediate to m/z 84. The simulations allowed for the determination of several reaction pathways that result in these products. A graph theory analysis enabled the elucidation of the significant structures that compose each peak. Methylation results in the preferential loss of the side-chain amide group and a reduction of cyclic structures within the m/z 84 peak population in simulations.

Determination of kinetic properties in unimolecular dissociation of complex systems from graph theory based analysis of an ensemble of reactive trajectories

In this paper, we report how graph theory can be used to analyze an ensemble of independent molecular trajectories, which can react during the simulation time-length, and obtain structural and kinetic information. This method is totally general and here is applied to the prototypical case of gas phase fragmentation of protonated cyclo-di-glycine. This methodology allows us to analyze the whole set of trajectories in an automatic computer-based way without the need of visual inspection but by getting all the needed information. In particular, we not only determine the appearance of different products and intermediates but also characterize the corresponding kinetics. The use of colored graph and canonical labeling allows for the correct characterization of the chemical species involved. In the present case, the simulations consist of an ensemble of unimolecular fragmentation trajectories at constant energy such that from the rate constants at different energies, the threshold energy can also be obtained for both global and specific pathways. This approach allows for the characterization of ion-molecule complexes, likely through a roaming mechanism, by properly taking into account the elusive nature of such species. Finally, it is possible to directly obtain the theoretical mass spectrum of the fragmenting species if the reacting system is an ion as in the specific example.

Coupling 193 nm Ultraviolet Photodissociation and Ion Mobility for Sequence Characterization of Conformationally-Selected Peptides

Ultraviolet photodissociation (UVPD) has emerged as a useful technique for characterizing peptide, protein, and protein complex primary and secondary structure. 193 nm UVPD, specifically, enables extensive covalent fragmentation of the peptide backbone without the requirement of a specific side chain chromophore and with no precursor charge state dependence. We have modified a commercial quadrupole-ion mobility-time-of-flight (Q-IM-TOF) mass spectrometer to include 193 nm UVPD following ion mobility. Ion mobility (IM) is a gas-phase separation technique that enables separation of ions by their size, shape, and charge, providing an orthogonal dimension of separation to mass analysis. Following instrument modifications, we characterized the performance of, and information that could be generated from, this new setup using the model peptides substance P, melittin, and insulin chain B. These experiments show extensive fragmentation across the peptide backbone and a variety of ion types as expected from 193 nm UVPD. Additionally, y-2 ions (along with complementary a+2 and b+2 ions) N-terminal to proline were observed. Combining the IM separation and mobility gating capabilities with UVPD, we demonstrate the ability to accomplish both mass- and mobility-selection of bradykinin des-Arg9 and des-Arg1 peptides followed by complete sequence characterization by UVPD. The new capabilities of this modified instrument demonstrate the utility of combining IM with UVPD because isobaric species cannot be independently selected with a traditional quadrupole alone.

Modeling the Effects of O-Sulfonation on the CID of Serine

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of O-sulfonation on the collision-induced dissociation for serine. Toward this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m/z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m/z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The m/z 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also observed that the RM1 semiempirical method was not able to obtain all of the major peaks seen in experimens for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.

Nonstatistical Reaction Dynamics

Nonstatistical dynamics is important for many chemical reactions. The Rice-Ramsperger-Kassel-Marcus (RRKM) theory of unimolecular kinetics assumes a reactant molecule maintains a statistical microcanonical ensemble of vibrational states during its dissociation so that its unimolecular dynamics are time independent. Such dynamics results when the reactant's atomic motion is chaotic or irregular. Intrinsic non-RRKM dynamics occurs when part of the reactant's phase space consists of quasiperiodic/regular motion and a bottleneck exists, so that the unimolecular rate constant is time dependent. Nonrandom excitation of a molecule may result in short-time apparent non-RRKM dynamics. For rotational activation, the 2J + 1 K levels for a particular J may be highly mixed, making K an active degree of freedom, or K may be a good quantum number and an adiabatic degree of freedom. Nonstatistical dynamics is often important for bimolecular reactions and their intermediates and for product-energy partitioning of bimolecular and unimolecular reactions. Post–transition state dynamics is often highly complex and nonstatistical.

A Tilted Surface and Ion Carpet Array for SID

The development of native mass spectrometry (MS) has provided structural biologists an additional tool to probe the structures of large macromolecular systems. Surface-induced dissociation (SID) is one activation method used within tandem MS experiments that has proven useful in interrogating the connectivity and topology of biologically-relevant protein complexes. We present here the use of a tilted surface and ion carpet array within a new SID device design, enabling decreased dimensions along the ion path and fewer lenses to tune. This device works well in fragmenting ions of both low (peptides) and high (protein complexes) m/z. Results show that the ion carpet array, while enabling simplification of the back-end of the device, has deficiencies in product collection and subsequently signal at higher SID energies when fragmenting protein complexes. However, the use of the tilted surface is advantageous as an effective way to shorten the device and reduce the number of independent voltages.

Role of Chemical Dynamics Simulations in Mass Spectrometry Studies of Collision-Induced Dissociation and Collisions of Biological Ions with Organic Surfaces

In this article, a perspective is given of chemical dynamics simulations of collisions of biological ions with surfaces and of collision-induced dissociation (CID) of ions. The simulations provide an atomic-level understanding of the collisions and, overall, are in quite good agreement with experiment. An integral component of ion/surface collisions is energy transfer to the internal degrees of freedom of both the ion and the surface. The simulations reveal how this energy transfer depends on the collision energy, incident angle, biological ion, and surface. With energy transfer to the ion's vibration fragmentation may occur, i.e. surface-induced dissociation (SID), and the simulations discovered a new fragmentation mechanism, called shattering, for which the ion fragments as it collides with the surface. The simulations also provide insight into the atomistic dynamics of soft-landing and reactive-landing of ions on surfaces. The CID simulations compared activation by multiple "soft" collisions, resulting in random excitation, versus high energy single collisions and nonrandom excitation. These two activation methods may result in different fragment ions. Simulations provide fragmentation products in agreement with experiments and, hence, can provide additional information regarding the reaction mechanisms taking place in experiment. Such studies paved the way on using simulations as an independent and predictive tool in increasing fundamental understanding of CID and related processes.

Collisional dynamics simulations revealing fragmentation properties of Zn(ii)-bound poly-peptide

Collisional simulations show how peptide fragmentation is modified by the presence of Zn(ii).

A Trajectory-Based Method to Explore Reaction Mechanisms

The tsscds method, recently developed in our group, discovers chemical reaction mechanisms with minimal human intervention. It employs accelerated molecular dynamics, spectral graph theory, statistical rate theory and stochastic simulations to uncover chemical reaction paths and to solve the kinetics at the experimental conditions. In the present review, its application to solve mechanistic/kinetics problems in different research areas will be presented. Examples will be given of reactions involved in photodissociation dynamics, mass spectrometry, combustion chemistry and organometallic catalysis. Some planned improvements will also be described.

Mass Spectrometry for Synthesis and Analysis

Mass spectrometry is the science and technology of ions. As such, it is concerned with generating ions, measuring their properties, following their reactions, isolating them, and using them to build and transform materials. Instrumentation is an essential element of these activities, and analytical applications are one driving force. Work from the Aston Laboratories at Purdue University's Department of Chemistry is described here, with an emphasis on accelerated reactions of ions in solution and small-scale synthesis; ion/surface collision processes, including surface-induced dissociation (SID) and ion soft landing; and applications to tissue imaging. Our special interest in chirality and the chemistry behind the origins of life is also featured together with the exciting area of tissue diagnostics.

Chemical dynamics simulations of CID of peptide ions: comparisons between TIK(H+)2 and TLK(H+)2 fragmentation dynamics, and with thermal simulations

Gas phase unimolecular fragmentation of the two model doubly protonated tripeptides threonine-isoleucine-lysine (TIK) and threonine-leucine-lysine (TLK) is studied using chemical dynamics simulations. Attention is focused on different aspects of collision induced dissociation (CID): fragmentation pathways, energy transfer, theoretical mass spectra, fragmentation mechanisms, and the possibility of distinguishing isoleucine (I) and leucine (L). Furthermore, discussion is given regarding the differences between single collision CID activation, which results from a localized impact between the ions and a colliding molecule N₂, and previous thermal activation simulation results; Z. Homayoon, S. Pratihar, E. Dratz, R. Snider, R. Spezia, G. L. Barnes, V. Macaluso, A. Martin-Somer and W. L. Hase, J. Phys. Chem. A, 2016, 120, 8211-8227. Upon thermal activation unimolecular fragmentation is statistical and in accord with RRKM unimolecular rate theory. Simulations show that in collisional activation some non-statistical fragmentation occurs, including shattering, which is not present when the ions dissociate statistically. Products formed by non-statistical shattering mechanisms may be related to characteristic mass spectrometry peaks which distinguish the two isomers I and L.

Threshold for shattering fragmentation in collision-induced dissociation of the doubly protonated tripeptide TIK(H+)2

Determination of shattering threshold for unimolecular dissociation of a model tripeptide.

Reactive Landing of Gramicidin S and Ubiquitin Ions onto Activated Self-Assembled Monolayer Surfaces

Using mass-selected ion deposition combined with in situ infrared reflection absorption spectroscopy (IRRAS), we examined the reactive landing of gramicidin S and ubiquitin ions onto activated self-assembled monolayer (SAM) surfaces terminated with N-hydroxysuccinimidyl ester (NHS-SAM) and acyl fluoride (COF-SAM) groups. Doubly protonated gramicidin S, [GS + 2H]²⁺, and two charge states of ubiquitin, [U + 5H]⁵⁺ and [U + 13H]¹³⁺, were used as model systems, allowing us to explore the effect of the number of free amino groups and the secondary structure on the efficiency of covalent bond formation between the projectile ion and the surface. For all projectile ions, ion deposition resulted in the depletion of IRRAS bands corresponding to the terminal groups on the SAM and the appearance of several new bands not associated with the deposited species. These new bands were assigned to the C=O stretching vibrations of COOH and COO^- groups formed on the surface as a result of ion deposition. The presence of these bands was attributed to an alternative reactive landing pathway that competes with covalent bond formation. This pathway with similar yields for both gramicidin S and ubiquitin ions is analogous to the hydrolysis of the NHS ester bond in solution. The covalent bond formation efficiency increased linearly with the number of free amino groups and was found to be lower for the more compact conformation of ubiquitin compared with the fully unfolded conformation. This observation was attributed to the limited availability of amino groups on the surface of the folded conformation. Our results have provided new insights on the efficiency and mechanism of reactive landing of peptides and proteins onto activated SAMs. Graphical Abstract ᅟ.