Charge-State Dependent Compaction and Dissociation of Protein Complexes: Insights from Ion Mobility and Molecular Dynamics

Revealing the Fates of Proteins in the Gas Phase

The ability to observe intact proteins by native mass spectrometry allows measurements of size, oligomeric state, numbers and types of ligands and post translational modifications bound, among many other characteristics. These studies have the potential to, and in some cases are, advancing our understanding of the role of structure in protein biology and biochemistry. However, there are some long-unresolved questions about to what extent solution-like structures persist without solvent in the vacuum of the mass spectrometer. Strong evidence from multiple sources over the years has demonstrated that well-folded proteins maintain native-like states if care is taken during sample preparation, ionization, and transmission through the gas phase. For partially unfolded states, dynamic and disordered proteins, and other important landmarks along the protein folding/unfolding pathway, caution has been urged in the interpretation of the results of native ion mobility/mass spectrometric data. New gas-phase tools allow us to provide insight into these questions with in situ, in vacuo labeling reactions delivered through ion/ion chemistry. This Young Scientist Perspective demonstrates the robustness of these tools in describing native-like structure as well as possible deviations from native-like structure during native ion mobility/mass spectrometry. This Perspective illustrates some of the changes in structure produced by the removal of solvent and details some of the challenges and potential of the field.

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.

Protein Complex Heterogeneity and Topology Revealed by Electron Capture Charge Reduction and Surface Induced Dissociation

We illustrate the utility of native mass spectrometry (nMS) combined with a fast, tunable gas-phase charge reduction, electron capture charge reduction (ECCR), for the characterization of protein complex topology and glycoprotein heterogeneity. ECCR efficiently reduces the charge states of tetradecameric GroEL, illustrating Orbitrap m/z measurements to greater than 100,000 m/z. For pentameric C-reactive protein and tetradecameric GroEL, our novel device combining ECCR with surface induced dissociation (SID) reduces the charge states and yields more topologically informative fragmentation. This is the first demonstration that ECCR yields more native-like SID fragmentation. ECCR also significantly improved mass and glycan heterogeneity measurements of heavily glycosylated SARS-CoV-2 spike protein trimer and thyroglobulin dimer. Protein glycosylation is important for structural and functional properties and plays essential roles in many biological processes. The immense heterogeneity in glycosylation sites and glycan structure poses significant analytical challenges that hinder a mechanistic understanding of the biological role of glycosylation. Without ECCR, average mass determination of glycoprotein complexes is available only through charge detection mass spectrometry or mass photometry. With narrow m/z selection windows followed by ECCR, multiple glycoform m/z values are apparent, providing quick global glycoform profiling and providing a future path for glycan localization on individual intact glycoforms.

Protein complex heterogeneity and topology revealed by electron capture charge reduction and surface induced dissociation

We illustrate the utility of native mass spectrometry (nMS) combined with a fast, tunable gas-phase charge reduction, electron capture charge reduction (ECCR), for the characterization of protein complex topology and glycoprotein heterogeneity. ECCR efficiently reduces the charge states of tetradecameric GroEL, illustrating Orbitrap m/z measurements to greater than 100,000 m/z. For pentameric C-reactive protein and tetradecameric GroEL, our novel device combining ECCR with surface induced dissociation (SID) reduces the charge states and yields more topologically informative fragmentation. This is the first demonstration that ECCR yields more native-like SID fragmentation. ECCR also significantly improved mass and glycan heterogeneity measurements of heavily glycosylated SARS-CoV-2 spike protein trimer and thyroglobulin dimer. Protein glycosylation is important for structural and functional properties and plays essential roles in many biological processes. The immense heterogeneity in glycosylation sites and glycan structure poses significant analytical challenges that hinder a mechanistic understanding of the biological role of glycosylation. Without ECCR, average mass determination of glycoprotein complexes is available only through charge detection mass spectrometry or mass photometry. With narrow m/z selection windows followed by ECCR, multiple glycoform m/z values are apparent, providing quick global glycoform profiling and providing a future path for glycan localization on individual intact glycoforms.

Study of the Conformational Dynamics of Prolyl Oligopeptidase by Mass Spectrometry: Lessons Learned

Ion mobility mass spectrometry (IM-MS) can be used to analyze native proteins according to their size and shape. By sampling individual molecules, it allows us to study mixtures of conformations, as long as they have different collision cross sections and maintain their native conformation after dehydration and vaporization in the mass spectrometer. Even though conformational heterogeneity of prolyl oligopeptidase has been demonstrated in solution, it is not detectable in IM-MS. Factors that affect the conformation in solution, binding of an active site ligand, the stabilizing Ser554Ala mutation, and acidification do not qualitatively affect the collision-induced unfolding pattern. However, measuring the protection of accessible cysteines upon ligand binding provides a principle for the development of MS-based ligand screening methods.

Ion Mobility-Mass Spectrometry and Collision-Induced Unfolding Rapidly Characterize the Structural Polydispersity and Stability of an Fc-Fusion Protein

Fc-fusion proteins are an emerging class of protein therapeutics that combine the properties of biological ligands with the unique properties of the fragment crystallizable (Fc) domain of an immunoglobulin G (IgG). Due to their diverse higher-order structures (HOSs), Fc-fusion proteins remain challenging characterization targets within biopharmaceutical pipelines. While high-resolution biophysical tools are available for HOS characterization, they frequently demand extended time frames and substantial quantities of purified samples, rendering them impractical for swiftly screening candidate molecules. Herein, we describe the development of ion mobility-mass spectrometry (IM-MS) and collision-induced unfolding (CIU) workflows that aim to fill this technology gap, where we focus on probing the HOS of a model Fc-Interleukin-10 (Fc-IL-10) fusion protein engineered using flexible glycine-serine linkers. We evaluate the ability of these techniques to probe the flexibility of Fc-IL-10 in the absence of bulk solvent relative to other proteins of similar size, as well as localize structural changes of low charge state Fc-IL-10 ions to specific Fc and IL-10 unfolding events during CIU. We subsequently apply these tools to probe the local effects of glycine-serine linkers on the HOS and stability of IL-10 homodimer, which is the biologically active form of IL-10. Our data reveals that Fc-IL-10 produces significantly more structural transitions during CIU and broader IM profiles when compared to a wide range of model proteins, indicative of its exceptional structural dynamism. Furthermore, we use a combination of enzymatic approaches to annotate these intricate CIU data and localize specific transitions to the unfolding of domains within Fc-IL-10. Finally, we detect a strong positive, quadratic relationship between average linker mass and fusion protein stability, suggesting a cooperative influence between glycine-serine linkers and overall fusion protein stability. This is the first reported study on the use of IM-MS and CIU to characterize HOS of Fc-fusion proteins, illustrating the practical applicability of this approach.

Top-down mass spectrometry of native proteoforms and their complexes: a community study

