Protein Structures under Electrospray Conditions

Characterization of Anticancer Drug Protomers Using Electrospray Ionization and Ion Mobility Spectrometry-Mass Spectrometry

We used electrospray ionization and ion mobility spectrometry-mass spectrometry to detect and characterize the three anticancer drugs palbociclib, copanlisib, and olaparib. Ion mobility-mass spectrometry and density functional theory revealed that these compounds generate isomers during ionization (protomers) due to the presence of multiple protonation sites within their chemical structures. Our work has implications for understanding the solution- and gas-phase chemistry of these molecules during spray-based ionization processes.

High-Performance Molecular Dynamics Simulations for Native Mass Spectrometry of Large Protein Complexes with the Fast Multipole Method

Native mass spectrometry (MS) is widely employed to study the structures and assemblies of proteins ranging from small monomers to megadalton complexes. Molecular dynamics (MD) simulation is a useful complement as it provides the spatial detail that native MS cannot offer. However, MD simulations performed in the gas phase have suffered from rapidly increasing computational costs with the system size. The primary bottleneck is the calculation of electrostatic forces, which are effective over long distances and must be explicitly computed for each atom pair, precluding efficient use of methods traditionally used to accelerate condensed-phase simulations. As a result, MD simulations have been unable to match the capacity of MS in probing large multimeric protein complexes. Here, we apply the fast multipole method (FMM) for computing the electrostatic forces, recently implemented by Kohnke et al. (J. Chem. Theory Comput., 2020, 16, 6938-6949), showing that it significantly enhances the performance of gas-phase simulations of large proteins. We assess how to achieve adequate accuracy and optimal performance with FMM, finding that it expands the accessible size range and time scales dramatically. Additionally, we simulate a 460 kDa ferritin complex over microsecond time scales, alongside complementary ion mobility (IM)-MS experiments, uncovering conformational changes that are not apparent from the IM-MS data alone.

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.

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.

X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals

The MD simulation package Amber offers an attractive platform to refine crystallographic structures of proteins: (i) state-of-the-art force fields help to regularize protein coordinates and reconstruct the poorly diffracting elements of the structure, such as flexible loops; (ii) MD simulations restrained by the experimental diffraction data provide an effective strategy to optimize structural models of protein crystals, including explicitly modeled interstitial solvent as well as crystal contacts. Here, we present the new crystallography module xray, released as a part of the Amber 2023 package. This module contains functions to calculate and scale structure factors (including the contributions from bulk solvent), evaluate the maximum-likelihood-type crystallographic potential, and compute its derivative forces. The X-ray functionality of Amber no longer relies on external dependencies so that the full advantage of GPU acceleration can be taken. This makes it possible to refine in a short time hundreds of crystal models, including supercell models comprised of multiple unit cells. The new automated Amber-based refinement procedure leads to an appreciable improvement in Rfree (in some cases, by as much as 0.067) as well as MolProbity scores.

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.

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.

Mass-selective and ice-free electron cryomicroscopy protein sample preparation via native electrospray ion-beam deposition

Despite tremendous advances in sample preparation and classification algorithms for electron cryomicroscopy (cryo-EM) and single-particle analysis (SPA), sample heterogeneity remains a major challenge and can prevent access to high-resolution structures. In addition, optimization of preparation conditions for a given sample can be time-consuming. In the current work, it is demonstrated that native electrospray ion-beam deposition (native ES-IBD) is an alternative, reliable approach for the preparation of extremely high-purity samples, based on mass selection in vacuum. Folded protein ions are generated by native electrospray ionization, separated from other proteins, contaminants, aggregates, and fragments, gently deposited on cryo-EM grids, frozen in liquid nitrogen, and subsequently imaged by cryo-EM. We demonstrate homogeneous coverage of ice-free cryo-EM grids with mass-selected protein complexes. SPA reveals that the complexes remain folded and assembled, but variations in secondary and tertiary structures are currently limiting information in 2D classes and 3D EM density maps. We identify and discuss challenges that need to be addressed to obtain a resolution comparable to that of the established cryo-EM workflow. Our results show the potential of native ES-IBD to increase the scope and throughput of cryo-EM for protein structure determination and provide an essential link between gas-phase and solution-phase protein structures.

