Proteins, Lipids, and Water in the Gas Phase

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.

Enhanced EMC-Advantages of partially known orientations in x-ray single particle imaging

Single particle imaging of proteins in the gas phase with x-ray free-electron lasers holds great potential to study fast protein dynamics, but is currently limited by weak and noisy data. A further challenge is to discover the proteins' orientation as each protein is randomly oriented when exposed to x-rays. Algorithms such as the expand, maximize, and compress (EMC) exist that can solve the orientation problem and reconstruct the three-dimensional diffraction intensity space, given sufficient measurements. If information about orientation were known, for example, by using an electric field to orient the particles, the reconstruction would benefit and potentially reach better results. We used simulated diffraction experiments to test how the reconstructions from EMC improve with particles' orientation to a preferred axis. Our reconstructions converged to correct maps of the three-dimensional diffraction space with fewer measurements if biased orientation information was considered. Even for a moderate bias, there was still significant improvement. Biased orientations also substantially improved the results in the case of missing central information, in particular in the case of small datasets. The effects were even more significant when adding a background with 50% the strength of the averaged diffraction signal photons to the diffraction patterns, sometimes reducing the data requirement for convergence by a factor of 10. This demonstrates the usefulness of having biased orientation information in single particle imaging experiments, even for a weaker bias than what was previously known. This could be a key component in overcoming the problems with background noise that currently plague these experiments.

Cryo-EM of soft-landed β-galactosidase: Gas-phase and native structures are remarkably similar

Native mass spectrometry (MS) has become widely accepted in structural biology, providing information on stoichiometry, interactions, homogeneity, and shape of protein complexes. Yet, the fundamental assumption that proteins inside the mass spectrometer retain a structure faithful to native proteins in solution remains a matter of intense debate. Here, we reveal the gas-phase structure of β-galactosidase using single-particle cryo-electron microscopy (cryo-EM) down to 2.6-Å resolution, enabled by soft landing of mass-selected protein complexes onto cold transmission electron microscopy (TEM) grids followed by in situ ice coating. We find that large parts of the secondary and tertiary structure are retained from the solution. Dehydration-driven subunit reorientation leads to consistent compaction in the gas phase. By providing a direct link between high-resolution imaging and the capability to handle and select protein complexes that behave problematically in conventional sample preparation, the approach has the potential to expand the scope of both native mass spectrometry and cryo-EM.

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.

Rehydration Post-orientation: Investigating Field-Induced Structural Changes via Computational Rehydration

Proteins can be oriented in the gas phase using strong electric fields, which brings advantages for structure determination using X-ray free electron lasers. Both the vacuum conditions and the electric-field exposure risk damaging the protein structures. Here, we employ molecular dynamics simulations to rehydrate and relax vacuum and electric-field exposed proteins in aqueous solution, which simulates a refinement of structure models derived from oriented gas-phase proteins. We find that the impact of the strong electric fields on the protein structures is of minor importance after rehydration, compared to that of vacuum exposure and ionization in electrospraying. The structures did not fully relax back to their native structure in solution on the simulated timescales of 200 ns, but they recover several features, including native-like intra-protein contacts, which suggests that the structures remain in a state from which the fully native structure is accessible. Our findings imply that the electric fields used in native mass spectrometry are well below a destructive level, and suggest that structures inferred from X-ray diffraction from gas-phase proteins are relevant for solution and in vivo conditions, at least after in silico rehydration.

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.

Protein shape sampled by ion mobility mass spectrometry consistently improves protein structure prediction

Ion mobility (IM) mass spectrometry provides structural information about protein shape and size in the form of an orientationally-averaged collision cross-section (CCS_IM). While IM data have been used with various computational methods, they have not yet been utilized to predict monomeric protein structure from sequence. Here, we show that IM data can significantly improve protein structure determination using the modelling suite Rosetta. We develop the Rosetta Projection Approximation using Rough Circular Shapes (PARCS) algorithm that allows for fast and accurate prediction of CCS_IM from structure. Following successful testing of the PARCS algorithm, we use an integrative modelling approach to utilize IM data for protein structure prediction. Additionally, we propose a confidence metric that identifies near native models in the absence of a known structure. The results of this study demonstrate the ability of IM data to consistently improve protein structure prediction.

Collision cross sections (CCS) from ion mobility mass spectrometry provide information about protein shape and size. Here, the authors develop an algorithm to predict CCS and integrate experimental ion mobility data into Rosetta-based molecular modelling to predict protein structures from sequence.

Exploring the Conformational Landscape and Stability of Aurora A Using Ion-Mobility Mass Spectrometry and Molecular Modeling

