Release of Native-like Gaseous Proteins from Electrospray Droplets via the Charged Residue Mechanism: Insights from Molecular Dynamics Simulations

Molecular Dynamics Simulations of Native Protein Charging in Electrosprayed Droplets with Experimentally Relevant Compositions

Electrospray ionization-mass spectrometry (ESI-MS) has been widely used to study proteins given its preservation of native-protein structure when transitioning to the gas-phase. Understanding the influence of experimental factors on ESI can provide insight into the resulting charge states, as well as the degree to which "native" structure is maintained. Experimentally, it is challenging to characterize nanometer-scale electrosprayed droplets; however, molecular dynamics (MD) simulations pose an attractive solution by providing a molecular perspective of protein-ion formation. By resolving many approximations used in past MD simulations of ESI, we demonstrate the capability of simulating electrosprayed droplets with experimentally relevant droplet compositions and behavior. This is accomplished by modeling proton transfers between all titratable molecules in simulated droplets under atmospheric conditions; thus, enabling simulated droplets containing ammonium acetate that form experimentally observed (de)protonated protein ions. Application of the proposed protocol to several native proteins in positive- and negative-ion mode ESI produced simulated weighted averages of charge-state distributions that differed by 0.14 compared to experimental values. Our simulations suggest that changes in residue basicity during the transition to the gas-phase play a significant role in moderating protein charging during native-ESI and can explain many experimentally observed trends. Ionic protein structures produced via the simulated model maintained, on average, 73% of their native contacts into the gas phase relative to solution-phase structures. While applied towards native proteins here, novel insights into effects of the transition to gas-phase enable a deeper understanding of the ESI process itself and thus, are informative regardless of analyte.

Characterizing Protein Solvent Accessible Surface Area in Solution by Dual Polarity Native Mass Spectrometry

Native mass spectrometry (nMS) is rapidly emerging as a pivotal technique for exploring protein conformations and protein-ligand interactions. Pioneering research has demonstrated that the charge state distribution (CSD) of proteins in native mass spectra can be indicative of their solvent accessible surface area (SASA). Moreover, beyond SASA, it is postulated that the abundance of acidic and basic amino acids on the protein surface may also impact the CSD. Specifically, basic amino acids tend to acquire positive charges during electrospray ionization (ESI), whereas acidic amino acids are prone to adopting negative charges. Consequently, this study investigates the CSDs of globular proteins in both positive and negative ion modes to provide a comprehensive characterization of protein SASA. Experiments were conducted under both native ESI and native nano-ESI conditions. By harnessing the average charges observed across dual polarity nMS data, we achieved significantly enhanced log linear correlations between protein SASA and its CSDs. The coefficient of determination (R2) improved from 0.9866 to 0.9888 under ESI conditions and from 0.9677 to 0.9902 under nano-ESI conditions when compared to models utilizing only positive ion mode data. These findings suggest that the SASA of globular proteins can be effectively characterized through the CSDs derived from dual polarity nMS analysis.

Effects of Surface Charge of Amphiphilic Peptides on Peptide-Lipid Interactions in the Gas Phase and in Solution

The interactions between peptides and lipids are fundamental for many biological processes. Therefore, exploring the noncovalent interactions that govern these interactions has become increasingly important. Native mass spectrometry is a valuable technique for the characterization of specific peptide-lipid interactions. However, native mass spectrometry requires the transfer of the analyte into the gas phase, and noncovalent interactions driven by the hydrophobic effect might be distorted. We, therefore, address the importance of electrostatic interactions for the formation of peptide-lipid interactions. For this, we make use of the amphipathic, antimicrobial peptide LL-37 as well as a positively and a negatively charged variant thereof and study binding of a variety of lipids by native mass spectrometry. We found that the surface charge of the peptides affects the transfer of stable peptide-lipid complexes into the gas phase and that the ionization mode is important to observe these interactions. We further compare our findings observed in the gas phase with interactions formed in solution between the peptides and lipid monolayers using a Langmuir film balance. The two approaches deliver comparable results and reveal a clear trend in the lipid preferences of all variants for those lipids with opposite charge. Notably, the unmodified wild-type peptide was more flexible in the formation of peptide-lipid interactions. We conclude that native mass spectrometry is indeed well-suited to explore the interactions between peptides and lipids and that electrostatic interactions as expressed by the surface charge of the peptides play an important role in the formation and stabilization of peptide-lipid interactions.

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.

Effects of anions on the electrospray ionization of proteins in strong acids

