Direct dynamics study of N-protonated diglycine surface-induced dissociation. Influence of collision energy

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.

The selective adsorption performance and mechanism of multiwall magnetic carbon nanotubes for heavy metals in wastewater

The safe treatment of heavy metals in wastewater is directly related to human health and social development. In this paper, a new type of recyclable adsorbent is synthesized through the oxidation of enhancer and modification with magnetic nanoparticles. The new adsorbent not only inherits the advantages of multiwall carbon nanotubes (6O-MWCNTs), but also exhibits a new magnetic property and further improved adsorption capacity, which is conducive to the magnetic separation and recovery of heavy metals. The adsorption results indicate that multiwall magnetic carbon nanotubes (6O-MWCNTs@Fe₃O₄) have a good performance for Pb(II) selective adsorption, with a maximum adsorption capacity of 215.05 mg/g, much higher than the existing adsorption capacity of the same type of adsorbents. Under the action of an external magnetic field, 6O-MWCNTs@Fe₃O₄ that adsorbed metal ions can quickly achieve good separation from the solution. The joint characterization results of FTIR and XPS show that under the action of both coordination and electrostatic attraction, the C=O bond in the -COOH group is induced to open by the metal ions and transforms into an ionic bond, and the metal ions are stably adsorbed on the surface of 6O-MWCNTs@Fe₃O₄. Pb(II) has a stronger attraction than Cu(II) and Cd(II) to the lone pair of electrons in oxygen atoms to form complexes, due to the covalent index of Pb (6.41) is more larger than that of Cu (2.98) and Cd (2.71).These data provide a new type of recyclable adsorbent for the efficient treatment of heavy metal ions in wastewater and enrich relevant theoretical knowledge.

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.

Quantum chemical calculation of electron ionization mass spectra for general organic and inorganic molecules

We introduce a fully stand-alone version of the Quantum Chemistry Electron Ionization Mass Spectra (QCEIMS) program [S. Grimme, Angew. Chem. Int. Ed., 2013, 52, 6306] allowing efficient simulations for molecules composed of elements with atomic numbers up to Z = 86. The recently developed extended tight-binding semi-empirical method GFN-xTB has been combined with QCEIMS, thereby eliminating dependencies on third-party electronic structure software. Furthermore, for reasonable calculations of ionization potentials, as required by the method, a second tight-binding variant, IPEA-xTB, is introduced here. This novel combination of methods allows the automatic, fast and reasonably accurate computation of electron ionization mass spectra for structurally different molecules across the periodic table. In order to validate and inspect the transferability of the method, we perform large-scale simulations for some representative organic, organometallic, and main-group inorganic systems. Theoretical spectra for 23 molecules are compared directly to experimental data taken from standard databases. For the first time, realistic quantum chemistry based EI-MS for organometallic systems like ferrocene or copper(ii)acetylacetonate are presented. Compared to previously used semiempirical methods, GFN-xTB is faster, more robust, and yields overall higher quality spectra. The partially analysed theoretical reaction and fragmentation mechanisms are chemically reasonable and reveal in unprecedented detail the extreme complexity of high energy gas phase ion chemistry including complicated rearrangement reactions prior to dissociation.

Classical Dynamics Simulations of Dissociation of Protonated Tryptophan in the Gas Phase

Gas phase decomposition of protonated amino acids are of great interest due to their role in understanding protein and peptide chemistry. Several experimental and theoretical studies have been reported in the literature on this subject. In the present work, decomposition of the aromatic amino acid protonated tryptophan was studied by on-the-fly classical chemical dynamics simulations using density functional theory. Mass spectrometry and electronic structure theory studies have shown multiple dissociation pathways for this biologically relevant molecule. Unlike aliphatic amino acids, protonated tryptophan dissociates via NH₃ elimination rather than the usual iminium ion formation by combined removal of H₂O and CO molecules. Also, a major fragmentation pathway in the present work involves C_α-C_β bond fission. Results of the chemical dynamics simulations reported here are in overall agreement with experiments, and detailed atomic level mechanisms are presented.

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