Protein kinase inhibitors are highly effective in treating diseases driven by aberrant kinase signaling and as chemical tools to help dissect the cellular roles of kinase signaling complexes. Evaluating the effects of binding of small molecule inhibitors on kinase conformational dynamics can assist in understanding both inhibition and resistance mechanisms. Using gas-phase ion-mobility mass spectrometry (IM-MS), we characterize changes in the conformational landscape and stability of the protein kinase Aurora A (Aur A) driven by binding of the physiological activator TPX2 or small molecule inhibition. Aided by molecular modeling, we establish three major conformations, the relative abundances of which were dependent on the Aur A activation status: one highly populated compact conformer similar to that observed in most crystal structures, a second highly populated conformer possessing a more open structure infrequently found in crystal structures, and an additional low-abundance conformer not currently represented in the protein databank. Notably, inhibitor binding induces more compact configurations of Aur A, as adopted by the unbound enzyme, with both IM-MS and modeling revealing inhibitor-mediated stabilization of active Aur A.

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.

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.

Lipid Head Group Adduction to Soluble Proteins Follows Gas-Phase Basicity Predictions: Dissociation Barriers and Charge Abstraction

Native mass spectrometry analysis of membrane proteins has yielded many useful insights in recent years with respect to membrane protein-lipid interactions, including identifying specific interactions and even measuring binding affinities based on observed abundances of lipid-bound ions after collision-induced dissociation (CID). However, the behavior of non-covalent complexes subjected to extensive CID can in principle be affected by numerous factors related to gas-phase chemistry, including gas-phase basicity (GB) and acidity, shared-proton bonds, and other factors. A recent report from our group showed that common lipids span a wide range of GB values. Notably, phosphatidylcholine (PC) and sphingomyelin lipids are more basic than arginine, suggesting they may strip charge upon dissociation in positive ion mode, while phosphoserine lipids are slightly less basic than arginine and may form especially strong shared-proton bonds. Here, we use CID to probe the strength of non-specific gas-phase interactions between lipid head groups and several soluble proteins, used to deliberately avoid possible physiological protein-lipid interactions. The strengths of the protein-head group interactions follow the trend predicted based solely on lipid and amino acid GBs: phosphoserine (PS) head group forms the strongest bonds with these proteins and out-competes the other head groups studied, while glycerophosphocholine (GPC) head groups form the weakest interactions and dissociate carrying away a positive charge. These results indicate that gas-phase thermochemistry can play an important role in determining which head groups remain bound to protein ions with native-like structures and charge states in positive ion mode upon extensive collisional activation.

Material and Charge Transport of Large Organic Salt Clusters and Nanoparticles in Electrospray Ion Beam Deposition

Electrospray ion beam deposition (ES-IBD) or ion soft landing has been demonstrated as a technique suitable for processing nonvolatile molecules in vacuum under perfectly controlled conditions, an approach also desirable for the deposition of nanoparticles. Here, we present results from several approaches to generate, characterize, and deposit nanoparticle ion beams in vacuum for deposition. We focus on cluster ion beams generated by ESI of organic salt solutions. Small cluster ions of the salts appear in the mass spectra as defined peaks. In addition, we find nanoparticle-sized aggregates, appearing as a low intensity background at high m/z-ratio, and show by IBD experiments that these clusters carry the major amount of material in the ion beam. This transition from clusters to nanoparticles, and their successful deposition, shows that ES-IBD can in principle handle ion beams of very heavy and highly charged nanoparticles. In related experiments, however, we found the deposition of nanoparticles from dispersions to be of low reproducibility, due to the lack of control by mass spectrometry.

Oligomeric complexes formed by Redβ single strand annealing protein in its different DNA bound states

Redβ is a single strand annealing protein from bacteriophage λ that binds loosely to ssDNA, not at all to pre-formed dsDNA, but tightly to a duplex intermediate of annealing. As viewed by electron microscopy, Redβ forms oligomeric rings on ssDNA substrate, and helical filaments on the annealed duplex intermediate. However, it is not clear if these are the functional forms of the protein in vivo. We have used size-exclusion chromatography coupled with multi-angle light scattering, analytical ultracentrifugation and native mass spectrometry (nMS) to characterize the size of the oligomers formed by Redβ in its different DNA-bound states. The nMS data, which resolve species with the highest resolution, reveal that Redβ forms an oligomer of 12 subunits in the absence of DNA, complexes ranging from 4 to 14 subunits on 38-mer ssDNA, and a much more distinct and stable complex of 11 subunits on 38-mer annealed duplex. We also measure the concentration of Redβ in cells active for recombination and find it to range from 7 to 27 μM. Collectively, these data provide new insights into the dynamic nature of the complex on ssDNA, and the more stable and defined complex on annealed duplex.

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.