The combination of native electrospray ionization with top-down fragmentation in mass spectrometry (MS) allows simultaneous determination of the stoichiometry of noncovalent complexes and identification of their component proteoforms and cofactors. Although this approach is powerful, both native MS and top-down MS are not yet well standardized, and only a limited number of laboratories regularly carry out this type of research. To address this challenge, the Consortium for Top-Down Proteomics initiated a study to develop and test protocols for native MS combined with top-down fragmentation of proteins and protein complexes across 11 instruments in nine laboratories. Here we report the summary of the outcomes to provide robust benchmarks and a valuable entry point for the scientific community.

Native Mass Spectrometry Dissects the Structural Dynamics of an Allosteric Heterodimer of SARS-CoV-2 Nonstructural Proteins

Structure-based drug design, which relies on precise understanding of the target protein and its interaction with the drug candidate, is dramatically expedited by advances in computational methods for candidate prediction. Yet, the accuracy needs to be improved with more structural data from high throughput experiments, which are challenging to generate, especially for dynamic and weak associations. Herein, we applied native mass spectrometry (native MS) to rapidly characterize ligand binding of an allosteric heterodimeric complex of SARS-CoV-2 nonstructural proteins (nsp) nsp10 and nsp16 (nsp10/16), a complex essential for virus survival in the host and thus a desirable drug target. Native MS showed that the dimer is in equilibrium with monomeric states in solution. Consistent with the literature, well characterized small cosubstrate, RNA substrate, and product bind with high specificity and affinity to the dimer but not the free monomers. Unsuccessfully designed ligands bind indiscriminately to all forms. Using neutral gas collision, the nsp16 monomer with bound cosubstrate can be released from the holo dimer complex, confirming the binding to nsp16 as revealed by the crystal structure. However, we observed an unusual migration of the endogenous zinc ions bound to nsp10 to nsp16 after collisional dissociation. The metal migration can be suppressed by using surface collision with reduced precursor charge states, which presumably resulted in minimal gas-phase structural rearrangement and highlighted the importance of complementary techniques. With minimal sample input (∼μg), native MS can rapidly detect ligand binding affinities and locations in dynamic multisubunit protein complexes, demonstrating the potential of an "all-in-one" native MS assay for rapid structural profiling of protein-to-AI-based compound systems to expedite drug discovery.

Collision induced unfolding and molecular dynamics simulations of norovirus capsid dimers reveal strain-specific stability profiles

Collision induced unfolding (CIU) is a method used with ion mobility mass spectrometry to examine protein structures and their stability. Such experiments yield information about higher order protein structures, yet are unable to provide details about the underlying processes. That information can however be provided using molecular dynamics simulations. Here, we investigate the gas-phase unfolding of norovirus capsid dimers from the Norwalk and Kawasaki strains by employing molecular dynamics simulations over a range of temperatures, representing different levels of activation, together with CIU experiments. The dimers have highly similar structures, but their CIU reveals different stability that can be explained by the different dynamics that arises in response to the activation seen in the simulations, including a part of the sequence with previously observed strain-specific dynamics in solution. Our findings show how similar protein variants can be examined using mass spectrometric techniques in conjunction with atomistic molecular dynamics simulations to reveal differences in stability as well as differences in how and where unfolding takes place upon activation.

Collision-Induced Unfolding Reveals Disease-Associated Stability Shifts in Mitochondrial Transfer Ribonucleic Acids

Ribonucleic acids (RNAs) remain challenging targets for structural biology, creating barriers to understanding their vast functions in cellular biology and fully realizing their applications in biotechnology. The inherent dynamism of RNAs creates numerous obstacles in capturing their biologically relevant higher-order structures (HOSs), and as a result, many RNA functions remain unknown. In this study, we describe the development of native ion mobility-mass spectrometry and collision-induced unfolding (CIU) for the structural characterization of a variety of RNAs. We evaluate the ability of these techniques to preserve native structural features in the gas phase across a wide range of functional RNAs. Finally, we apply these tools to study the elusive mitochondrial encephalopathy, lactic acidosis, and stroke-like episodes-associated A3243G mutation. Our data demonstrate that our experimentally determined conditions preserve some solution-state memory of RNAs via the correlated complexity of CIU fingerprints and RNA HOS, the observation of predicted stability shifts in the control RNA samples, and the retention of predicted magnesium binding events in gas-phase RNA ions. Significant differences in collision cross section and stability are observed as a function of the A3243G mutation across a subset of the mitochondrial tRNA maturation pathway. We conclude by discussing the potential application of CIU for the development of RNA-based biotherapeutics and, more broadly, transcriptomic characterization.

Transition Metal Ion FRET-Based Probe to Study Cu(II)-Mediated Amyloid- Ligand Binding

Recent therapeutic strategies suggest that small peptides can act as aggregation inhibitors of monomeric amyloid-β (Αβ) by inducing structural rearrangements upon complexation. However, characterizing the binding events in such dynamic and transient noncovalent complexes, especially in the presence of natively occurring metal ions, remains a challenge. Here, we deploy a combined transition metal ion Förster resonance energy transfer (tmFRET) and native ion mobility-mass spectrometry (IM-MS) approach to characterize the structure of mass- and charge-selected Aβ complexes with Cu(II) ions (a quencher) and a potential aggregation inhibitor, a small neuropeptide named leucine enkephalin (LE). We show conformational changes of monomeric Αβ species upon Cu(II)-binding, indicating an uncoiled N-terminus and a close interaction between the C-terminus and the central hydrophobic region. Furthermore, we introduce LE labeled at the N-terminus with a metal-chelating agent, nitrilotriacetic acid (NTA). This allows us to employ tmFRET to probe the binding even in low-abundance and transient Aβ-inhibitor-metal ion complexes. Complementary intramolecular distance and global shape information from tmFRET and native IM-MS, respectively, confirmed Cu(II) displacement toward the N-terminus of Αβ, which discloses the binding region and the inhibitor's orientation.

Quantitatively Differentiating Antibodies Using Charge-State Manipulation, Collisional Activation, and Ion Mobility-Mass Spectrometry

Antibody-based therapeutics continue to expand both in the number of products and in their use in patients. These heterogeneous proteins challenge traditional drug characterization strategies, but ion mobility (IM) and mass spectrometry (MS) approaches have eased the challenge of higher-order structural characterization. Energy-dependent IM-MS, e.g., collision-induced unfolding (CIU), has been demonstrated to be sensitive to subtle differences in structure. In this study, we combine a charge-reduction method, cation-to-anion proton-transfer reactions (CAPTR), with energy-dependent IM-MS and varied solution conditions to probe their combined effects on the gas-phase structures of IgG1κ and IgG4κ from human myeloma. CAPTR paired with MS-only analysis improves the confidence of charge-state assignments and the resolution of the interfering protein species. Collision cross-section distributions were determined for each of the charge-reduced products. Similarity scoring was used to quantitatively compare distributions determined from matched experiments analyzing samples of the two antibodies. Relative to workflows using energy-dependent IM-MS without charge-state manipulation, combining CAPTR and energy-dependent IM-MS enhanced the differentiation of these antibodies. Combined, these results indicate that CAPTR can benefit many aspects of antibody characterization and differentiation.

Collision-Induced Unfolding, Tandem MS, Bottom-up Proteomics, and Interactomics for Identification of Protein Complexes in Native Surface Mass Spectrometry

