Effects of Multidentate Metal Interactions on the Structure of Collisionally Activated Proteins: Insights from Ion Mobility Spectrometry and Molecular Dynamics Simulations

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.

Effects of Hydrogen/Deuterium Exchange on Protein Stability in Solution and in the Gas Phase

Mass spectrometry (MS)-based techniques are widely used for probing protein structure and dynamics in solution. H/D exchange (HDX)-MS is one of the most common approaches in this context. HDX is often considered to be a "benign" labeling method, in that it does not perturb protein behavior in solution. However, several studies have reported that D₂O pushes unfolding equilibria toward the native state. The origin, and even the existence of this protein stabilization remain controversial. Here we conducted thermal unfolding assays in solution to confirm that deuterated proteins in D₂O are more stable, with 2-4 K higher melting temperatures than unlabeled proteins in H₂O. Previous studies tentatively attributed this phenomenon to strengthened H-bonds after deuteration, an effect that may arise from the lower zero-point vibrational energy of the deuterated species. Specifically, it was proposed that strengthened water-water bonds (W···W) in D₂O lower the solubility of nonpolar side chains. The current work takes a broader view by noting that protein stability in solution also depends on water-protein (W···P) and protein-protein (P···P) H-bonds. To help unravel these contributions, we performed collision-induced unfolding (CIU) experiments on gaseous proteins generated by native electrospray ionization. CIU profiles of deuterated and unlabeled proteins were indistinguishable, implying that P···P contacts are insensitive to deuteration. Thus, protein stabilization in D₂O is attributable to solvent effects, rather than alterations of intraprotein H-bonds. Strengthening of W···W contacts represents one possible explanation, but the stabilizing effect of D₂O can also originate from weakened W···P bonds. Future work will be required to elucidate which of these two scenarios is correct, or if both contribute to protein stabilization in D₂O. In any case, the often-repeated adage that "D-bonds are more stable than H-bonds" does not apply to intramolecular contacts in native proteins.

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.

ATP Kinetically Modulates Pathogenic Tau Fibrillations

Advanced understanding of Alzheimer's disease (AD) and several tauopathies over the past decades indicates the pathological importance of tau aggregation in these diseases. Herein, we demonstrated that adenosine triphosphate (ATP), a highly charged anionic molecule found abundantly in the cytosol of cells, catalyzes fibrillation of tau as well as human islet amyloid polypeptide, a representative of basic intrinsically disordered proteins. Our results showed that ATP attracts multiple lysine residues of the four-repeat domain of tau (K18) via supramolecular complexation, thereby forming dimers that are converted to nuclei and accelerate fibril elongation. However, ATP was not directly incorporated into the K18 fibrils, suggesting that ATP plays the role of a catalyst, rather than a reactant, during K18 fibrillation. We also characterized the correlation between ATP dyshomeostasis and tau aggregation in the cellular environment. Our multiple biophysical approaches, including native mass spectrometry (MS), small-angle X-ray scattering (SAXS), and molecular dynamics (MD) simulation, provided insights into the molecular-level influence of ATP on the structural changes and fibrillation of tau.

Compaction of RNA Hairpins and Their Kissing Complexes in Native Electrospray Mass Spectrometry

When electrosprayed from typical native MS solution conditions, RNA hairpins and kissing complexes acquire charge states at which they get significantly more compact in the gas phase than their initial structure in solution. Here, we also show the limits of using force field molecular dynamics to interpret the structures of nucleic acid complexes in the gas phase, as the predicted CCS distributions do not fully match the experimental ones. We suggest that higher level calculation levels should be used in the future.

Enhancing Protein Electrospray Charge States by Multivalent Metal Ions: Mechanistic Insights from MD Simulations and Mass Spectrometry Experiments

The structure and reactivity of electrosprayed protein ions is governed by their net charge. Native proteins in non-denaturing aqueous solutions produce low charge states. More highly charged ions are formed when electrospraying proteins that are unfolded and/or exposed to organic supercharging agents. Numerous studies have explored the electrospray process under these various conditions. One phenomenon that has received surprisingly little attention is the charge enhancement caused by multivalent metal ions such as La³⁺ when electrospraying proteins out of non-denaturing solutions. Here, we conducted mass spectrometry and ion mobility spectrometry experiments, in combination with molecular dynamics (MD) simulations, to uncover the mechanistic basis of this charge enhancement. MD simulations of aqueous ESI droplets reproduced the experimental observation that La³⁺ boosts protein charge states relative to monovalent metals (e.g., Na⁺). The simulations showed that gaseous proteins were released by solvent evaporation to dryness, consistent with the charged residue model. Metal ion ejection kept the shrinking droplets close to the Rayleigh limit until ∼99% of the solvent had left. For droplets charged with Na⁺, metal adduction during the final stage of solvent evaporation produced low protein charge states. Droplets containing La³⁺ showed a very different behavior. The trivalent nature of La³⁺ favored adduction to the protein at a very early stage, when most of the solvent had not evaporated yet. This irreversible binding via multidentate contacts suppressed La³⁺ ejection from the vanishing droplets, such that the resulting gaseous proteins carried significantly more charge. Our results illustrate that MD simulations are suitable for uncovering intricate aspects of electrospray mechanisms, paving the way toward an atomistic understanding of mass spectrometry based analytical workflows.

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.

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.

Cation-Dependent Conformations in 25-Hydroxyvitamin D3-Cation Adducts Measured by Ion Mobility-Mass Spectrometry and Theoretical Modeling

