A systematic and efficient method to estimate the vibrational frequencies of linear peptide and protein ions with any amino acid sequence for the calculation of Rice-Ramsperger-Kassel-Marcus rate constant

Enantioselective Reduction of Noncovalent Complexes of Amino Acids with CuII via Resonant Collision-Induced Dissociation: Collision Energy, Activation Duration Effects, and RRKM Modeling

Formation of noncovalent complexes is one of the approaches to perform chiral analysis with mass spectrometry. Enantiomeric distinction of amino acids (AAs) based on the relative rate constants of competitive fragmentations of quaternary copper complexes is an efficient method for chiral differentiation. Here, we studied the complex [CuII,(Phe,PhG,Pro-H)]+ (m/z 493) under resonant collision-induced dissociation conditions while varying the activation time. The precursor ion can yield two main fragments through the loss of the non-natural AA phenylglycine (PhG): the expected product ion [CuII,(Phe,Pro-H)]+ (m/z 342) and the reduced product ion [CuI,(Phe,Pro)]+ (m/z 343). Enantioselective reduction describes the difference in relative abundance of these ions, which depends on the chirality of the precursor ion: the formation of the reduced ion m/z 343 is favored in homochiral complexes (DDD) compared to heterochiral complexes (such as LDD). Energy-resolved mass spectrometry data show that reduction, which arises from rearrangement, is favored at a low collision energy (CE) and long activation time (ActT), whereas direct cleavage preferentially occurs at a high CE and short ActT. These results were confirmed with kinetic modeling based on RRKM theory. For this modeling, it was necessary to set a pre-exponential factor as a reference, so that the E0 values obtained are relative values. Interestingly, these simulations showed that the critical energy E0 required to form the reduced ion is comparable in both homochiral and heterochiral complexes. However, the formation of product ion m/z 342 through direct cleavage is associated with a lower E0 in heterochiral complexes. Consequently, enantioselectivity would not be caused by enhanced reduction in homochiral complexes but rather by direct cleavage being favored in heterochiral complexes.

Benchmarking higher energy collision dissociation (HCD) by investigation of binding energies of gas-phase host-guest complexes of hemicryptophane cages

Synthesis of host molecules that feature well-defined characteristics for molecular recognition of guest molecules is often a major aim of synthetic host-guest (H-G) chemistry. A key consideration in evaluating the selectivity of hosts and the affinities of guests is the measurement of binding energies of obtained H-G complexes. In contrast to nuclear magnetic resonance (NMR) or fluorescence measurements that are capable of measuring binding strengths in solution, mass spectrometry offers the opportunity to measure gas-phase binding energies. Presented in this article is a higher energy collision dissociation (HCD) approach for determining critical energies of dissociation of H-G complexes. Experiments were performed on electrospray ionization (ESI)-generated H-G pairs in an LTQ-XL/Orbitrap hybrid instrument. The presented HCD approach requires preliminary calibration of the internal energy distribution of generated ions that was achieved by the use of activation parameters that were known from previous low-energy collision-induced dissociation (low-energy CID) experiments. Internal energy deposition was modeled based on a truncated Maxwell-Boltzmann distribution and characteristic temperature (T_char ). Using this method, critical energies of dissociation were determined for 10 H-G biologically relevant complexes of the heteroditopic hemicryptophane cage host (Host). Obtained results are compared with those found previously by low-energy CID. The use of this HCD technique is relatively straightforward, although its implementation does require knowledge (or a presumption) about the Arrhenius pre-exponential factor of the complexes to obtain their critical energies of dissociation.

TUTORIAL: ION ACTIVATION IN TANDEM MASS SPECTROMETRY USING ULTRA-HIGH RESOLUTION INSTRUMENTATION

Tandem mass spectrometry involves isolation of specific precursor ions and their subsequent excitation through collision-, photon-, or electron-mediated activation techniques in order to induce unimolecular dissociation leading to formation of fragment ions. These powerful ion activation techniques, typically used in between mass selection and mass analysis steps for structural elucidation, have not only found a wide variety of analytical applications in chemistry and biology, but they have also been used to study the fundamental properties of ions in the gas phase. In this tutorial paper, a brief overview is presented of the theories that have been used to describe the activation of ions and their subsequent unimolecular dissociation. Acronyms of the presented techniques include CID, PQD, HCD, SORI, SID, BIRD, IRMPD, UVPD, EPD, ECD, EDD, ETD, and EID. The fundamental principles of these techniques are discussed in the context of their implementation on ultra-high resolution tandem mass spectrometers. © 2020 John Wiley & Sons Ltd. Mass Spec Rev.

Sonification based de novo protein design using artificial intelligence, structure prediction, and analysis using molecular modeling

