On the zwitterionic nature of gas-phase peptides and protein ions

High-Performance Molecular Dynamics Simulations for Native Mass Spectrometry of Large Protein Complexes with the Fast Multipole Method

Native mass spectrometry (MS) is widely employed to study the structures and assemblies of proteins ranging from small monomers to megadalton complexes. Molecular dynamics (MD) simulation is a useful complement as it provides the spatial detail that native MS cannot offer. However, MD simulations performed in the gas phase have suffered from rapidly increasing computational costs with the system size. The primary bottleneck is the calculation of electrostatic forces, which are effective over long distances and must be explicitly computed for each atom pair, precluding efficient use of methods traditionally used to accelerate condensed-phase simulations. As a result, MD simulations have been unable to match the capacity of MS in probing large multimeric protein complexes. Here, we apply the fast multipole method (FMM) for computing the electrostatic forces, recently implemented by Kohnke et al. (J. Chem. Theory Comput., 2020, 16, 6938-6949), showing that it significantly enhances the performance of gas-phase simulations of large proteins. We assess how to achieve adequate accuracy and optimal performance with FMM, finding that it expands the accessible size range and time scales dramatically. Additionally, we simulate a 460 kDa ferritin complex over microsecond time scales, alongside complementary ion mobility (IM)-MS experiments, uncovering conformational changes that are not apparent from the IM-MS data alone.

Observation of Solvent-Dependence in the Mechanism of Neutral-Catalyzed Isomerization of para-Aminobenzoic Acid Protomers

Proton-transfer reactions are commonplace during electrospray ionization (ESI) mass spectrometry experiments and are often responsible for imparting charge to analyte molecules. Multiple protonation-site isomers (protomers) can arise for polyfunctional molecules and these isomers can interconvert via solvent-mediated proton transfer reactions during various stages of the ESI process. Studying the populations and interconversion of protonation isomers provides key insight into the ESI process, ion-molecule interactions, and ion dissociation mechanisms. An archetype molecule to study protomer interconversion fundamentals in this context is para-aminobenzoic acid (pABA), where both the amino and carboxylic acid protomers are typically formed under ESI and the mechanisms for interconversion are still under refinement. Using ion-trap mass spectrometry reaction kinetics (2.5 mTorr, 300 K), this study examines gas-phase interconversion catalysis of pABA protomers by seven neutral species, which are commen solvents and additives used for ESI: water, formic acid, methanol, ethanol, propanol, ammonia, and acetonitrile. Three distinct reaction cases are reported: (i) formic acid, methanol, ethanol, propanol, and ammonia each catalyze the interconversion between the amino and carboxylic acid protomers via a n = 1 solvent-molecule vehicle mechanism; (ii) for water, however, a n = 6 adduct complex is detected and this suggests that the observed protomer interconversion occurs through a Grotthuss mechanism, in accord with literature reports; (iii) acetonitrile inhibits proton transfer by the formation of particularly stable n = 1 and 2 adduct complexes. The second-order rate constants for the protomer interconversion are observed to increase in the following order: H2O < HCO2H < MeOH < EtOH < PrOH < NH3. Potential energy schemes are reported for all neutral-catalyzed proton transfer reactions using the DSD-PBEP86-D3(BJ)/aug-cc-pVDZ level of theory. A central transition state, which connects the protonation site adducts, is shown to be the key rate-limiting step. The energy of this transition state is sensitive to the proton affinity of the neutral solvent, and this is supported by the correlation between the reaction rate and the solvent proton affinity.

Negative-Ion Electron Capture Dissociation of MALDI-Generated Peptide Anions

Negative-ion electron capture dissociation (niECD) is an anion MS/MS technique that provides fragmentation analogous to conventional ECD, including high peptide sequence coverage and retention of labile post-translational modifications (PTMs). niECD has been proposed to be the most efficient for salt-bridged zwitterionic precursor ion structures. Several important PTMs, e.g., sulfation and phosphorylation, are acidic and can, therefore, be challenging to characterize in the positive-ion mode. Furthermore, PTM-friendly techniques, such as ECD, require multiple precursor ion-positive charges. By contrast, singly charged ions, refractory to ECD, are most compatible with niECD. Because electrospray ionization (ESI) typically yields multiply charged ions, we sought to explore matrix-assisted laser desorption/ionization (MALDI) in combination with niECD. However, the requirement for zwitterionic gaseous structures may preclude efficient niECD of MALDI-generated anions. Unexpectedly, we found that niECD of anions from MALDI is not only possible but proceeds with similar or higher efficiency compared with ESI-generated anions. Matrix selection did not appear to have a major effect. With MALDI, niECD is demonstrated up to m/z ∼4300. For such larger analytes, multiple electron captures are observed, resulting in triply charged fragments from singly charged precursor ions. Such charge-increased fragments show improved detectability. Furthermore, significantly improved (∼20-fold signal-to-noise increase) niECD spectral quality is achieved with equivalent sample amounts from MALDI vs ESI. Overall, the reported combination with MALDI significantly boosts the analytical utility of niECD.

