Molecular Dynamics Simulations on Gas-Phase Proteins with Mobile Protons: Inclusion of All-Atom Charge Solvation

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.

Molecular Dynamics Simulations of Native Protein Charging via Proton Transfer during Electrospray Ionization with Grotthuss Diffuse H3O. Anal Chem 2024;96:4146-4153. [PMID: 38427846 DOI: 10.1021/acs.analchem.3c05089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Unraveling the mechanism by which native proteins are charged through electrospray ionization (ESI) has been the focus of considerable research because observable charge states can be correlated to biophysical characteristics, such as protein folding and, thus, solution conformation. Difficulties in characterizing electrosprayed droplets have catalyzed the use of molecular dynamics (MD) to provide insights into the mechanisms by which proteins are charged and transferred to the gas phase. However, prior MD studies have utilized metal ions, primarily Na+, as charge carriers, even though proteins are primarily detected as protonated ions in the mass spectra. Here, we propose a modified MD protocol for simulating discrete Grotthuss diffuse H3O+ that is capable of dynamically altering amino-acid protonation states to model electrospray charging and gaseous ion formation of model proteins, ubiquitin, and myoglobin. Application of the protocol to the evaporation of acidic droplets enables a molecular perspective of H3O+ coordination and proton transfer to/from proteins, which is unfeasible with the metal charge carriers used in previous MD studies of ESI. Our protocol recreates experimentally observed charge-state distributions and supports the charge residue model (CRM) as the dominant mechanism of native protein ionization during ESI. Additionally, our results suggest that protonation is highly specific to individual residues and is correlated to the formation of localized hydrated regions on the protein surface as droplets desolvate. Considering the use of discrete H3O+ instead of Na+, the developed protocol is a necessary step toward developing a more comprehensive model of protein ionization during ESI.

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.

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.

Diserinol Isophthalamide: A Novel Reagent for Complexation with Biomolecular Anions in Electrospray Ionization Mass Spectrometry

Transferring biomolecules from solution to vacuum facilitates a detailed analysis of molecular structure and dynamics by isolating molecules of interest from a complex environment. However, inherent in the ion desolvation process is the loss of solvent hydrogen bonding partners, which are critical for the stability of a condensed-phase structure. Thus, transfer of ions to vacuum can favor structural rearrangement, especially near solvent-accessible charge sites, which tend to adopt intramolecular hydrogen bonding motifs in the absence of solvent. Complexation of monoalkylammonium moieties (e.g., lysine side chains) with crown ethers such as 18-crown-6 can disfavor structural rearrangement of protonated sites, but no equivalent ligand has been investigated for deprotonated groups. Herein we describe diserinol isophthalamide (DIP), a novel reagent for the gas-phase complexation of anionic moieties within biomolecules. Complexation is observed to the C-terminus or side chains of the small model peptides GD, GE, GG, DF-OMe, VYV, YGGFL, and EYMPME in electrospray ionization mass spectrometry (ESI-MS) studies. In addition, complexation is observed with the phosphate and carboxylate moieities of phosphoserine and phosphotyrosine. DIP performs favorably in comparison to an existing anion recognition reagent, 1,1'-(1,2-phenylene)bis(3-phenylurea), that exhibits moderate carboxylate binding in organic solvent. This improved performance in ESI-MS experiments is attributed to reduced steric constraints to complexation with carboxylate groups of larger molecules. Overall, diserinol isophthalamide is an effective complexation reagent that can be applied in future work to study retention of solution-phase structure, investigate intrinsic molecular properties, and examine solvation effects.

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.

Understanding of protomers/deprotomers by combining mass spectrometry and computation

Multifunctional compounds may form different prototropic isomers under different conditions, which are known as protomers/deprotomers. In biological systems, these protomer/deprotomer isomers affect the interaction modes and conformational landscape between compounds and enzymes and thus present different biological activities. Study on protomers/deprotomers is essentially the study on the acidity/basicity of each intramolecular functional group and its effect on molecular structure. In recent years, the combination of mass spectrometry (MS) and computational chemistry has been proven to be a powerful and effective means to study prototropic isomers. MS-based technologies are developed to discriminate and characterize protomers/deprotomers to provide structural information and monitor transformations, showing great superiority than other experimental methods. Computational chemistry is used to predict the thermodynamic stability of protomers/deprotomers, provide the simulated MS/MS spectra, infrared spectra, and calculate collision cross-section values. By comparing the theoretical data with the corresponding experimental results, the researchers can not only determine the protomer/deprotomer structure, but also investigate the structure-activity relationship in a given system. This review covers various MS methods and theoretical calculations and their devotion to isomer discrimination, structure identification, conformational transformation, and phase transition investigation of protomers/deprotomers.

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.

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.

•

Presented extensive molecular dynamics (MD) simulation data investigating protein vacuum exposure and rehydration dynamics.

•

Demonstrated that the majority of the protein structure recovers their initial solution conformation after vacuum exposure.

•

Explored the potential gain for structural biology of using MD simulation to refine gas-phase determined protein structures.

Grotthuss Molecular Dynamics Simulations for Modeling Proton Hopping in Electrosprayed Water Droplets

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H₃O⁺ + H₂O → H₂O + H₃O⁺ transitions. In accordance with ab initio MD data, these Grotthuss hopping events are executed in "bursts" with quasi-instantaneous hopping across one or more waters. The bursts are separated by regular MD periods during which H₃O⁺ ions undergo Brownian diffusion. The resulting proton diffusion coefficient agrees with the literature value. We apply this Grotthuss MD technique to highly charged water droplets that are in a size regime encountered during electrospray ionization (5 nm radius, ∼17,000 H₂O). The droplets undergo rapid solvent evaporation and occasional H₃O⁺ ejection, keeping them at ca. 81% of the Rayleigh limit. The simulated behavior is consistent with phase Doppler anemometry data. The Grotthuss MD technique developed here should be useful for modeling the behavior of various proton-containing systems that are too large for high-level computational approaches. In particular, we envision future applications related to electrospray processes, where earlier simulations used metal cations while in reality excess protons dominate.

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.

Protection by desiccation-tolerance proteins probed at the residue level

Extremotolerant organisms from all domains of life produce protective intrinsically disordered proteins (IDPs) in response to desiccation stress. In vitro, many of these IDPs protect enzymes from dehydration stress better than U.S. Food and Drug Administration-approved excipients. However, as with most excipients, their protective mechanism is poorly understood. Here, we apply thermogravimetric analysis, differential scanning calorimetry, and liquid-observed vapor exchange (LOVE) NMR to study the protection of two model globular proteins (the B1 domain of staphylococcal protein G [GB1] and chymotrypsin inhibitor 2 [CI2]) by two desiccation-tolerance proteins (CAHS D from tardigrades and PvLEA4 from an anhydrobiotic midge), as well as by disordered and globular protein controls. We find that all protein samples retain similar amounts of water and possess similar glass transition temperatures, suggesting that neither enhanced water retention nor vitrification is responsible for protection. LOVE NMR reveals that IDPs protect against dehydration-induced unfolding better than the globular protein control, generally protect the same regions of GB1 and CI2, and protect GB1 better than CI2. These observations suggest that electrostatic interactions, charge patterning, and expanded conformations are key to protection. Further application of LOVE NMR to additional client proteins and protectants will deepen our understanding of dehydration protection, enabling the streamlined production of dehydrated proteins for expanded use in the medical, biotechnology, and chemical industries.

Charge Engineering Reveals the Roles of Ionizable Side Chains in Electrospray Ionization Mass Spectrometry

In solution, the charge of a protein is intricately linked to its stability, but electrospray ionization distorts this connection, potentially limiting the ability of native mass spectrometry to inform about protein structure and dynamics. How the behavior of intact proteins in the gas phase depends on the presence and distribution of ionizable surface residues has been difficult to answer because multiple chargeable sites are present in virtually all proteins. Turning to protein engineering, we show that ionizable side chains are completely dispensable for charging under native conditions, but if present, they are preferential protonation sites. The absence of ionizable side chains results in identical charge state distributions under native-like and denaturing conditions, while coexisting conformers can be distinguished using ion mobility separation. An excess of ionizable side chains, on the other hand, effectively modulates protein ion stability. In fact, moving a single ionizable group can dramatically alter the gas-phase conformation of a protein ion. We conclude that although the sum of the charges is governed solely by Coulombic terms, their locations affect the stability of the protein in the gas phase.

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.

Higher-order structural characterisation of native proteins and complexes by top-down mass spectrometry

In biology, it can be argued that if the genome contains the script for a cell's life cycle, then the proteome constitutes an ensemble cast of actors that brings these instructions to life. Their interactions with each other, co-factors, ligands, substrates, and so on, are key to understanding nearly any biological process. Mass spectrometry is well established as the method of choice to determine protein primary structure and location of post-translational modifications. In recent years, top-down fragmentation of intact proteins has been increasingly combined with ionisation of noncovalent assemblies under non-denaturing conditions, i.e., native mass spectrometry. Sequence, post-translational modifications, ligand/metal binding, protein folding, and complex stoichiometry can thus all be probed directly. Here, we review recent developments in this new and exciting field of research. While this work is written primarily from a mass spectrometry perspective, it is targeted to all bioanalytical scientists who are interested in applying these methods to their own biochemistry and chemical biology research.

Charge Movement and Structural Changes in the Gas-Phase Unfolding of Multimeric Protein Complexes Captured by Native Top-Down Mass Spectrometry

The extent to which noncovalent protein complexes retain native structure in the gas phase is highly dependent on experimental conditions. Energetic collisions with background gas can cause structural changes ranging from unfolding to subunit dissociation. Additionally, recent studies have highlighted the role of charge in such structural changes, but the mechanism is not completely understood. In this study, native top down (native TD) mass spectrometry was used to probe gas-phase structural changes of alcohol dehydrogenase (ADH, 4mer) under varying degrees of in-source activation. Changes in covalent backbone fragments produced by electron capture dissociation (ECD) or 193 nm ultraviolet photodissociation (UVPD) were attributed to structural changes of the ADH 4mer. ECD fragments indicated unfolding started at the N-terminus, and the charge states of UVPD fragments enabled monitoring of charge migration to the unfolded regions. Interestingly, UVPD fragments also indicated that the charge at the "unfolding" N-terminus of ADH decreased at high in-source activation energies after the initial increase. We proposed a possible "refolding-after-unfolding" mechanism, as further supported by monitoring hydrogen elimination from radical a-ions produced by UVPD at the N-terminus of ADH. However, "refolding-after-unfolding" with increasing in-source activation was not observed for charge-reduced ADH, which likely adopted compact structures that are resistant to both charge migration and unfolding. When combined, these results support a charge-directed unfolding mechanism for protein complexes. Overall, an experimental framework was outlined for utilizing native TD to generate structure-informative mass spectral signatures for protein complexes that complement other structure characterization techniques, such as ion mobility and computational modeling.

Collision-Induced Unfolding Is Sensitive to the Polarity of Proteins and Protein Complexes

Collision-induced unfolding (CIU) uses ion mobility to probe the structures of ions of proteins and noncovalent complexes as a function of the extent of gas-phase activation prior to analysis. CIU can be sensitive to domain structures, isoform identities, and binding partners, which makes it appealing for many applications. Almost all previous applications of CIU have probed cations. Here, we evaluate the application of CIU to anions and compare the results for anions with those for cations. Towards that end, we developed a "similarity score" that we used to quantify the differences between the results of different CIU experiments and evaluate the significance of those differences relative to the variance of the underlying measurements. Many of the differences between anions and cations that were identified can be attributed to the lower absolute charge states of anions. For example, the extents of the increase in collision cross section over the full range of energies depended strongly on absolute charge state. However, over intermediate energies, there are significant difference between anions and cations with the same absolute charge state. Therefore, CIU is sensitive to the polarity of protein ions. Based on these results, we propose that the utility of CIU to differentiate similar proteins or noncovalent complexes may also depend on polarity. More generally, these results indicate that the relationship between the structures and dynamics of native-like cations and anions deserve further attention and that future studies may benefit from integrating results from ions of both polarities.

Computational Insights into Compaction of Gas-Phase Protein and Protein Complex Ions in Native Ion Mobility-Mass Spectrometry

Native ion mobility-mass spectrometry (IM-MS) is a rapidly growing field for studying the composition and structure of biomolecules and biomolecular complexes using gas-phase methods. Typically, ions are formed in native IM-MS using gentle nanoelectrospray ionization conditions, which in many cases can preserve condensed-phase stoichiometry. Although much evidence shows that large-scale condensed-phase structure, such as quaternary structure and topology, can also be preserved, it is less clear to what extent smaller-scale structure is preserved in native IM-MS. This review surveys computational and experimental efforts aimed at characterizing compaction and structural rearrangements of protein and protein complex ions upon transfer to the gas phase. A brief summary of gas-phase compaction results from molecular dynamics simulations using multiple common force fields and a wide variety of protein ions is presented and compared to literature IM-MS data.

Role of hydrophobic residues for the gaseous formation of helical motifs

The secondary structure content of proteins and their complexes may change significantly on passing from aqueous solution to the gas phase (as in mass spectrometry experiments). In this work, we investigate the impact of hydrophobic residues on the formation of the secondary structure of a real protein complex in the gas phase. We focus on a well-studied protein complex, the amyloid-β (1-40) dimer (2Aβ). Molecular dynamics simulations reproduce the results of ion mobility-mass spectrometry experiments. In addition, a helix (not present in the solution) is identified involving 19FFAED23, consistent with infrared spectroscopy data on an Aβ segment. Our simulations further point to the role of hydrophobic residues in the formation of helical motifs - hydrophobic sidechains "shield" helices from being approached by residues that carry hydrogen bond sites. In particular, two hydrophobic phenylalanine residues, F19 and F20, play an important role for the helix, which is induced in the gas phase in spite of the presence of two carboxyl-containing residues.

A Mass-Spectrometry-Based Modelling Workflow for Accurate Prediction of IgG Antibody Conformations in the Gas Phase

Immunoglobulins are biomolecules involved in defence against foreign substances. Flexibility is key to their functional properties in relation to antigen binding and receptor interactions. We have developed an integrative strategy combining ion mobility mass spectrometry (IM-MS) with molecular modelling to study the conformational dynamics of human IgG antibodies. Predictive models of all four human IgG subclasses were assembled and their dynamics sampled in the transition from extended to collapsed state during IM-MS. Our data imply that this collapse of IgG antibodies is related to their intrinsic structural features, including Fab arm flexibility, collapse towards the Fc region, and the length of their hinge regions. The workflow presented here provides an accurate structural representation in good agreement with the observed collision cross section for these flexible IgG molecules. These results have implications for studying other nonglobular flexible proteins.

Replacing H+ by Na+ or K+ in phosphopeptide anions and cations prevents electron capture dissociation

By successively replacing H+ by Na+ or K+ in phosphopeptide anions and cations, we show that the efficiency of fragmentation into c and z˙ or c˙ and z fragments from N-Cα backbone bond cleavage by negative ion electron capture dissociation (niECD) and electron capture dissociation (ECD) substantially decreases with increasing number of alkali ions attached. In proton-deficient phosphopeptide ions with a net charge of 2-, we observed an exponential decrease in electron capture efficiency with increasing number of Na+ or K+ ions attached, suggesting that electrons are preferentially captured at protonated sites. In proton-abundant phosphopeptide ions with a net charge of 3+, the electron capture efficiency was not affected by replacing up to four H+ ions with Na+ or K+ ions, but the yield of c, z˙ and c˙, z fragments from N-Cα backbone bond cleavage generally decreased next to Na+ or K+ binding sites. We interpret the site-specific decrease in fragmentation efficiency as Na+ or K+ binding to backbone amide oxygen in competition with interactions of protonated sites that would otherwise lead to backbone cleavage into c, z˙ or c˙, z fragments. Our findings seriously challenge the hypothesis that the positive charge responsible for ECD into c, z˙ or c˙, z fragments can generally be a sodium or other metal ion instead of a proton.