In this Perspective, mass spectrometry experiments and chemical dynamics simulations are described that have explored the atomistic dynamics of protonated peptide ions, peptide-H(+), colliding with organic surfaces. These studies have investigated the energy transfer and fragmentation dynamics for peptide-H(+) surface-induced dissociation (SID), peptide-H(+) physisorption on the surface, soft landing (SL), and peptide-H(+) reaction with the surface, reactive landing (RL). SID provides primary structures of biological ions and information regarding their fragmentation pathways and energetics. Two SID mechanisms are found for peptide-H(+) fragmentation. A traditional mechanism in which peptide-H(+) is vibrationally excited by its collision with the surface, rebounds off the surface and then dissociates in accord with the statistical, RRKM unimolecular rate theory. The other, shattering, is a nonstatistical mechanism in which peptide-H(+) fragments as it collides with the surface, dissociating via many pathways and forming many product ions. Shattering is important for collisions with diamond and perfluorinated self-assembled monolayer (F-SAM) surfaces, increasing in importance with the peptide-H(+) collision energy. Chemical dynamics simulations also provide important mechanistic insights on SL and RL of biological ions on surfaces. The simulations indicate that SL occurs via multiple mechanisms consisting of sequences of peptide-H(+) physisorption on and penetration in the surface. SL and RL have a broad range of important applications including preparation of protein or peptide microarrays, development of biocompatible substrates and biosensors, and preparation of novel synthetic materials, including nanomaterials. An important RL mechanism is intact deposition of peptide-H(+) on the surface.

Chemical dynamics simulations of energy transfer, surface-induced dissociation, soft-landing, and reactive-landing in collisions of protonated peptide ions with organic surfaces

Different simulation approaches like MM, QM + MM, and QM/MM, were used to study surface-induced dissociation, soft-landing, and reactive-landing for the peptide-H⁺ + surface collisions.

Dynamics of energy transfer and soft-landing in collisions of protonated dialanine with perfluorinated self-assembled monolayer surfaces

Chemical dynamics simulations are reported which provide atomistic details of collisions of protonated dialanine, ala2-H(+), with a perfluorinated octanethiolate self-assembled monolayer (F-SAM) surface. The simulations are performed at collision energies Ei of 5.0, 13.5, 22.5, 30.00, and 70 eV, and incident angles 0° (normal) and 45° (grazing). Excellent agreement with experiment (J. Am. Chem. Soc., 2000, 122, 9703-9714) is found for both the average fraction and distribution of the collision energy transferred to the ala2-H(+) internal degrees of freedom. The dominant pathway for this energy transfer is to ala2-H(+) vibration, but for Ei = 5.0 eV ∼20% of the energy transfer is to ala2-H(+) rotation. Energy transfer to ala2-H(+) rotation decreases with increase in Ei and becomes negligible at high Ei. Three types of collisions are observed in the simulations: i.e. those for which ala2-H(+) (1) directly scatters off the F-SAM surface; (2) sticks/physisorbs on/in the surface, but desorbs within the 10 ps numerical integration of the simulations; and (3) remains trapped (i.e. soft-landed) on/in the surface when the simulations are terminated. Penetration of the F-SAM by ala2-H(+) is important for the latter two types of events. The trapped trajectories are expected to have relatively long residence times on the surface, since a previous molecular dynamics simulation (J. Phys. Chem. B, 2014, 118, 5577-5588) shows that thermally accommodated ala2-H(+) ions have an binding energy with the F-SAM surface of at least ∼15 kcal mol(-1).

Intermolecular potential for binding of protonated peptide ions with perfluorinated hydrocarbon surfaces