The effect of anions on the positive electrospray ionization (ESI) of proteins in different strong acids with varying pH values from 3 to 1 is studied using high-pressure ESI. Reducing the pH from ∼2 to 1 caused a drastic shift in charge state from a high-charge-state distribution (HCSD) to a narrow low-charge-state distribution (LCSD). The shift in charge state was consistent with the circular dichroism result that showed a conformational change due to the "acid-induced folding" of proteins from an unfolding state to a compact molten globule state. Acids of different anions produced noticeable differences in the average charge for HCSD and LCSD. For HCSD, the average charge was lower than the value typically observed using formic and acetic acids. As for LCSD, the average charge was lower than the "native" charge. The high abundance of acid anion that induces the protein compaction was believed to play a role in charge reduction. The effectiveness of anions to "refold" a highly unfolded protein to a compact state and the propensity to reduce the charge of HCSD for proteins appeared to follow the selectivity series of anions towards the stationary phase in ion chromatography. However, the propensity of anions to reduce the charge for LCSD follows quite an opposite trend. The presence of ammonium salt in the acidic solution was found to increase the charge of LCSD. The simple mass spectrum with a narrow distribution of charge state obtained with perchloric acid at pH 1 was demonstrated to facilitate the counting of basic sites.

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 PMCID: PMC11337394 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.

Triboelectric Nanogenerators: Low-Cost Power Supplies for Improved Electrospray Ionization

Electrospray ionization (ESI) is one of the most popular methods to generate ions for mass spectrometry (MS). When compared with other ionization techniques, it can generate ions from liquid-phase samples without additives, retaining covalent and non-covalent interactions of the molecules of interest. When hyphenated to liquid chromatography, it greatly expands the versatility of MS analysis of complex mixtures. However, despite the extensive growth in the application of ESI, the technique still suffers from some drawbacks when powered by direct current (DC) power supplies. Triboelectric nanogenerators promise to be a new power source for the generation of ions by ESI, improving on the analytical capabilities of traditional DC ESI. In this review we highlight the fundamentals of ESI driven by DC power supplies, its contrasting qualities to triboelectric nanogenerator power supplies, and its applications to three distinct fields of research: forensics, metabolomics, and protein structure analysis.

Native Mass Spectrometry: Insights and Opportunities for Targeted Protein Degradation

Native mass spectrometry (nMS) is one of the most powerful biophysical methods for the direct observation of noncovalent protein interactions with both small molecules and other proteins. With the advent of targeted protein degradation (TPD), nMS is now emerging as a compelling approach to characterize the multiple fundamental interactions that underpin the TPD mechanism. Specifically, nMS enables the simultaneous observation of the multiple binary and ternary complexes [i.e., all combinations of E3 ligase, target protein of interest, and small molecule proximity-inducing reagents (such as PROteolysis TArgeting Chimeras (PROTACs) and molecular glues)], formed as part of the TPD equilibrium; this is not possible with any other biophysical method. In this paper we overview the proof-of-concept applications of nMS within the field of TPD and demonstrate how it is providing researchers with critical insight into the systems under study. We also provide an outlook on the scope and future opportunities offered by nMS as a core and agnostic biophysical tool for advancing research developments in TPD.

Evidence of H/D Exchange within Metal-Adducted Carbohydrates after Ion/Ion-Dissociation Reactions

Tandem mass spectrometry (MS/MS) using fragmentation has become one of the most effective methods for gaining sequence and structural information on biomolecules. Ion/ion reactions are competitive reactions, where either proton transfer (PT) or electron transfer (ET) can occur from interactions between multiply charged cations and singly charged anions. Utilizing ion/ion reactions with fluoranthene has offered a unique method of fragment formation for the structural elucidation of biomolecules. Fluoranthene is considered an ideal anion reagent because it selectively causes electron-transfer dissociation (ETD) and minimizes PT when interacting with peptides. However, limited investigations have sought to understand how fluoranthene─the primary, commercially available anion reagent─interacts with other biomolecules. Here, we apply deuterium labeling to investigate ion/ion reaction mechanisms between fluoranthene and divalent, metal-adducted carbohydrates (Ca2+, Mg2+, Co2+, and Ni2+). Deuterium labeling of carbohydrates allowed us to observe evidence of hydrogen/deuterium exchange (HDX) occurring after ion/ion dissociation reactions. The extent of deuterium loss is dependent on several factors, including the physical properties of the metal ion and the fragment structure. Based on the deuterium labeling data, we have proposed ETD, PTD, and intermolecular PT─also described as HDX─mechanisms. This research provides a fundamental perspective of ion/ion and ion/molecule reaction mechanisms and illustrates properties that impact ion/ion and ion/molecule reactions for carbohydrates. Together, this could improve the capability to distinguish complex and heterogeneous biomolecules, such as carbohydrates.