Tracking local and global structural changes in a protein by cold ion spectroscopy

Characterization of native structures of proteins in the gas phase remains challenging due to the unpredictable conformational changes the molecules undergo during desolvation and ionization. We spectroscopically studied cryogenically cooled protonated protein ubiquitin and its microhydrated complexes prepared in the gas phase in a range of charge states under different ionization conditions. The UV spectra appear vibrationally resolved for the unfolded protein, but become redshifted and smooth for the native-like structures of ubiquitin. This spectroscopic change results from the H-bonding of the hydroxyl of Tyr to the amide group of Glu-51 in the compact structures; the minimum length of this bond was estimated to be ∼1.7 Å. IR spectroscopy reflects the global structural change by observing redshifts of free NH/OH-stretch vibrational transitions. Evaporative cooling of microhydrated complexes of ubiquitin keeps the protein chilly during ionization, enabling native-like conformers with up to eight protons to survive in the gas phase.

Modeling a unit cell: crystallographic refinement procedure using the biomolecular MD simulation platform Amber

A procedure has been developed for the refinement of crystallographic protein structures based on the biomolecular simulation program Amber. The procedure constructs a model representing a crystal unit cell, which generally contains multiple protein molecules and is fully hydrated with TIP3P water. Periodic boundary conditions are applied to the cell in order to emulate the crystal lattice. The refinement is conducted in the form of a specially designed short molecular-dynamics run controlled by the Amber ff14SB force field and the maximum-likelihood potential that encodes the structure-factor-based restraints. The new Amber-based refinement procedure has been tested on a set of 84 protein structures. In most cases, the new procedure led to appreciably lower R free values compared with those reported in the original PDB depositions or obtained by means of the industry-standard phenix.refine program. In particular, the new method has the edge in refining low-accuracy scrambled models. It has also been successful in refining a number of molecular-replacement models, including one with an r.m.s.d. of 2.15 Å. In addition, Amber-refined structures consistently show superior MolProbity scores. The new approach offers a highly realistic representation of protein-protein interactions in the crystal, as well as of protein-water interactions. It also offers a realistic representation of protein crystal dynamics (akin to ensemble-refinement schemes). Importantly, the method fully utilizes the information from the available diffraction data, while relying on state-of-the-art molecular-dynamics modeling to assist with those elements of the structure that do not diffract well (for example mobile loops or side chains). Finally, it should be noted that the protocol employs no tunable parameters, and the calculations can be conducted in a matter of several hours on desktop computers equipped with graphical processing units or using a designated web service.

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.”

Interrogating Membrane Protein Structure and Lipid Interactions by Native Mass Spectrometry

Native mass spectrometry and native ion mobility mass spectrometry are now established techniques in structural biology, with recent work developing these methods for the study of integral membrane proteins reconstituted in both lipid bilayer and detergent environments. Here we show how native mass spectrometry can be used to interrogate integral membrane proteins, providing insights into conformation, oligomerization, subunit composition/stoichiometry, and interactions with detergents/lipids/drugs. Furthermore, we discuss the sample requirements and experimental considerations unique to integral membrane protein native mass spectrometry research.

Sulfolane-Induced Supercharging of Electrosprayed Salt Clusters: An Experimental/Computational Perspective