Endogenously occurring salts and nonvolatile matrix components in untreated biological surfaces can suppress protein ionization and promote adduct formation, challenging protein identification. Characterization of labile proteins within biological specimens is particularly demanding because additional purification or sample treatment steps can be time-intensive and can disrupt noncovalent interactions. It is demonstrated that the combined use of collision-induced unfolding, tandem mass spectrometry, and bottom-up proteomics improves protein characterization in native surface mass spectrometry (NSMS). This multiprong analysis is achieved by acquiring NSMS, MS/MS, ion mobility (IM), and bottom-up proteomics data from a single surface extracted sample. The validity of this multiprong approach was confirmed by the successful characterization of nine surface-deposited proteins, with molecular weights ranging from 8 to 147 kDa, in two separate mixtures. Bottom-up proteomics provided a list of proteins to match against observed proteins in NSMS and their detected subunits in tandem MS. The method was applied to characterize endogenous proteins from untreated chicken liver samples. The subcapsular liver sampling for NSMS analysis allowed for the detection of endogenous proteins with molecular weights of up to ∼220 kDa. Moreover, using IM-MS, collision cross sections and collision-induced unfolding pathways of enzymatic proteins and protein complexes of up to 145 kDa were obtained.

Giving the Green Light to Photochemical Uncaging of Large Biomolecules in High Vacuum

The isolation of biomolecules in a high vacuum enables experiments on fragile species in the absence of a perturbing environment. Since many molecular properties are influenced by local electric fields, here we seek to gain control over the number of charges on a biopolymer by photochemical uncaging. We present the design, modeling, and synthesis of photoactive molecular tags, their labeling to peptides and proteins as well as their photochemical validation in solution and in the gas phase. The tailored tags can be selectively cleaved off at a well-defined time and without the need for any external charge-transferring agents. The energy of a single or two green photons can already trigger the process, and it is soft enough to ensure the integrity of the released biomolecular cargo. We exploit differences in the cleavage pathways in solution and in vacuum and observe a surprising robustness in upscaling the approach from a model system to genuine proteins. The interaction wavelength of 532 nm is compatible with various biomolecular entities, such as oligonucleotides or oligosaccharides.

Evaluation for Ion Heating of H2A-H2B Dimer in Ion Mobility Spectrometry-Mass Spectrometry

Ion mobility spectrometry-mass spectrometry (IMS-MS) provides m/z values and collision cross sections (CCSs) of gas-phase ions. In our previous study, an intrinsically disordered protein, the H2A-H2B dimer, was analyzed using IMS-MS, resulting in two conformational populations of CCS. Based on experimental and theoretical approaches, this resulted from a structural diversity of intrinsically disordered regions. We predicted that this phenomenon is related to ion heating in the IMS-MS instrument. In this study, to reveal the effect of ion heating from parameters in the IMS-MS instrument on the conformational population of the H2A-H2B dimer, we investigated the arrival time distributions of the H2A-H2B dimer by changing values of three instrumental parameters, namely, cone voltage located in the first vacuum chamber, trap collision energy (trap CE) for tandem mass spectrometry, and trap bias voltage for the entrance of IMS. These results revealed that the two populations observed for the H2A-H2B dimer were due to the trap bias voltage. Furthermore, to evaluate the internal energies of the analyte ions with respect to each parameter, benzylpyridinium derivatives were used as temperature-sensitive probes. The results showed that the trap CE voltage imparts greater internal energy to the ions than the trap bias voltage. In addition, this slight change in the internal energy caused by the trap bias voltage resulted in the structural diversity of the H2A-H2B dimer. Therefore, the trap bias voltage should be set with attention to the properties of the analytes, even if the effect of the trap bias voltage on the internal energy is negligible.

Elucidating Structures of Protein Complexes by Collision-Induced Dissociation at Elevated Gas Pressures

Ion activation methods carried out at gas pressures compatible with ion mobility separations are not yet widely established. This limits the analytical utility of emerging tandem-ion mobility spectrometers that conduct multiple ion mobility separations in series. The present work investigates the applicability of collision-induced dissociation (CID) at 1 to 3 mbar in a tandem-trapped ion mobility spectrometer (tandem-TIMS) to study the architecture of protein complexes. We show that CID of the homotetrameric protein complexes streptavidin (53 kDa), neutravidin (60 kDa), and concanavalin A (110 kDa) provides access to all subunits of the investigated protein complexes, including structurally informative dimers. We report on an "atypical" dissociation pathway, which for concanavalin A proceeds via symmetric partitioning of the precursor charges and produces dimers with the same charge states that were previously reported from surface induced dissociation. Our data suggest a correlation between the formation of subunits by CID in tandem-TIMS/MS, their binding strengths in the native tetramer structures, and the applied activation voltage. Ion mobility spectra of in situ-generated subunits reveal a marked structural heterogeneity inconsistent with annealing into their most stable gas phase structures. Structural transitions are observed for in situ-generated subunits that resemble the transitions reported from collision-induced unfolding of natively folded proteins. These observations indicate that some aspects of the native precursor structure is preserved in the subunits generated from disassembly of the precursor complex. We rationalize our observations by an approximately 100-fold shorter activation time scale in comparison to traditional CID in a collision cell. Finally, the approach discussed here to conduct CID at elevated pressures appears generally applicable also for other types of tandem-ion mobility spectrometers.

Oxidized and Reduced Dimeric Protein Complexes Illustrate Contrasting CID and SID Charge Partitioning

Charge partitioning during the dissociation of protein complexes in the gas phase is influenced by many factors, such as interfacial interactions, protein flexibility, protein conformation, and dissociation methods. In the present work, two cysteine-containing homodimer proteins, β-lactoglobulin and α-lactalbumin, with the disulfide bonds intact and reduced, were used to gain insight into the charge partitioning behaviors of collision-induced dissociation (CID) and surface-induced dissociation (SID) processes. For these proteins, we find that restructuring dominates with CID and dissociation with symmetric charge partitioning dominates with SID, regardless of whether intramolecular disulfide bonds are oxidized or reduced. CID of the charge-reduced dimeric protein complex leads to a precursor with a slightly smaller collision cross section (CCS), greater stability, and more symmetrically distributed charges than the significantly expanded form produced by CID of the higher charged dimer. Collision-induced unfolding plots demonstrate that the unfolding-restructuring of the protein complexes initiates the charge migration of higher charge-state precursors. Overall, gas collisions reveal the charge-dependent restructuring/unfolding properties of the protein precursor, while surface collisions lead predominantly to more charge-symmetric monomer separation. CID's multiple low-energy collisions sequentially reorganize intra- and intermolecular bonds, while SID's large-step energy jump cleaves intermolecular interfacial bonds in preference to reorganizing intramolecular bonds. The activated population of precursors that have taken on energy without dissociating (populated in CID over a wide range of collision energies, populated in SID for only a narrow distribution of collision energies near the onset of dissociation) is expected to be restructured, regardless of the activation method.

Top-down mass spectrometry of native proteoforms and their complexes: A community study