SARS-CoV-2 external structures interacting with nanospheres using docking and molecular dynamics

Coronavirus is caused by the SARS-CoV-2 virus has shown rapid proliferation and scarcity of treatments with proven effectiveness. In this way, we simulated the hospitalization of carbon nanospheres, with external active sites of the SARS-CoV-2 virus (M-Pro, S-Gly and E-Pro), which can be adsorbed or inactivated when interacting with the nanospheres. The computational procedures performed in this work were developed with the SwissDock server for molecular docking and the GROMACS software for molecular dynamics, making it possible to extract relevant data on affinity energy, distance between molecules, free Gibbs energy and mean square deviation of atomic positions, surface area accessible to solvents. Molecular docking indicates that all ligands have an affinity for the receptor's active sites. The nanospheres interact favorably with all proteins, showing promising results, especially C60, which presented the best affinity energy and RMSD values ​​for all protein macromolecules investigated. The C60 with E-Pro exhibited the highest affinity energy of -9.361 kcal/mol, demonstrating stability in both molecular docking and molecular dynamics simulations. Our RMSD calculations indicated that the nanospheres remained predominantly stable, fluctuating within a range of 2 to 3 Å. Additionally, the analysis of other structures yielded promising results that hold potential for application in other proteases.Communicated by Ramaswamy H. Sarma.

Gas-phase structure of polymer ions: Tying together theoretical approaches and ion mobility spectrometry

An increasing number of studies take advantage of ion mobility spectrometry (IMS) coupled to mass spectrometry (IMS-MS) to investigate the spatial structure of gaseous ions. Synthetic polymers occupy a unique place in the field of IMS-MS. Indeed, due to their intrinsic dispersity, they offer a broad range of homologous ions with different lengths. To help rationalize experimental data, various theoretical approaches have been described. First, the study of trend lines is proposed to derive physicochemical and structural parameters. However, the evaluation of data fitting reflects the overall behavior of the ions without reflecting specific information on their conformation. Atomistic simulations constitute another approach that provide accurate information about the ion shape. The overall scope of this review is dedicated to the synergy between IMS-MS and theoretical approaches, including computational chemistry, demonstrating the essential role they play to fully understand/interpret IMS-MS data.

Inter-Domain Repulsion of Dumbbell-Shaped Calmodulin during Electrospray Ionization Revealed by Molecular Dynamics Simulations

The mechanisms whereby protein ions are released from nanodroplets at the liquid-gas interface have continued to be controversial since electrospray ionization (ESI) mass spectrometry was widely applied in biomolecular structure analysis in solution. Several viable pathways have been proposed and verified for single-domain proteins. However, the ESI mechanism of multi-domain proteins with more complicated and flexible structures remains unclear. Herein, dumbbell-shaped calmodulin was chosen as a multi-domain protein model to perform molecular dynamics simulations to investigate the structural evolution during the ESI process. For [Ca₄CAM], the protein followed the classical charge residue model. As the inter-domain electrostatic repulsion increased, the droplet was found to split into two sub-droplets, while stronger-repulsive apo-calmodulin unfolded during the early evaporation stage. We designated this novel ESI mechanism as the domain repulsion model, which provides new mechanistic insights into further exploration of proteins containing more domains. Our results suggest that greater attention should be paid to the effect of domain-domain interactions on structure retention during liquid-gas interface transfer when mass spectrometry is used as the developing technique in gas phase structural biology.

Molecular dynamics study on evaporation of metal nitrate-containing nanodroplets in flame spray pyrolysis

Flame spray pyrolysis (FSP) provides an advantageous synthetic route for LiNi_1-x-yCo_xMn_yO₂ (NCM) materials, which are one of the most practical and promising cathode materials for Li-ion batteries. However, a detailed understanding of the NCM nanoparticle formation mechanisms through FSP is lacking. To shed light on the evaporation of NCM precursor droplets in FSP, in this work, we employ classical molecular dynamics (MD) simulations to explore the dynamic evaporation process of nanodroplets composed of metal nitrates (including LiNO₃, Ni(NO₃)₂, Co(NO₃)₂, and Mn(NO₃)₂ as solutes) and water (as solvent) from a microscopic point of view. Quantitative analysis on the evaporation process has been performed by tracking the temporal evolution of key features including the radial distribution of mass density, the radial distribution of number density of metal ions, droplet diameter, and coordination number (CN) of metal ions with oxygen atoms. Our MD simulation results show that during the evaporation of an MNO₃-containing (M = Li, Ni, Co, or Mn) nanodroplet, Ni²⁺, Co²⁺, and Mn²⁺ will precipitate on the droplet surface, forming a solvent-core-solute-shell structure; whereas the distribution of Li⁺ within the evaporating LiNO₃-containing droplet is more even due to the high diffusivity of Li⁺ compared with other metal ions. For the evaporation of a Ni(NO₃)₂- or Co(NO₃)₂-containing nanodroplet, the temporal evolution of the CN of M-OW (M = Ni or Co; OW represents O atoms from water) suggests a "free H₂O" evaporation stage, during which both CN of M-OW and CN of M-ON are unchanged with time. Evaporation rate constants at various conditions are extracted by making analogy to the classical D² law for droplet evaporation. Unlike Ni or Co, CN of Mn-OW keeps changing with time, yet the temporal evolution of the squared droplet diameter indicates the evaporation rate for a Ni(NO₃)₂-, Co(NO₃)₂-, or Mn(NO₃)₂-containing droplet is hardly affected by the different types of the metal ions.

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.