A GPU-Accelerated Fast Multipole Method for GROMACS: Performance and Accuracy

An important and computationally demanding part of molecular dynamics simulations is the calculation of long-range electrostatic interactions. Today, the prevalent method to compute these interactions is particle mesh Ewald (PME). The PME implementation in the GROMACS molecular dynamics package is extremely fast on individual GPU nodes. However, for large scale multinode parallel simulations, PME becomes the main scaling bottleneck as it requires all-to-all communication between the nodes; as a consequence, the number of exchanged messages scales quadratically with the number of involved nodes in that communication step. To enable efficient and scalable biomolecular simulations on future exascale supercomputers, clearly a method with a better scaling property is required. The fast multipole method (FMM) is such a method. As a first step on the path to exascale, we have implemented a performance-optimized, highly efficient GPU FMM and integrated it into GROMACS as an alternative to PME. For a fair performance comparison between FMM and PME, we first assessed the accuracies of the methods for various sets of input parameters. With parameters yielding similar accuracies for both methods, we determined the performance of GROMACS with FMM and compared it to PME for exemplary benchmark systems. We found that FMM with a multipole order of 8 yields electrostatic forces that are as accurate as PME with standard parameters. Further, for typical mixed-precision simulation settings, FMM does not lead to an increased energy drift with multipole orders of 8 or larger. Whereas an ≈50 000 atom simulation system with our FMM reaches only about a third of the performance with PME, for systems with large dimensions and inhomogeneous particle distribution, e.g., aerosol systems with water droplets floating in a vacuum, FMM substantially outperforms PME already on a single node.

Structural Heterogeneity in Single Particle Imaging Using X-ray Lasers

One of the challenges facing single particle imaging with ultrafast X-ray pulses is the structural heterogeneity of the sample to be imaged. For the method to succeed with weakly scattering samples, the diffracted images from a large number of individual proteins need to be averaged. The more the individual proteins differ in structure, the lower the achievable resolution in the final reconstructed image. We use molecular dynamics to simulate two globular proteins in vacuum, fully desolvated as well as with two different solvation layers, at various temperatures. We calculate the diffraction patterns based on the simulations and evaluate the noise in the averaged patterns arising from the structural differences and the surrounding water. Our simulations show that the presence of a minimal water coverage with an average 3 Å thickness will stabilize the protein, reducing the noise associated with structural heterogeneity, whereas additional water will generate more background noise.

Comparative Structural Analysis of 20S Proteasome Ortholog Protein Complexes by Native Mass Spectrometry

Ortholog protein complexes are responsible for equivalent functions in different organisms. However, during evolution, each organism adapts to meet its physiological needs and the environmental challenges imposed by its niche. This selection pressure leads to structural diversity in protein complexes, which are often difficult to specify, especially in the absence of high-resolution structures. Here, we describe a multilevel experimental approach based on native mass spectrometry (MS) tools for elucidating the structural preservation and variations among highly related protein complexes. The 20S proteasome, an essential protein degradation machinery, served as our model system, wherein we examined five complexes isolated from different organisms. We show that throughout evolution, from the Thermoplasma acidophilum archaeal prokaryotic complex to the eukaryotic 20S proteasomes in yeast (Saccharomyces cerevisiae) and mammals (rat - Rattus norvegicus, rabbit - Oryctolagus cuniculus and human - HEK293 cells), the proteasome increased both in size and stability. Native MS structural signatures of the rat and rabbit 20S proteasomes, which heretofore lacked high-resolution, three-dimensional structures, highly resembled that of the human complex. Using cryoelectron microscopy single-particle analysis, we were able to obtain a high-resolution structure of the rat 20S proteasome, allowing us to validate the MS-based results. Our study also revealed that the yeast complex, and not those in mammals, was the largest in size and displayed the greatest degree of kinetic stability. Moreover, we also identified a new proteoform of the PSMA7 subunit that resides within the rat and rabbit complexes, which to our knowledge have not been previously described. Altogether, our strategy enables elucidation of the unique structural properties of protein complexes that are highly similar to one another, a framework that is valid not only to ortholog protein complexes, but also for other highly related protein assemblies.

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.

Weighing-up protein dynamics: the combination of native mass spectrometry and molecular dynamics simulations

Structural dynamics underpin biological function at the molecular level, yet many biophysical and structural biology approaches give only a static or averaged view of proteins. Native mass spectrometry yields spectra of the many states and interactions in the structural ensemble, but its spatial resolution is limited. Conversely, molecular dynamics simulations are innately high-resolution, but have a limited capacity for exploring all structural possibilities. The two techniques hence differ fundamentally in the information they provide, returning data that reflect different length scales and time scales, making them natural bedfellows. Here we discuss how the combination of native mass spectrometry with molecular dynamics simulations is enabling unprecedented insights into a range of biological questions by interrogating the motions of proteins, their assemblies, and interactions.