The combination of native electrospray ionisation with top-down fragmentation in mass spectrometry allows simultaneous determination of the stoichiometry of noncovalent complexes and identification of their component proteoforms and co-factors. While this approach is powerful, both native mass spectrometry and top-down mass spectrometry are not yet well standardised, and only a limited number of laboratories regularly carry out this type of research. To address this challenge, the Consortium for Top-Down Proteomics (CTDP) initiated a study to develop and test protocols for native mass spectrometry combined with top-down fragmentation of proteins and protein complexes across eleven instruments in nine laboratories. The outcomes are summarised in this report to provide robust benchmarks and a valuable entry point for the scientific community.

Structure Relaxation Approximation (SRA) for Elucidation of Protein Structures from Ion Mobility Measurements (II). Protein Complexes

Characterizing structures of protein complexes and their disease-related aberrations is essential to understanding molecular mechanisms of many biological processes. Electrospray ionization coupled with hybrid ion mobility/mass spectrometry (ESI-IM/MS) methods offer sufficient sensitivity, sample throughput, and dynamic range to enable systematic structural characterization of proteomes. However, because ESI-IM/MS characterizes ionized protein systems in the gas phase, it generally remains unclear to what extent the protein ions characterized by IM/MS have retained their solution structures. Here, we discuss the first application of our computational structure relaxation approximation [Bleiholder, C.; et al. J. Phys. Chem. B 2019, 123 (13), 2756-2769] to assign structures of protein complexes in the range from ∼16 to ∼60 kDa from their "native" IM/MS spectra. Our analysis shows that the computed IM/MS spectra agree with the experimental spectra within the errors of the methods. The structure relaxation approximation (SRA) indicates that native backbone contacts appear largely retained in the absence of solvent for the investigated protein complexes and charge states. Native contacts between polypeptide chains of the protein complex appear to be retained to a comparable extent as contacts within a folded polypeptide chain. Our computations also indicate that the hallmark "compaction" often observed for protein systems in native IM/MS measurements appears to be a poor indicator of the extent to which native residue-residue interactions are lost in the absence of solvent. Further, the SRA indicates that structural reorganization of the protein systems in IM/MS measurements appears driven largely by remodeling of the protein surface that increases its hydrophobic content by approximately 10%. For the systems studied here, this remodeling of the protein surface appears to occur mainly by structural reorganization of surface-associated hydrophilic amino acid residues not associated with β-strand secondary structure elements. Properties related to the internal protein structure, as assessed by void volume or packing density, appear unaffected by remodeling of the surface. Taken together, the structural reorganization of the protein surface appears to be generic in nature and to sufficiently stabilize protein structures to render them metastable on the time scale of IM/MS measurements.

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.

Ion Mobility Mass Spectrometry (IM-MS) for Structural Biology: Insights Gained by Measuring Mass, Charge, and Collision Cross Section

The investigation of macromolecular biomolecules with ion mobility mass spectrometry (IM-MS) techniques has provided substantial insights into the field of structural biology over the past two decades. An IM-MS workflow applied to a given target analyte provides mass, charge, and conformation, and all three of these can be used to discern structural information. While mass and charge are determined in mass spectrometry (MS), it is the addition of ion mobility that enables the separation of isomeric and isobaric ions and the direct elucidation of conformation, which has reaped huge benefits for structural biology. In this review, where we focus on the analysis of proteins and their complexes, we outline the typical features of an IM-MS experiment from the preparation of samples, the creation of ions, and their separation in different mobility and mass spectrometers. We describe the interpretation of ion mobility data in terms of protein conformation and how the data can be compared with data from other sources with the use of computational tools. The benefit of coupling mobility analysis to activation via collisions with gas or surfaces or photons photoactivation is detailed with reference to recent examples. And finally, we focus on insights afforded by IM-MS experiments when applied to the study of conformationally dynamic and intrinsically disordered proteins.

Simulation toolkits at the molecular scale for trans-scale thermal signaling

Thermogenesis is a physiological activity of releasing heat that originates from intracellular biochemical reactions. Recent experimental studies discovered that externally applied heat changes intracellular signaling locally, resulting in global changes in cell morphology and signaling. Therefore, we hypothesize an inevitable contribution of thermogenesis in modulating biological system functions throughout the spatial scales from molecules to individual organisms. One key issue examining the hypothesis, namely, the "trans-scale thermal signaling," resides at the molecular scale on the amount of heat released via individual reactions and by which mechanism the heat is employed for cellular function operations. This review introduces atomistic simulation tool kits for studying the mechanisms of thermal signaling processes at the molecular scale that even state-of-the-art experimental methodologies of today are hardly accessible. We consider biological processes and biomolecules as potential heat sources in cells, such as ATP/GTP hydrolysis and multiple biopolymer complex formation and disassembly. Microscopic heat release could be related to mesoscopic processes via thermal conductivity and thermal conductance. Additionally, theoretical simulations to estimate these thermal properties in biological membranes and proteins are introduced. Finally, we envisage the future direction of this research field.

Gas-Phase Unfolding of Protein Complexes Distinguishes Conformational Isomers

Proteins can adopt different conformational states that are important for their biological function and, in some cases, can be responsible for their dysfunction. The essential roles that proteins play in biological systems make distinguishing the structural differences between these conformational states both fundamentally and practically important. Here, we demonstrate that collision-induced unfolding (CIU), in combination with ion mobility-mass spectrometry (IM-MS) measurements, distinguish subtly different conformational states for protein complexes. Using the open and closed states of the β-lactoglobulin (βLG) dimer as a model, we show that these two conformational isomers unfold during collisional activation to generate distinct states that are readily separated by IM-MS. Extensive molecular modeling of the CIU process reproduces the distinct unfolding intermediates and identifies the molecular details that explain why the two conformational states unfold in distinct ways. Strikingly, the open conformational state forms new electrostatic interactions upon collisional heating, while the closed state does not. These newly formed electrostatic interactions involve residues on the loop differentially positioned in the two βLG conformational isomers, highlighting that gas-phase unfolding pathways reflect aspects of solution structure. This combination of experiment and theory provides a path forward for distinguishing subtly different conformational isomers for protein complexes via gas-phase unfolding experiments. Our results also have implications for understanding how protein complexes dissociate in the gas phase, indicating that current models need to be refined to explain protein complex dissociation.

Native Mass Spectrometry and Collision-Induced Unfolding of Laser-Ablated Proteins

Infrared laser ablation sample transfer (LAST) was used to collect samples from solid surfaces for mass spectrometry under native spray conditions. Native mass spectrometry was utilized to probe the charge states and collision-induced unfolding (CIU) characteristics of bovine serum albumin (BSA), bovine hemoglobin (BHb), and jack-bean concanavalin A (ConA) via direct injection electrospray, after liquid extraction surface sampling, and after LAST. Each protein was deposited from solution on solid surfaces and laser-ablated for off-line analysis or sampled for online analysis. It was found that the protein ion gas-phase charge-state distributions were comparable for direct infusion, liquid extraction, and laser ablation experiments. Moreover, calculated average collision cross section (CCS) values from direct injection, liquid extraction, and laser ablation experiments were consistent with previously reported literature values. Additionally, an equivalent number of mobility features and conformational turnovers were identified from unfolding pathways from all three methods for all charge states of each protein analyzed in this work. The presented work suggests that laser ablation yields intact proteins (BSA, BHb, and ConA), is compatible with native mass spectrometry, and could be suitable for spatially resolved interrogation of unfolding pathways of proteins.

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.