Recent advances in mass spectrometry analysis of neuropeptides

Due to their involvement in numerous biochemical pathways, neuropeptides have been the focus of many recent research studies. Unfortunately, classic analytical methods, such as western blots and enzyme-linked immunosorbent assays, are extremely limited in terms of global investigations, leading researchers to search for more advanced techniques capable of probing the entire neuropeptidome of an organism. With recent technological advances, mass spectrometry (MS) has provided methodology to gain global knowledge of a neuropeptidome on a spatial, temporal, and quantitative level. This review will cover key considerations for the analysis of neuropeptides by MS, including sample preparation strategies, instrumental advances for identification, structural characterization, and imaging; insightful functional studies; and newly developed absolute and relative quantitation strategies. While many discoveries have been made with MS, the methodology is still in its infancy. Many of the current challenges and areas that need development will also be highlighted in this review.

Release of nanodiscs from charged nano-droplets in the electrospray ionization revealed by molecular dynamics simulations

Electrospray ionization (ESI) is essential for application of mass spectrometry in biological systems, as it prevents the analyte being split into fragments. However, due to lack of a clear understanding of the mechanism of ESI, the interpretation of mass spectra is often ambiguous. This is a particular challenge for complex biological systems. Here, we focus on systems that include nanodiscs as membrane environment, which are essential for membrane proteins. We performed microsecond atomistic molecular dynamics simulations to study the release of nanodiscs from highly charged nano-droplets into the gas phase, the late stage of ESI. We observed two distinct major scenarios, highlighting the diversity of morphologies of gaseous product ions. Our simulations are in reasonable agreement with experimental results. Our work provides a detailed atomistic view of the ESI process of a heterogeneous system (lipid nanodisc), which may give insights into the interpretation of mass spectra of all lipid-protein systems.

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.

Structural Characterization of Dendriplexes In Vacuo: A Joint Ion Mobility/Molecular Dynamics Investigation

The combination between ion mobility mass spectrometry and molecular dynamics simulations is demonstrated for the first time to afford valuable information on structural changes undergone by dendriplexes containing ds-DNA and low-generation dendrimers when transferred from the solution to the gas phase. Dendriplex ions presenting 1:1 and 2:1 stoichiometries are identified using mass spectrometry experiments, and the collision cross sections (CCS) of the 1:1 ions are measured using drift time ion mobility experiments. Structural predictions using Molecular Dynamics (MD) simulations showed that gas-phase relevant structures, i.e., with a good match between the experimental and theoretical CCS, are generated when the global electrospray process is simulated, including the solvent molecule evaporation, rather than abruptly transferring the ions from the solution to the gas phase. The progressive migration of ammonium groups (either NH₄⁺ from the buffer or protonated amines of the dendrimer) into the minor and major grooves of DNA all along the evaporation processes is shown to compact the DNA structure by electrostatic and hydrogen-bond interactions. The subsequent proton transfer from the ammonium (NH₄⁺ or protonated amino groups) to the DNA phosphate groups allows creation of protonated phosphate/phosphate hydrogen bonds within the compact structures. MD simulations showed major structural differences between the dendriplexes in solution and in the gas phase, not only due to the loss of the solvent but also due to the proton transfers and the huge difference between the solution and gas-phase charge states.

Comparing Selected-Ion Collision Induced Unfolding with All Ion Unfolding Methods for Comprehensive Protein Conformational Characterization