It is well-known that supercharging agents (SCAs) such as sulfolane enhance the electrospray ionization (ESI) charge states of proteins, although the mechanistic origins of this effect remain contentious. Only very few studies have explored SCA effects on analytes other than proteins or peptides. This work examines how sulfolane affects electrosprayed NaI salt clusters. Such alkali metal halide clusters have played a key role for earlier ESI mechanistic studies, making them interesting targets for supercharging investigations. ESI of aqueous NaI solutions predominantly generated singly charged [Na_nI_(n-1)]⁺ clusters. The addition of sulfolane resulted in abundant doubly charged [Na_nI_(n-2)Sulfolane_s]²⁺ species. These experimental data for the first time demonstrate that electrosprayed salt clusters can undergo supercharging. Molecular dynamics (MD) simulations of aqueous ESI nanodroplets containing Na⁺/I^- with and without sulfolane were conducted to obtain atomistic insights into the supercharging mechanism. The simulations produced [Na_nI_i]^z+ and [Na_nI_iSulfolane_s]^z+ clusters similar to those observed experimentally. The MD trajectories demonstrated that these clusters were released into the gas phase upon droplet evaporation to dryness, in line with the charged residue model. Sulfolane was found to evaporate much more slowly than water. This slow evaporation, in conjunction with the large dipole moment of sulfolane, resulted in electrostatic stabilization of the shrinking ESI droplets and the final clusters. Hence, charge-dipole stabilization causes the sulfolane-containing droplets and clusters to retain more charge, thereby providing the mechanistic foundation of salt cluster supercharging.

Reproducibility in the unfolding process of protein induced by an external electric field

The dynamics of proteins are crucial for their function. However, commonly used techniques for studying protein structures are limited in monitoring time-resolved dynamics at high resolution. Combining electric fields with existing techniques to study gas-phase proteins, such as single particle imaging using free-electron lasers and gas-phase small angle X-ray scattering, has the potential to open up a new era in time-resolved studies of gas-phase protein dynamics. Using molecular dynamics simulations, we identify well-defined unfolding pathways of a protein, induced by experimentally achievable external electric fields. Our simulations show that strong electric fields in conjunction with short-pulsed X-ray sources such as free-electron lasers can be a new path for imaging dynamics of gas-phase proteins at high spatial and temporal resolution.

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.

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

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

Is radiation damage the limiting factor in high-resolution single particle imaging with X-ray free-electron lasers?

The prospect of single particle imaging with atomic resolution is one of the scientific drivers for the development of X-ray free-electron lasers. The assumption since the beginning has been that damage to the sample caused by intense X-ray pulses is one of the limiting factors for achieving subnanometer X-ray imaging of single particles and that X-ray pulses need to be as short as possible. Based on the molecular dynamics simulations of proteins in X-ray fields of various durations (5 fs, 25 fs, and 50 fs), we show that the noise in the diffracted signal caused by radiation damage is less than what can be expected from other sources, such as sample inhomogeneity and X-ray shot-to-shot variations. These findings show a different aspect of the feasibility of high-resolution single particle imaging using free-electron lasers, where employing X-ray pulses of longer durations could still provide a useful diffraction signal above the noise due to the Coulomb explosion.

Preservation of Protein Zwitterionic States in the Transition from Solution to Gas Phase Revealed by Sodium Adduction Mass Spectrometry

The structural characterization of proteins and their interaction network mapping in the gas phase highlights the need to preserve their most nativelike conformers in the transition from the solution to gas phase. Zwitterionic interactions in a protein are weak bonds between oppositely charged residues, which make an important contribution to protein stability. However, it is still not clear whether the native zwitterionic states of proteins can be retained or not when it is transferred from the solution to gas phase. Using the nonspecific Na⁺ adduction as a novel signature, here we show that the zwitterionic states of proteins can be preserved when a moderated droplet desolvation condition (temperature <30 °C) is used in native electrospray ionization mass spectrometry. The very low-level nonspecific metal adduction to proteins under such conditions also enables rapid and direct determination of the binding states of metal-binding proteins and sensitive detection of proteins from solutions containing highly concentrated involatile salts (e.g., 50 mM NaCl). We believe that our findings can be instructive for performing mass spectrometric analysis of proteins and useful for protein ions desalting which simply involves altering the temperature and flow rate of drying gas in the desolvation region.

The effects of solution additives and gas-phase modifiers on the molecular environment and conformational space of common heme proteins

RATIONALE

The molecular environment is known to impact the secondary and tertiary structures of biomolecules both in solution and in the gas phase, shifting the equilibrium between different conformational and oligomerization states. However, there is a lack of studies monitoring the impacts of solution additives and gas-phase modifiers on biomolecules characterized using ion mobility techniques.