We report the use of a deep learning model to design de novo proteins, based on the interplay of elementary building blocks via hierarchical patterns. The deep neural network model is based on translating protein sequences and structural information into a musical score that features different pitches for each of the amino acids, and variations in note length and note volume reflecting secondary structure information and information about the chain length and distinct protein molecules. We train a deep learning model whose architecture is composed of several long short-term memory units from data consisting of musical representations of proteins classified by certain features, focused here on alpha-helix rich proteins. Using the deep learning model, we then generate de novo musical scores and translate the pitch information and chain lengths into sequences of amino acids. We use a Basic Local Alignment Search Tool to compare the predicted amino acid sequences against known proteins, and estimate folded protein structures using the Optimized protein fold RecognitION method (ORION) and MODELLER. We find that the method proposed here can be used to design de novo proteins that do not exist yet, and that the designed proteins fold into specified secondary structures. We validate the newly predicted protein by molecular dynamics equilibration in explicit water and subsequent characterization using a normal mode analysis. The method provides a tool to design novel protein materials that could find useful applications as materials in biology, medicine, and engineering.

A Self-Consistent Sonification Method to Translate Amino Acid Sequences into Musical Compositions and Application in Protein Design Using Artificial Intelligence

We report a self-consistent method to translate amino acid sequences into audible sound, use the representation in the musical space to train a neural network, and then apply it to generate protein designs using artificial intelligence (AI). The sonification method proposed here uses the normal mode vibrations of the amino acid building blocks of proteins to compute an audible representation of each of the 20 natural amino acids, which is fully defined by the overlay of its respective natural vibrations. The vibrational frequencies are transposed to the audible spectrum following the musical concept of transpositional equivalence, playing or writing music in a way that makes it sound higher or lower in pitch while retaining the relationships between tones or chords played. This transposition method ensures that the relative values of the vibrational frequencies within each amino acid and among different amino acids are retained. The characteristic frequency spectrum and sound associated with each of the amino acids represents a type of musical scale that consists of 20 tones, the "amino acid scale". To create a playable instrument, each tone associated with the amino acids is assigned to a specific key on a piano roll, which allows us to map the sequence of amino acids in proteins into a musical score. To reflect higher-order structural details of proteins, the volume and duration of the notes associated with each amino acid are defined by the secondary structure of proteins, computed using DSSP and thereby introducing musical rhythm. We then train a recurrent neural network based on a large set of musical scores generated by this sonification method and use AI to generate musical compositions, capturing the innate relationships between amino acid sequence and protein structure. We then translate the de novo musical data generated by AI into protein sequences, thereby obtaining de novo protein designs that feature specific design characteristics. We illustrate the approach in several examples that reflect the sonification of protein sequences, including multihour audible representations of natural proteins and protein-based musical compositions solely generated by AI. The approach proposed here may provide an avenue for understanding sequence patterns, variations, and mutations and offers an outreach mechanism to explain the significance of protein sequences. The method may also offer insight into protein folding and understanding the context of the amino acid sequence in defining the secondary and higher-order folded structure of proteins and could hence be used to detect the effects of mutations through sound.

Investigation of Hemicryptophane Host-Guest Binding Energies Using High-Pressure Collision-Induced Dissociation in Combination with RRKM Modeling

In advancing host-guest (H-G) chemistry, considerable effort has been spent to synthesize host molecules with specific and well-defined molecular recognition characteristics including selectivity and adjustable affinity. An important step in the process is the characterization of binding strengths of the H-G complexes that is typically performed in solution using NMR or fluorescence. Here, we present a mass spectrometry-based multimodal approach to obtain critical energies of dissociation for two hemicryptophane cages with three biologically relevant guest molecules. A combination of blackbody infrared radiative dissociation (BIRD) and high-pressure collision-induced dissociation (high-pressure CID), along with RRKM modeling, was employed for this purpose. For the two tested hemicryptophane hosts, the cage containing naphthyl linkages exhibited stronger interactions than the cage bearing phenyl linkages. For both cages, the order of guest stability is choline > acetylcholine > betaine. The information obtained by these types of mass spectrometric studies can provide new insight into the structural features that most influence the stability of H-G pairs, thereby providing guidance for future syntheses. Graphical Abstract.

Formation of gas-phase peptide ions and their dissociation in MALDI: insights from kinetic and ion yield studies

