Proteins in the gas phase

Structural Elucidation of Ubiquitin via Gas-Phase Ion/Ion Cross-Linking Reactions Using Sodium-Cationized Reagents Coupled with Infrared Multiphoton Dissociation

Accurate structural determination of proteins is critical to understanding their biological functions and the impact of structural disruption on disease progression. Gas-phase cross-linking mass spectrometry (XL-MS) via ion/ion reactions between multiply charged protein cations and singly charged cross-linker anions has previously been developed to obtain low-resolution structural information on proteins. This method significantly shortens experimental time relative to conventional solution-phase XL-MS but has several technical limitations: (1) the singly deprotonated N-hydroxysulfosuccinimide (sulfo-NHS)-based cross-linker anions are restricted to attachment at neutral amine groups of basic amino acid residues and (2) analyzing terminal cross-linked fragment ions is insufficient to unambiguously localize sites of linker attachment. Herein, we demonstrate enhanced structural information for alcohol-denatured A-state ubiquitin obtained from an alternative gas-phase XL-MS approach. Briefly, singly sodiated ethylene glycol bis(sulfosuccinimidyl succinate) (sulfo-EGS) cross-linker anions enable covalent cross-linking at both ammonium and amine groups. Additionally, covalently modified internal fragment ions, along with terminal b-/y-type counterparts, improve the determination of linker attachment sites. Molecular dynamics simulations validate experimentally obtained gas-phase conformations of denatured ubiquitin. This method has identified four cross-linking sites across 8+ ubiquitin, including two new sites in the N-terminal region of the protein that were originally inaccessible in prior gas-phase XL approaches. The two N-terminal cross-linking sites suggest that the N-terminal half of ubiquitin is more compact in gas-phase conformations. By comparison, the two C-terminal linker sites indicate the signature transformation of this region of the protein from a native to a denatured conformation. Overall, the results suggest that the solution-phase secondary structures of the A-state ubiquitin are conserved in the gas phase. This method also provides sufficient sensitivity to differentiate between two gas-phase conformers of the same charge state with subtle structural variations.

Giving the Green Light to Photochemical Uncaging of Large Biomolecules in High Vacuum

The isolation of biomolecules in a high vacuum enables experiments on fragile species in the absence of a perturbing environment. Since many molecular properties are influenced by local electric fields, here we seek to gain control over the number of charges on a biopolymer by photochemical uncaging. We present the design, modeling, and synthesis of photoactive molecular tags, their labeling to peptides and proteins as well as their photochemical validation in solution and in the gas phase. The tailored tags can be selectively cleaved off at a well-defined time and without the need for any external charge-transferring agents. The energy of a single or two green photons can already trigger the process, and it is soft enough to ensure the integrity of the released biomolecular cargo. We exploit differences in the cleavage pathways in solution and in vacuum and observe a surprising robustness in upscaling the approach from a model system to genuine proteins. The interaction wavelength of 532 nm is compatible with various biomolecular entities, such as oligonucleotides or oligosaccharides.

Elucidating Structures of Protein Complexes by Collision-Induced Dissociation at Elevated Gas Pressures

Ion activation methods carried out at gas pressures compatible with ion mobility separations are not yet widely established. This limits the analytical utility of emerging tandem-ion mobility spectrometers that conduct multiple ion mobility separations in series. The present work investigates the applicability of collision-induced dissociation (CID) at 1 to 3 mbar in a tandem-trapped ion mobility spectrometer (tandem-TIMS) to study the architecture of protein complexes. We show that CID of the homotetrameric protein complexes streptavidin (53 kDa), neutravidin (60 kDa), and concanavalin A (110 kDa) provides access to all subunits of the investigated protein complexes, including structurally informative dimers. We report on an "atypical" dissociation pathway, which for concanavalin A proceeds via symmetric partitioning of the precursor charges and produces dimers with the same charge states that were previously reported from surface induced dissociation. Our data suggest a correlation between the formation of subunits by CID in tandem-TIMS/MS, their binding strengths in the native tetramer structures, and the applied activation voltage. Ion mobility spectra of in situ-generated subunits reveal a marked structural heterogeneity inconsistent with annealing into their most stable gas phase structures. Structural transitions are observed for in situ-generated subunits that resemble the transitions reported from collision-induced unfolding of natively folded proteins. These observations indicate that some aspects of the native precursor structure is preserved in the subunits generated from disassembly of the precursor complex. We rationalize our observations by an approximately 100-fold shorter activation time scale in comparison to traditional CID in a collision cell. Finally, the approach discussed here to conduct CID at elevated pressures appears generally applicable also for other types of tandem-ion mobility spectrometers.