METHODS

The effect of solution additives and gas-phase modifiers on the molecular environment of two common heme proteins, bovine cytochrome c and equine myoglobin, is investigated as a function of the time after desolvation (e.g., 100-500 ms) using nanoelectrospray ionization coupled to trapped ion mobility spectrometry with detection by time-of-flight mass spectrometry. Organic compounds used as additives/modifiers (methanol, acetonitrile, acetone) were either added to the aqueous protein solution before ionization or added to the ion mobility bath gas by nebulization.

RESULTS

Changes in the mobility profiles are observed depending on the starting solution composition (i.e., in aqueous solution at neutral pH or in the presence of organic content: methanol, acetone, or acetonitrile) and the protein. In the presence of gas-phase modifiers (i.e., N₂ doped with methanol, acetone, or acetonitrile), a shift in the mobility profiles driven by the gas-modifier mass and size and changes in the relative abundances and number of IMS bands are observed.

CONCLUSIONS

We attribute the observed changes in the mobility profiles in the presence of gas-phase modifiers to a clustering/declustering mechanism by which organic molecules adsorb to the protein ion surface and lower energetic barriers for interconversion between conformational states, thus redefining the free energy landscape and equilibria between conformers. These structural biology experiments open new avenues for manipulation and interrogation of biomolecules in the gas phase with the potential to emulate a large suite of solution conditions, ultimately including conditions that more accurately reflect a variety of intracellular environments.

Ion Mobility Mass Spectrometry Uncovers the Impact of the Patterning of Oppositely Charged Residues on the Conformational Distributions of Intrinsically Disordered Proteins

The global dimensions and amplitudes of conformational fluctuations of intrinsically disordered proteins are governed, in part, by the linear segregation versus clustering of oppositely charged residues within the primary sequence. Ion mobility-mass spectrometry (IM-MS) affords unique advantages for probing the conformational consequences of the linear patterning of oppositely charged residues because it measures and separates proteins electrosprayed from solution on the basis of charge and shape. Here, we use IM-MS to measure the conformational consequences of charge patterning on the C-terminal intrinsically disordered region (p27 IDR) of the cell cycle inhibitory protein p27^Kip1. We report the range of charge states and accompanying collisional cross section distributions for wild-type p27 IDR and two variants with identical amino acid compositions, κ14 and κ56, distinguished by the extent of linear mixing versus segregation of oppositely charged residues. Wild-type p27 IDR (κ31) and κ14, where the oppositely charged residues are more evenly distributed, exhibit a broad distribution of charge states. This is concordant with high degrees of conformational heterogeneity in solution. By contrast, κ56 with linear segregation of oppositely charged residues leads to limited conformational heterogeneity and a narrow distribution of charged states. Gas-phase molecular dynamics simulations demonstrate that the interplay between chain solvation and intrachain interactions (self-solvation) leads to conformational distributions that are modulated by salt concentration, with the wild-type sequence showing the most sensitivity to changes in salt concentration. These results suggest that the charge patterning within the wild-type p27 IDR may be optimized to sample both highly solvated and self-solvated conformational states.

Charging and supercharging of proteins for mass spectrometry: recent insights into the mechanisms of electrospray ionization

Molecular dynamics simulations have uncovered mechanistic details of the protein ESI process under various experimental conditions.

Reproducibility of single protein explosions induced by X-ray lasers

