“Denaturation” and Refolding of Cytochrome c in Vacuo

A Modified Traveling Wave Ion Mobility Mass Spectrometer as a Versatile Platform for Gas-Phase Ion-Molecule Reactions

Taken individually, chemical labeling and mass spectrometry are two well-established tools for the structural characterization of biomolecular complexes. A way to combine their respective advantages is to perform gas-phase ion-molecule reactions (IMRs) inside the mass spectrometer. This is, however, not so well developed because of the limited range of usable chemicals and the lack of commercially available IMR devices. Here, we modified a traveling wave ion mobility mass spectrometer to enable IMRs in the trapping region of the instrument. Only one minor hardware modification is needed to allow vapors of a variety of liquid reagents to be leaked into the trap traveling wave ion guide of the instrument. A diverse set of IMRs can then readily be performed without any loss in instrument performance. We demonstrate the advantages of implementing IMR capabilities in general, and to this quadrupole-ion mobility-time-of-flight (Q-IM-TOF) mass spectrometer in particular, by exploiting the full functionality of the instrument, including mass selection, ion mobility separation, and post-mobility fragmentation. The potential to carry out gas-phase IMR kinetics experiments is also illustrated. We demonstrate the versatility of the setup using gas-phase IMRs of established utility for biological mass spectrometry, including hydrogen-deuterium exchange, ion-molecule proton transfer reactions, and covalent modification of DNA anions using trimethylsilyl chloride.

Mass spectrometric characterization of protein structures and protein complexes in condensed and gas phase

Proteins are essential for almost all physiological processes of life. They serve a myriad of functions which are as varied as their unique amino acid sequences and their corresponding three-dimensional structures. To fulfill their tasks, most proteins depend on stable physical associations, in the form of protein complexes that evolved between themselves and other proteins. In solution (condensed phase), proteins and/or protein complexes are in constant energy exchange with the surrounding solvent. Albeit methods to describe in-solution thermodynamic properties of proteins and of protein complexes are well established and broadly applied, they do not provide a broad enough access to life-science experimentalists to study all their proteins' properties at leisure. This leaves great desire to add novel methods to the analytical biochemist's toolbox. The development of electrospray ionization created the opportunity to characterize protein higher order structures and protein complexes rather elegantly by simultaneously lessening the need of sophisticated sample preparation steps. Electrospray mass spectrometry enabled us to translate proteins and protein complexes very efficiently into the gas phase under mild conditions, retaining both, intact protein complexes, and gross protein structures upon phase transition. Moreover, in the environment of the mass spectrometer (gas phase, in vacuo), analyte molecules are free of interactions with surrounding solvent molecules and, therefore, the energy of inter- and intramolecular forces can be studied independently from interference of the solvating environment. Provided that gas phase methods can give information which is relevant for understanding in-solution processes, gas phase protein structure studies and/or investigations on the characterization of protein complexes has rapidly gained more and more attention from the bioanalytical scientific community. Recent reports have shown that electrospray mass spectrometry provides direct access to six prime protein complex properties: stabilities, compositions, binding surfaces (epitopes), disassembly processes, stoichiometries, and thermodynamic parameters.

Benchmark Comparison for a Multi-Processing Ion Mobility Calculator in the Free Molecular Regime

A benchmark comparison between two ion mobility and collision cross-section (CCS) calculators, MOBCAL and IMoS, is presented here as a standard to test the efficiency and performance of both programs. Utilizing 47 organic ions, results are in excellent agreement between IMoS and MOBCAL in He and N_2, when both programs use identical input parameters. Due to a more efficiently written algorithm and to its parallelization, IMoS is able to calculate the same CCS (within 1%) with a speed around two orders of magnitude faster than its MOBCAL counterpart when seven cores are used. Due to the high computational cost of MOBCAL in N₂, reaching tens of thousands of seconds even for small ions, the comparison between IMoS and MOBCAL is stopped at 70 atoms. Large biomolecules (>10000 atoms) remain computationally expensive when IMoS is used in N₂ (even when employing 16 cores). Approximations such as diffuse trajectory methods (DHSS, TDHSS) with and without partial charges and projected area approximation corrections can be used to reduce the total computational time by several folds without hurting the accuracy of the solution. These latter methods can in principle be used with coarse-grained model structures and should yield acceptable CCS results. Graphical Abstract ᅟ.

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.