Ion mobility-mass spectrometry is a useful tool in separation of biological isomers, including clinically relevant analytes such as 25-hydroxyvitamin D3 (25OHD3) and its epimer, 3-epi-25-hydroxyvitamin D3 (epi25OHD3). Previous research indicates that these epimers adopt different gas-phase sodiated monomer structures, either the "open" or "closed" conformer, which allow 25OHD3 to be readily resolved in mixtures. In the current work, alternative metal cation adducts are investigated for their relative effects on the ratio of "open" and "closed conformers. Alkali and alkaline earth metal adducts caused changes in the 25OHD3 conformer ratio, where the proportion of the "open" conformer generally increases with the size of the metal cation in a given group. As such, the ratio of the "open" conformer, which is unique to 25OHD3 and absent for its epimer, can be increased from approximately 1:1 for the sodiated monomer to greater than 8:1 for the barium adduct. Molecular modeling and energy calculations agree with the experimental results, indicating that the Gibbs free energy of conversion from the "closed" to the "open" conformation decreased with increasing cation size, correlating with the variation in ratio between the conformers. This work demonstrates the effect of cation adducts on gas-phase conformations of small, flexible molecules and offers an additional strategy for resolution of clinically relevant epimers.

Native Top-Down Mass Spectrometry and Ion Mobility MS for Characterizing the Cobalt and Manganese Metal Binding of α-Synuclein Protein

Structural characterization of intrinsically disordered proteins (IDPs) has been a major challenge in the field of protein science due to limited capabilities to obtain full-length high-resolution structures. Native ESI-MS with top-down MS was utilized to obtain structural features of protein-ligand binding for the Parkinson's disease-related protein, α-synuclein (αSyn), which is natively unstructured. Binding of heavy metals has been implicated in the accelerated formation of αSyn aggregation. Using high-resolution Fourier transform ion cyclotron resonance (FT-ICR) mass spectrometry, native top-down MS with various fragmentation methods, including electron capture dissociation (ECD), collisional activated dissociation (CAD), and multistage tandem MS (MS³), deduced the binding sites of cobalt and manganese to the C-terminal region of the protein. Ion mobility MS (IM-MS) revealed a collapse toward compacted states of αSyn upon metal binding. The combination of native top-down MS and IM-MS provides structural information of protein-ligand interactions for intrinsically disordered proteins. Graphical Abstract ᅟ.

Combining Structural Probes in the Gas Phase - Ion Mobility-Resolved Action-FRET

In the context of native mass spectrometry, the development of gas-phase structural probes sensitive to the different levels of structuration of biomolecular assemblies is necessary to push forward conformational studies. In this paper, we provide the first example of the combination of ion mobility (IM) and Förster resonance energy transfer (FRET) measurements within the same experimental setup. The possibility to obtain mass- and mobility-resolved FRET measurements is demonstrated on a model peptide and applied to monitor the collision-induced unfolding of ubiquitin. Graphical Abstract ᅟ.

Effects of Detergent Micelles on Lipid Binding to Proteins in Electrospray Ionization Mass Spectrometry

A wide variety of biological processes rely upon interactions between proteins and lipids, ranging from molecular transport to the organization of the cell membrane. It was recently established that electrospray ionization mass spectrometry (ESI-MS) is capable of capturing transient interactions between membrane proteins and their lipid environment, and a detailed understanding of the underlying processes is therefore of high importance. Here, we apply ESI-MS to investigate the factors that govern complex formation in solution and gas phases by comparing nonselective lipid binding with soluble and membrane proteins. We find that exogenously added lipids did not bind to soluble proteins, suggesting that lipids have a low propensity to form electrospray ionization adducts. The presence of detergents at increasing micelle concentrations, on the other hand, resulted in moderate lipid binding to soluble proteins. A direct ESI-MS comparison of lipid binding to the soluble protein serum albumin and to the integral membrane protein NapA shows that soluble proteins acquire fewer lipid adducts. Our results suggest that protein–lipid complexes form via contacts between proteins and mixed lipid/detergent micelles. For soluble proteins, these complexes arise from nonspecific contacts between the protein and detergent/lipid micelles in the electrospray droplet. For membrane proteins, lipids are incorporated into the surrounding micelle in solution, and complex formation occurs independently of the ESI process. We conclude that the lipids in the resulting complexes interact predominantly with sites located in the transmembrane segments, resulting in nativelike complexes that can be interrogated by MS.

Selective cysteine modification of metal-free human metallothionein 1a and its isolated domain fragments: Solution structural properties revealed via ESI-MS

Human metallothionein 1a, a protein with two cysteine-rich metal-binding domains (α with 11 Cys and β with 9), was analyzed in its metal-free form by selective, covalent Cys modification coupled with ESI-MS. The modification profiles of the isolated β- and α-fragments reacted with p-benzoquinone (Bq), N-ethylmalemide (NEM) and iodoacetamide (IAM) were compared with the full length protein using ESI-mass spectral data to follow the reaction pathway. Under denaturing conditions at low pH, the reaction profile with each modifier followed pathways that resulted in stochastic, Normal distributions of species whose maxima was equal to the mol. eq. of modifier added. Our interpretation of modification at this pH is that reaction with the cysteines is unimpeded when the full protein or those of its isolated domains are denatured. At neutral pH, where the protein is expected to be folded in a more compact structure, there is a difference in the larger Bq and NEM modification, whose reaction profiles indicate a cooperative pattern. The reaction profile with IAM under native conditions follows a similar stochastic distribution as at low pH, suggesting that this modifier is small enough to access the cysteines unimpeded by the compact structure. The data emphasize the utility of residue modification coupled with electrospray ionization mass spectrometry for the study of protein structure.