An integrated native mass spectrometry and top-down proteomics method that connects sequence to structure and function of macromolecular complexes

Mass spectrometry (MS) has become a crucial technique for the analysis of protein complexes. Native MS has traditionally examined protein subunit arrangements, while proteomics MS has focused on sequence identification. These two techniques are usually performed separately without taking advantage of the synergies between them. Here we describe the development of an integrated native MS and top-down proteomics method using Fourier-transform ion cyclotron resonance (FTICR) to analyse macromolecular protein complexes in a single experiment. We address previous concerns of employing FTICR MS to measure large macromolecular complexes by demonstrating the detection of complexes up to 1.8 MDa, and we demonstrate the efficacy of this technique for direct acquirement of sequence to higher-order structural information with several large complexes. We then summarize the unique functionalities of different activation/dissociation techniques. The platform expands the ability of MS to integrate proteomics and structural biology to provide insights into protein structure, function and regulation.

Combining Mass Spectrometry and X-Ray Crystallography for Analyzing Native-Like Membrane Protein Lipid Complexes

Membrane proteins represent a challenging family of macromolecules, particularly related to the methodology aimed at characterizing their three-dimensional structure. This is mostly due to their amphipathic nature as well as requirements of ligand bindings to stabilize or control their function. Recently, Mass Spectrometry (MS) has become an important tool to identify the overall stoichiometry of native-like membrane proteins complexed to ligand bindings as well as to provide insights into the transport mechanism across the membrane, with complementary information coming from X-ray crystallography. This perspective article emphasizes MS findings coupled with X-ray crystallography in several membrane protein lipid complexes, in particular transporters, ion channels and molecular machines, with an overview of techniques that allows a more thorough structural interpretation of the results, which can help us to unravel hidden mysteries on the membrane protein function.

Micelles of the chiral biocompatible surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethylammonium bromide (DMEB): molecular dynamics and fragmentation patterns in the gas phase

RATIONALE

The study of self-assembly processes of surfactant molecules in the gas phase is of great interest for several theoretical and technological reasons related to their possible exploitation as drug carriers, protein shields and cleaning agents in the gas phase.

METHODS

The stability and fragmentation patterns of singly and multiply charged (either positively or negatively) aggregates of the surfactant (1R,2S)-dodecyl(2-hydroxy-1-methyl-2-phenylethyl)dimethyl ammonium bromide (DMEB) in the gas phase have been studied by ion mobility mass spectrometry and tandem mass spectrometry. Molecular dynamics (MD) simulations of positively and negatively singly and multiply charged DMEB aggregates have been performed to obtain structural and energetics information. Finally, in order to ascertain some clues on the DMEB growth mechanism, quantum mechanics calculations were carried out.

RESULTS

It has been evidenced that positively and negatively singly charged aggregates at low collision energy decompose preferentially by loss of only one DMEB molecule. Increasing the collision energy, the loss of neutrals becomes increasingly abundant. Multiply charged DMEB aggregates are unstable and decompose forming singly charged monomers or dimers. MD simulations show reverse micelle-like structures with polar heads somewhat segregated into the aggregate interior. Finally, a good correlation between experimental and calculated collisional cross sections (CCS) was found.

CONCLUSIONS

The fragmentation pathways of DMEB charged species evidenced for singly charged aggregates exhibit features similar to that of other detergent aggregates, but multiply charged aggregates showed a system-specific behavior. QM calculations on the optimized structures (2¹⁺ , 3¹⁺ , 1^1- and 2^1- ) indicate that the most determinant interactions are due to an OH---Br hydrogen bonding that is also involved in the link between monomeric DMEB units. The MD models gave CCS values in good agreement with experimental ones, evidenced by a less strict reverse micelle-like structure and a reasonably spread bromine anion distribution Copyright © 2017 John Wiley & Sons, Ltd.

Micelles, Rods, Liposomes, and Other Supramolecular Surfactant Aggregates: Computational Approaches