Application of Multiple Length Cross-linkers to the Characterization of Gaseous Protein Structure

The speed, sensitivity, and tolerance of heterogeneity, as well as the kinetic trapping of solution-like states during electrospray, make native mass spectrometry an attractive method to study protein structure. Increases in the resolution of ion mobility measurements and in mass resolving power and range are leading to the increase of the information content of intact protein measurements and an expanded role of mass spectrometry in structural biology. Herein, a suite of different length noncovalent (sulfonate to positively charged side chain) cross-linkers was introduced via gas-phase ion/ion chemistry and used to determine distance restraints of kinetically trapped gas-phase structures of native-like cytochrome c ions. Electron capture dissociation allowed for the identification of cross-linked sites. Different length linkers resulted in distinct pairs of side chains being linked, supporting the ability of gas-phase cross-linking to be structurally specific. The gas-phase lengths of the cross-linkers were determined by conformational searches and density functional theory, allowing for the interpretation of the cross-links as distance restraints. These distance restraints were used to model gas-phase structures with molecular dynamics simulations, revealing a mixture of structures with similar overall shape/size but distinct features, thereby illustrating the kinetic trapping of multiple native-like solution structures in the gas phase.

Step toward Probing the Nonannular Belt of Membrane Proteins

Integral membrane proteins are embedded in the biological membrane, where they carry out numerous biological processes. Although lipids present in the membrane are crucial for membrane protein function, it remains difficult to characterize many lipid binding events to membrane proteins, such as those that form the annular belt. Here, we use native mass spectrometry along with the charge-reducing properties of trimethylamine N-oxide (TMAO) to characterize a large number of lipid binding events to the bacterial ammonia channel (AmtB). In the absence of TMAO, significant peak overlap between neighboring charge states is observed, resulting in erroneous abundances for different molecular species. With the addition of TMAO, the weighted average charge state (Z_avg) was decreased. In addition, the increased spacing between nearby charge states enabled a higher number of lipid binding species to be observed while minimizing mass spectral peak overlap. These conditions helped us to determine the equilibrium binding constants (K_d) for up to 16 lipid binding events. The binding constants for the first few lipid binding events display the highest affinity, whereas the binding constants for higher lipid binding events converge to a similar value. These findings suggest a transition from nonannular to annular lipid binding to AmtB.

Mass Spectrometry and Machine Learning Reveal Determinants of Client Recognition by Antiamyloid Chaperones

The assembly of proteins and peptides into amyloid fibrils is causally linked to serious disorders such as Alzheimer's disease. Multiple proteins have been shown to prevent amyloid formation in vitro and in vivo, ranging from highly specific chaperone-client pairs to completely nonspecific binding of aggregation-prone peptides. The underlying interactions remain elusive. Here, we turn to the machine learning-based structure prediction algorithm AlphaFold2 to obtain models for the nonspecific interactions of β-lactoglobulin, transthyretin, or thioredoxin 80 with the model amyloid peptide amyloid β and the highly specific complex between the BRICHOS chaperone domain of C-terminal region of lung surfactant protein C and its polyvaline target. Using a combination of native mass spectrometry (MS) and ion mobility MS, we show that nonspecific chaperoning is driven predominantly by hydrophobic interactions of amyloid β with hydrophobic surfaces in β-lactoglobulin, transthyretin, and thioredoxin 80, and in part regulated by oligomer stability. For C-terminal region of lung surfactant protein C, native MS and hydrogen-deuterium exchange MS reveal that a disordered region recognizes the polyvaline target by forming a complementary β-strand. Hence, we show that AlphaFold2 and MS can yield atomistic models of hard-to-capture protein interactions that reveal different chaperoning mechanisms based on separate ligand properties and may provide possible clues for specific therapeutic intervention.

Simulation of Energy-Resolved Mass Spectrometry Distributions from Surface-Induced Dissociation

Understanding the relationship between protein structure and experimental data is crucial for utilizing experiments to solve biochemical problems and optimizing the use of sparse experimental data for structural interpretation. Tandem mass spectrometry (MS/MS) can be used with a variety of methods to collect structural data for proteins. One example is surface-induced dissociation (SID), which is used to break apart protein complexes (via a surface collision) into intact subcomplexes and can be performed at multiple laboratory frame SID collision energies. These energy-resolved MS/MS experiments have shown that the profile of the breakages depends on the acceleration energy of the collision. It is possible to extract an appearance energy (AE) from energy-resolved mass spectrometry (ERMS) data, which shows the relative intensity of each type of subcomplex as a function of SID acceleration energy. We previously determined that these AE values for specific interfaces correlated with structural features related to interface strength. In this study, we further examined the structural relationships by developing a method to predict the full ERMS plot from the structure, rather than extracting a single value. First, we noted that for proteins with multiple interface types, we could reproduce the correct shapes of breakdown curves, further confirming previous structural hypotheses. Next, we demonstrated that interface size and energy density (measured using Rosetta) correlated with data derived from the ERMS plot (R2 = 0.71). Furthermore, based on this trend, we used native crystal structures to predict ERMS. The majority of predictions resulted in good agreement, and the average root-mean-square error was 0.20 for the 20 complexes in our data set. We also show that if additional information on cleavage as a function of collision energy could be obtained, the accuracy of predictions improved further. Finally, we demonstrated that ERMS prediction results were better for the native than for inaccurate models in 17/20 cases. An application to run this simulation has been developed in Rosetta, which is freely available for use.

Surface Activity of Amines Provides Evidence for the Combined ESI Mechanism of Charge Reduction for Protein Complexes

Charge reduction reactions are important for native mass spectrometry (nMS) because lower charge states help retain native-like conformations and preserve noncovalent interactions of protein complexes. While mechanisms of charge reduction reactions are not well understood, they are generally achieved through the addition of small molecules, such as polyamines, to traditional nMS buffers. Here, we present new evidence that surface-active, charge reducing reagents carry away excess charge from the droplet after being emitted due to Coulombic repulsion, thereby reducing the overall charge of the droplet. Furthermore, these processes are directly linked to two mechanisms for electrospray ionization, specifically the charge residue and ion evaporation models (CRM and IEM). Selected protein complexes were analyzed in solutions containing ammonium acetate and selected trialkylamines or diaminoalkanes of increasing alkyl chain lengths. Results show that amines with higher surface activity have increased propensities for promoting charge reduction of the protein ions. The electrospray ionization (ESI) emitter potential was also found to be a major contributing parameter to the prevalence of charge reduction; higher emitter potentials consistently coincided with lower average charge states among all protein complexes analyzed. These results offer experimental evidence for the mechanism of charge reduction in ESI and also provide insight into the final stages of the ESI and their impact on biological ions.

Cold Denaturation of Proteins in the Absence of Solvent: Implications for Protein Storage