Uncovering Distinct Peptide Charging Behaviors in Electrospray Ionization Mass Spectrometry Using a Large-Scale Dataset

Electrospray ionization is a powerful and prevalent technique used to ionize analytes in mass spectrometry. The distribution of charges that an analyte receives (charge state distribution, CSD) is an important consideration for interpreting mass spectra. However, due to an incomplete understanding of the ionization mechanism, the analyte properties that influence CSDs are not fully understood. Here, we employ a machine learning-based approach and analyze CSDs of hundreds of thousands of peptides. Interestingly, half of the peptides exhibit charges that differ from what one would naively expect (the number of basic sites). We find that these peptides can be classified into two regimes (undercharging and overcharging) and that these two regimes display markedly different charging characteristics. Notably, peptides in the overcharging regime show minimal dependence on basic site count, and more generally, the two regimes exhibit distinct sequence determinants. These findings highlight the rich ionization behavior of peptides and the potential of CSDs for enhancing peptide identification.

Evaluation of Hydrogen-Deuterium Exchange during Transient Vapor Binding of MeOD with Model Peptide Systems Angiotensin II and Bradykinin

The advancement of hybrid mass spectrometric tools as an indirect probe of molecular structure and dynamics relies heavily upon a clear understanding between gas-phase ion reactivity and ion structural characteristics. This work provides new insights into gas-phase ion-neutral reactions of the model peptides (i.e., angiotensin II and bradykinin) on a per-residue basis by integrating hydrogen/deuterium exchange, ion mobility, tandem mass spectrometry, selective vapor binding, and molecular dynamics simulations. By comparing fragmentation patterns with simulated probabilities of vapor uptake, a clear link between gas-phase hydrogen/deuterium exchange and the probabilities of localized vapor association is established. The observed molecular dynamics trends related to the sites and duration of vapor binding track closely with experimental observation. Additionally, the influence of additional charges and structural characteristics on exchange kinetics and ion-neutral cluster formation is examined. These data provide a foundation for the analysis of solvation dynamics of larger, native-like conformations of proteins in the gas phase.

Insight on physicochemical properties governing peptide MS1 response in HPLC-ESI-MS/MS: A deep learning approach

Accurate and absolute quantification of peptides in complex mixtures using quantitative mass spectrometry (MS)-based methods requires foreground knowledge and isotopically labeled standards, thereby increasing analytical expenses, time consumption, and labor, thus limiting the number of peptides that can be accurately quantified. This originates from differential ionization efficiency between peptides and thus, understanding the physicochemical properties that influence the ionization and response in MS analysis is essential for developing less restrictive label-free quantitative methods. Here, we used equimolar peptide pool repository data to develop a deep learning model capable of identifying amino acids influencing the MS1 response. By using an encoder-decoder with an attention mechanism and correlating attention weights with amino acid physicochemical properties, we obtain insight on properties governing the peptide-level MS1 response within the datasets. While the problem cannot be described by one single set of amino acids and properties, distinct patterns were reproducibly obtained. Properties are grouped in three main categories related to peptide hydrophobicity, charge, and structural propensities. Moreover, our model can predict MS1 intensity output under defined conditions based solely on peptide sequence input. Using a refined training dataset, the model predicted log-transformed peptide MS1 intensities with an average error of 9.7 ± 0.5% based on 5-fold cross validation, and outperformed random forest and ridge regression models on both log-transformed and real scale data. This work demonstrates how deep learning can facilitate identification of physicochemical properties influencing peptide MS1 responses, but also illustrates how sequence-based response prediction and label-free peptide-level quantification may impact future workflows within quantitative proteomics.

Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins

Current experiments only provide low-resolution information on gaseous protein ions generated by electrospray ionization (ESI). Molecular dynamics (MD) simulations can yield complementary insights. Unfortunately, conventional MD does not capture the mobile nature of protons in gaseous proteins. Mobile-proton MD (MPMD) overcomes this limitation. Earlier MPMD data at 300 K indicated that protein ions generated by "native" ESI retain solution-like structures with a hydrophobic core and zwitterionic exterior [Bakhtiari, M.; Konermann, L. J. Phys. Chem. B 2019, 123, 1784-1796]. MPMD redistributes protons using electrostatic and proton affinity calculations. The robustness of this approach has never been scrutinized. Here, we close this gap by benchmarking MPMD against density functional theory (DFT) at the B3LYP/6-31G* level, which is well suited for predicting proton affinities. The computational cost of DFT necessitated the use of small peptides. The MPMD energetic ranking of proton configurations was found to be consistent with DFT single-point energies, implying that MPMD can reliably identify favorable protonation sites. Peptide MPMD runs converged to DFT-optimized structures only when applying 300-500 K temperature cycling, which was necessary to prevent trapping in local minima. Temperature cycling MPMD was then applied to gaseous protein ions. Native ubiquitin converted to slightly expanded structures with a zwitterionic core and a nonpolar exterior. Our data suggest that such inside-out protein structures are intrinsically preferred in the gas phase, and that they form in ESI experiments after moderate collisional excitation. This is in contrast to native ESI (with minimal collisional excitation, simulated by MPMD at 300 K), where kinetic trapping promotes the survival of solution-like structures. In summary, this work validates the MPMD approach for simulations on gaseous peptides and proteins.