Sample pH Can Drift during Native Mass Spectrometry Experiments: Results from Ratiometric Fluorescence Imaging

The ability of nanoelectrospray ionization (nanoESI) to generate a continuous flow of charged droplets relies on the electrolytic nature of the process. This electrochemistry can lead to the accumulation of redox products in the sample solution. This consequence can have significant implications for native mass spectrometry (MS), which aims to probe the structures and interactions of biomolecules in solution. Here, ratiometric fluorescence imaging and a pH-sensitive, fluorescent probe are used to quantify changes in solution pH during nanoESI under conditions relevant to native MS. Results show that the extent and rate of change in sample pH depends on several experimental parameters. There is a strong correlation between the extent and rate of change in solution pH and the magnitude of both the nanoESI current and electrolyte concentration. Smaller changes in solution pH are observed during experiments when a negative potential is applied than for those when a positive potential is applied. Finally, we make specific recommendations for designing native MS experiments that control for these effects.

Reconstructing the infrared spectrum of a peptide from representative conformers of the full canonical ensemble

Leucine enkephalin (LeuEnk), a biologically active endogenous opioid pentapeptide, has been under intense investigation because it is small enough to allow efficient use of sophisticated computational methods and large enough to provide insights into low-lying minima of its conformational space. Here, we reproduce and interpret experimental infrared (IR) spectra of this model peptide in gas phase using a combination of replica-exchange molecular dynamics simulations, machine learning, and ab initio calculations. In particular, we evaluate the possibility of averaging representative structural contributions to obtain an accurate computed spectrum that accounts for the corresponding canonical ensemble of the real experimental situation. Representative conformers are identified by partitioning the conformational phase space into subensembles of similar conformers. The IR contribution of each representative conformer is calculated from ab initio and weighted according to the population of each cluster. Convergence of the averaged IR signal is rationalized by merging contributions in a hierarchical clustering and the comparison to IR multiple photon dissociation experiments. The improvements achieved by decomposing clusters containing similar conformations into even smaller subensembles is strong evidence that a thorough assessment of the conformational landscape and the associated hydrogen bonding is a prerequisite for deciphering important fingerprints in experimental spectroscopic data.

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.

Tuning Electrospray Ionization with Low-Frequency Sound

Electrospray ionization (ESI) mass spectrometry (MS) is one of the key techniques used in biomolecular analysis nowadays. It relies on formation of polydisperse microdroplets, which undergo desolvation and liberate ions to the gas phase. Here we demonstrate low-frequency-sound-modulated ESI for analysis of biomolecules. By using a low-frequency (50-350 Hz) sound, it is possible to deflect electrospray microdroplets toward the mass spectrometer's orifice. Microdroplets of different sizes are deflected to a different extent leading to a partial size segregation. This effect leads to either an increase or decrease of MS signal intensity as well as signal-to-noise ratio. It also affects the selectivity of the ESI-MS analysis. The observations are rationalized by taking into account different pathways of ion formation and the likelihood of deflecting microdroplets of certain size. The online ESI-MS observations are supported with offline shadowgraphs obtained at varied sound frequencies, signal amplitudes, and phase shifts.

Interfacial Water in the SARS Spike Protein: Investigating the Interaction with Human ACE2 Receptor and In Vitro Uptake in A549 Cells