Folding of Protein Ions in the Gas Phase after Cation-to-Anion Proton-Transfer Reactions

The structure and folding of a protein in solution depends on noncovalent interactions within the protein and those with surrounding ions and molecules. Decoupling these interactions in solution is challenging, which has hindered the development of accurate physics-based models for structure prediction. Investigations of proteins in the gas phase can be used to selectively decouple factors affecting the structures of proteins. Here, we use cation-to-anion proton-transfer reactions (CAPTR) to reduce the charge states of denatured ubiquitin ions in the gas phase, and ion mobility to probe their structures. In CAPTR, a precursor charge state is selected (P) and reacted with monoanions to generate charge-reduced product ions (C). Following each CAPTR event, denatured ubiquitin ions (13+ to 6+) yield products that rapidly isomerize to structures that have smaller collision cross sections (Ω). The Ω values of CAPTR product ions depend strongly on C and very weakly on P. Pre- and post-CAPTR activation was then used to probe the potential-energy surfaces of the precursor and product ions, respectively. Post-CAPTR activation showed that ions of different P fold differently and populate different regions of the potential-energy surface of that ion. Finally, pre-CAPTR activation showed that the structures of protein ions can be indirectly investigated using ion mobility of their CAPTR product ions, even for subtle structural differences that are not apparent from ion mobility characterization of the activated precursor ions. More generally, these results show that CAPTR strongly complements existing techniques for characterizing the structures and dynamics of biological molecules in the gas phase.

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.

Insights into the conformations of three structurally diverse proteins: cytochrome c, p53, and MDM2, provided by variable-temperature ion mobility mass spectrometry

Thermally induced conformational transitions of three proteins of increasing intrinsic disorder-cytochrome c, the tumor suppressor protein p53 DNA binding domain (p53 DBD), and the N-terminus of the oncoprotein murine double minute 2 (NT-MDM2)-have been studied by native mass spectrometry and variable-temperature drift time ion mobility mass spectrometry (VT-DT-IM-MS). Ion mobility measurements were carried out at temperatures ranging from 200 to 571 K. Multiple conformations are observable over several charge states for all three monomeric proteins, and for cytochrome c, dimers of significant intensity are also observed. Cytochrome c [M + 5H](5+) ions present in one conformer of CCS ∼1200 Å(2), undergoing compaction in line with the reported Tmelt = 360.15 K before slight unfolding at 571 K. The more extended [M + 7H](7+) cytochrome c monomer presents as two conformers undergoing similar compaction and structural rearrangements, prior to thermally induced unfolding. The [D + 11H](11+) dimer presents as two conformers, which undergo slight structural compaction or annealing before dissociation. p53 DBD follows a trend of structural collapse before an increase in the observed collision cross section (CCS), akin to that observed for cytochrome c but proceeding more smoothly. At 300 K, the monomeric charge states present in two conformational families, which compact to one conformer of CCS ∼1750 Å(2) at 365 K, in line with the low solution Tmelt = 315-317 K. The protein then extends to produce either a broad unresolved CCS distribution or, for z > 9, two conformers. NT-MDM2 exhibits a greater number of structural rearrangements, displaying charge-state-dependent unfolding pathways. DT-IM-MS experiments at 200 K resolve multiple conformers. Low charge state species of NT-MDM2 present as a single compact conformational family centered on CCS ∼1250 Å(2) at 300 K. This undergoes conformational tightening in line with the solution Tmelt = 348 K before unfolding at the highest temperatures. The more extended charge states present in two or more conformers at room temperature, undergoing thermally induced unfolding before significant structural collapse or annealing at high temperatures. Variable-temperature IM-MS is here shown to be an exciting approach to discern protein unfolding pathways for conformationally diverse proteins.

Time-resolved pulsed hydrogen/deuterium exchange mass spectrometry probes gaseous proteins structural kinetics