A high-throughput approach reveals distinct peptide charging behaviors in electrospray ionization mass spectrometry

Electrospray ionization is a powerful and prevalent technique used to ionize analytes in mass spectrometry. The distribution of charges that an analyte receives (charge state distribution, CSD) is an important consideration for interpreting mass spectra. However, due to an incomplete understanding of the ionization mechanism, the analyte properties that influence CSDs are not fully understood. Here, we employ a machine learning-based high-throughput approach and analyze CSDs of hundreds of thousands of peptides. Interestingly, half of the peptides exhibit charges that differ from what one would naively expect (number of basic sites). We find that these peptides can be classified into two regimes-undercharging and overcharging-and that these two regimes display markedly different charging characteristics. Strikingly, peptides in the overcharging regime show minimal dependence on basic site count, and more generally, the two regimes exhibit distinct sequence determinants. These findings highlight the rich ionization behavior of peptides and the potential of CSDs for enhancing peptide identification.

Stability and conformational memory of electrosprayed and rehydrated bacteriophage MS2 virus coat proteins

Proteins are innately dynamic, which is important for their functions, but which also poses significant challenges when studying their structures. Gas-phase techniques can utilise separation and a range of sample manipulations to transcend some of the limitations of conventional techniques for structural biology in crystalline or solution phase, and isolate different states for separate interrogation. However, the transfer from solution to the gas phase risks affecting the structures, and it is unclear to what extent different conformations remain distinct in the gas phase, and if resolution in silico can recover the native conformations and their differences. Here, we use extensive molecular dynamics simulations to study the two distinct conformations of dimeric capsid protein of the MS2 bacteriophage. The protein undergoes notable restructuring of its peripheral parts in the gas phase, but subsequent simulation in solvent largely recovers the native structure. Our results suggest that despite some structural loss due to the experimental conditions, gas-phase structural biology techniques provide meaningful data that inform not only about the structures but also conformational dynamics of proteins.

Native mass spectrometry for the investigation of protein structural (dis)order

A central challenge in structural biology is represented by dynamic and heterogeneous systems, as typically represented by proteins in solution, with the extreme case of intrinsically disordered proteins (IDPs) [1-3]. These proteins lack a specific three-dimensional structure and have poorly organized secondary structure. For these reasons, they escape structural characterization by conventional biophysical methods. The investigation of these systems requires description of conformational ensembles, rather than of unique, defined structures or bundles of largely superimposable structures. Mass spectrometry (MS) has become a central tool in this field, offering a variety of complementary approaches to generate structural information on either folded or disordered proteins [4-6]. Two main categories of methods can be recognized. On one side, conformation-dependent reactions (such as cross-linking, covalent labeling, H/D exchange) are exploited to label molecules in solution, followed by the characterization of the labeling products by denaturing MS [7-11]. On the other side, non-denaturing ("native") MS can be used to directly explore the different conformational components in terms of geometry and structural compactness [12-16]. All these approaches have in common the capability to conjugate protein structure investigation with the peculiar analytical power of MS measurements, offering the possibility of assessing species distributions for folding and binding equilibria and the combination of both. These methods can be combined with characterization of noncovalent complexes [17, 18] and post-translational modifications [19-23]. This review focuses on the application of native MS to protein structure and dynamics investigation, with a general methodological section, followed by examples on specific proteins from our laboratory.

Protonation of Serine in Gas and Condensed and Microsolvated States in Aqueous Solution

Identification of molecules and elucidation of their chemical structure are ubiquitous problems in chemistry. Mass spectrometry (MS) can be used due to its sensitivity and versatility. For detection to occur, analytes must be ionized and transferred to the gas phase. Soft ionization processes such as electrospray ionization are popular; however, resulting microsolvated phases can alter the chemistry of analytes and therefore detection and identification. To understand these processes, we use computational methods to probe the ionization propensity of serine in the gas phase, aqueous microsolvated clusters, and aqueous solution. We show that the tautomeric form of serine is altered by the presence of water, as five water molecules can stabilize the zwitterionic tautomer. Inclusion of cosolutes such as ions can stabilize the zwitterion with as few as one or two water molecules present. We demonstrate that ionization propensity, as measured by gas phase bacisity, can increase by over 100 kJ/mol when placed in a small water-serine cluster, showing the sensitivity of the chemistry of microsolvated analytes. Finally, detailed analysis reveals that small droplets (less than seven water molecules) are extremely sensitive to addition of further water molecules. Beyond this limit, structural and electronic properties change little with droplet size.

Atomistic Insights into the Formation of Nonspecific Protein Complexes during Electrospray Ionization