The effect of temperature on the stability of proteins is well explored above 298 K, but harder to track experimentally below 273 K. Variable-temperature ion mobility mass spectrometry (VT IM-MS) allows us to measure the structure of molecules at sub-ambient temperatures. Here we monitor conformational changes that occur to two isotypes of monoclonal antibodies (mAbs) on cooling by measuring their collision cross sections (CCS) at discrete drift gas temperatures from 295 to 160 K. The CCS at 250 K is larger than predicted from collisional theory and experimental data at 295 K. This restructure is attributed to change in the strength of stabilizing intermolecular interactions. Below 250 K the CCS of the mAbs increases in line with prediction implying no rearrangement. Comparing data from isotypes suggest disulfide bridging influences thermal structural rearrangement. These findings indicate that in vacuo deep-freezing minimizes denaturation and maintains the native fold and VT IM-MS measurements at sub ambient temperatures provide new insights to the phenomenon of cold denaturation.

Cold Denaturation of Proteins in the Absence of Solvent: Implications for Protein Storage

The effect of temperature on the stability of proteins is well explored above 298 K, but harder to track experimentally below 273 K. Variable-temperature ion mobility mass spectrometry (VT IM-MS) allows us to measure the structure of molecules at sub-ambient temperatures. Here we monitor conformational changes that occur to two isotypes of monoclonal antibodies (mAbs) on cooling by measuring their collision cross sections (CCS) at discrete drift gas temperatures from 295 to 160 K. The CCS at 250 K is larger than predicted from collisional theory and experimental data at 295 K. This restructure is attributed to change in the strength of stabilizing intermolecular interactions. Below 250 K the CCS of the mAbs increases in line with prediction implying no rearrangement. Comparing data from isotypes suggest disulfide bridging influences thermal structural rearrangement. These findings indicate that in vacuo deep-freezing minimizes denaturation and maintains the native fold and VT IM-MS measurements at sub ambient temperatures provide new insights to the phenomenon of cold denaturation.

Native Mass Spectrometry: Recent Progress and Remaining Challenges

Native mass spectrometry (nMS) has emerged as an important tool in studying the structure and function of macromolecules and their complexes in the gas phase. In this review, we cover recent advances in nMS and related techniques including sample preparation, instrumentation, activation methods, and data analysis software. These advances have enabled nMS-based techniques to address a variety of challenging questions in structural biology. The second half of this review highlights recent applications of these technologies and surveys the classes of complexes that can be studied with nMS. Complementarity of nMS to existing structural biology techniques and current challenges in nMS are also addressed.

Approaches to Heterogeneity in Native Mass Spectrometry

Native mass spectrometry (MS) is aimed at preserving and determining the native structure, composition, and stoichiometry of biomolecules and their complexes from solution after they are transferred into the gas phase. Major improvements in native MS instrumentation and experimental methods over the past few decades have led to a concomitant increase in the complexity and heterogeneity of samples that can be analyzed, including protein-ligand complexes, protein complexes with multiple coexisting stoichiometries, and membrane protein-lipid assemblies. Heterogeneous features of these biomolecular samples can be important for understanding structure and function. However, sample heterogeneity can make assignment of ion mass, charge, composition, and structure very challenging due to the overlap of tens or even hundreds of peaks in the mass spectrum. In this review, we cover data analysis, experimental, and instrumental advances and strategies aimed at solving this problem, with an in-depth discussion of theoretical and practical aspects of the use of available deconvolution algorithms and tools. We also reflect upon current challenges and provide a view of the future of this exciting field.

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.

Electrospray ionization of native membrane proteins proceeds via a charge equilibration step

Electrospray ionization mass spectrometry is increasingly applied to study the structures and interactions of membrane protein complexes. However, the charging mechanism is complicated by the presence of detergent micelles during ionization. Here, we show that the final charge of membrane proteins can be predicted by their molecular weight when released from the non-charge reducing saccharide detergents. Our data indicate that PEG detergents lower the charge depending on the number of detergent molecules in the surrounding micelle, whereas fos-choline detergents may additionally participate in ion–ion reactions after desolvation. The supercharging reagent sulfolane, on the other hand, has no discernible effect on the charge of detergent-free membrane proteins. Taking our observations into the context of protein-detergent interactions in the gas phase, we propose a charge equilibration model for the generation of native-like membrane protein ions. During ionization of the protein-detergent complex, the ESI charges are distributed between detergent and protein according to proton affinity of the detergent, number of detergent molecules, and surface area of the protein. Charge equilibration influenced by detergents determines the final charge state of membrane proteins. This process likely contributes to maintaining a native-like fold after detergent release and can be harnessed to stabilize particularly labile membrane protein complexes in the gas phase.

The electrospray ionization mechanism contributes to preserving the structures and interactions of membrane protein complexes in native mass spectrometry.

Investigation of Charge-State-Dependent Compaction of Protein Ions with Native Ion Mobility-Mass Spectrometry and Theory

The precise relationship between native gas-phase protein ion structure, charge, desolvation, and activation remains elusive. Much evidence supports the Charge Residue Model for native protein ions formed by electrospray ionization, but scaling laws derived from it relate only to overall ion size. Closer examination of drift tube CCSs across individual native protein ion charge state distributions (CSDs) reveals deviations from global trends. To investigate whether this is due to structure variation across CSDs or contributions of long-range charge-dipole interactions, we performed in vacuo force field molecular dynamics (MD) simulations of multiple charge conformers of three proteins representing a variety of physical and structural features: β-lactoglobulin, concanavalin A, and glutamate dehydrogenase. Results from these simulated ions indicate subtle structure variation across their native CSDs, although effects of these structural differences and long-range charge-dependent interactions on CCS are small. The structure and CCS of smaller proteins may be more sensitive to charge due to their low surface-to-volume ratios and reduced capacity to compact. Secondary and higher order structure from condensed-phase structures is largely retained in these simulations, supporting the use of the term "native-like" to describe results from native ion mobility-mass spectrometry experiments, although, notably, the most compact structure can be the most different from the condensed-phase structure. Collapse of surface side chains to self-solvate through formation of new hydrogen bonds is a major feature of gas-phase compaction and likely occurs during the desolvation process. Results from these MD simulations provide new insight into the relationship of gas-phase protein ion structure, charge, and CCS.

Surface-Induced Dissociation for Protein Complex Characterization

Native mass spectrometry (nMS) enables intact non-covalent complexes to be studied in the gas phase. nMS can provide information on composition, stoichiometry, topology, and, when coupled with surface-induced dissociation (SID), subunit connectivity. Here we describe the characterization of protein complexes by nMS and SID. Substructural information obtained using this method is consistent with the solved complex structure, when a structure exists. This provides confidence that the method can also be used to obtain substructural information for unknowns, providing insight into subunit connectivity and arrangements. High-energy SID can also provide information on proteoforms present. Previously SID has been limited to a few in-house modified instruments and here we focus on SID implemented within an in-house-modified Q Exactive UHMR. However, SID is currently commercially available within the Waters Select Series Cyclic IMS instrument. Projects are underway that involve the NIH-funded native MS resource (nativems.osu.edu), instrument vendors, and third-party vendors, with the hope of bringing the technology to more platforms and labs in the near future. Currently, nMS resource staff can perform SID experiments for interested research groups.