The severity of global pandemic due to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has engaged the researchers and clinicians to find the key features triggering the viral infection to lung cells. By utilizing such crucial information, researchers and scientists try to combat the spread of the virus. Here, in this work, we performed in silico analysis of the protein-protein interactions between the receptor-binding domain (RBD) of the viral spike protein and the human angiotensin-converting enzyme 2 (hACE2) receptor to highlight the key alteration that happened from SARS-CoV to SARS-CoV-2. We analyzed and compared the molecular differences between spike proteins of the two viruses using various computational approaches such as binding affinity calculations, computational alanine, and molecular dynamics simulations. The binding affinity calculations showed that SARS-CoV-2 binds a little more firmly to the hACE2 receptor than SARS-CoV. The major finding obtained from molecular dynamics simulations was that the RBD-ACE2 interface is populated with water molecules and interacts strongly with both RBD and ACE2 interfacial residues during the simulation periods. The water-mediated hydrogen bond by the bridge water molecules is crucial for stabilizing the RBD and ACE2 domains. Near-ambient pressure X-ray photoelectron spectroscopy (NAP-XPS) confirmed the presence of vapor and molecular water phases in the protein-protein interfacial domain, further validating the computationally predicted interfacial water molecules. In addition, we examined the role of interfacial water molecules in virus uptake by lung cell A549 by binding and maintaining the RBD/hACE2 complex at varying temperatures using nanourchins coated with spike proteins as pseudoviruses and fluorescence-activated cell sorting (FACS) as a quantitative approach. The structural and dynamical features presented here may serve as a guide for developing new drug molecules, vaccines, or antibodies to combat the COVID-19 pandemic.

Protein orientation in time-dependent electric fields: orientation before destruction

Proteins often have nonzero electric dipole moments, making them interact with external electric fields and offering a means for controlling their orientation. One application that is known to benefit from orientation control is single-particle imaging with x-ray free-electron lasers, in which diffraction is recorded from proteins in the gas phase to determine their structures. To this point, theoretical investigations into this phenomenon have assumed that the field experienced by the proteins is constant or a perfect step function, whereas any real-world pulse will be smooth. Here, we explore the possibility of orienting gas-phase proteins using time-dependent electric fields. We performed ab initio simulations to estimate the field strength required to break protein bonds, with 45 V/nm as a breaking point value. We then simulated ubiquitin in time-dependent electric fields using classical molecular dynamics. The minimal field strength required for orientation within 10 ns was on the order of 0.5 V/nm. Although high fields can be destructive for the structure, the structures in our simulations were preserved until orientation was achieved regardless of field strength, a principle we denote “orientation before destruction.”

Side-Chain Polarity Modulates the Intrinsic Conformational Landscape of Model Dipeptides

The intrinsic conformational preferences of small peptides may provide additional insight into the thermodynamics and kinetics of protein folding. In this study, we explore the underlying energy landscapes of two model peptides, namely, Ac-Ala-NH₂ and Ac-Ser-NH₂, using geometry-optimization-based tools developed within the context of energy landscape theory. We analyze not only how side-chain polarity influences the structural preferences of the dipeptides, but also other emergent properties of the landscape, including heat capacity profiles, and kinetics of conformational rearrangements. The contrasting topographies of the free energy landscape agree with recent results from Fourier transform microwave spectroscopy experiments, where Ac-Ala-NH₂ was found to exist as a mixture of two conformers, while Ac-Ser-NH₂ remained structurally locked, despite exhibiting an apparently rich conformational landscape.

Water-Mediated Interactions Determine Helix Formation of Peptides in Open Nanotubes

Water-mediated interactions (WMIs) play diverse roles in molecular biology. They are particularly relevant in geometrically confined spaces such as the interior of the chaperonin, at the interface between ligands and their binding partners, and in the ribosome tunnel. Inspired in part by the geometry of the ribosome tunnel, we consider confinement effects on the stability of peptides. We describe results from replica exchange molecular dynamics simulations of a system containing a 23-alanine or 23-serine polypeptide confined to nonpolar and polar nanotubes in the gas phase and when open to a water reservoir. We quantify the effect of water in determining the preferred conformational states of these polypeptides by calculating the difference in the solvation free energy for the helix and coil states in the open nanotube in the two phases. Our simulations reveal several possibilities. We find that nanoscopic confinement preferentially stabilizes the helical state of polypeptides with hydrophobic side chains, which is explained by the entropic stabilization mechanism proposed on the basis of polymer physics. Polypeptide chains with hydrophilic side chains can adopt helical structures within nanotubes, but helix formation is sensitive to the nature of the nanotube due to WMIs. We elaborate on the potential implications of our findings to the stability of peptides in the ribosome tunnel.

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.