Structural analysis by native ion mobility-mass spectrometry provides a direct means to characterize protein interactions, stability, and other biophysical properties of disease-associated biomolecules. Such information is often extracted from collision-induced unfolding (CIU) experiments, performed by ramping a voltage used to accelerate ions entering a trap cell prior to an ion mobility separator. Traditionally, to simplify data analysis and achieve confident ion identification, precursor ion selection with a quadrupole is performed prior to collisional activation. Only one charge state can be selected at one time, leading to an imbalance between the total time required to survey CIU data across all protein charge states and the resulting structural analysis efficiency. Furthermore, the arbitrary selection of a single charge state can inherently bias CIU analyses. We herein aim to compare two conformation sampling methods for protein gas-phase unfolding: (1) traditional quadrupole selection-based CIU and (2) nontargeted, charge selection-free and shotgun workflow, all ion unfolding (AIU). Additionally, we provide a new data interpretation method that integrates across all charge states to project collisional cross section (CCS) data acquired over a range of activation voltages to produce a single unfolding fingerprint, regardless of charge state distributions. We find that AIU in combination with CCS accumulation across all charges offers an opportunity to maximize protein conformational information with minimal time cost, where additional benefits include (1) an improved signal-to-noise ratios for unfolding fingerprints and (2) a higher tolerance to charge state shifts induced by either operating parameters or other factors that affect protein ionization efficiency.

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.

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.

Formation of Carbohydrate-Metal Adducts from Solvent Mixtures during Electrospray: A Molecular Dynamics and ESI-MS Study

Electrospray ionization (ESI) is frequently used to produce gas-phase ions for mass spectrometry (MS)-based techniques. The composition of solvents used in ESI-MS is often manipulated to enhance analyte ionization, including for carbohydrates. Moreover, to characterize analyte structures, ESI has been coupled to hydrogen/deuterium exchange, ion mobility, and tandem MS. Therefore, it is important to understand how solvent composition affects the structure of carbohydrates during and after ESI. In this work, we use molecular dynamics to simulate the desolvation of ESI droplets containing a model carbohydrate and observe the formation of carbohydrate adducts with metal ions. Molecular-level details on the effects of formulating mixtures of water, methanol, and acetonitrile to achieve enhanced ionization are presented. We complement our simulations with ESI-MS experiments. We report that when sprayed from aqueous mixtures containing volatile solvents, carbohydrates ionize to form metal-ion adducts rapidly due to rapid solvent evaporation rather than changes in the ionization mechanism. We find that when sprayed from solvent mixtures, carbohydrates are primarily solvated by water due to the migration of more volatile solvents to the surface of the droplet. Ultimately, the structure of the carbohydrate varies depending on its solvent environment, as inter- and intramolecular interactions are affected. We propose that solvents with 25% or more water may be used to enhance the ionization of carbohydrates with minimal effect on the structure during and after ESI.

p-Aminobenzoic acid protonation dynamics in an evaporating droplet by ab initio molecular dynamics

Protonation equilibria are known to vary from the bulk to microdroplet conditions, which could induce many chemical and physical phenomena. Protonated p-aminobenzoic acid (PABA + H⁺) can be considered a model for probing the protonation dynamics in an evaporating droplet, as its protonation equilibrium is highly dependent on the formation conditions from solution via atmospheric pressure ionization sources. Experiments using diverse experimental techniques have shown that protic solvents allow formation of the O-protomer (PABA protonated in the carboxylic acid group) stable in the gas phase, while aprotic solvents yield the N-protomer (protonated in the amino group) that is the most stable protomer in solution. In this work, we explore the protonation equilibrium of PABA solvated by different numbers of water molecules (n = 0 to 32) using ab initio molecular dynamics. For n = 8-32, the protonation is either at the NH₂ group or in the solvent network. The solvent network interacts with the carboxylic acid group, but there is no complete proton transfer to form the O-protomer. For smaller clusters, however, solvent-mediated proton transfers to the carboxylic acid were observed, both via the Grotthuss mechanism and the vehicle or shuttle mechanism (for n = 1 and 2). Thermodynamic considerations allowed a description of the origins of the kinetic trapping effect, which explains the observation of the solution structure in the gas phase. This effect likely occurs in the final evaporation steps, which are outside the droplet size range covered by previous classical molecular dynamics simulations of charged droplets. These results may be considered relevant in determining the nature of the species observed in the ubiquitous ESI based mass spectrometry analysis, and in general for droplet chemistry, explaining how protonation equilibria are drastically changed from bulk to microdroplet conditions.

Experimental Determination of Activation Energies for Covalent Bond Formation via Ion/Ion Reactions and Competing Processes

