Proton Dynamics in Protein Mass Spectrometry

Sialylation-induced stabilization of dynamic glycoprotein conformations unveiled by time-aligned parallel unfolding and glycan release mass spectrometry

Sialylation, a critical post-translational modification, regulates glycoprotein structure and function by tuning their molecular heterogeneity. However, characterizing its subtle and dynamic conformational effects at the intact glycoprotein level remains challenging. We introduce a glycoform-resolved unfolding approach based on a high-throughput ion mobility-mass spectrometry (IM-MS) platform. This method integrates high-throughput unfolding with parallel fragmentation, enabling simultaneous analysis of sialylation patterns, stoichiometries, and their impact on conformational stability. Applying this approach to fetuin, we identified distinct sialylation patterns and their differential influence on protein conformation, namely sialylation-induced stabilization during early unfolding and increased flexibility in later unfolding stages. IM-MS-guided molecular dynamics simulations revealed that increased sialylation enhances the initial conformational stability, likely through enhanced electrostatic interactions and hydrogen bonding. These findings highlight the complex interplay between sialylation and protein dynamics and establish glycoform-resolved unfolding IM-MS as a powerful tool for characterizing glycoprotein conformational landscapes.

Development, Optimization, and Clinical Relevance of Lactoferrin Delivery Systems: A Focus on Ocular Delivery

Lactoferrin (Lf), a multifunctional protein found abundantly in secretions, including tears, plays a crucial role in ocular health through its antimicrobial, immunoregulatory, anti-inflammatory, and antioxidant activities. Advanced delivery systems are desirable to fully leverage its therapeutic potential in treating ocular diseases. The process of Lf quantification for diagnostic purposes underscores the importance of developing reliable, cost-effective detection methods, ranging from conventional techniques to advanced nano-based sensors. Despite the ease and non-invasiveness of topical administration for ocular surface diseases, challenges such as rapid drug elimination necessitate innovations, such as Lf-loaded contact lenses and biodegradable polymeric nanocapsules, to enhance drug stability and bioavailability. Furthermore, overcoming ocular barriers for the treatment of posterior segment disease calls for nano-formulations. The scope of this review is to underline the advancements in nanotechnology-based Lf delivery methods, emphasizing the pivotal role of multidisciplinary approaches and cross-field strategies in improving ocular drug delivery and achieving better therapeutic outcomes for a wide spectrum of eye conditions.

Structural studies of supramolecular complexes and assemblies by ion mobility mass spectrometry

Recent advances in instrumentation and development of computational strategies for ion mobility mass spectrometry (IM-MS) studies have contributed to an extensive growth in the application of this analytical technique to comprehensive structural description of supramolecular systems. Apart from the benefits of IM-MS for interrogation of intrinsic properties of noncovalent aggregates in the experimental gas-phase environment, its merits for the description of native structural aspects, under the premises of having maintained the noncovalent interactions innate upon the ionization process, have attracted even more attention and gained increasing interest in the scientific community. Thus, various types of supramolecular complexes and assemblies relevant for biological, medical, material, and environmental sciences have been characterized so far by IM-MS supported by computational chemistry. This review covers the state-of-the-art in this field and discusses experimental methods and accompanying computational approaches for assessing the reliable three-dimensional structural elucidation of supramolecular complexes and assemblies by IM-MS.

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.

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.

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

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

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.

Changes in Protonation States of In-Pathway Residues can Alter Ligand Binding Pathways Obtained From Spontaneous Binding Molecular Dynamics Simulations

Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations.

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.

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

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

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.

HIV-1 Gag protein with or without p6 specifically dimerizes on the viral RNA packaging signal

The HIV-1 Gag protein is responsible for genomic RNA (gRNA) packaging and immature viral particle assembly. Although the presence of gRNA in virions is required for viral infectivity, in its absence, Gag can assemble around cellular RNAs and form particles resembling gRNA-containing particles. When gRNA is expressed, it is selectively packaged despite the presence of excess host RNA, but how it is selectively packaged is not understood. Specific recognition of a gRNA packaging signal (Psi) has been proposed to stimulate the efficient nucleation of viral assembly. However, the heterogeneity of Gag-RNA interactions renders capturing this transient nucleation complex using traditional structural biology approaches challenging. Here, we used native MS to investigate RNA binding of wild-type (WT) Gag and Gag lacking the p6 domain (GagΔp6). Both proteins bind to Psi RNA primarily as dimers, but to a control RNA primarily as monomers. The dimeric complexes on Psi RNA require an intact dimer interface within Gag. GagΔp6 binds to Psi RNA with high specificity in vitro and also selectively packages gRNA in particles produced in mammalian cells. These studies provide direct support for the idea that Gag binding to Psi specifically promotes nucleation of Gag-Gag interactions at the early stages of immature viral particle assembly in a p6-independent manner.

Secondary structure effects on internal proton transfer in poly-peptides

A pump-probe approach was designed to determine the internal proton transfer (PT) rate in a series of poly-peptide radical cations containing both histidine and tryptophan. The proton transfer is driven by the gas-phase basicity difference between residues. The fragmentation scheme indicates that the gas-phase basicity of histidine is lower than that of radical tryptophan so that histidine is always pulling the proton away from tryptophan. However, the proton transfer requires the two basic sites to be in close proximity, which is rate limited by the peptide conformational dynamics. PT rate measurements were used to probe and explore the peptide conformational dynamics in several poly-glycines/prolines/alanines. For small and unstructured peptides, the PT rate decreases with the size, as expected from a statistical point of view in a flat conformational space. Conversely, if structured conformations are accessible, the structural flexibility of the peptide is decreased. This slows down the occurrence of conformations favorable to proton transfer. A dramatic decrease in the PT rates was observed for peptides HAnW, when n changes from 5 to 6. This is attributed to the onset of a stable helix for n = 6. No such discontinuity is observed for poly-glycines or poly-prolines. In HAnW, the gas-phase basicity and helix propensity compete for the position of the charge. Interestingly, in this competition between PT and helix formation in HA6W, the energy gain associated with helix formation is large enough to slow down the PT beyond experimental time but does not ultimately prevail over the proton preference for histidine.

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

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

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.

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.

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.