Native electrospray ionization (ESI)-mass spectrometry (MS) is widely used for the detection and characterization of multi-protein complexes. A well-known problem with this approach is the possible occurrence of nonspecific protein clustering in the ESI plume. This effect can distort the results of binding affinity measurements, and it can even generate gas-phase complexes from proteins that are strictly monomeric in bulk solution. By combining experiments and molecular dynamics (MD) simulations, the current work for the first time provides detailed insights into the ESI clustering of proteins. Using ubiquitin as a model system, we demonstrate how the entrapment of more than one protein molecule in an ESI droplet can generate nonspecific clusters (e.g., dimers or trimers) via solvent evaporation to dryness. These events are in line with earlier proposals, according to which protein clustering is associated with the charged residue model (CRM). MD simulations on cytochrome c (which carries a large intrinsic positive charge) confirmed the viability of this CRM avenue. In addition, the cytochrome c data uncovered an alternative mechanism where protein-protein contacts were formed early within ESI droplets, followed by cluster ejection from the droplet surface. This second pathway is consistent with the ion evaporation model (IEM). The observation of these IEM events for large protein clusters is unexpected because the IEM has been thought to be associated primarily with low-molecular-weight analytes. In all cases, our MD simulations produced protein clusters that were stabilized by intermolecular salt bridges. The MD-generated charge states agreed with experiments. Overall, this work reveals that ESI-induced protein clustering does not follow a tightly orchestrated pathway but can proceed along different avenues.

Mobile Protons Limit the Stability of Salt Bridges in the Gas Phase: Implications for the Structures of Electrosprayed Protein Ions

Electrosprayed protein ions can retain native-like conformations. The intramolecular contacts that stabilize these compact gas-phase structures remain poorly understood. Recent work has uncovered abundant salt bridges in electrosprayed proteins. Salt bridges are zwitterionic BH⁺/A^- contacts. The low dielectric constant in the vacuum strengthens electrostatic interactions, suggesting that salt bridges could be a key contributor to the retention of compact protein structures. A problem with this assertion is that H⁺ are mobile, such that H⁺ transfer can convert salt bridges into neutral B⁰/HA⁰ contacts. This possible salt bridge annihilation puts into question the role of zwitterionic motifs in the gas phase, and it calls for a detailed analysis of BH⁺/A^- versus B⁰/HA⁰ interactions. Here, we investigate this issue using molecular dynamics (MD) simulations and electrospray experiments. MD data for short model peptides revealed that salt bridges with static H⁺ have dissociation energies around 700 kJ mol^-1. The corresponding B⁰/HA⁰ contacts are 1 order of magnitude weaker. When considering the effects of mobile H⁺, BH⁺/A^- bond energies were found to be between these two extremes, confirming that H⁺ migration can significantly weaken salt bridges. Next, we examined the protein ubiquitin under collision-induced unfolding (CIU) conditions. CIU simulations were conducted using three different MD models: (i) Positive-only runs with static H⁺ did not allow for salt bridge formation and produced highly expanded CIU structures. (ii) Zwitterionic runs with static H⁺ resulted in abundant salt bridges, culminating in much more compact CIU structures. (iii) Mobile H⁺ simulations allowed for the dynamic formation/annihilation of salt bridges, generating CIU structures intermediate between scenarios (i) and (ii). Our results uncover that mobile H⁺ limit the stabilizing effects of salt bridges in the gas phase. Failure to consider the effects of mobile H⁺ in MD simulations will result in unrealistic outcomes under CIU conditions.

Cryogenic "Iodide-Tagging" Photoelectron Spectroscopy: A Sensitive Probe for Specific Binding Sites of Amino Acids

This work showcases cryogenic and temperature-dependent "iodide-tagging" photoelectron spectroscopy to probe specific binding sites of amino acids using the glycine-iodide complex (Gly·I^-) as a case study. Multiple Gly·I^- isomers were generated from ambient electrospray ionization and kinetically isolated in a cryogenic ion trap. These structures were characterized with temperature-dependent "iodide-tagging" negative ion photoelectron spectroscopy (NIPES), where iodide was used as the "messenger" to interpret electronic energetics and structural information of various Gly·I^- isomers. Accompanied by theoretical computations and Franck-Condon simulations, a total of five cluster structures have been identified along with their various binding motifs. This work demonstrates that "iodide-tagging" NIPES is a powerful general means for probing specific binding interactions in biological molecules of interest.

Mechanistic Study of Enhanced Protonation by Chromium(III) in Electrospray Ionization: A Superacid Bound to a Peptide

Addition of trivalent chromium, Cr(III), to solutions undergoing electrospray ionization (ESI) enhances protonation and leads to formation of [M + 2H]²⁺ for peptides that normally produce [M + H]⁺. This effect is explored using electronic structure calculations at the density functional theory (DFT) level to predict the energetics of various species that are potentially important to the mechanism. Gas- and solution-phase reaction free energies for glycine and its anion reacting with [Cr(III)(H₂O)₆]³⁺ and for dehydration of these species have been predicted, where glycine is used as a simple model for a peptide. For comparison, calculations were also performed with Fe(III), Al(III), Sc(III), Y(III), and La(III). Removal of water from these complexes, as would occur during the ESI desolvation process, results in species that are highly acidic. The calculated pK_a of Cr(III) with a single solvation shell is -10.8, making [Cr(III)(H₂O)₆]³⁺ a superacid that is more acidic than sulfuric acid (pK_a = -8.8). Binding to glycine requires removal of two aqua ligands, which gives [Cr(III)(H₂O)₄]³⁺ that has an extremely acidic pK_a of -28.8. Removal of additional water further enhances acidity, reaching a pK_a of -84.7 for [Cr(III)(H₂O)]³⁺. A mechanism for enhanced protonation is proposed that incorporates computational and experiment results, as well as information on the known chemistry of Cr(III), which is substitutionally inert. The initial step involves binding of [Cr(III)(H₂O)₄]³⁺ to the deprotonated C-terminus of a peptide. As the drying process during ESI strips water from the complex, the resulting superacid transfers protons to the bound peptide, eventually leading to formation of [M + 2H]²⁺.