The combination of ion/ion chemistry with commercially available ion mobility/mass spectrometry systems has allowed rich structural information to be obtained for gaseous protein ions. Recently, the simple modification of such an instrument with an electrospray reagent source has allowed three-dimensional gas-phase interrogation of protein structures through covalent and noncovalent interactions coupled with collision cross section measurements. However, the energetics of these processes have not yet been studied quantitatively. In this work, previously developed Monte Carlo simulations of ion temperatures inside traveling wave ion guides are used to characterize the energetics of the transition state of activated ubiquitin cation/sulfo-benzoyl-HOAt reagent anion long-lived complexes formed via ion/ion reactions. The ΔH^‡ and ΔS^‡ of major processes observed from collisional activation of long-lived gas-phase ion/ion complexes, namely collision induced unfolding (CIU), covalent bond formation, or neutral loss of the anionic reagent via intramolecular proton transfer, were determined. Covalent bond formation via ion/ion complexes was found to be significantly lower energy compared to unfolding and bond cleavage. The ΔG^‡ values of activation of all three processes lie between 55 and 75 kJ/mol, easily accessible with moderate collisional activation. Bond formation is favored over reagent loss at lower activation energies, whereas reagent loss becomes competitive at higher collision energies. Though the ΔG^‡ values between CIU of a precursor ion and covalent bond formation of its ion/ion product complex are comparable, our data suggest covalent bond formation does not require extensive isomerization.

Protons Are Fast and Smart; Proteins Are Slow and Dumb: On the Relationship of Electrospray Ionization Charge States and Conformations

We present simple considerations of how differences in time scales of motions of protons, the lightest and fastest chemical moiety, and the much longer time scales associated with the dynamics of proteins, among the heaviest and slowest analytes, may allow many protein conformations from solution to be kinetically trapped during the process of electrospraying protein solutions into the gas phase. In solution, the quantum nature of protons leads them to change locations by tunneling, an instantaneous process; moreover, the Grotthuss mechanism suggests that these small particles can respond nearly instantaneously to the dynamic motions of proteins that occur on much longer time scales. A conformational change is accompanied by favorable or unfavorable variations in the free energy of the system, providing the impetus for solvent ↔ protein proton exchange. Thus, as thermal distributions of protein conformations interconvert, protonation states rapidly respond, as specific acidic and basic sites are exposed or protected. In the vacuum of the mass spectrometer, protons become immobilized in locations that are specific to the protein conformations from which they were incorporated. In this way, conformational states from solution are preserved upon electrospraying them into the gas phase. These ideas are consistent with the exquisite sensitivity of electrospray mass spectra to small changes of the local environment that alter protein structure in solution. We might remember this approximation for the protonation of proteins in solution with the colloquial expression-protons are fast and smart; proteins are slow and dumb.

Native mass spectrometry and gas-phase fragmentation provide rapid and in-depth topological characterization of a PROTAC ternary complex

Proteolysis-targeting chimeras (PROTACs) represent a new direction in small-molecule therapeutics whereby a heterobifunctional linker to a protein of interest (POI) induces its ubiquitination-based proteolysis by recruiting an E3 ligase. Here, we show that charge reduction, native mass spectrometry, and gas-phase activation methods combine for an in-depth analysis of a PROTAC-linked ternary complex. Electron capture dissociation (ECD) of the intact POI-PROTAC-VCB complex (a trimeric subunit of an E3 ubiquitin ligase) promotes POI dissociation. Collision-induced dissociation (CID) causes elimination of the nonperipheral PROTAC, producing an intact VCB-POI complex not seen in solution but consistent with PROTAC-induced protein-protein interactions. In addition, we used ion mobility spectrometry (IMS) and collisional activation to identify the source of this unexpected dissociation. Together, the evidence shows that this integrated approach can be used to screen for ternary complex formation and PROTAC-protein contacts and may report on PROTAC-induced protein-protein interactions, a characteristic correlated with PROTAC selectivity and efficacy.

Toward In Silico Prediction of CO2 Diffusion in Champagne Wines

Carbon dioxide diffusion is the main physical process behind the formation and growth of bubbles in sparkling wines, especially champagne wines. By approximating brut-labeled champagnes as carbonated hydroalcoholic solutions, molecular dynamics (MD) simulations are carried out with six rigid water models and three CO₂ models to evaluate CO₂ diffusion coefficients. MD simulations are little sensitive to the CO₂ model but proper water modeling is essential to reproduce experimental measurements. A satisfactory agreement with nuclear magnetic resonance (NMR) data is only reached at all temperatures for simulations based on the OPC and TIP4P/2005 water models; the similar efficiency of these two models is attributed to their common properties such as low mixture enthalpy, same number of hydrogen bonds, alike water tetrahedrality, and multipole values. Correcting CO₂ diffusion coefficients to take into account their system-size dependence does not significantly alter the quality of the results. Estimates of viscosities deduced from the Stokes-Einstein formula are found in excellent agreement with viscometry on brut-labeled champagnes, while theoretical densities tend to underestimate experimental values. OPC and TIP4P/2005 water models appear to be choice water models to investigate CO₂ solvation and transport properties in carbonated hydroalcoholic mixtures and should be the best candidates for any MD simulations concerning wines, spirits, or multicomponent mixtures with alike chemical composition.