Insights on mechanisms for the generation of gas-phase peptide ions and their dissociation in matrix-assisted laser desorption ionization (MALDI) gained from the kinetic and ion yield studies are presented. Even though the time-resolved photodissociation technique was initially used to determine the dissociation kinetics of peptide ions and their effective temperature, it was replaced by a simpler method utilizing dissociation yields from in-source decay (ISD) and post-source decay (PSD). The ion yields for a matrix and a peptide were measured by repeatedly irradiating a region on a sample and collecting ion signals until the sample in the region was completely depleted. Matrix- and peptide-derived gas-phase cations were found to be generated by pre-formed ion emission or by ion-pair emission followed by anion loss, but not by laser-induced ionization. The total number of ions, that is, matrix plus peptide, was found to be equal to the number of ions emitted from a pure matrix. A matrix plume was found to cool as it expanded, from around 800-1,000 K to 400-500 K. Dissociation of peptide ions along b/y channels was found to occur statistically, that is, following RRKM behavior. Small critical energy (E0  = 0.6-0.7 eV) and highly negative critical entropy (ΔS(‡)  = -30 to -25 eu) suggested that the transition structure was stabilized by multiple intramolecular interactions.

Internal energy distribution in electrospray ionization: towards the evaluation of a thermal-like distribution from the multiple-collision model

RATIONALE

The internal energy deposition in ions that cross the desolvation region of an electrospray ionization (ESI) source affects the mass spectra that are obtained using in-source collision-induced dissociation (CID) or in tandem mass spectrometry (MS/MS) mode. It is thus important to evaluate the internal energy distributions of the ions in different parts of an ESI mass spectrometer.

METHODS