An arsenal of methods for the experimental characterization of intrinsically disordered proteins - How to choose and combine them?

In this review, we detail the most common experimental approaches to assess and characterize protein intrinsic structural disorder, with the notable exception of NMR and EPR spectroscopy, two ideally suited approaches that will be described in depth in two other reviews within this special issue. We discuss the advantages, the limitations, as well as the caveats of the various methods. We also describe less common and more demanding approaches that enable achieving further insights into the conformational properties of IDPs. Finally, we present recent developments that have enabled assessment of structural disorder in living cells, and discuss the currently available methods to model IDPs as conformational ensembles.

Preservation of Protein Zwitterionic States in the Transition from Solution to Gas Phase Revealed by Sodium Adduction Mass Spectrometry

The structural characterization of proteins and their interaction network mapping in the gas phase highlights the need to preserve their most nativelike conformers in the transition from the solution to gas phase. Zwitterionic interactions in a protein are weak bonds between oppositely charged residues, which make an important contribution to protein stability. However, it is still not clear whether the native zwitterionic states of proteins can be retained or not when it is transferred from the solution to gas phase. Using the nonspecific Na⁺ adduction as a novel signature, here we show that the zwitterionic states of proteins can be preserved when a moderated droplet desolvation condition (temperature <30 °C) is used in native electrospray ionization mass spectrometry. The very low-level nonspecific metal adduction to proteins under such conditions also enables rapid and direct determination of the binding states of metal-binding proteins and sensitive detection of proteins from solutions containing highly concentrated involatile salts (e.g., 50 mM NaCl). We believe that our findings can be instructive for performing mass spectrometric analysis of proteins and useful for protein ions desalting which simply involves altering the temperature and flow rate of drying gas in the desolvation region.

The Fragment Molecular Orbital Method Based on Long-Range Corrected Density-Functional Tight-Binding

The presently available linear scaling approaches to density-functional tight-binding (DFTB) based on the fragment molecular orbital (FMO) method are severely impacted by the problem of artificial charge transfer due to the self-interaction error (SIE), which hampers the simulation of zwitterionic systems such as biopolymers or ionic liquids. Here we report an extension of FMO-DFTB where we included a long-range corrected (LC) functional designed to mitigate the DFTB SIE, called the FMO-LC-DFTB method, resulting in a robust method which succeeds in simulating zwitterionic systems. Both energy and analytic gradient are developed for the gas phase and the polarizable continuum model of solvation. The scaling of FMO-LC-DFTB with system size N is shown to be almost linear, O( N^1.13-1.28), and its numerical accuracy is established for a variety of representative systems including neutral and charged polypeptides. It is shown that pair interaction energies between fragments for two mini-proteins are in excellent agreement with results from long-range corrected density functional theory. The new method was employed in long time scale (1 ns) molecular dynamics simulations of the tryptophan cage protein (PDB: 1L2Y ) in the gas phase for four different protonation states and in stochastic global minimum structure searches for 1-ethyl-3-methylimidazolium nitrate ionic liquid clusters containing up to 2300 atoms.

Gas-Phase Ion/Ion Chemistry as a Probe for the Presence of Carboxylate Groups in Polypeptide Cations

The reactivity of 1-hydroxybenzoyl triazole (HOBt) esters with the carboxylate functionality present in peptides is demonstrated in the gas phase with a doubly deprotonated dianion. The reaction forms an anhydride linkage at the carboxylate site. Upon ion trap collisional-induced dissociation (CID) of the modified peptide, the resulting spectrum shows a nominal loss of the mass of the reagent and a water molecule. Analogous phenomenology was also noted for model peptide cations that likely contain zwitterionic/salt-bridged motifs in reactions with a negatively charged HOBt ester. Control experiments indicate that a carboxylate group is the likely reactive site, rather than other possible nucleophilic sites present in the peptide. These observations suggest that HOBt ester chemistry may be used as a chemical probe for the presence and location of carboxylate groups in net positively charged polypeptide ions. As an illustration, deprotonated sulfobenzoyl HOBt was reacted with the [M+7H]⁷⁺ ion of ubiquitin. The ion was shown to react with the reagent and CID of the covalent reaction product yielded an abundant [M+6H-H₂O]⁶⁺ ion. Comparison of the CID product ion spectrum of this ion with that of the water loss product generated from CID of the unmodified [M+6H]⁶⁺ ion revealed the glutamic acid at residue 64 as a reactive site, suggesting that it is present in the deprotonated form. Graphical Abstract ᅟ.