Comparing the Conformational Stability of Pyruvate Kinase in the Gas Phase and in Solution

Collision induced unfolding (CIU) is increasingly used to characterize protein complexes in the gas phase and is often employed to detect ligand binding-induced conformational stabilization. However, the extent to which gas-phase conformational stabilities measured by CIU reflect analogous parameters in solution is not yet clear, particularly for systems where conformational and protein complex stability are modulated by point mutation. Here, we compare CIU-derived relative stabilities of four point mutants of the homotetramer pyruvate kinase to solution stabilities measured by differential scanning fluorimetry (DSF) and solution conformational dynamics measured by time-resolved electrospray ionization hydrogen-deuterium exchange (TRESI-HDX). Our results demonstrate that both destabilization of the tetrameric state and generally reduced conformational stability of the monomer in solution are well correlated to lower onset energies for specific unfolding transitions observed in CIU. However, this correlation not fully retained when comparing CIU to HDX data, where the latter measurement is strongly impacted by conformational dynamics within the tetramer.

The Fragment Molecular Orbital Method Based on Long-Range Corrected Density-Functional Tight-Binding

The presently available linear scaling approaches to density-functional tight-binding (DFTB) based on the fragment molecular orbital (FMO) method are severely impacted by the problem of artificial charge transfer due to the self-interaction error (SIE), which hampers the simulation of zwitterionic systems such as biopolymers or ionic liquids. Here we report an extension of FMO-DFTB where we included a long-range corrected (LC) functional designed to mitigate the DFTB SIE, called the FMO-LC-DFTB method, resulting in a robust method which succeeds in simulating zwitterionic systems. Both energy and analytic gradient are developed for the gas phase and the polarizable continuum model of solvation. The scaling of FMO-LC-DFTB with system size N is shown to be almost linear, O( N^1.13-1.28), and its numerical accuracy is established for a variety of representative systems including neutral and charged polypeptides. It is shown that pair interaction energies between fragments for two mini-proteins are in excellent agreement with results from long-range corrected density functional theory. The new method was employed in long time scale (1 ns) molecular dynamics simulations of the tryptophan cage protein (PDB: 1L2Y ) in the gas phase for four different protonation states and in stochastic global minimum structure searches for 1-ethyl-3-methylimidazolium nitrate ionic liquid clusters containing up to 2300 atoms.

Initial Protein Unfolding Events in Ubiquitin, Cytochrome c and Myoglobin Are Revealed with the Use of 213 nm UVPD Coupled to IM-MS

The initial stages of protein unfolding may reflect the stability of the entire fold and can also reveal which parts of a protein can be perturbed, without restructuring the rest. In this work, we couple UVPD with activated ion mobility mass spectrometry to measure how three model proteins start to unfold. Ubiquitin, cytochrome c and myoglobin ions produced via nESI from salty solutions are subjected to UV irradiation pre-mobility separation; experiments are conducted with a range of source conditions which alter the conformation of the precursor ion as shown by the drift time profiles. For all three proteins, the compact structures result in less fragmentation than more extended structures which emerge following progressive in-source activation. Cleavage sites are found to differ between conformational ensembles, for example, for the dominant charge state of cytochrome c [M + 7H]⁷⁺, cleavage at Phe10, Thr19 and Val20 was only observed in activating conditions whilst cleavage at Ala43 is dramatically enhanced. Mapping the photo-cleaved fragments onto crystallographic structures provides insight into the local structural changes that occur as protein unfolding progresses, which is coupled to global restructuring observed in the drift time profiles. Graphical Abstract.

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

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

Analysis of Proteins, Protein Complexes, and Organellar Proteomes Using Sheathless Capillary Zone Electrophoresis - Native Mass Spectrometry

