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

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.

A suicidal and extensively disordered luciferase with a bright luminescence

Gaussia luciferase (GLuc) is one of the most luminescent luciferases known and is widely used as a reporter in biochemistry and cell biology. During catalysis, GLuc undergoes inactivation by irreversible covalent modification. The mechanism by which GLuc generates luminescence and how it becomes inactivated are however not known. Here, we show that GLuc unlike other enzymes has an extensively disordered structure with a minimal hydrophobic core and no apparent binding pocket for the main substrate, coelenterazine. From an alanine scan, we identified two Arg residues required for light production. These residues separated with an average of about 22 Å and a major structural rearrangement is required if they are to interact with the substrate simultaneously. We furthermore show that in addition to coelenterazine, GLuc also can oxidize furimazine, however, in this case without production of light. Both substrates result in the formation of adducts with the enzyme, which eventually leads to enzyme inactivation. Our results demonstrate that a rigid protein structure and substrate-binding site are no prerequisites for high enzymatic activity and specificity. In addition to the increased understanding of enzymes in general, the findings will facilitate future improvement of GLuc as a reporter luciferase.

Characterization of Mass, Diameter, Density, and Surface Properties of Colloidal Nanoparticles Enabled by Charge Detection Mass Spectrometry

A variety of scattering-based, microscopy-based, and mobility-based methods are frequently used to probe the size distributions of colloidal nanoparticles with transmission electron microscopy (TEM) often considered to be the "gold standard". Charge detection mass spectrometry (CDMS) is an alternative method for nanoparticle characterization that can rapidly measure the mass and charge of individual nanoparticle ions with high accuracy. Two low polydispersity, ∼100 nm diameter nanoparticle size standards with different compositions (polymethyl methacrylate/polystyrene copolymer and 100% polystyrene) were characterized using both TEM and CDMS to explore the merits and complementary aspects of both methods. Mass and diameter distributions are rapidly obtained from CDMS measurements of thousands of individual ions of known spherical shape, requiring less time than TEM sample preparation and image analysis. TEM image-to-image variations resulted in a ∼1-2 nm range in the determined mean diameters whereas the CDMS mass precision of ∼1% in these experiments leads to a diameter uncertainty of just 0.3 nm. For the 100% polystyrene nanoparticles with known density, the CDMS and TEM particle diameter distributions were in excellent agreement. For the copolymer nanoparticles with unknown density, the diameter from TEM measurements combined with the mass from CDMS measurements enabled an accurate measurement of nanoparticle density. Differing extents of charging for the two nanoparticle standards measured by CDMS show that charging is sensitive to nanoparticle surface properties. A mixture of the two samples was separated based on their different extents of charging despite having overlapping mass distributions centered at 341.5 and 331.0 MDa.

Metastability of Protein Solution Structures in the Absence of a Solvent: Rugged Energy Landscape and Glass-like Behavior

Native ion mobility/mass spectrometry is well-poised to structurally screen proteomes but characterizes protein structures in the absence of a solvent. This raises long-standing unanswered questions about the biological significance of protein structures identified through ion mobility/mass spectrometry. Using newly developed computational and experimental ion mobility/ion mobility/mass spectrometry methods, we investigate the unfolding of the protein ubiquitin in a solvent-free environment. Our data suggest that the folded, solvent-free ubiquitin observed by ion mobility/mass spectrometry exists in a largely native fold with an intact β-grasp motif and α-helix. The ensemble of folded, solvent-free ubiquitin ions can be partitioned into kinetically stable subpopulations that appear to correspond to the structural heterogeneity of ubiquitin in solution. Time-resolved ion mobility/ion mobility/mass spectrometry measurements show that folded, solvent-free ubiquitin exhibits a strongly stretched-exponential time dependence, which simulations trace to a rugged energy landscape with kinetic traps. Unfolding rate constants are estimated to be approximately 800 to 20,000 times smaller than in the presence of water, effectively quenching the unfolding process on the time scale of typical ion mobility/mass spectrometry measurements. Our proposed unfolding pathway of solvent-free ubiquitin shares substantial characteristics with that established for the presence of solvent, including a polarized transition state with significant native content in the N-terminal β-hairpin and α-helix. Our experimental and computational data suggest that (1) the energy landscape governing the motions of folded, solvent-free proteins is rugged in analogy to that of glassy systems; (2) large-scale protein motions may at least partially be determined by the amino acid sequence of a polypeptide chain; and (3) solvent facilitates, rather than controls, protein motions.