An analytic potential energy function was developed to model both short-range and long-range interactions between protonated peptide ions and perfluorinated hydrocarbon chains. The potential function is defined as a sum of two-body potentials of the Buckingham form. The parameters of the two-body potentials were obtained by fits to intermolecular potential energy curves (IPECs) calculated for CF4, which represents the F and C atoms of a perfluoroalkane chain, interacting with small molecules chosen as representatives of the main functional groups and atoms present in protonated peptide ions: specifically, CH4, NH3, NH4(+), and HCOOH. The IPECs were calculated at the MP2/aug-cc-pVTZ level of theory, with basis set superposition error (BSSE) corrections. Good fits were obtained for an energy range extending up to about 400 kcal/mol. It is shown that the pair potentials derived from the NH3/CF4 and HCOOH/CF4 fits reproduce acceptably well the intermolecular interactions in HCONH2/CF4, which indicates that the parameters obtained for the amine and carbonyl atoms may be transferable to the corresponding atoms of the amide group. The derived potential energy function may be used in chemical dynamics simulations of collisions of peptide-H(+) ions with perfluorinated hydrocarbon surfaces.

Systematic fragmentation method and the effective fragment potential: an efficient method for capturing molecular energies

The systematic fragmentation method fragments a large molecular system into smaller pieces, in such a way as to greatly reduce the computational cost while retaining nearly the accuracy of the parent ab initio electronic structure method. In order to attain the desired (sub-kcal/mol) accuracy, one must properly account for the nonbonded interactions between the separated fragments. Since, for a large molecular species, there can be a great many fragments and therefore a great many nonbonded interactions, computations of the nonbonded interactions can be very time-consuming. The present work explores the efficacy of employing the effective fragment potential (EFP) method to obtain the nonbonded interactions since the EFP method has been shown previously to capture nonbonded interactions with an accuracy that is often comparable to that of second-order perturbation theory. It is demonstrated that for nonbonded interactions that are not high on the repulsive wall (generally >2.7 A), the EFP method appears to be a viable approach for evaluating the nonbonded interactions. The efficacy of the EFP method for this purpose is illustrated by comparing the method to ab initio methods for small water clusters, the ZOVGAS molecule, retinal, and the alpha-helix. Using SFM with EFP for nonbonded interactions yields an error of 0.2 kcal/mol for the retinal cis-trans isomerization and a mean error of 1.0 kcal/mol for the isomerization energies of five small (120-170 atoms) alpha-helices.

Importance of shattering fragmentation in the surface-induced dissociation of protonated octaglycine

A QM + MM direct chemical dynamics simulation was performed to study collisions of protonated octaglycine, gly(8)-H(+), with the diamond {111} surface at an initial collision energy E(i) of 100 eV and incident angle theta(i) of 0 degrees and 45 degrees. The semiempirical model AM1 was used for the gly(8)-H(+) intramolecular potential, so that its fragmentation could be studied. Shattering dominates gly(8)-H(+) fragmentation at theta(i) = 0 degrees, with 78% of the ions dissociating in this way. At theta(i) = 45 degrees shattering is much less important. For theta(i) = 0 degrees there are 304 different pathways, many related by their backbone cleavage patterns. For the theta(i) = 0 degrees fragmentations, 59% resulted from both a-x and b-y cleavages, while for theta(i) = 45 degrees 70% of the fragmentations occurred with only a-x cleavage. For theta(i) = 0 degrees, the average percentage energy transfers to the internal degrees of freedom of the ion and the surface, and the energy remaining in ion translation are 45%, 26%, and 29%. For 45 degrees these percentages are 26%, 12%, and 62%. The percentage energy-transfer to DeltaE(int) for theta(i) = 0 degrees is larger than that reported in previous experiments for collisions of des-Arg(1)-bradykinin with a diamond surface at the same theta(i). This difference is discussed in terms of differences between the model diamond surface used in the simulations and the diamond surface prepared for the experiments.

Surface-induced dissociation of small molecules, peptides, and non-covalent protein complexes

This article provides a perspective on collisions of ions with surfaces, including surface-induced dissociation (SID) and reactive ion scattering spectrometry (RISS). The content is organized into sections on surface-induced dissociation of small ions, surface characterization of organic thin films by collision of well-characterized ions into surfaces, the use of SID to probe peptide fragmentation, and the dissociation of large non-covalent complexes by SID. Examples are given from the literature with a focus on experiments from the authors' laboratory. The article is not a comprehensive review but is designed to provide the reader with an overview of the types of results possible by collisions of ions into surfaces.