Fluorescent Janus emulsions for biosensing of Listeria monocytogenes

Here we report a sensing method for Listeria monocytogenes based on the agglutination of all-liquid Janus emulsions. This two-dye assay enables the rapid detection of trace Listeria in less than 2 h via an emissive signal produced in response to Listeria binding. The biorecognition interface between the Janus emulsions is assembled by attaching antibodies to a functional surfactant polymer with a tetrazine/transcyclooctene click reaction. The strong binding between Listeria and the Listeria antibody located at the hydrocarbon surface of the emulsions results in the tilting of the Janus structure from its equilibrium position to produce emission that would ordinarily be obscured by a blocking dye. This method provides rapid and inexpensive Listeria detection with high sensitivity (<100 CFU/mL in 2 h) that can be paired with many antibody or related recognition elements to create a new class of biosensors.

Gas-phase conformations of intrinsically disordered proteins and their complexes with ligands: Kinetically trapped states during transfer from solution to the gas phase

Flexible structures of intrinsically disordered proteins (IDPs) are crucial for versatile functions in living organisms, which involve interaction with diverse partners. Electrospray ionization ion mobility mass spectrometry (ESI-IM-MS) has been widely applied for structural characterization of apo-state and ligand-associated IDPs via two-dimensional separation in the gas phase. Gas-phase IDP structures have been regarded as kinetically trapped states originated from conformational features in solution. However, an implication of the states remains elusive in the structural characterization of IDPs, because it is unclear what structural property of IDPs is preserved. Recent studies have indicated that the conformational features of IDPs in solution are not fully reproduced in the gas phase. Nevertheless, the molecular interactions captured in the gas phase amplify the structural differences between IDP conformers. Therefore, an IDP conformational change that is not observed in solution is observable in the gas-phase structures obtained by ESI-IM-MS. Herein, we have presented up-to-date researches on the key implications of kinetically trapped states in the gas phase with a brief summary of the structural dynamics of IDPs in ESI-IM-MS.

Structural mass spectrometry goes viral

Over the last 20 years, mass spectrometry (MS), with its ability to analyze small sample amounts with high speed and sensitivity, has more and more entered the field of structural virology, aiming to investigate the structure and dynamics of viral proteins as close to their native environment as possible. The use of non-perturbing labels in hydrogen-deuterium exchange MS allows for the analysis of interactions between viral proteins and host cell factors as well as their dynamic responses to the environment. Cross-linking MS, on the other hand, can analyze interactions in viral protein complexes and identify virus-host interactions in cells. Native MS allows transferring viral proteins, complexes and capsids into the gas phase and has broken boundaries to overcome size limitations, so that now even the analysis of intact virions is possible. Different MS approaches not only inform about size, stability, interactions and dynamics of virus assemblies, but also bridge the gap to other biophysical techniques, providing valuable constraints for integrative structural modeling of viral complex assemblies that are often inaccessible by single technique approaches. In this review, recent advances are highlighted, clearly showing that structural MS approaches in virology are moving towards systems biology and ever more experiments are performed on cellular level.

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.

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

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

Charging and supercharging of proteins for mass spectrometry: recent insights into the mechanisms of electrospray ionization

Molecular dynamics simulations have uncovered mechanistic details of the protein ESI process under various experimental conditions.

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.

Resolving the Discrepancies Between Empirical and Rayleigh Charge Limiting Models for Globular Proteins

Starting with the Rayleigh charge limiting model, a slightly different approach is used to account for the well-known discrepancy that exists between the said model and experimental ESI MS data for globular proteins. It is shown using published datasets that for globular proteins, the mass density ρ exhibits a weak second-order dependence on its mass M, according to ρ(M)∝ M^-α, α ~ 0.14. A direct equivalence established between ESI MS and x-ray techniques suggests a minimum but critical surface tension of 15.6 ± 5.2 mN/m for the droplet at the liquid-to-gas phase transition point. The packing density factor η for globular proteins is believed to lie between 1 (very tightly packed) and 4.6 (less tight, natively packed). While the Rayleigh charge limiting model has been linked historically to the CRM (J. Chem. Phys. 49:2240-2249, 1968; Anal. Chim. Acta 406:93-104, 2000), this paper does not expressly seek to justify the CRM, but rather uses empirical data and existing knowledge across subfields to help build a consistent picture of ESI MS phenomena that might be difficult to explain otherwise. These results would be useful in molecular dynamics (MD) simulations, understanding liquid-to-gas phase transitions and in opening up new routes for cross-calibration between ESI MS, IM MS, NMR and x-ray crystallography studies. Graphical Abstract ᅟ.