Size-Dependent Hydrogen Atom Attachment to Gas-Phase Hydrogen-Deficient Polypeptide Radical Cations

Despite significant affinity to carbonyl oxygens, thermal hydrogen atoms attach to unmodified polypeptides at a very low rate, while the hydrogen-hydrogen exchange rate is high. Here, using the novel omnitrap setup, we found that attachment to polypeptides is much more facile when radical site is already present, but the rate decreases for larger radical ions. The likely explanation is the intramolecular hydrogen atom rearrangement in hydrogen-deficient radicals to a more stable or less accessible site.

Photoelectron Transfer Dissociation Reveals Surprising Favorability of Zwitterionic States in Large Gaseous Peptides and Proteins

Structural characterization of proteins in the gas phase is becoming increasingly popular, highlighting the need for a greater understanding of how proteins behave in the absence of solvent. It is clear that charged residues exert significant influence over structures in the gas phase due to strong Coulombic and hydrogen-bonding interactions. The net charge for a gaseous ion is easily identified by mass spectrometry, but the presence of zwitterionic pairs or salt bridges has previously been more difficult to detect. We show that these sites can be revealed by photoinduced electron transfer dissociation, which produces characteristic c and z ions only if zwitterionic species are present. Although previous work on small molecules has shown that zwitterionic pairs are rarely stable in the gas phase, we now demonstrate that charge-separated states are favored in larger molecules. Indeed, we have detected zwitterionic pairs in peptides and proteins where the net charge equals the number of basic sites, requiring additional protonation at nonbasic residues. For example, the small protein ubiquitin can sustain a zwitterionic conformer for all charge states up to 14+, despite having only 13 basic sites. Virtually all of the peptides/proteins examined herein contain zwitterionic sites if both acidic and basic residues are present and the overall charge density is low. This bias in favor of charge-separated states has important consequences for efforts to model gaseous proteins via computational analysis, which should consider not only charge state isomers that include salt bridges but also protonation at nonbasic residues.

Proton Dynamics in Protein Mass Spectrometry

Native electrospray ionization/ion mobility-mass spectrometry (ESI/IM-MS) allows an accurate determination of low-resolution structural features of proteins. Yet, the presence of proton dynamics, observed already by us for DNA in the gas phase, and its impact on protein structural determinants, have not been investigated so far. Here, we address this issue by a multistep simulation strategy on a pharmacologically relevant peptide, the N-terminal residues of amyloid-β peptide (Aβ(1-16)). Our calculations reproduce the experimental maximum charge state from ESI-MS and are also in fair agreement with collision cross section (CCS) data measured here by ESI/IM-MS. Although the main structural features are preserved, subtle conformational changes do take place in the first ∼0.1 ms of dynamics. In addition, intramolecular proton dynamics processes occur on the picosecond-time scale in the gas phase as emerging from quantum mechanics/molecular mechanics (QM/MM) simulations at the B3LYP level of theory. We conclude that proton transfer phenomena do occur frequently during fly time in ESI-MS experiments (typically on the millisecond time scale). However, the structural changes associated with the process do not significantly affect the structural determinants.

Are Charge-State Distributions a Reliable Tool Describing Molecular Ensembles of Intrinsically Disordered Proteins by Native MS?

Native mass spectrometry (MS) has become a central tool of structural proteomics, but its applicability to the peculiar class of intrinsically disordered proteins (IDPs) is still object of debate. IDPs lack an ordered tridimensional structure and are characterized by high conformational plasticity. Since they represent valuable targets for cancer and neurodegeneration research, there is an urgent need of methodological advances for description of the conformational ensembles populated by these proteins in solution. However, structural rearrangements during electrospray-ionization (ESI) or after the transfer to the gas phase could affect data obtained by native ESI-MS. In particular, charge-state distributions (CSDs) are affected by protein conformation inside ESI droplets, while ion mobility (IM) reflects protein conformation in the gas phase. This review focuses on the available evidence relating IDP solution ensembles with CSDs, trying to summarize cases of apparent consistency or discrepancy. The protein-specificity of ionization patterns and their responses to ligands and buffer conditions suggests that CSDs are imprinted to protein structural features also in the case of IDPs. Nevertheless, it seems that these proteins are more easily affected by electrospray conditions, leading in some cases to rearrangements of the conformational ensembles. Graphical Abstract ᅟ.

Salt Bridge Rearrangement (SaBRe) Explains the Dissociation Behavior of Noncovalent Complexes