Advances in mass spectrometry to unravel the structure and function of protein condensates

Membrane-less organelles assemble through liquid-liquid phase separation (LLPS) of partially disordered proteins into highly specialized microenvironments. Currently, it is challenging to obtain a clear understanding of the relationship between the structure and function of phase-separated protein assemblies, owing to their size, dynamics and heterogeneity. In this Perspective, we discuss recent advances in mass spectrometry (MS) that offer several promising approaches for the study of protein LLPS. We survey MS tools that have provided valuable insights into other insoluble protein systems, such as amyloids, and describe how they can also be applied to study proteins that undergo LLPS. On the basis of these recent advances, we propose to integrate MS into the experimental workflow for LLPS studies. We identify specific challenges and future opportunities for the analysis of protein condensate structure and function by MS.

Exploiting ion-mobility mass spectrometry for unraveling proteome complexity

Ion mobility spectrometry-mass spectrometry (IMS-MS) is experiencing rapid growth in proteomic studies, driven by its enhancements in dynamic range and throughput, increasing the quantitation precision, and the depth of proteome coverage. The core principle of ion mobility spectrometry is to separate ions in an inert gas under the influence of an electric field based on differences in drift time. This minireview provides an introduction to IMS operation modes and a description of advantages and limitations is presented. Moreover, the principles of trapped IMS-MS (TIMS-MS), including parallel accumulation-serial fragmentation are discussed. Finally, emerging applications linked to TIMS focusing on sample throughput (in clinical proteomics) and sensitivity (single-cell proteomics) are reviewed, and the possibilities of intact protein analysis are discussed.

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.

Ultraviolet Photodissociation Reveals the Molecular Mechanism of Crown Ether Microsolvation Effect on the Gas-Phase Native-like Protein Structure

Maintaining the protein high-order structures and interactions during the transition from aqueous solution to gas phase is essential to the structural analysis of native mass spectrometry (nMS). Herein, we systematically interrogate the effects of charge state and crown ether (CE) complexation on the gas-phase native-like protein structure by integrating nMS with 193 nm ultraviolet photodissociation (UVPD). The alterations of photofragmentation yields of protein residues and the charge site distribution of fragment ions reveal the specific sites and sequence regions where charge and CE take effect. Our results exhibit the CE complexation on protonated residues can largely alleviate the structure disruption induced by the intramolecular solvation of charged side chains. The influences of CE complexation and positive charge on gas-phase protein structure exhibit generally opposite trends because the CE microsolvation avoids the hydrogen-bonding formation between the charged side chains with backbone carbonyls. Thus, CE complexation leads to a more stable and native-like protein structure in the gas phase.

A "grappling hook" interaction connects self-assembly and chaperone activity of Nucleophosmin 1