Single particle imaging (SPI) using X-ray pulses has become increasingly attainable with the advent of high-intensity free electron lasers. Eliminating the need for crystallized samples enables structural studies of molecules previously inaccessible by conventional crystallography. While this emerging technique already demonstrates substantial promise, some obstacles need to be overcome before SPI can reach its full potential. One such problem is determining the spatial orientation of the sample at the time of X-ray interaction. Existing solutions rely on diffraction data and are computationally demanding and sensitive to noise. In this in silico study, we explore the possibility of aiding these methods by mapping the ion distribution as the sample undergoes a Coulomb explosion following the intense ionization. By detecting the ions ejected from the fragmented sample, the orientation of the original sample should be possible to determine. Knowledge of the orientation has been shown earlier to be of substantial advantage in the reconstruction of the original structure. 150 explosions of each of twelve separate systems - four polypeptides with different amounts of surface bound water - were simulated with molecular dynamics (MD) and the average angular distribution of carbon and sulfur ions was investigated independently. The results show that the explosion maps are reproducible in both cases, supporting the idea that orientation information is preserved. Additional water seems to restrict the carbon ion trajectories further through a shielding mechanism, making the maps more distinct. For sulfurs, water has no significant impact on the trajectories, likely due to their higher mass and greater ionization cross section, indicating that they could be of particular interest. Based on these findings, we conclude that explosion data can aid spatial orientation in SPI experiments and could substantially improve the capabilities of the novel technique.

Charging and Release Mechanisms of Flexible Macromolecules in Droplets

We study systematically the charging and release mechanisms of a flexible macromolecule, modeled by poly(ethylene glycol) (PEG), in a droplet by using molecular dynamics simulations. We compare how PEG is solvated and charged by sodium Na⁺ ions in a droplet of water (H₂O), acetonitrile (MeCN), and their mixtures. Initially, we examine the location and the conformation of the macromolecule in a droplet bearing no net charge. It is revealed that the presence of charge carriers do not affect the location of PEG in aqueous and MeCN droplets compared with that in the neutral droplets, but the location of the macromolecule and the droplet size do affect the PEG conformation. PEG is charged on the surface of a sodiated aqueous droplet that is found close to the Rayleigh limit. Its charging is coupled to the extrusion mechanism, where PEG segments leave the droplet once they coordinate a Na⁺ ion or in a correlated motion with Na⁺ ions. In contrast, as PEG resides in the interior of a MeCN droplet, it is sodiated inside the droplet. The compact macro-ion transitions through partially unwound states to an extended conformation, a process occurring during the final stage of desolvation and in the presence of only a handful of MeCN molecules. For charged H₂O/MeCN droplets, the sodiation of PEG is determined by the H₂O component, reflecting its slower evaporation and preference over MeCN for solvating Na⁺ ions. We use the simulation data to construct an analytical model that suggests that the droplet surface electric field may play a role in the macro-ion-droplet interactions that lead to the extrusion of the macro-ion. This study provides the first evidence of the effect of the surface electric field by using atomistic simulations. Graphical Abstract ᅟ.

Controlling Protein Orientation in Vacuum Using Electric Fields

Single-particle imaging using X-ray free-electron lasers is an emerging technique that could provide high-resolution structures of macromolecules in the gas phase. One of the largest difficulties in realizing this goal is the unknown orientation of the individual sample molecules at the time of exposure. Preorientation of the molecules has been identified as a possible solution to this problem. Using molecular dynamics simulations, we identify a range of electric field strengths where proteins become oriented without losing their structure. For a number of experimentally relevant cases we show that structure determination is possible only when orientation information is included in the orientation-recovery process. We conclude that nondestructive field orientation of intact proteins is feasible and that it enables a range of new structural investigations with single-particle imaging.

Critical Evaluation of Native Electrospray Ionization Mass Spectrometry for Fragment-Based Screening

Fragment-based screening presents a promising alternative to high-throughput screening and has gained great attention in recent years. So far, only a few studies have discussed mass spectrometry as a screening technology for fragments. Herein, we report the application of native electrospray ionization mass spectrometry (MS) for screening defined sets of fragments against four different target proteins. Fragments were selected from a primary screening conducted with a thermal shift assay (TSA) and represented different binding categories. Our data indicated that, beside specific complex formation, many fragments show extensive multiple binding and also charge-state shifts. Both of these factors complicate automated data analysis and decrease the attractiveness of native MS as a primary screening tool for fragments. A comparison of the hits identified by native MS and TSA showed good agreement for two of the proteins. Furthermore, we discuss general challenges, including the determination of an optimal fragment concentration and the question of how to rank fragment hits according to their affinity. In conclusion, we consider native MS to be a highly valuable tool for the validation and deeper investigation of promising fragment hits rather than a method for primary screening.