Native electrospray ionization-mass spectrometry, with gas-phase activation and solution compositions that partially release subcomplexes, can elucidate topologies of macromolecular assemblies. That so much complexity can be preserved in gas-phase assemblies is remarkable, although a long-standing conundrum has been the differences between their gas- and solution-phase decompositions. Collision-induced dissociation of multimeric noncovalent complexes typically distributes products asymmetrically (i.e., by ejecting a single subunit bearing a large percentage of the excess charge). That unexpected behavior has been rationalized as one subunit "unfolding" to depart with more charge. We present an alternative explanation based on heterolytic ion-pair scission and rearrangement, a mechanism that inherently partitions charge asymmetrically. Excessive barriers to dissociation are circumvented in this manner, when local charge rearrangements access a lower-barrier surface. An implication of this ion pair consideration is that stability differences between high- and low-charge state ions usually attributed to Coulomb repulsion may, alternatively, be conveyed by attractive forces from ion pairs (salt bridges) stabilizing low-charge state ions. Should the number of ion pairs be roughly inversely related to charge, symmetric dissociations would be favored from highly charged complexes, as observed. Correlations between a gas-phase protein's size and charge reflect the quantity of restraining ion pairs. Collisionally-facilitated salt bridge rearrangement (SaBRe) may explain unusual size "contractions" seen for some activated, low charge state complexes. That some low-charged multimers preferentially cleave covalent bonds or shed small ions to disrupting noncovalent associations is also explained by greater ion pairing in low charge state complexes. Graphical Abstract ᅟ.

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.

Conformational effects in protein electrospray-ionization mass spectrometry

Electrospray-ionization mass spectrometry (ESI-MS) is a key tool of structural biology, complementing the information delivered by conventional biochemical and biophysical methods. Yet, the mechanism behind the conformational effects in protein ESI-MS is an object of debate. Two parameters-solvent-accessible surface area (As) and apparent gas-phase basicity (GBapp)-are thought to play a role in controlling the extent of protein ionization during ESI-MS experiments. This review focuses on recent experimental and theoretical investigations concerning the influence of these parameters on ESI-MS results and the structural information that can be derived. The available evidence supports a unified model for the ionization mechanism of folded and unfolded proteins. These data indicate that charge-state distribution (CSD) analysis can provide valuable structural information on normally folded, as well as disordered structures.

The use of chromium(III) to supercharge peptides by protonation at low basicity sites

The addition of chromium(III) nitrate to solutions of peptides with seven or more residues greatly increases the formation of doubly protonated peptides, [M + 2H](2+), by electrospray ionization. The test compound heptaalanine has only one highly basic site (the N-terminal amino group) and undergoes almost exclusive single protonation using standard solvents. When Cr(III) is added to the solution, abundant [M + 2H](2+) forms, which involves protonation of the peptide backbone or the C-terminus. Salts of Al(III), Mn(II), Fe(III), Fe(II), Cu(II), Zn (II), Rh(III), La(III), Ce(IV), and Eu(III) were also studied. Although several metal ions slightly enhance protonation, Cr(III) has by far the greatest ability to generate [M + 2H](2+). Cr(III) does not supercharge peptide methyl esters, which suggests that the mechanism involves interaction of Cr(III) with a carboxylic acid group. Other factors may include the high acidity of hexa-aquochromium(III) and the resistance of Cr(III) to reduction. Nitrate salts enhance protonation more than chloride salts and a molar ratio of 10:1 Cr(III):peptide produces the most intense [M + 2H](2+). Cr(III) also supercharges numerous other small peptides, including highly acidic species. For basic peptides, Cr(III) increases the charge state (2+ versus 1+) and causes the number of peptide molecules being protonated to double or triple. Chromium(III) does not supercharge the proteins cytochrome c and myoglobin. The ability of Cr(III) to enhance [M + 2H](2+) intensity may prove useful in tandem mass spectrometry because of the resulting overall increase in signal-to-noise ratio, the fact that [M + 2H](2+) generally dissociate more readily than [M + H](+), and the ability to produce [M + 2H](2+) precursors for electron-based dissociation techniques.

What protein charging (and supercharging) reveal about the mechanism of electrospray ionization

Understanding the charging mechanism of electrospray ionization is central to overcoming shortcomings such as ion suppression or limited dynamic range, and explaining phenomena such as supercharging. Towards that end, we explore what accumulated observations reveal about the mechanism of electrospray. We introduce the idea of an intermediate region for electrospray ionization (and other ionization methods) to account for the facts that solution charge state distributions (CSDs) do not correlate with those observed by ESI-MS (the latter bear more charge) and that gas phase reactions can reduce, but not increase, the extent of charging. This region incorporates properties (e.g., basicities) intermediate between solution and gas phase. Assuming that droplet species polarize within the high electric field leads to equations describing ion emission resembling those from the equilibrium partitioning model. The equations predict many trends successfully, including CSD shifts to higher m/z for concentrated analytes and shifts to lower m/z for sprays employing smaller emitter opening diameters. From this view, a single mechanism can be formulated to explain how reagents that promote analyte charging ("supercharging") such as m-NBA, sulfolane, and 3-nitrobenzonitrile increase analyte charge from "denaturing" and "native" solvent systems. It is suggested that additives' Brønsted basicities are inversely correlated to their ability to shift CSDs to lower m/z in positive ESI, as are Brønsted acidities for negative ESI. Because supercharging agents reduce an analyte's solution ionization, excess spray charge is bestowed on evaporating ions carrying fewer opposing charges. Brønsted basicity (or acidity) determines how much ESI charge is lost to the agent (unavailable to evaporating analyte).