Suppression of Protein Structural Perturbations in Native Electrospray Ionization during the Final Evaporation Stages Revealed by Molecular Dynamics Simulations

Native electrospray ionization was known to preserve the protein structure in solution, which overcame the uncontrollable acidification of droplets during transfer from solution into the gas phase in conventional electrospray ionization. However, detailed experimental studies on when and how could native electrospray ionization minimize structural perturbations remain quite unclear. Herein, we conducted molecular dynamics simulations to investigate the protein structure evolution during electrospray ionization. At a neutral droplet pH, the protein structure in solution could be retained after evaporation, which was in accordance with previous reports. As the droplet pH deviated from neutral, we have found that the compact protein structure would not unfold until the last 10 ns prior to the final desolvation, which demonstrated that the role of native electrospray ionization in preserving the protein structure was mainly reflected on the final evaporation stages. The present study might provide new insights into studying the microscopic biomolecular events occurring during the liquid-gas interface transition and their influence on solution-structure retention.

Influence of Primary Structure on Fragmentation of Native-Like Proteins by Ultraviolet Photodissociation

Analysis of native-like protein structures in the gas phase via native mass spectrometry and auxiliary techniques has become a powerful tool for structural biology applications. In combination with ultraviolet photodissociation (UVPD), native top-down mass spectrometry informs backbone flexibility, topology, hydrogen bonding networks, and conformational changes in protein structure. Although it is known that the primary structure affects dissociation of peptides and proteins in the gas phase, its effect on the types and locations of backbone cleavages promoted by UVPD and concomitant influence on structural characterization of native-like proteins is not well understood. Here, trends in the fragmentation of native-like proteins were evaluated by tracking the propensity of 10 fragment types (a, a+1, b, c, x, x+1, y, y-1, Y, and z) in relation to primary structure in a native-top down UVPD data set encompassing >9600 fragment ions. Differing fragmentation trends are reported for the production of distinct fragment types, attributed to a combination of both direct dissociation pathways from excited electronic states and those surmised to involve intramolecular vibrational energy redistribution after internal conversion. The latter pathways were systematically evaluated to evince the role of proton mobility in the generation of "CID-like" fragments through UVPD, providing pertinent insight into the characterization of native-like proteins. Fragmentation trends presented here are envisioned to enhance analysis of the protein higher-order structure or augment scoring algorithms in the high-throughput analysis of intact proteins.

Rapid structural discrimination of IgG antibodies by multicharge-state collision-induced unfolding

Immunoglobulin G (IgG) antibodies are an important class of biotherapeutics that target various diseases, such as cancers, neurodegenerative disorders, and autoimmune diseases, yet rapid discrimination of IgG antibodies remains a great challenge due to heterogeneity, flexibility, and large size. Herein, we demonstrate a simplified multicharge-state collision-induced unfolding (CIU) method for rapid differentiation of four IgG isotypes that differ in terms of the numbers and patterns of disulfide bonds, bypassing tedious single charge-state selection in advance. The results presented herein reveal that gas-phase unfolding behaviors have a strong dependence on charge states outside IgG surfaces; therefore, multicharge-state CIU analysis of IgG subtypes could offer a great opportunity to gain deeper insights into their gas-phase structural differentiation. The full discrimination of IgG antibody isoforms that possess different disulfide bond numbers and even subtle disulfide bonding patterns can be achieved based on their charge-dependent gas-phase unfolding behaviors and root-mean square deviation in CIU difference spectra. Taken together, the incorporation of all charge states observed in a native ion mobility-mass spectrometry (IM-MS) experiment to CIU analysis could make this strategy sensitive to more subtle structural discrepancies, facilitating the rapid discrimination and evaluation of innovative structurally similar biotherapeutic candidates with unexplored functions.

Structural Characterization of Human Histone H4.1 by Tandem Nonlinear and Linear Ion Mobility Spectrometry Complemented with Molecular Dynamics Simulations

Extracellular histone H4 is an attractive drug target owing to its roles in organ failure in sepsis and other diseases. To identify inhibitors using in silico methods, information on histone H4 structural dynamics and three-dimensional (3D) structural coordinates is required. Here, DNA-free histone H4 type 1 (H4.1) was characterized by utilizing tandem nonlinear and linear ion mobility spectrometry (FAIMS-TIMS) coupled to mass spectrometry (MS) complemented with molecular dynamics (MD) simulations. The gas-phase structures of H4.1 are dependent on the starting solution conditions, evidenced by differences in charge state distributions, mobility distributions, and collision-induced unfolding (CIU) pathways. The experimental results show that H4.1 adopts diverse conformational types from compact (C) to partially folded (P) and subsequently elongated (E) structures. Molecular dynamics simulations provided candidate structures for the histone H4.1 monomer in solution and for the gas-phase structures observed using FAIMS-IMS-TOF MS as a function of the charge state and mobility distribution. A combination of the FAIMS-TIMS experimental results with theoretical dipole calculations reveals the important role of charge distribution in the dipole alignment of H4.1 elongated structures at high electric fields. A comparison of the secondary and primary structures of DNA-free H2A.1 and H4.1 is made based on the experimental IMS-MS and MD findings.

Surface-Induced Dissociation of Anionic vs Cationic Native-Like Protein Complexes

Characterizing protein-protein interactions, stoichiometries, and subunit connectivity is key to understanding how subunits assemble into biologically relevant, multisubunit protein complexes. Native mass spectrometry (nMS) has emerged as a powerful tool to study protein complexes due to its low sample consumption and tolerance for heterogeneity. In nMS, positive mode ionization is routinely used and charge reduction, through the addition of solution additives, is often used, as the resulting lower charge states are often considered more native-like. When fragmented by surface-induced dissociation (SID), charge reduced complexes often give increased structural information over their "normal-charged" counterparts. A disadvantage of solution phase charge reduction is that increased adduction, and hence peak broadening, is often observed. Previous studies have shown that protein complexes ionized using negative mode generally form lower charge states relative to positive mode. Here we demonstrate that the lower charged protein complex anions activated by surface collisions fragment in a manner consistent with their solved structures, hence providing substructural information. Negative mode ionization in ammonium acetate offers the advantage of charge reduction without the peak broadening associated with solution phase charge reduction additives and provides direct structural information when coupled with SID. SID of 20S human proteasome (a 28-mer comprised of four stacked heptamer rings in an αββα formation), for example, provides information on both substructure (e.g., splitting into a 7α ring and the corresponding ββα 21-mer, and into α dimers and trimers to provide connectivity around the 7 α ring) and proteoform information on monomers.

Structural Heterogeneity of Human Histone H2A.1