Proton Dynamics in Protein Mass Spectrometry

Native electrospray ionization/ion mobility-mass spectrometry (ESI/IM-MS) allows an accurate determination of low-resolution structural features of proteins. Yet, the presence of proton dynamics, observed already by us for DNA in the gas phase, and its impact on protein structural determinants, have not been investigated so far. Here, we address this issue by a multistep simulation strategy on a pharmacologically relevant peptide, the N-terminal residues of amyloid-β peptide (Aβ(1-16)). Our calculations reproduce the experimental maximum charge state from ESI-MS and are also in fair agreement with collision cross section (CCS) data measured here by ESI/IM-MS. Although the main structural features are preserved, subtle conformational changes do take place in the first ∼0.1 ms of dynamics. In addition, intramolecular proton dynamics processes occur on the picosecond-time scale in the gas phase as emerging from quantum mechanics/molecular mechanics (QM/MM) simulations at the B3LYP level of theory. We conclude that proton transfer phenomena do occur frequently during fly time in ESI-MS experiments (typically on the millisecond time scale). However, the structural changes associated with the process do not significantly affect the structural determinants.

Terahertz Frequency Spectroscopy to Determine Cold Shock Protein Stability upon Solvation and Evaporation - A Molecular Dynamics Study

Infrared (IR) and Terahertz (THz) spectroscopy simulations were carried out using CHARMM35b2 to determine protein stability. The stabilities of three bacterial cold shock proteins (Csps) originating from mesophiles, thermophiles and hyper- thermophiles respectively were investigated in this study. The three different Csps were investigated by Normal-Mode analysis and Molecular Dynamics simulation of THz spectra using the Hessian matrix for solvated systems, interpreted in the harmonic approximation at optimum near-melting temperatures of each homologue, by incorporating differences in the hydrous and anhydrous states of the Csps. The results show slight variations in the large scale protein motion. However, the IR spectra of Csps observed at the low frequency saddle surface region, clearly distinguishes the thermophilic and mesophilic proteins based on their stability. Further studies on protein stability employing low-frequency collective modes have the potential to reveal functionally important conformational changes that are biologically significant.

The fate of proteins in outer space

It is well established that any properly conducted biophysical studies of proteins must take appropriate account of solvent. For water-soluble proteins it has been an article of faith that water is largely responsible for stabilizing the fold, a notion that has recently come under increasing scrutiny. Further, there are some instances when proteins are studied experimentally in the absence of solvent, as in matrix-assisted laser desorption/ionization or electrospray mass spectrometry, for example, or in organic solvents for protein engineering purposes. Apart from these considerations, there is considerable speculation as to whether there is life on planets other than Earth, where conditions including the presence of water (both in liquid or vapor form and indeed ice), temperature and pressure may be vastly different from those prevailing on Earth. Mars, for example, has only 0.6% of Earth's mean atmospheric pressure which presents profound problems to protein structures, as this paper and a large corpus of experimental work demonstrate. Similar objections will most likely apply in the case of most exoplanets and other bodies such as comets whose chemistry and climate are still largely unknown. This poses the question, how do proteins survive in these different environments? In order to cast some light on these issues we have conducted a series of molecular dynamics simulations on protein dehydration under a variety of conditions. We find that, while proteins undergoing dehydration can retain their integrity for a short duration they ultimately become disordered, and we further show that the disordering can be retarded if superficial water is kept in place on the surface. These findings are compared with other published results on protein solvation in an astrobiological and astrochemical setting. Inter alia, our results suggest that there are limits as to what to expect in terms of the existence of possible extraterrestrial forms as well to what can be achieved in experimental investigations on living systems despatched from Earth. This finding may appear to undermine currently held hopes that life will be found on nearby planets, but it is important to be aware that the presence of ice and water are by themselves not sufficient; there has to be an atmosphere which includes water vapor at a sufficiently high partial pressure for proteins to be active. A possible scenario in which there has been a history of adequate water vapor pressure which allowed organisms to prepare for a future dessicated state by forming suitable protective capsules cannot of course be ruled out.