Molecular dynamics and metadynamics simulations of electrosprayed water nanodroplets including sodium bis(2-ethylhexyl)sulfosuccinate micelles

The behavior of aqueous solutions of sodium bis(2-ethylhexyl)sulfosuccinate (AOTNa) under conditions of electrospray ionization (ESI) has been investigated by molecular dynamics (MD) and well-tempered metadynamics (WTM) simulations at 300 K and 400 K. We have examined water droplets with initial fixed numbers of water molecules (1000) and AOT^- anions (100), and with sodium cations in the range of 70-130. At 300 K, all charged droplets show the water evaporation rate increasing with the absolute value of the initial droplet charge state (Z), accompanied by ejection of an increasing number of solvated sodium ions or by expulsion of AOT^- anions depending on the sign of Z and by fragmentation in the case of high |Z|. At 400 K, the water evaporation becomes more rapid and the fission process more extensive. In all cases, the AOTNa molecules, arranged as a direct micelle inside the aqueous system, undergo a rapid inversion in vacuo so that the hydrophilic heads and sodium ions surrounded by water molecules move toward the droplet interior. At the end of the 100-ns MD simulations, some water molecules remain within the aggregates at both temperatures. The subsequent metadynamics simulations accelerate the droplet evolution and show that all systems become anhydrous, in agreement with the experimental results of ESI mass spectrometry. This complete water loss is accompanied by sodium counterion emission for positively charged aggregates at 300 K. The analysis shows how the temperature and droplet charge state affect the populations of the generated surfactant aggregates, providing information potentially useful in designing future ESI experimental conditions.

Charge-State-Dependent Variation of Signal Intensity Ratio between Unbound Protein and Protein-Ligand Complex in Electrospray Ionization Mass Spectrometry: The Role of Solvent-Accessible Surface Area

Native electrospray ionization mass spectrometry (ESI-MS) is nowadays widely used for the direct and sensitive determination of protein complex stoichiometry and binding affinity constants ( K_a). A common yet poorly understood phenomenon in native ESI-MS is the difference between the charge-state distributions (CSDs) of the bound protein-ligand complex (PL) and unbound protein (P) signals. This phenomenon is typically attributed to experimental artifacts such as nonspecific binding or in-source dissociation and is considered highly undesirable, because the determined K_a values display strong variation with charge state. This situation raises serious concerns regarding the reliability of ESI-MS for the analysis of protein complexes. Here we demonstrate that, contrary to the common belief, the CSD difference between P and PL ions can occur without any loss of complex integrity, simply due to a change in the solvent-accessible surface area (ΔSASA) of the protein upon ligand binding in solution. The experimental CSD shifts for PL and P ions in ESI-MS are explained in relation to the magnitude of ΔSASA for diverse protein-ligand systems using a simple model based on the charged residue mechanism. Our analysis shows that the revealed ΔSASA factor should be considered rather general and be given attention for the correct spectral interpretation of protein complexes.

Addressing a Common Misconception: Ammonium Acetate as Neutral pH "Buffer" for Native Electrospray Mass Spectrometry

Native ESI-MS involves the transfer of intact proteins and biomolecular complexes from solution into the gas phase. One potential pitfall is the occurrence of pH-induced changes that can affect the analyte while it is still surrounded by solvent. Most native ESI-MS studies employ neutral aqueous ammonium acetate solutions. It is a widely perpetuated misconception that ammonium acetate buffers the analyte solution at neutral pH. By definition, a buffer consists of a weak acid and its conjugate weak base. The buffering range covers the weak acid pK_a ± 1 pH unit. NH₄⁺ and CH₃-COO^- are not a conjugate acid/base pair, which means that they do not constitute a buffer at pH 7. Dissolution of ammonium acetate salt in water results in pH 7, but this pH is highly labile. Ammonium acetate does provide buffering around pH 4.75 (the pK_a of acetic acid) and around pH 9.25 (the pK_a of ammonium). This implies that neutral ammonium acetate solutions electrosprayed in positive ion mode will likely undergo acidification down to pH 4.75 ± 1 in the ESI plume. Ammonium acetate nonetheless remains a useful additive for native ESI-MS. It is a volatile electrolyte that can mimic the solvation properties experienced by proteins under physiological conditions. Also, a drop from pH 7 to around pH 4.75 is less dramatic than the acidification that would take place in pure water. It is hoped that the habit of referring to pH 7 solutions as ammonium acetate "buffer" will disappear from the literature. Ammonium acetate "solution" should be used instead. Graphical Abstract ᅟ.