Surfactants are an interesting class of compounds characterized by the segregation of polar and apolar domains in the same molecule. This peculiarity makes possible a whole series of microscopic and macroscopic effects. Among their features, their ability to segregate particles (fluids or entire domains) and to reduce the surface/interfacial tension is the utmost important. The interest in the chemistry of surfactants never weakened; instead, waves of increasing interest have occurred every time a new field of application of these molecules has been discovered. All these special characteristics depend largely on the ability of surfactants to self-assemble and constitute supramolecular structures where their chemical properties are amplified. The possibility to obtain structural and energy information and, above all, the possibility of forecast the self-organizing mechanisms of surfactants have had a significant boost via computational chemistry. The molecular dynamics models, initially coarse-grained and subsequently (with the increasing computer power) using more accurate models, allowed, over the years, to better understand different aspects of the processes of dispersion, self-assembly, segregation of surfactant. Moreover, several other aspects have been investigated as the effect of the counterions of many ionic surfactants in defining the final supramolecular structures, the mobility of side chains, and the capacity of some surfactant to envelope entire proteins. This review constitutes a perspective/prospective view of these results. On the other hand, some comparison of in silico results with experimental information recently acquired through innovative analytical techniques such as ion mobility mass spectrometry which have been introduced.

Surface Accessibility and Dynamics of Macromolecular Assemblies Probed by Covalent Labeling Mass Spectrometry and Integrative Modeling

Mass spectrometry (MS) has become an indispensable tool for investigating the architectures and dynamics of macromolecular assemblies. Here we show that covalent labeling of solvent accessible residues followed by their MS-based identification yields modeling restraints that allow mapping the location and orientation of subunits within protein assemblies. Together with complementary restraints derived from cross-linking and native MS, we built native-like models of four heterocomplexes with known subunit structures and compared them with available X-ray crystal structures. The results demonstrated that covalent labeling followed by MS markedly increased the predictive power of the integrative modeling strategy enabling more accurate protein assembly models. We applied this strategy to the F-type ATP synthase from spinach chloroplasts (cATPase) providing a structural basis for its function as a nanomotor. By subjecting the models generated by our restraint-based strategy to molecular dynamics (MD) simulations, we revealed the conformational states of the peripheral stalk and assigned flexible regions in the enzyme. Our strategy can readily incorporate complementary chemical labeling strategies and we anticipate that it will be applicable to many other systems providing new insights into the structure and function of protein complexes.

Integrating mass spectrometry with MD simulations reveals the role of lipids in Na+/H+ antiporters

Na+/H+ antiporters are found in all kingdoms of life and exhibit catalysis rates that are among the fastest of all known secondary-active transporters. Here we combine ion mobility mass spectrometry and molecular dynamics simulations to study the conformational stability and lipid-binding properties of the Na+/H+ exchanger NapA from Thermus thermophilus and compare this to the prototypical antiporter NhaA from Escherichia coli and the human homologue NHA2. We find that NapA and NHA2, but not NhaA, form stable dimers and do not selectively retain membrane lipids. By comparing wild-type NapA with engineered variants, we show that the unfolding of the protein in the gas phase involves the disruption of inter-domain contacts. Lipids around the domain interface protect the native fold in the gas phase by mediating contacts between the mobile protein segments. We speculate that elevator-type antiporters such as NapA, and likely NHA2, use a subset of annular lipids as structural support to facilitate large-scale conformational changes within the membrane.

Are Charge-State Distributions a Reliable Tool Describing Molecular Ensembles of Intrinsically Disordered Proteins by Native MS?

Native mass spectrometry (MS) has become a central tool of structural proteomics, but its applicability to the peculiar class of intrinsically disordered proteins (IDPs) is still object of debate. IDPs lack an ordered tridimensional structure and are characterized by high conformational plasticity. Since they represent valuable targets for cancer and neurodegeneration research, there is an urgent need of methodological advances for description of the conformational ensembles populated by these proteins in solution. However, structural rearrangements during electrospray-ionization (ESI) or after the transfer to the gas phase could affect data obtained by native ESI-MS. In particular, charge-state distributions (CSDs) are affected by protein conformation inside ESI droplets, while ion mobility (IM) reflects protein conformation in the gas phase. This review focuses on the available evidence relating IDP solution ensembles with CSDs, trying to summarize cases of apparent consistency or discrepancy. The protein-specificity of ionization patterns and their responses to ligands and buffer conditions suggests that CSDs are imprinted to protein structural features also in the case of IDPs. Nevertheless, it seems that these proteins are more easily affected by electrospray conditions, leading in some cases to rearrangements of the conformational ensembles. Graphical Abstract ᅟ.

A data set from flash X-ray imaging of carboxysomes