Native Mass Spectrometry for the Characterization of Structure and Interactions of Membrane Proteins

Over the past years, native mass spectrometry and ion mobility have grown into techniques that are widely applicable to the study of aspects of protein structure. More recently, it has become apparent that this approach provides a very promising avenue for the investigation of integral membrane proteins in lipid or detergent environments.In this chapter, we discuss applications of native mass spectrometry and ion mobility in membrane protein research-what is important to take into consideration when working with membrane proteins, and what the requirements are for sample preparation for native mass spectrometry. Furthermore, we will discuss the types of information provided by the measurements, including the oligomeric state, subunit composition and stoichiometry, interactions with detergents or lipids, conformational transitions, and the binding and structural effect of ligands and drugs.

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.

Stepwise Hydration of 2-Aminooxazole: Theoretical Insight into the Structure, Finite Temperature Behavior and Proton-Induced Charge Transfer

It was recently suggested that 2-aminooxazole (AO) could contribute to the formation of RNA nucleotides on primitive earth. In this article we have considered by means of computational modeling the influence of microhydration on the structural and spectral properties of this potential prebiotic molecule. The stable structures of AO(H2O)n were obtained first by sampling the potential energy landscapes of clusters containing up to n = 20 water molecules, using a simple but reasonably accurate force field and replica-exchange molecular dynamics simulations. Through reoptimization using an explicit description of electronic structure at the level of density functional theory with the M06-2X functional, the formation energies, ionization energies and electron affinities were determined in the vertical and adiabatic treatments, as well as vibrational and optical spectra covering the far-IR, mid-IR, and lower part of the UV ranges. The results generally show a clear segregation between the aminooxazole solute and the water molecules, a water cluster being formed near the nitrogen and amino group side leaving the hydrocarbon side dry even at temperatures corresponding to the liquid state. The spectral signatures generally concur and show distinct contributions of the solute and solvent, spectral shifts to lower energies being in agreement with earlier calculations in bulk solvent. We have also investigated the importance of microhydration on the charge transfer cross section upon collision with a proton, thereby extending an earlier investigation on the bare AO molecule. The presence of water molecules generally reduces the propensity for charge transfer at small sizes, but the influence of the solvent steadily decreases in larger droplets.

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.

Resolution of Stepwise Cooperativities of Copper Binding by the Homotetrameric Copper-Sensitive Operon Repressor (CsoR): Impact on Structure and Stability

The cooperativity of ligand binding is central to biological regulation and new approaches are needed to quantify these allosteric relationships. Herein, we exploit a suite of mass spectrometry (MS) experiments to provide novel insights into homotropic Cu-binding cooperativity, gas-phase stabilities and conformational ensembles of the D2 -symmetric, homotetrameric copper-sensitive operon repressor (CsoR) as a function of Cu(I) ligation state. Cu(I) binding is overall positively cooperative, but is characterized by distinct ligation state-specific cooperativities. Structural transitions occur upon binding the first and fourth Cu(I) , with the latter occurring with significantly higher cooperativity than previous steps; this results in the formation of a holo-tetramer that is markedly more resistant than apo-, and partially ligated CsoR tetramers toward surface-induced dissociation (SID).

Extending native mass spectrometry approaches to integral membrane proteins

Recent developments in native mass spectrometry and ion mobility have made it possible to analyze the composition and structure of membrane protein complexes in the gas-phase. In this short review we discuss the experimental strategies that allow to elucidate aspects of the dynamic structure of these important drug targets, such as the structural effects of lipid binding or detection of co-populated conformational and assembly states during gating on an ion channel. As native mass spectrometry relies on nano-electrospray of natively reconstituted proteins, a number of commonly used lipid- and detergent-based reconstitution systems have been evaluated for their compatibility with this approach, and parameters for the release of intact, native-like folded membrane proteins studied in the gas-phase. The strategy thus developed can be employed for the investigation of the subunit composition and stoichiometry, oligomeric state, conformational changes, and lipid and drug binding of integral membrane proteins.

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.