What protein charging (and supercharging) reveal about the mechanism of electrospray ionization

Understanding the charging mechanism of electrospray ionization is central to overcoming shortcomings such as ion suppression or limited dynamic range, and explaining phenomena such as supercharging. Towards that end, we explore what accumulated observations reveal about the mechanism of electrospray. We introduce the idea of an intermediate region for electrospray ionization (and other ionization methods) to account for the facts that solution charge state distributions (CSDs) do not correlate with those observed by ESI-MS (the latter bear more charge) and that gas phase reactions can reduce, but not increase, the extent of charging. This region incorporates properties (e.g., basicities) intermediate between solution and gas phase. Assuming that droplet species polarize within the high electric field leads to equations describing ion emission resembling those from the equilibrium partitioning model. The equations predict many trends successfully, including CSD shifts to higher m/z for concentrated analytes and shifts to lower m/z for sprays employing smaller emitter opening diameters. From this view, a single mechanism can be formulated to explain how reagents that promote analyte charging ("supercharging") such as m-NBA, sulfolane, and 3-nitrobenzonitrile increase analyte charge from "denaturing" and "native" solvent systems. It is suggested that additives' Brønsted basicities are inversely correlated to their ability to shift CSDs to lower m/z in positive ESI, as are Brønsted acidities for negative ESI. Because supercharging agents reduce an analyte's solution ionization, excess spray charge is bestowed on evaporating ions carrying fewer opposing charges. Brønsted basicity (or acidity) determines how much ESI charge is lost to the agent (unavailable to evaporating analyte).

Molecular simulation-based structural prediction of protein complexes in mass spectrometry: the human insulin dimer

Protein electrospray ionization (ESI) mass spectrometry (MS)-based techniques are widely used to provide insight into structural proteomics under the assumption that non-covalent protein complexes being transferred into the gas phase preserve basically the same intermolecular interactions as in solution. Here we investigate the applicability of this assumption by extending our previous structural prediction protocol for single proteins in ESI-MS to protein complexes. We apply our protocol to the human insulin dimer (hIns2) as a test case. Our calculations reproduce the main charge and the collision cross section (CCS) measured in ESI-MS experiments. Molecular dynamics simulations for 0.075 ms show that the complex maximizes intermolecular non-bonded interactions relative to the structure in water, without affecting the cross section. The overall gas-phase structure of hIns2 does exhibit differences with the one in aqueous solution, not inferable from a comparison with calculated CCS. Hence, care should be exerted when interpreting ESI-MS proteomics data based solely on NMR and/or X-ray structural information.

Extracting structural information from charge-state distributions of intrinsically disordered proteins by non-denaturing electrospray-ionization mass spectrometry

Intrinsically disordered proteins (IDPs) exert key biological functions but tend to escape identification and characterization due to their high structural dynamics and heterogeneity. The possibility to dissect conformational ensembles by electrospray-ionization mass spectrometry (ESI-MS) offers an attracting possibility to develop a signature for this class of proteins based on their peculiar ionization behavior. This review summarizes available data on charge-state distributions (CSDs) obtained for IDPs by non-denaturing ESI-MS, with reference to globular or chemically denatured proteins. The results illustrate the contributions that direct ESI-MS analysis can give to the identification of new putative IDPs and to their conformational investigation.

A computational model for protein ionization by electrospray based on gas-phase basicity

Identifying the key factor(s) governing the overall protein charge is crucial for the interpretation of electrospray-ionization mass spectrometry data. Current hypotheses invoke different principles for folded and unfolded proteins. Here, first we investigate the gas-phase structure and energetics of several proteins of variable size and different folds. The conformer and protomer space of these proteins ions is explored exhaustively by hybrid Monte-Carlo/molecular dynamics calculations, allowing for zwitterionic states. From these calculations, the apparent gas-phase basicity of desolvated protein ions turns out to be the unifying trait dictating protein ionization by electrospray. Next, we develop a simple, general, adjustable-parameter-free model for the potential energy function of proteins. The model is capable to predict with remarkable accuracy the experimental charge of folded proteins and its well-known correlation with the square root of protein mass.

Computational and statistical analysis of protein mass spectrometry data

High-throughput proteomics experiments involving tandem mass spectrometry produce large volumes of complex data that require sophisticated computational analyses. As such, the field offers many challenges for computational biologists. In this article, we briefly introduce some of the core computational and statistical problems in the field and then describe a variety of outstanding problems that readers of PLoS Computational Biology might be able to help solve.

Separation of a set of peptide sequence isomers using differential ion mobility spectrometry

Protein identification in bottom-up proteomics requires disentangling isomers of proteolytic peptides, a major class of which are sequence inversions. Their separation using ion mobility spectrometry (IMS) has been limited to isomeric pairs. Here we demonstrate baseline separation of all seven 8-mer tryptic peptide isomers using differential IMS. Evaluation of peak capacity implies that even larger libraries should be resolved for heavier peptides with higher charge states.