Native mass spectrometry (MS) is a rapidly advancing field in the analysis of proteins, protein complexes, and macromolecular species of various types. The majority of native MS experiments reported to-date has been conducted using direct infusion of purified analytes into a mass spectrometer. In this study, capillary zone electrophoresis (CZE) was coupled online to Orbitrap mass spectrometers using a commercial sheathless interface to enable high-performance separation, identification, and structural characterization of limited amounts of purified proteins and protein complexes, the latter with preserved non-covalent associations under native conditions. The performance of both bare-fused silica and polyacrylamide-coated capillaries was assessed using mixtures of protein standards known to form non-covalent protein-protein and protein-ligand complexes. High-efficiency separation of native complexes is demonstrated using both capillary types, while the polyacrylamide neutral-coated capillary showed better reproducibility and higher efficiency for more complex samples. The platform was then evaluated for the determination of monoclonal antibody aggregation and for analysis of proteomes of limited complexity using a ribosomal isolate from E. coli. Native CZE-MS, using accurate single stage and tandem-MS measurements, enabled identification of proteoforms and non-covalent complexes at femtomole levels. This study demonstrates that native CZE-MS can serve as an orthogonal and complementary technique to conventional native MS methodologies with the advantages of low sample consumption, minimal sample processing and losses, and high throughput and sensitivity. This study presents a novel platform for analysis of ribosomes and other macromolecular complexes and organelles, with the potential for discovery of novel structural features defining cellular phenotypes (e.g., specialized ribosomes). Graphical Abstract ᅟ.

Quantum-Assisted Metrology of Neutral Vitamins in the Gas Phase

It has recently been shown that matter-wave interferometry can be used to imprint a periodic nanostructure onto a molecular beam, which provides a highly sensitive tool for beam displacement measurements. Herein, we used this feature to measure electronic properties of provitamin A, vitamin E, and vitamin K1 in the gas phase for the first time. The shift of the matter-wave fringes in a static electric field encodes the molecular susceptibility and the time-averaged dynamic electric dipole moment. The dependence of the fringe pattern on the intensity of the central light-wave diffraction grating was used to determine the molecular optical polarizability. Comparison of our experimental findings with molecular dynamics simulations and density functional theory provides a rich picture of the electronic structures and dynamics of these biomolecules in the gas phase with β-carotene as a particularly interesting example.

The Solution Assembly of Biological Molecules Using Ion Mobility Methods: From Amino Acids to Amyloid β-Protein

Ion mobility spectrometry-mass spectrometry (IMS-MS) methods are increasingly used to study noncovalent assemblies of peptides and proteins. This review focuses on the noncovalent self-assembly of amino acids and peptides, systems at the heart of the amyloid process that play a central role in a number of devastating diseases. Three different systems are discussed in detail: the 42-residue peptide amyloid-β42 implicated in the etiology of Alzheimer's disease, several amyloid-forming peptides with 6-11 residues, and the assembly of individual amino acids. We also discuss from a more fundamental perspective the processes that determine how quickly proteins and their assemblies denature when the analyte ion has been stripped of its solvent in an IMS-MS measurement and how to soften the measurement so that biologically meaningful data can be recorded.

Effects of Solution Structure on the Folding of Lysozyme Ions in the Gas Phase

The fidelity between the structures of proteins in solution and protein ions in the gas phase is critical to experiments that use gas-phase measurements to infer structures in solution. Here we generate ions of lysozyme, a 129-residue protein whose native tertiary structure contains four internal disulfide bonds, from three solutions that preserve varying extents of the original native structure. We then use cation-to-anion proton-transfer reactions (CAPTR) to reduce the charge states of those ions in the gas phase and ion mobility to probe their structures. The collision cross section (Ω) distributions of each CAPTR product depends to varying extents on the original solution, the charge state of the product, and the charge state of the precursor. For example, the Ω distributions of the 6+ ions depend strongly on the original solutions conditions and to a lesser extent on the charge state of the precursor. Energy-dependent experiments suggest that very different structures are accessible to disulfide-reduced and disulfide-intact ions, but similar Ω distributions are formed at high energy for disulfide-intact ions from denaturing and from aqueous conditions. The Ω distributions of the 3+ ions are all similar but exhibit subtle differences that depend more strongly on the original solutions conditions than other factors. More generally, these results suggest that specific CAPTR products may be especially sensitive to specific elements of structure in solution.