Ultra-intense femtosecond X-ray pulses from X-ray lasers permit structural studies on single particles and biomolecules without crystals. We present a large data set on inherently heterogeneous, polyhedral carboxysome particles. Carboxysomes are cell organelles that vary in size and facilitate up to 40% of Earth's carbon fixation by cyanobacteria and certain proteobacteria. Variation in size hinders crystallization. Carboxysomes appear icosahedral in the electron microscope. A protein shell encapsulates a large number of Rubisco molecules in paracrystalline arrays inside the organelle. We used carboxysomes with a mean diameter of 115±26 nm from Halothiobacillus neapolitanus. A new aerosol sample-injector allowed us to record 70,000 low-noise diffraction patterns in 12 min. Every diffraction pattern is a unique structure measurement and high-throughput imaging allows sampling the space of structural variability. The different structures can be separated and phased directly from the diffraction data and open a way for accurate, high-throughput studies on structures and structural heterogeneity in biology and elsewhere.

Role of single-point mutations and deletions on transition temperatures in ideal proteinogenic heteropolymer chains in the gas phase

A coarse-grained statistical mechanics-based model for ideal heteropolymer proteinogenic chains of non-interacting residues is presented in terms of the size K of the chain and the set of helical propensities [Formula: see text] associated with each residue j along the chain. For this model, we provide an algorithm to compute the degeneracy tensor [Formula: see text] associated with energy level [Formula: see text] where [Formula: see text] is the number of residues with a native contact in a given conformation. From these results, we calculate the equilibrium partition function [Formula: see text] and characteristic temperature [Formula: see text] at which a transition from a low to a high entropy states is observed. The formalism is applied to analyze the effect on characteristic temperatures [Formula: see text] of single-point mutations and deletions of specific amino acids [Formula: see text] along the chain. Two probe systems are considered. First, we address the case of a random heteropolymer of size K and given helical propensities [Formula: see text] on a conformational phase space. Second, we focus our attention to a particular set of neuropentapeptides, [Met-5] and [Leu-5] enkephalins whose thermodynamic stability is a key feature on their coupling to [Formula: see text] and [Formula: see text] receptors and the triggering of biochemical responses.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Two-dimensional honeycomb network through sequence-controlled self-assembly of oligopeptides

The sequence of a peptide programs its self-assembly and hence the expression of specific properties through non-covalent interactions. A large variety of peptide nanostructures has been designed employing different aspects of these non-covalent interactions, such as dispersive interactions, hydrogen bonding or ionic interactions. Here we demonstrate the sequence-controlled fabrication of molecular nanostructures using peptides as bio-organic building blocks for two-dimensional (2D) self-assembly. Scanning tunnelling microscopy reveals changes from compact or linear assemblies (angiotensin I) to long-range ordered, chiral honeycomb networks (angiotensin II) as a result of removal of steric hindrance by sequence modification. Guided by our observations, molecular dynamic simulations yield atomistic models for the elucidation of interpeptide-binding motifs. This new approach to 2D self-assembly on surfaces grants insight at the atomic level that will enable the use of oligo- and polypeptides as large, multi-functional bio-organic building blocks, and opens a new route towards rationally designed, bio-inspired surfaces.

Peptide nanostructures are currently arousing interest thanks to their potential applications in medicine, electronics and coatings. Here, through experiment and theory, the authors demonstrate exquisite control over surface peptide assembly behaviour through manipulation of amino acid sequence.

Conformational effects in protein electrospray-ionization mass spectrometry

Electrospray-ionization mass spectrometry (ESI-MS) is a key tool of structural biology, complementing the information delivered by conventional biochemical and biophysical methods. Yet, the mechanism behind the conformational effects in protein ESI-MS is an object of debate. Two parameters-solvent-accessible surface area (As) and apparent gas-phase basicity (GBapp)-are thought to play a role in controlling the extent of protein ionization during ESI-MS experiments. This review focuses on recent experimental and theoretical investigations concerning the influence of these parameters on ESI-MS results and the structural information that can be derived. The available evidence supports a unified model for the ionization mechanism of folded and unfolded proteins. These data indicate that charge-state distribution (CSD) analysis can provide valuable structural information on normally folded, as well as disordered structures.

Charged supramolecular assemblies of surfactant molecules in gas phase

The aim of this review is to critically analyze recent literature on charged supramolecular assemblies formed by surfactant molecules in gas phase. Apart our specific interest on this research area, the stimuli to undertake the task arise from the widespread theoretical and applicative benefits emerging from a comprehensive view of this topic. In fact, the study of the formation, stability, and physicochemical peculiarities of non-covalent assemblies of surfactant molecules in gas phase allows to unveil interesting aspects such as the role of attractive, repulsive, and steric intermolecular interactions as driving force of supramolecular organization in absence of interactions with surrounding medium and the size and charge state dependence of aggregate structural and dynamical properties. Other interesting aspects worth to be investigated are joined to the ability of these assemblies to incorporate selected solubilizates molecules as well as to give rise to chemical reactions within a single organized structure. In particular, the incorporation of large molecules such as proteins has been of recent interest with the objective to protect their structure and functionality during the transition from solution to gas phase. Exciting fall-out of the study of gas phase surfactant aggregates includes mass and energy transport in the atmosphere, origin of life and simulation of supramolecular aggregation in the interstellar space. Moreover, supramolecular assemblies of amphiphilic molecules in gas phase could find remarkable applications as atmospheric cleaning agents, nanosolvents and nanoreactors for specialized chemical processes in confined space. Mass spectrometry techniques have proven to be particularly suitable to generate these assemblies and to furnish useful information on their size, size polydispersity, stability, and structural organization. On the other hand molecular dynamics simulations have been very useful to rationalize many experimental findings and to furnish a vivid picture of the structural and dynamic features of these aggregates. Thus, in this review, we will focus on the most important achievements gained in recent years by both these investigative tools.