A pulsed hydrogen/deuterium exchange (HDX) method has been developed for rapid monitoring of the exchange kinetics of protein ions with D2O a few milliseconds after electrospray ionization (ESI). The stepwise gradual evolution of HDX of multiply charged protein ions was monitored using the pulsed HDX mass spectrometry technique. Upon introducing a very short pulse of D2O (in the μs to ms time scale) into the linear ion trap (LIT) of a time-of-flight (TOF) mass spectrometer, bimodal distributions were detected for the ions of cytochrome c and ubiquitin. Mechanistic details of HDX reactions for ubiquitin and cytochrome c in the gas phase were uncovered and the structural transitions were followed by analyzing the kinetics of HDX.

Developing liquid chromatography ion mobility mass spectometry techniques

When a packet of ions in a buffer gas is exposed to a weak electric field, the ions will separate according to differences in their mobilities through the gas. This separation forms the basis of the analytical method known as ion mobility spectroscopy and is highly efficient, in that it can be carried out in a very short time frame (micro- to milliseconds). Recently, efforts have been made to couple the approach with liquid-phase separations and mass spectrometry in order to create a high-throughput and high-coverage approach for analyzing complex mixtures. This article reviews recent work to develop this approach for proteomics analyses. The instrumentation is described briefly. Several multidimensional data sets obtained upon analyzing complex mixtures are shown in order to illustrate the approach as well as provide a view of the limitations and required future work.

Critical examination of gas-phase protein conformation/multimer ion formation by electrospray ion mobility-mass spectrometry

The generally accepted view of protein structure in the gas-phase is that protein ions produced by electrospray ionization (ESI) exist in a number of different states, and the resulting charge state distribution (CSD) and ion mobility spectrum is interpreted as evidence for protein ions retaining some memory of solution-phase conformation. Even with the inclusion of ion mobility information, reports of protein ion structure in the gas-phase are oftentimes in disagreement not only within the discipline but also as interpreted by other gas-phase techniques. The focus of this work will be to correctly distinguish truly different ion conformations formed by ESI versus homomultimeric complexes with the same m/z. The concentration of cytochrome c in solution was varied over a wide range, and the multiply charged multimers (MCMs) present in the ion mobility/mass spectrum were unambiguously assigned by m/z selection and dissociation prior to ion mobility/mass spectrometry analysis. The results revealed false negatives for protein oligomer formation and false positives for protein conformational states and no evidence that gas-phase cytochrome c ions retain memory of solution-phase conformation, characteristics of great importance for structural biology. The results also suggest that the total IM-MS distribution for a protein is the complex result of individual MCMs either surviving until detection (undissociated) or dissociating into lower order multimers or a number of product ions for each m/z.

A close look at proteins: submolecular resolution of two- and three-dimensionally folded cytochrome c at surfaces

Imaging of individual protein molecules at the single amino acid level has so far not been possible due to the incompatibility of proteins with the vacuum environment necessary for high-resolution scanning probe microscopy. Here we demonstrate electrospray ion beam deposition of selectively folded and unfolded cytochrome c protein ions on atomically defined solid surfaces in ultrahigh vacuum (10(-10) mbar) and achieve unprecedented resolution with scanning tunneling microscopy. On the surface folded proteins are found to retain their three-dimensional structure. Unfolded proteins are observed as extended polymer strands displaying submolecular features with resolution at the amino acid level. On weakly interacting surfaces, unfolded proteins refold into flat, irregular patches composed of individual molecules. This suggests the possibility of two-dimensionally confined folding of peptides of an appropriate sequence into regular two-dimensional structures as a new approach toward functional molecular surface coatings.

Alkali metal-cationized serine clusters studied by sonic spray ionization tandem mass spectrometry