Histones are highly basic chromatin proteins that tightly package and order eukaryotic DNA into nucleosomes. While the atomic structure of the nucleosomes has been determined, the three-dimensional structure of DNA-free histones remains unresolved. Here, we combine tandem nonlinear and linear ion mobility spectrometry (FAIMS-TIMS) coupled to mass spectrometry in parallel with molecular modeling to study the conformational space of a DNA-free histone H2A type 1 (H2A.1). Experimental results showed the dependence of the gas-phase structures on the starting solution conditions, characterized by charge state distributions, mobility distributions, and collision-induced-unfolding pathways. The measured H2A.1 gas-phase structures showed a high diversity of structural features ranging from compact (C) to partially folded (P) and then highly elongated (E) conformations. Molecular dynamics simulations provided candidate structures for the solution H2A.1 native conformation with folded N- and C-terminal tails, as well as gas-phase candidate structures associated with the mobility trends. Complementary collision cross section and dipole calculations showed that the charge distribution in the case of elongated gas-phase structures, where basic and acidic residues are mostly exposed (e.g., z > 15+), is sufficient to induce differences in the dipole alignment at high electric fields, in good agreement with the trends observed during the FAIMS-TIMS experiments.

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.

Reaction Pathway Sampling and Free-Energy Analyses for Multimeric Protein Complex Disassembly by Employing Hybrid Configuration Bias Monte Carlo/Molecular Dynamics Simulation

Physicochemical characterization of multimeric biomacromolecule assembly and disassembly processes is a milestone to understand the mechanisms for biological phenomena at the molecular level. Mass spectroscopy (MS) and structural bioinformatics (SB) approaches have become feasible to identify subcomplexes involved in assembly and disassembly, while they cannot provide atomic information sufficient for free-energy calculation to characterize transition mechanism between two different sets of subcomplexes. To combine observations derived from MS and SB approaches with conventional free-energy calculation protocols, we here designed a new reaction pathway sampling method by employing hybrid configuration bias Monte Carlo/molecular dynamics (hcbMC/MD) scheme and applied it to simulate the disassembly process of serum amyloid P component (SAP) pentamer. The results we obtained are consistent with those of the earlier MS and SB studies with respect to SAP subcomplex species and the initial stage of SAP disassembly processes. Furthermore, we observed a novel dissociation event, ring-opening reaction of SAP pentamer. Employing free-energy calculation combined with the hcbMC/MD reaction pathway trajectories, we moreover obtained experimentally testable observations on (1) reaction time of the ring-opening reaction and (2) importance of Asp42 and Lys117 for stable formation of SAP oligomer.

Sulfolane-Induced Supercharging of Electrosprayed Salt Clusters: An Experimental/Computational Perspective

It is well-known that supercharging agents (SCAs) such as sulfolane enhance the electrospray ionization (ESI) charge states of proteins, although the mechanistic origins of this effect remain contentious. Only very few studies have explored SCA effects on analytes other than proteins or peptides. This work examines how sulfolane affects electrosprayed NaI salt clusters. Such alkali metal halide clusters have played a key role for earlier ESI mechanistic studies, making them interesting targets for supercharging investigations. ESI of aqueous NaI solutions predominantly generated singly charged [Na_nI_(n-1)]⁺ clusters. The addition of sulfolane resulted in abundant doubly charged [Na_nI_(n-2)Sulfolane_s]²⁺ species. These experimental data for the first time demonstrate that electrosprayed salt clusters can undergo supercharging. Molecular dynamics (MD) simulations of aqueous ESI nanodroplets containing Na⁺/I^- with and without sulfolane were conducted to obtain atomistic insights into the supercharging mechanism. The simulations produced [Na_nI_i]^z+ and [Na_nI_iSulfolane_s]^z+ clusters similar to those observed experimentally. The MD trajectories demonstrated that these clusters were released into the gas phase upon droplet evaporation to dryness, in line with the charged residue model. Sulfolane was found to evaporate much more slowly than water. This slow evaporation, in conjunction with the large dipole moment of sulfolane, resulted in electrostatic stabilization of the shrinking ESI droplets and the final clusters. Hence, charge-dipole stabilization causes the sulfolane-containing droplets and clusters to retain more charge, thereby providing the mechanistic foundation of salt cluster supercharging.

Interrogating the Quaternary Structure of Noncanonical Hemoglobin Complexes by Electrospray Mass Spectrometry and Collision-Induced Dissociation

Various activation methods are available for the fragmentation of gaseous protein complexes produced by electrospray ionization (ESI). Such experiments can potentially yield insights into quaternary structure. Collision-induced dissociation (CID) is the most widely used fragmentation technique. Unfortunately, CID of protein complexes is dominated by the ejection of highly charged monomers, a process that does not yield any structural insights. Using hemoglobin (Hb) as a model system, this work examines under what conditions CID generates structurally informative subcomplexes. Native ESI mainly produced tetrameric Hb ions. In addition, "noncanonical" hexameric and octameric complexes were observed. CID of all these species [(αβ)₂, (αβ)₃, and (αβ)₄] predominantly generated highly charged monomers. In addition, we observed hexamer → tetramer + dimer dissociation, implying that hexamers have a tetramer··dimer architecture. Similarly, the observation of octamer → two tetramer dissociation revealed that octamers have a tetramer··tetramer composition. Gas-phase candidate structures of Hb assemblies were produced by molecular dynamics (MD) simulations. Ion mobility spectrometry was used to identify the most likely candidates. Our data reveal that the capability of CID to produce structurally informative subcomplexes depends on the fate of protein-protein interfaces after transfer into the gas phase. Collapse of low affinity interfaces conjoins the corresponding subunits and favors CID via monomer ejection. Structurally informative subcomplexes are formed only if low affinity interfaces do not undergo a major collapse. However, even in these favorable cases CID is still dominated by monomer ejection, requiring careful analysis of the experimental data for the identification of structurally informative subcomplexes.

Digital ion trap mass analysis of high mass protein complexes using IR activation coupled with ion/ion reactions

Native mass spectrometry (MS) focuses on measuring the masses of large biomolecular complexes and probing their structures. Large biomolecular complexes are readily introduced into mass spectrometers as gas-phase ions using electrospray ionization (ESI); however, the ions tend to be heavily adducted with solvent and salts, which leads to mass measurement errors. Various solution clean-up approaches can reduce the degree of adduction prior to introduction to the mass spectrometer. Gas-phase activation of trapped ions can provide additional adduct reduction, and charge reduction ion/ion reactions increase charge state separation. Together, gas-phase activation and charge reduction can combine to yield spectra of well separated charge states for improved mass measurements. A simple gas-phase collisional activation technique is to apply a dipolar DC (DDC) field to opposing electrodes in an ion trap. DDC activation loses its efficacy when ions are trapped at low q values, which is true of the high m/z ions generated by charge reduction ion/ion reactions. Digital ion trapping (DIT) readily traps high m/z ions at higher q values by varying trapping frequency rather than amplitude, but the low frequencies used to trap high m/z ions also decreases the efficacy of DDC activation. We demonstrate here using ions derived from GroEL that IR activation of ions shows no discrimination against high m/z ions trapped with DIT, because they can be focused equally well to the trap center to interact with the IR laser beam. Following pump out of excess background gas, IR activation can also induce efficient dissociation of the GroEL complex. This work demonstrates that IR activation is an effective approach for ion heating in native MS over the unusually wide range of charge states accessible via gas-phase ion/ion reactions.