Implementation and benchmark of a long-range corrected functional in the density functional based tight-binding method

Bridging the gap between first principles methods and empirical schemes, the density functional based tight-binding method (DFTB) has become a versatile tool in predictive atomistic simulations over the past years. One of the major restrictions of this method is the limitation to local or gradient corrected exchange-correlation functionals. This excludes the important class of hybrid or long-range corrected functionals, which are advantageous in thermochemistry, as well as in the computation of vibrational, photoelectron, and optical spectra. The present work provides a detailed account of the implementation of DFTB for a long-range corrected functional in generalized Kohn-Sham theory. We apply the method to a set of organic molecules and compare ionization potentials and electron affinities with the original DFTB method and higher level theory. The new scheme cures the significant overpolarization in electric fields found for local DFTB, which parallels the functional dependence in first principles density functional theory (DFT). At the same time, the computational savings with respect to full DFT calculations are not compromised as evidenced by numerical benchmark data.

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.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

SDS-PAGE analysis of Aβ oligomers is disserving research into Alzheimer´s disease: appealing for ESI-IM-MS

The characterization of amyloid-beta peptide (Aβ) oligomer forms and structures is crucial to the advancement in the field of Alzheimer´s disease (AD). Here we report a critical evaluation of two methods used for this purpose, namely sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), extensively used in the field, and ion mobility coupled to electrospray ionization mass spectrometry (ESI-IM-MS), an emerging technique with great potential for oligomer characterization. To evaluate their performance, we first obtained pure cross-linked Aβ40 and Aβ42 oligomers of well-defined order. Analysis of these samples by SDS-PAGE revealed that SDS affects the oligomerization state of Aβ42 oligomers, thus providing flawed information on their order and distribution. In contrast, ESI-IM-MS provided accurate information, while also reported on the chemical nature and on the structure of the oligomers. Our findings have important implications as they challenge scientific paradigms in the AD field built upon SDS-PAGE characterization of Aβ oligomer samples.

Linking molecular models with ion mobility experiments. Illustration with a rigid nucleic acid structure

Ion mobility spectrometry experiments allow the mass spectrometrist to determine an ion's rotationally averaged collision cross section Ω(EXP). Molecular modelling is used to visualize what ion three-dimensional structure(s) is(are) compatible with the experiment. The collision cross sections of candidate molecular models have to be calculated, and the resulting Ω(CALC) are compared with the experimental data. Researchers who want to apply this strategy to a new type of molecule face many questions: (1) What experimental error is associated with Ω(EXP) determination, and how to estimate it (in particular when using a calibration for traveling wave ion guides)? (2) How to generate plausible 3D models in the gas phase? (3) Different collision cross section calculation models exist, which have been developed for other analytes than mine. Which one(s) can I apply to my systems? To apply ion mobility spectrometry to nucleic acid structural characterization, we explored each of these questions using a rigid structure which we know is preserved in the gas phase: the tetramolecular G-quadruplex [dTGGGGT]4, and we will present these detailed investigation in this tutorial.

Microheterogeneity within conformational states of ubiquitin revealed by high resolution trapped ion mobility spectrometry

The present work employs trapped ion mobility spectrometry (TIMS) for the analysis of ubiquitin ions known to display a multitude of previously unresolved interchangeable conformations upon electrospray ionization.

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.

Proteins with highly similar native folds can show vastly dissimilar folding behavior when desolvated

Proteins can be exposed to vastly different environments such as the cytosol or membranes, but the delicate balance between external factors and intrinsic determinants of protein structure, stability, and folding is only poorly understood. Here we used electron capture dissociation to study horse and tuna heart Cytochromes c in the complete absence of solvent. The significantly different stability of their highly similar native folds after transfer into the gas phase, and their strikingly different folding behavior in the gas phase, can be rationalized on the basis of electrostatic interactions such as salt bridges. In the absence of hydrophobic bonding, protein folding is far slower and more complex than in solution.