Protein Structural Studies by Traveling Wave Ion Mobility Spectrometry: A Critical Look at Electrospray Sources and Calibration Issues

The question whether electrosprayed protein ions retain solution-like conformations continues to be a matter of debate. One way to address this issue involves comparisons of collision cross sections (Ω) measured by ion mobility spectrometry (IMS) with Ω values calculated for candidate structures. Many investigations in this area employ traveling wave IMS (TWIMS). It is often implied that nanoESI is more conducive for the retention of solution structure than regular ESI. Focusing on ubiquitin, cytochrome c, myoglobin, and hemoglobin, we demonstrate that Ω values and collisional unfolding profiles are virtually indistinguishable under both conditions. These findings suggest that gas-phase structures and ion internal energies are independent of the type of electrospray source. We also note that TWIMS calibration can be challenging because differences in the extent of collisional activation relative to drift tube reference data may lead to ambiguous peak assignments. It is demonstrated that this problem can be circumvented by employing collisionally heated calibrant ions. Overall, our data are consistent with the view that exposure of native proteins to electrospray conditions can generate kinetically trapped ions that retain solution-like structures on the millisecond time scale of TWIMS experiments. ᅟ

Ion Mobility-Mass Spectrometry Differentiates Protein Quaternary Structures Formed in Solution and in Electrospray Droplets

Electrospray ionization coupled to mass spectrometry is a key technology for determining the stoichiometries of multiprotein complexes. Despite highly accurate results for many assemblies, challenging samples can generate signals for artifact protein-protein binding born of the crowding forces present within drying electrospray droplets. Here, for the first time, we study the formation of preferred protein quaternary structures within such rapidly evaporating nanodroplets. We use ion mobility and tandem mass spectrometry to investigate glutamate dehydrogenase dodecamers and serum amyloid P decamers as a function of protein concentration, along with control experiments using carefully chosen protein analogues, to both establish the formation of operative mechanisms and assign the bimodal conformer populations observed. Further, we identify an unprecedented symmetric collision-induced dissociation pathway that we link directly to the quaternary structures of the precursor ions selected.

Structural organization of surfactant aggregates in vacuo: a molecular dynamics and well-tempered metadynamics study

MD and well tempered metadynamics indicate that the structural organization of large surfactant aggregates in the gas phase and in the condensed apolar phase are different.

Spectroscopic studies of kinetically trapped conformations in the gas phase: the case of triply protonated bradykinin

We explore conformational space of triply protonated bradykinin. Three conformational families are mobility-separated and spectroscopically characterized. Kinetically trapped structures are identified via annealing.

Active conformation control of unfolded proteins by hyperthermal collision with a metal surface

The physical and chemical properties of macromolecules like proteins are strongly dependent on their conformation. The degrees of freedom of their chemical bonds generate a huge conformational space, of which, however, only a small fraction is accessible in thermal equilibrium. Here we show that soft-landing electrospray ion beam deposition (ES-IBD) of unfolded proteins allows to control their conformation. The dynamics and result of the deposition process can be actively steered by selecting the molecular ion beam's charge state or tuning the incident energy. Using these parameters, protein conformations ranging from fully extended to completely compact can be prepared selectively on a surface, as evidenced on the subnanometer/amino acid resolution level by scanning tunneling microscopy (STM). Supported by molecular dynamics (MD) simulations, our results demonstrate that the final conformation on the surface is reached through a mechanical deformation during the hyperthermal ion surface collision. Our experimental results independently confirm the findings of ion mobility spectrometry (IMS) studies of protein gas phase conformations. Moreover, we establish a new route for the processing of macromolecular materials, with the potential to reach conformations that would be inaccessible otherwise.

Molecular simulation-based structural prediction of protein complexes in mass spectrometry: the human insulin dimer