Serine solutions containing salts of alkali metals yield magic number clusters of the type (Ser(4)+C)(+), (Ser(8)+C)(+), (Ser(12)+C)(+), and (Ser(17)+2C)(+2) (where C = Li(+), Na(+), K(+), Rb(+), or Cs(+)), in relative abundances which are strongly dependent on the cation size. Strong selectivity for homochirality is involved in the formation of serine tetramers cationized by K(+), Rb(+), and Cs(+). This is also the case for the octamers cationized by the smaller alkalis but there is a strong preference for heterochirality in the octamers cationized by the larger alkali cations. Tandem mass spectrometry shows that the octamers and dodecamers cationized by K(+), Rb(+), and Cs(+) dissociate mainly by the loss of Ser(4) units, suggesting that the neutral tetramers are the stable building blocks of the observed larger aggregates, (Ser(8)+C)(+) and (Ser(12)+C)(+). Remarkably, although the Ser(4) units are formed with a strong preference for homochirality, they aggregate further regardless of their handedness and, therefore, with a preference for the nominally racemic 4D:4L structure and an overall strong heterochiral preference. The octamers cationized by K(+), Rb(+), or Cs(+) therefore represent a new type of cluster ion that is homochiral in its internal subunits, which then assemble in a random fashion to form octamers. We tentatively interpret the homochirality of these tetramers as a consequence of assembly of the serine molecules around a central metal ion. The data provide additional evidence that the neutral serine octamer is homochiral and is readily cationized by smaller ions.

Resolution and structural transitions of elongated states of ubiquitin

Electrospray ionization, combined with two-dimensional ion mobility spectrometry and mass spectrometry, is used to produce, select, and activate distributions of elongated ions, [M + 11H]11+ to [M + 13H]13+, of ubiquitin. The analysis makes it possible to examine state-to-state transitions for structural types, and transition diagrams associated with the efficiencies of structural changes are presented. The +11 and +12 charge states can form four resolvable states while only one state is formed for [M + 13H]13+. Some conformations, which appear to belong to the same family based on mobility analysis of different charge states, undergo similar transitions, others do not. Activation of ions that exist in low-abundance conformations, having mobilities that fall in between sharp peaks associated with higher abundances species, shows that the low-abundance forms undergo efficient (approximately 90 to 100%) conversion into states associated with well-defined peaks. This efficiency is significantly higher than the approximately 10 to 60% efficiency of transitions of structures associated with well-defined peaks. The formation of sharp features from a range of low-intensity species with different cross sections indicates that large regions of conformation space must be unfavorable or inaccessible in the gas phase. These results are compared with several previous IMS measurements of this system as well as information about gas-phase structure provided by other techniques.

Collision-induced alignment of H2O+ drifting in helium

The collision-induced alignment of H(2)O(+) drifting in helium is studied with a molecular dynamics method that has been extended to treat nonlinear rigid ions. Rotational distribution functions and averaged quantities are presented in terms of the rho formalism [M. Thachuk, Phys. Rev. A 72, 032722 (2005)], and it is shown that this description gives a very good agreement with simulation results. In addition to velocity and angular velocity autocorrelation functions, a velocity-angular velocity cross correlation function is introduced. This cross correlation function provides insight into the dynamical nature of the alignment mechanism.

Characterizing the Structures and Folding of Free Proteins Using 2-D Gas-Phase Separations: Observation of Multiple Unfolded Conformers

Understanding the 3-D structure and dynamics of proteins and other biological macromolecules in various environments is among the central challenges of chemistry. Electrospray ionization can often transfer ions from solution to gas phase with only limited structural distortion, allowing their profiling using mass spectrometry and other gas-phase approaches. Ion mobility spectrometry (IMS) can separate and characterize macroion conformations with high sensitivity and speed. However, IMS separation power is generally insufficient for full resolution of major structural variants of protein ions and elucidation of their interconversion dynamics. Here we report characterization of macromolecular conformations using field asymmetric waveform IMS (FAIMS) coupled to conventional IMS in conjunction with mass spectrometry. The collisional heating of ions in the electrodynamic funnel trap between FAIMS and IMS stages enables investigating the structural evolution of particular isomeric precursors as a function of the intensity and duration of activation that can be varied over large ranges. These new capabilities are demonstrated for ubiquitin and cytochrome c, two common model proteins for structure and folding studies. For nearly all charge states, two-dimensional FAIMS/IMS separations distinguish many more conformations than either FAIMS or IMS alone, including some with very low abundance. For cytochrome c in high charge states, we find several abundant "unfolded" isomer series not distinguishable by IMS, possibly corresponding to different "string of beads" geometries. The unfolding of specific ubiquitin conformers selected by FAIMS has been studied by employing their heating in the FAIMS/IMS interface.