Classical Trajectory Simulations of Photoionization Dynamics of Tryptophan: Intramolecular Energy Flow, Hydrogen-Transfer Processes and Conformational Transitions

One-photon and two-photon ionization dynamics of tryptophan is studied by classical trajectory simulations using the semiempirical parametric method number 3 (PM3) potential surface in "on the fly" calculations. The tryptophan conformer is assumed to be in the vibrational ground state prior to ionization. Initial conditions for the trajectories are weighted according to the Wigner distribution function computed for that state. Vertical ionization in the spirit of the classical Franck-Condon principle is assumed. For the two-photon ionization process the ionization is assumed to go resonantively through the first excited state. Most trajectories are computed, and the analysis is carried out for the first 10 ps. A range of interesting effects are observed. The main findings are as follows: (1) Multiple conformational transitions are observed in most of the trajectories within the ultrafast duration of 10 ps. (2) Hydrogen transfer from the carboxyl group to the amino group and back has been observed. A zwitterion is formed as a transient state. (3) Two new isomers are formed during the dynamics, which have apparently not been previously observed. (4) Fast energy flow between the ring modes and the amino acid backbone is observed for both one- and two-photon ionization. However, the effective vibrational temperatures only approach the same value after 90 ps. The conformation transition dynamics, the proton-transfer processes and the vibrational energy flow are discussed and analyzed.

Photoionization dynamics of glycine adsorbed on a silicon cluster: “On-the-fly” simulations

Dynamics of glycine chemisorbed on the surface of a silicon cluster is studied for a process that involves single-photon ionization, followed by recombination with the electron after a selected time delay. The process is studied by "on-the-fly" molecular dynamics simulations, using the semiempirical parametric method number 3 (PM3) potential energy surface. The system is taken to be in the ground state prior to photoionization, and time delays from 5 to 50 fs before the recombination are considered. The time evolution is computed over 10 ps. The main findings are (1) the positive charge after ionization is initially mostly distributed on the silicon cluster. (2) After ionization the major structural changes are on the silicon cluster. These include Si-Si bond breaking and formation and hydrogen transfer between different silicon atoms. (3) The transient ionization event gives rise to dynamical behavior that depends sensitively on the ion state lifetime. Subsequent to 45 fs evolution in the charged state, the glycine molecule starts to rotate on the silicon cluster. Implications of the results to various processes that are induced by transient transition to a charged state are discussed. These include inelastic tunneling in molecular devices, photochemistry on conducting surfaces, and electron-molecule scattering.

Intermolecular Potential To Represent Collisions of Protonated Peptide Ions with Fluorinated Alkane Surfaces

The MP2/6-311++G(2df,2pd) level of theory was used to calculate intermolecular potential curves between CF(4), as a model for the C and F atoms of a fluorinated alkane surface, and CH(4), NH(3), NH(4)(+), H(2)CO, and H(2)O as models for different types of atoms and functional groups comprising protonated peptide ions. This level of theory was tested by comparisons with the MP2/aug-cc-pVTZ and CCSD(T)/aug-cc-pVTZ theories. Explicit-atom (EA) analytic potential energy functions were then derived by fitting these potential energy curves with two-body potentials between the atoms of the two interacting molecules. An intermolecular potential for the interaction of a protonated peptide ion with a fluorinated alkane surface may be constructed from these two-body potentials. Intermolecular potentials, for which CF(4) is treated as a united atom (UA), were developed by isotropically averaging the CF(4) orientation for each of the EA potential energy curves. The intermolecular potential energy curves calculated for CF(4) are compared with curves calculated previously for CH(4) interacting with the same molecules, to consider the relative efficiency of energy transfer for protonated peptide ion collisions with hydrogenated and fluorinated alkane surfaces.