Protein electrospray ionization (ESI) mass spectrometry (MS)-based techniques are widely used to provide insight into structural proteomics under the assumption that non-covalent protein complexes being transferred into the gas phase preserve basically the same intermolecular interactions as in solution. Here we investigate the applicability of this assumption by extending our previous structural prediction protocol for single proteins in ESI-MS to protein complexes. We apply our protocol to the human insulin dimer (hIns2) as a test case. Our calculations reproduce the main charge and the collision cross section (CCS) measured in ESI-MS experiments. Molecular dynamics simulations for 0.075 ms show that the complex maximizes intermolecular non-bonded interactions relative to the structure in water, without affecting the cross section. The overall gas-phase structure of hIns2 does exhibit differences with the one in aqueous solution, not inferable from a comparison with calculated CCS. Hence, care should be exerted when interpreting ESI-MS proteomics data based solely on NMR and/or X-ray structural information.

Insights into the mechanism of protein electrospray ionization from salt adduction measurements

The mechanisms whereby protein ions are liberated from charged droplets during electrospray ionization (ESI) remain under investigation. Compact conformers electrosprayed from aqueous solution in positive ion mode likely follow the charged residue model (CRM), which envisions analyte release after solvent evaporation to dryness. The concentration of nonvolatile salts such as NaCl increases sharply within vanishing CRM droplets, promoting nonspecific pairing of Cl(-) and Na(+) with charged groups on the protein surface. For unfolded proteins, it has been proposed that ion formation occurs via the chain ejection model (CEM). During the CEM proteins are expelled from the droplet long before complete solvent evaporation has taken place. Here we examine whether salt adduction levels support the view that folded and unfolded proteins follow different ESI mechanisms. Solvent evaporation during the CEM is expected to be less extensive and, hence, the salt concentration at the point of protein release should be substantially lower than for the CRM. CEM ions should therefore exhibit lower adduction levels than CRM species. We explore the adduction behavior of several proteins that were chosen to allow comparative studies on folded and unfolded structures in the same solution. In-source activation eliminates chloride adducts via HCl release, generating protein ions that are heterogeneously charged because of sodiation and protonation. Sodiation levels measured under such conditions provide estimates of the salt adduction behavior experienced by the "nascent" analyte ions. Sodiation levels are significantly reduced for unfolded proteins, supporting the view that these species are indeed formed via the CEM.

On the Modeling of Polar Component of Solvation Energy using Smooth Gaussian-Based Dielectric Function

Traditional implicit methods for modeling electrostatics in biomolecules use a two-dielectric approach: a biomolecule is assigned low dielectric constant while the water phase is considered as a high dielectric constant medium. However, such an approach treats the biomolecule-water interface as a sharp dielectric border between two homogeneous dielectric media and does not account for inhomogeneous dielectric properties of the macromolecule as well. Recently we reported a new development, a smooth Gaussian-based dielectric function which treats the entire system, the solute and the water phase, as inhomogeneous dielectric medium (J Chem Theory Comput. 2013 Apr 9; 9(4): 2126-2136.). Here we examine various aspects of the modeling of polar solvation energy in such inhomogeneous systems in terms of the solute-water boundary and the inhomogeneity of the solute in the absence of water surrounding. The smooth Gaussian-based dielectric function is implemented in the DelPhi finite-difference program, and therefore the sensitivity of the results with respect to the grid parameters is investigated, and it is shown that the calculated polar solvation energy is almost grid independent. Furthermore, the results are compared with the standard two-media model and it is demonstrated that on average, the standard method overestimates the magnitude of the polar solvation energy by a factor 2.5. Lastly, the possibility of the solute to have local dielectric constant larger than of a bulk water is investigated in a benchmarking test against experimentally determined set of pKa's and it is speculated that side chain rearrangements could result in local dielectric constant larger than 80.

Twenty years of gas phase structural biology

Over the past two decades, mass spectrometry (MS) of protein complexes from their native state has made inroads into structural biology. To coincide with the 20(th) anniversary of Structure, and given that it is now approximately 20 years since the first mass spectra of noncovalent protein complexes were reported, it is timely to consider progress of MS as a structural biology tool. Early reports focused on soluble complexes, contributing to ligand binding studies, subunit interaction maps, and topological models. Recent discoveries have enabled delivery of membrane complexes, encapsulated in detergent micelles, prompting new opportunities. By maintaining interactions between membrane and cytoplasmic subunits in the gas phase, it is now possible to investigate the effects of lipids, nucleotides, and drugs on intact membrane assemblies. These investigations reveal allosteric and synergistic effects of small molecule binding and expose the consequences of posttranslational modifications. In this review, we consider recent progress in the study of protein complexes, focusing particularly on complexes extracted from membranes, and outline future prospects for gas phase structural biology.