Blackbody infrared radiative dissociation of nonspecific protein-carbohydrate complexes produced by nanoelectrospray ionization: the nature of the noncovalent interactions

Gas-phase thermal dissociation experiments, implemented with blackbody infrared radiative dissociation (BIRD) and Fourier transform ion cyclotron resonance mass spectrometry, have been performed on a series of protonated and deprotonated 1:1 and protonated 1:2 protein-carbohydrate complexes formed by nonspecific interactions during the nanoflow electrospray (nanoES) ionization process. Nonspecific interactions between the proteins bovine carbonic anhydrase II (CA), bovine ubiquitin (Ubq), and bovine pancreatic trypsin inhibitor and several carbohydrates, ranging in size from mono- to tetrasaccharides, have been investigated. Over the range of temperatures studied (60-190 degrees C), BIRD of the protonated and deprotonated complexes proceeds exclusively by the loss of the carbohydrate in its neutral form. The rates of dissociation of the 1:1 complexes containing a mono- or disaccharide decrease with reaction time, suggesting the presence of two or more kinetically distinct structures produced during nanoES or by gas-phase processes. In contrast, the 1:1 complexes of the tri- and tetrasaccharides exhibit simple first-order dissociation kinetics, a result that, on its own, is suggestive of a single preferred carbohydrate binding site or multiple equivalent sites in the gas phase. A comparative analysis of the dissociation kinetics measured for protonated 1:1 and 1:2 complexes of Ubq with alphaTal[alphaAbe]alphaMan further supports the presence of a single preferred binding site. However, a similar analysis performed on the complexes of CA and alphaTal[alphaAbe]alphaMan suggests that equivalent but dependent carbohydrate binding sites exist in the gas phase. Analysis of the Arrhenius activation parameters (E(a) and A) determined for the dissociation of 1:1 complexes of CA with structurally related trisaccharides provides evidence that neutral intermolecular hydrogen bonds contribute, at least in part, to the stability of the gaseous complexes. Surprisingly, the E(a) values for the complexes of the same charge state are not sensitive to the structure (primary or higher order) of the protein, suggesting that the carbohydrates are able to form energetically equivalent interactions with the various functional groups presented by the protein. For a given protein-carbohydrate complex, the dissociation E(a) is sensitive to charge state, initially increasing and then decreasing with increasing charge. It is proposed that both ionic and neutral hydrogen bonds stabilize the nonspecific protein-carbohydrate complexes in the gas phase and that the relative contribution of the neutral and ionic interactions is strongly influenced by the charge state of the complex, with neutral interactions dominating at low charge states and ionic interactions dominating at high charge states.

Nanoflow LC/IMS-MS and LC/IMS-CID/MS of protein mixtures

A simple ion trap/ion mobility/time-of-flight (TOF) mass spectrometer has been coupled with nanoflow liquid chromatography to examine the feasibility of analyzing mixtures of intact proteins. In this approach proteins are separated using reversed-phase chromatography. As components elute from the column, they are electrosprayed into the gas phase and separated again in a drift tube prior to being dispersed and analyzed in a TOF mass spectrometer. The mobilities of ions through a buffer gas depend upon their collision cross sections and charge states; separation based on these gas-phase parameters provides a new means of simplifying mass spectra and characterizing mixtures. Additionally it is possible to induce dissociation at the exit of the drift tube and examine the fragmentation patterns of specific protein ion charge states and conformations. The approach is demonstrated by examining a simple three-component mixture containing ubiquitin, cytochrome c, and myoglobin and several larger prepared protein mixtures. The potential of this approach for use in proteomic applications is considered.

A Split-Field Drift Tube for Separation and Efficient Fragmentation of Biomolecular Ions