How the self-assembly of partially disordered proteins generates functional compartments in the cytoplasm and particularly in the nucleus is poorly understood. Nucleophosmin 1 (NPM1) is an abundant nucleolar protein that forms large oligomers and undergoes liquid-liquid phase separation by binding RNA or ribosomal proteins. It provides the scaffold for ribosome assembly but also prevents protein aggregation as part of the cellular stress response. Here, we use aggregation assays and native mass spectrometry (MS) to examine the relationship between the self-assembly and chaperone activity of NPM1. We find that oligomerization of full-length NPM1 modulates its ability to retard amyloid formation in vitro. Machine learning-based structure prediction and cryo-electron microscopy reveal fuzzy interactions between the acidic disordered region and the C-terminal nucleotide-binding domain, which cross-link NPM1 pentamers into partially disordered oligomers. The addition of basic peptides results in a tighter association within the oligomers, reducing their capacity to prevent amyloid formation. Together, our findings show that NPM1 uses a "grappling hook" mechanism to form a network-like structure that traps aggregation-prone proteins. Nucleolar proteins and RNAs simultaneously modulate the association strength and chaperone activity, suggesting a mechanism by which nucleolar composition regulates the chaperone activity of NPM1.

Linking Gas-Phase and Solution-Phase Protein Unfolding via Mobile Proton Simulations

Native mass spectrometry coupled to ion mobility (IM-MS) combined with collisional activation (CA) of ions in the gas phase (in vacuo) is an important method for the study of protein unfolding. It has advantages over classical biophysical and structural techniques as it can be used to analyze small volumes of low-concentration heterogeneous mixtures while maintaining solution-like behavior and does not require labeling with fluorescent or other probes. It is unclear, however, whether the unfolding observed during collision activation experiments mirrors solution-phase unfolding. To bridge the gap between in vacuo and in-solution behavior, we use unbiased molecular dynamics (MD) to create in silico models of in vacuo unfolding of a well-studied protein, the N-terminal domain of ribosomal L9 (NTL9) protein. We utilize a mobile proton algorithm (MPA) to create 100 thermally unfolded and coulombically unfolded in silico models for observed charge states of NTL9. The unfolding behavior in silico replicates the behavior in-solution and is in line with the in vacuo observations; however, the theoretical collision cross section (CCS) of the in silico models was lower compared to that of the in vacuo data, which may reflect reduced sampling.

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.

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.

Molecular Dynamics-Assisted Interpretation of Experimentally Determined Intrinsically Disordered Protein Conformational Components: The Case of Human α-Synuclein

Mass spectrometry and single molecule force microscopy are two experimental approaches able to provide structural information on intrinsically disordered proteins (IDPs). These techniques allow the dissection of conformational ensembles in their main components, although at a low-resolution level. In this work, we interpret the results emerging from these experimental approaches on human alpha synuclein (AS) by analyzing a previously published 73 μs-long molecular-dynamics (MD) simulation of the protein in explicit solvent. We further compare MD-based predictions of single molecule Förster resonance energy transfer (smFRET) data of AS in solution with experimental data. The combined theoretical and experimental data provide a description of AS main conformational ensemble, shedding light into its intramolecular interactions and overall structural compactness. This analysis could be easily transferred to other IDPs.

Condensable Vapor Sorption by Low Charge State Protein Ions

Measurement of the gas-phase ion mobility of proteins provides a means to quantitatively assess the relative sizes of charged proteins. However, protein ion mobility measurements are typically singular values. Here, we apply tandem mobility analysis to low charge state protein ions (+1 and +2 ions) introduced into the gas phase by nanodroplet nebulization. We first determine protein ion mobilities in dry air and subsequently examine shifts in mobilities brought about by the clustering of vapor molecules. Tandem mobility analysis yields mobility-vapor concentration curves for each protein ion, expanding the information obtained from mobility analysis. This experimental procedure and analysis is extended to bovine serum albumin, transferrin, immunoglobulin G, and apoferritin with water, 1-butanol, and nonane. All protein ions appear to adsorb vapor molecules, with mobility "diameter" shifts of up to 6-7% at conditions just below vapor saturation. We parametrize results using κ-Köhler theory, where the term κ quantifies the extent of uptake beyond Köhler model expectations. For 1-butanol and nonane, κ decreases with increasing protein ion size, while it increases with increasing protein ion size for water. For the systems probed, the extent of mobility shift for the organic vapors is unaffected by the nebulized solution pH, while shifts with water are sensitive to pH.