The desolvation region is considered as a collision zone and a partially elastic multiple-collision model is used to account for the accumulation of internal energy in the ions. The ion survival yields (SY(Theo) of the theoretical mass spectra calculated by MassKinetics software are fitted with the experimental ion survival yields (SY(Exp)) of the substituted benzylpyridinium cations that have been obtained with an ESI source interfaced with a quadrupole mass spectrometer. The theoretical parameters used for fitting the calculation data with the experimental results are the center-of-mass collision energy (Ecom ) of the colliding ions and a term related to the pressure of the desolvation area of the ESI interface.

RESULTS

In the proposed model, an average number of 'effective' collisions of close to 30 in the desolvation area is employed. The voltages applied to the orifice of this interface are correlated to a theoretical initial kinetic energy (E(init,Kin)) in the laboratory frame of the ions. In the present case, these theoretical initial kinetic energies range from 5.5 to 9 eV. The internal energy distributions evaluated from this model resemble the thermal distributions of ions having 'characteristic temperatures' between 1020 and 1550 K, and the results of calculations show that the mean internal energy of the ions increases linearly with the orifice voltage.

CONCLUSIONS

The model used in this study can account for the energy build-up of the ions in an ESI interface and allows the change in the internal energy distribution of the electrosprayed ions in different regions of a mass spectrometer to be evaluated.

Dissociation kinetics of singly protonated leucine enkephalin investigated by time-resolved photodissociation tandem mass spectrometry

The yields of post-source decay (PSD) and time-resolved photodissociation (PD) at 193 and 266 nm were measured for singly protonated leucine enkephalin ([YGGFL + H](+)), a benchmark in the study of peptide ion dissociation, by using tandem time-of-flight mass spectrometry. The peptide ion was generated by matrix-assisted laser desorption ionization (MALDI) using 2,5-dihydroxybenzoic acid as the matrix. The critical energy (E(0)) and entropy (DeltaS(++) at 1000 K) for the dissociation were determined by Rice-Ramsperger-Kassel-Marcus fit of the experimental data. MALDI was done for a mixture of YGGFL and Y(6) and the plume temperature determined by the kinetic analysis of [Y(6) + H](+) data were used to improve the precision of E(0) and DeltaS(++) for [YGGFL + H](+). E(0) and DeltaS(++) thus determined (E(0) = 0.67 +/- 0.08 eV, DeltaS(++) = -24.4 +/- 3.2 eu with 1 eu = 4.184 J K(-1)mol(-1)) were significantly different from those determined by blackbody infrared radiative dissociation (BIRD) (E(0) = 1.10 eV, DeltaS(++) = -14.9 eu), and by surface-induced dissociation (SID) (E(0) = 1.13 eV, DeltaS(double dagger) = -10.3 eu). Analysis of the present experimental data with the SID kinetics (and BIRD kinetics also) led to an unrealistic situation where not only PSD and PD but also MALDI-TOF signals could not be detected. As an explanation for the discrepancy, it was suggested that transition-state switching occurs from an energy bottleneck (SID/BIRD) to an entropy bottleneck (PSD/PD) as the internal energy increases.

Influence of basic residues on dissociation kinetics and dynamics of singly protonated peptides: time-resolved photodissociation study

Product ion yields in postsource decay and time-resolved photodissociation at 193 and 266 nm were measured for some peptide ions with lysine ([KF6 + H]+, [F6K + H]+, and [F3KF3 + H]+) formed by matrix-assisted laser desorption ionization. The critical energy (E0) and entropy (DeltaS(double dagger)) were determined by RRKM fitting of the data. The results were similar to those found previously for peptide ions with histidine. To summarize, the presence of a basic residue, histidine or lysine, inside a peptide ion retarded its dissociation by lowering DeltaS(double dagger). On the basis of highly negative DeltaS(double dagger), presence of intramolecular interaction involving a basic group in the transition structure was proposed.

Time-resolved photodissociation study of singly protonated peptides with a histidine residue generated by matrix-assisted laser desorption ionization: dissociation rate constant and internal temperature

Product ion yields in post-source decay and time-resolved photodissociation at 193 and 266 nm were measured for some peptide ions with a histidine residue ([HF(6) + H](+), [F(6)H + H](+), and [F(3)HF(3) + H](+)) formed by matrix-assisted laser desorption ionization (MALDI). Compared with similar data for peptide ions without any basic residue reported previously, significant reduction in dissociation efficiency was observed. Internal temperatures (T) of the peptide ions and their dissociation kinetic parameters-the critical energy (E(0)) and entropy (DeltaS(double dagger))-were determined by the method reported previously. Slight decreases in E(0), DeltaS(double dagger), and T were responsible for the histidine effect-reduction in dissociation rate constant. Regardless of the presence of the residue, DeltaS(double dagger) was far more negative than previous quantum chemical results. Based on this, we propose the existence of transition structures in which the nitrogen atoms in the histidine residue or at the N-terminus coordinate to the reaction centers. Reduction in T in the presence of a histidine residue could not be explained based on popular models for ion formation in MALDI, such as the gas-phase proton transfer model.

Analysis of caged xanthones from the resin of Garcinia hanburyi using ultra-performance liquid chromatography/electrospray ionization quadrupole time-of-flight tandem mass spectrometry

On-line ultra high-performance liquid chromatography (UHPLC) coupled with electrospray quadrupole time-of-flight tandem mass spectrometry (ESI-QTOF-MS/MS/MS) has been developed for the analysis of a series of caged xanthones in the resin of Garcinia hanburyi. The fragmentation of protonated molecular ions for 12 known cadged xanthones was carried out using low-energy collision-induced electrospray ionization tandem mass spectrometry. It was found that Retro-Diels-Alder rearrangement occurred in the CID processes and produced the characteristic fragment ions, which are especially valuable for the identification of this class of xanthones. The fragmentation differential between some cis-, trans-isomers was uncovered. Computation methods were utilized to rationalize the observed MS behaviors. On-line UHPLC-ESI-MS/MS/MS method has proved to be rapid and efficient in that within 6min, 15 caged scaffold xanthones, including three pairs of epimers and four pairs of isomers in gamboges, were effectively separated and identified. Among them, two known, namely isogambogenin (13) and isomorellinol (14) and one likely new caged Garcinia xanthones from the Garcinia hanburyi were tentatively characterized based on the tandem mass spectra of known ones.

Time-resolved photodissociation of singly protonated peptides with an arginine at the N-terminus: a statistical interpretation

Time-evolution of product ion signals in ultraviolet photodissociation (UV-PD) of singly protonated peptides with an arginine at the N-terminus was investigated by using a tandem time-of-flight mass spectrometer equipped with a cell floated at high voltage. Observation of different time-evolution patterns for different product ion types--an apparently nonstatistical behavior--could be explained within the statistical framework by invoking consecutive formation of some product ions and broad internal energy distributions for precursor ions. a(n) + 1 and b(n) ions were taken as the primary product ions from this type of peptide ions. Spectral characteristics in post-source decay, UV-PD, and collisionally activated dissociation at low and high kinetic energies could be explained via rough statistical calculation of rate constants. Specifically, the striking characteristics in high-energy CAD and UV-PD--dominance of a(n) and d(n) formed via a(n) + 1--were not due to the peculiarity of the excitation processes themselves, but due to quenching of the b(n) channels caused by the presence of arginine.

Efficient and reliable calculation of Rice-Ramsperger-Kassel-Marcus unimolecular reaction rate constants for biopolymers: modification of Beyer-Swinehart algorithm for degenerate vibrations

The Beyer-Swinehart (BS) algorithm, which calculates vibrational state density and sum, was modified for simultaneous treatment of degenerate vibrations. The modified algorithm was used in the grouped-frequency mode of the Rice-Ramsperger-Kassel-Marcus (RRKM) unimolecular reaction rate constant calculation for proteins with relative molecular mass as large as 100,000. Compared to the original BS method, reduction in computation time by a factor of around 3000 was achieved. Even though large systematic errors arising from frequency grouping were observed for state densities and sums, they more or less canceled each other, thus enabling reliable rate constant calculation. The present method is thought to be adequate for efficient and reliable RRKM calculations for any macromolecule in the gas phase such as the molecular ions of proteins, nucleic acids, and carbohydrates generated inside a mass spectrometer. The algorithm can also be used to calculate the internal energy distribution of a macromolecule at thermal equilibrium.