A new ion mobility instrument that incorporates a low-field region for ion separation and a high-field region for collisional activation is described. In this approach, mixtures of ions are separated based on differences in their mobilities in a approximately 20-cm-long low-field ( approximately 5 V cm(-)(1)) region of a drift tube. As the ions approach the drift tube exit, they are exposed to a large focusing potential drop; at high fields, ions are efficiently collisionally activated and dissociate as they exit the drift tube. We have demonstrated this approach by examining the fragmentation pattern for electrosprayed bradykinin ions. These studies show that the activation process is highly tunable; it is possible to modulate the field such that precursor ion mass spectra as well as several high-field collision-induced fragmentation patterns can be obtained. The approach also appears to be a simple means of activating protein ions, as demonstrated by examining electrosprayed myoglobin ions.

Temperature-dependent H/D exchange of compact and elongated cytochrome c ions in the gas phase

Isotopic exchange reactions of compact and elongated conformations of gaseous cytochrome c ions (+5 and +9 states) with D2O have been measured as a function of temperature (from 300 to approximately 440 K) using ion mobility techniques. Rate constants for those sites that exchange at high temperatures (>400 K) are about an order of magnitude smaller than rate constants for sites that exchange at 300 K. Although the exchange rates decrease, the maximum exchange levels for rapidly exchanging sites increase with temperature. At 300 K, exchange levels of 53 +/- 3 and 63 +/- 3 are measured for the compact and elongated states, respectively. From 300 to 335 K, the exchange levels increase slightly to approximately 60 to 70 hydrogens. Above 335 K, the levels increase to a value of approximately 200 for the +5 state and approximately 190 for the +9 state, near the maximum possible levels, 200 and 204 for these respective charge states. Molecular dynamics simulations have been carried out on structures having calculated cross sections that are near the experimental values in order to explore the exchange process. Overall, it appears that charge site and exchange site proximities are important factors in the exchange profiles for the elongated +9 ion and the compact +5 ion.

Structural transitions in neutral and charged proteins in vacuo

In vacuo proteins provide a simple laboratory to explore the roles of sequence, temperature, charge state, and initial configuration in protein folding. Moreover, by the very absence of solvent, the study of anhydrous proteins in vacuo will also help us to understand specific environmental effects. From the experimental viewpoint, these systems are now beginning to be characterized at low resolution. Molecular dynamics (MD) simulations, in combination with tools for protein shape analysis, can complement experiments and provide further insights on the folding-unfolding transitions of these proteins. We review some aspects of this issue by using the results from a detailed MD study of hen egg-white lysozyme. For lysozyme ions, unfolding can be triggered by Coulombic repulsion. In neutral lysozyme, unfolding can be induced by centrifugal forces and also by weakening the monomer-monomer interaction. In both cases, the resulting unfolded transients can be used as initial configurations for relaxation dynamics. All trajectories are analyzed in terms of global molecular shape features of the backbone, including its anisometry and chain entanglement complexity. This strategy allows us to quantify separately the degree of polymer collapse and the evolution of large-scale folding features. Using these last two notions, we discuss some basic questions regarding the nature of the accessible paths associated with unfolding from, and refolding into, compact conformers.

Proteins in vacuo: denaturing and folding mechanisms studied with computer-simulated molecular dynamics

Mounting evidence from experiments suggests that the native fold in solution is metastable in dehydrated proteins. Results from a number of experiments that use mass spectrometry indicate also that folding-unfolding transitions take place in protein ions even in the absence of water. These observations on anhydrous proteins call for a re-evaluation of our understanding of the folding transition. In this context, computer-assisted simulations are an important complementary tool. Here, we provide an overview of recent progress on the simulation of proteins in vacuo. In particular, we discuss the response of proteins and protein ions to perturbations that trigger unfolding and re-folding transitions. By comparing the general patterns emerging from theory and experiment, we propose a series of new measurements that could help to validate, and improve, current simulation models.

Peptides and proteins in the vapor phase

This article provides a review of recent studies of the properties of unsolvated (and partially solvated) peptides and proteins. The methods used to produce vapor-phase peptide and protein ions are described along with some of the techniques used to study them, such as H/D exchange, blackbody infrared radiative dissociation, and ion mobility measurements. Studies of unsolvated peptides and proteins provide information about their intrinsic intramolecular interactions. The topics covered include the role of zwitterions and salt bridges in the vapor phase, Coulomb interactions in multiply charged ions, the unfolding and refolding of vapor-phase proteins, and the stability of unsolvated helices and sheets. Finally, dehydration and rehydration studies of proteins in the vapor phase are described. These can provide exquisitely detailed information about hydration interactions, such as the enthalpy and entropy changes associated with adsorbing individual water molecules.

A database of 660 peptide ion cross sections: use of intrinsic size parameters for bona fide predictions of cross sections

An ion trap/ion mobility/time-of-flight mass spectrometry technique has been used to measure collision cross sections for 660 peptide ions generated by tryptic digestion of 34 common proteins. Measured cross sections have been compiled into a database that contains peptide molecular weight and sequence information. The database is used to generate average intrinsic contributions to cross section (size parameters) for different amino acid residues by solving systems of equations that relate the unknown contributions of individual residues to the sequences and cross sections of database peptides. Size parameters are combined with information about amino acid composition to calculate cross sections for database peptides. Bona fide cross section predictions (made prior to measurement) for peptides observed in tryptic digests of sperm whale myoglobin and yeast enolase are made. Eight of 10 predicted cross sections are within 2% of the experimental values and all 10 are within 3.2%. The utility of size parameters for cross section prediction is explored and discussed.

Effects of cysteic acid groups on the gas-phase reactivity and dissociation of [M + 4H]4+ ions from insulin chain B

Gas-phase ion/molecule reactions and collision-induced dissociation (CID) were conducted on [M + 4H]4+ of insulin chain B. This Fourier transform mass spectrometry work involved ions from the oxidized peptide (with two cysteic acid residues) and its reduced form (with two cysteine residues). Kinetic behavior during deprotonation and hydrogen/deuterium exchange reactions indicates that insulin B (ox) ions have two distinct structural types. In contrast, insulin B (red) ions have only one major reacting population, which has a more compact structure than the oxidized ions. No significant differences in fragmentation patterns for the two insulin B (ox) populations were observed when CID was performed as a function of deprotonating reaction time. However, markedly different fragmentation was found between [M + 4H]4+ of insulin B (ox) and (red). Therefore, the presence of cysteic acid groups in insulin B (ox) significantly impacts dissociation and presumably structure. This suggests that some insulin B (ox) ions are zwitterionic, with the five basic sites protonated and one cysteic acid group deprotonated. Molecular dynamics calculations revealed several viable structures that are consistent with the experimental results. For example, the most stable form of the reduced ion, which is unprotonated at the His10, is very compact and has lost the alpha-helix of native insulin. Low energy structures for the oxidized ions include a zwitterion with an intraionic interaction between anionic Cyx7 and cationic His10, as well as a nonzwitterionic conformer that lacks a proton at Phe1; both structures retain the alpha-helix. These structures may account for the two experimentally observed isomers, although others are possible. In addition, experiments on oxidized insulin B were conducted from methanolic solution, which may denature the conformation, and pure aqueous solution, which may leave a native conformation. These differences in solvent composition had no effect on the gas-phase results.

Characterization of fold diversity among proteins with the same number of amino acid residues

Chain entanglements provide a simple and global measure of folding in a macromolecule. The complexity of these entanglements can be expressed by the pattern of projected bond-bond crossings, or "overcrossings", associated with the molecular backbone. In this work, we use this approach to characterize quantitatively the range of tertiary folds observed in proteins with a given chain length. To discriminate among folding features, we use two shape descriptors derived from the probability distribution of overcrossings: the mean overcrossing number, N, and the most probable overcrossing number, N*. The values of N and N* relate to the content of secondary structure in a protein as well as its global three-dimensional organization. We propose a measure of folding diversity based on the properties of these descriptors. In addition, we discuss the application of our method to study how tertiary folds evolve during protein dynamics.