Molecular Dynamics Simulations of the Charge-Induced Unfolding and Refolding of Unsolvated Cytochrome c

Oxidized and Reduced Dimeric Protein Complexes Illustrate Contrasting CID and SID Charge Partitioning

Charge partitioning during the dissociation of protein complexes in the gas phase is influenced by many factors, such as interfacial interactions, protein flexibility, protein conformation, and dissociation methods. In the present work, two cysteine-containing homodimer proteins, β-lactoglobulin and α-lactalbumin, with the disulfide bonds intact and reduced, were used to gain insight into the charge partitioning behaviors of collision-induced dissociation (CID) and surface-induced dissociation (SID) processes. For these proteins, we find that restructuring dominates with CID and dissociation with symmetric charge partitioning dominates with SID, regardless of whether intramolecular disulfide bonds are oxidized or reduced. CID of the charge-reduced dimeric protein complex leads to a precursor with a slightly smaller collision cross section (CCS), greater stability, and more symmetrically distributed charges than the significantly expanded form produced by CID of the higher charged dimer. Collision-induced unfolding plots demonstrate that the unfolding-restructuring of the protein complexes initiates the charge migration of higher charge-state precursors. Overall, gas collisions reveal the charge-dependent restructuring/unfolding properties of the protein precursor, while surface collisions lead predominantly to more charge-symmetric monomer separation. CID's multiple low-energy collisions sequentially reorganize intra- and intermolecular bonds, while SID's large-step energy jump cleaves intermolecular interfacial bonds in preference to reorganizing intramolecular bonds. The activated population of precursors that have taken on energy without dissociating (populated in CID over a wide range of collision energies, populated in SID for only a narrow distribution of collision energies near the onset of dissociation) is expected to be restructured, regardless of the activation method.

Using Density Functional Theory for Testing the Robustness of Mobile-Proton Molecular Dynamics Simulations on Electrosprayed Ions: Structural Implications for Gaseous Proteins

Current experiments only provide low-resolution information on gaseous protein ions generated by electrospray ionization (ESI). Molecular dynamics (MD) simulations can yield complementary insights. Unfortunately, conventional MD does not capture the mobile nature of protons in gaseous proteins. Mobile-proton MD (MPMD) overcomes this limitation. Earlier MPMD data at 300 K indicated that protein ions generated by "native" ESI retain solution-like structures with a hydrophobic core and zwitterionic exterior [Bakhtiari, M.; Konermann, L. J. Phys. Chem. B 2019, 123, 1784-1796]. MPMD redistributes protons using electrostatic and proton affinity calculations. The robustness of this approach has never been scrutinized. Here, we close this gap by benchmarking MPMD against density functional theory (DFT) at the B3LYP/6-31G* level, which is well suited for predicting proton affinities. The computational cost of DFT necessitated the use of small peptides. The MPMD energetic ranking of proton configurations was found to be consistent with DFT single-point energies, implying that MPMD can reliably identify favorable protonation sites. Peptide MPMD runs converged to DFT-optimized structures only when applying 300-500 K temperature cycling, which was necessary to prevent trapping in local minima. Temperature cycling MPMD was then applied to gaseous protein ions. Native ubiquitin converted to slightly expanded structures with a zwitterionic core and a nonpolar exterior. Our data suggest that such inside-out protein structures are intrinsically preferred in the gas phase, and that they form in ESI experiments after moderate collisional excitation. This is in contrast to native ESI (with minimal collisional excitation, simulated by MPMD at 300 K), where kinetic trapping promotes the survival of solution-like structures. In summary, this work validates the MPMD approach for simulations on gaseous peptides and proteins.

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

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

Suppression of Protein Structural Perturbations in Native Electrospray Ionization during the Final Evaporation Stages Revealed by Molecular Dynamics Simulations

Native electrospray ionization was known to preserve the protein structure in solution, which overcame the uncontrollable acidification of droplets during transfer from solution into the gas phase in conventional electrospray ionization. However, detailed experimental studies on when and how could native electrospray ionization minimize structural perturbations remain quite unclear. Herein, we conducted molecular dynamics simulations to investigate the protein structure evolution during electrospray ionization. At a neutral droplet pH, the protein structure in solution could be retained after evaporation, which was in accordance with previous reports. As the droplet pH deviated from neutral, we have found that the compact protein structure would not unfold until the last 10 ns prior to the final desolvation, which demonstrated that the role of native electrospray ionization in preserving the protein structure was mainly reflected on the final evaporation stages. The present study might provide new insights into studying the microscopic biomolecular events occurring during the liquid-gas interface transition and their influence on solution-structure retention.

Interrogating the Quaternary Structure of Noncanonical Hemoglobin Complexes by Electrospray Mass Spectrometry and Collision-Induced Dissociation

Various activation methods are available for the fragmentation of gaseous protein complexes produced by electrospray ionization (ESI). Such experiments can potentially yield insights into quaternary structure. Collision-induced dissociation (CID) is the most widely used fragmentation technique. Unfortunately, CID of protein complexes is dominated by the ejection of highly charged monomers, a process that does not yield any structural insights. Using hemoglobin (Hb) as a model system, this work examines under what conditions CID generates structurally informative subcomplexes. Native ESI mainly produced tetrameric Hb ions. In addition, "noncanonical" hexameric and octameric complexes were observed. CID of all these species [(αβ)₂, (αβ)₃, and (αβ)₄] predominantly generated highly charged monomers. In addition, we observed hexamer → tetramer + dimer dissociation, implying that hexamers have a tetramer··dimer architecture. Similarly, the observation of octamer → two tetramer dissociation revealed that octamers have a tetramer··tetramer composition. Gas-phase candidate structures of Hb assemblies were produced by molecular dynamics (MD) simulations. Ion mobility spectrometry was used to identify the most likely candidates. Our data reveal that the capability of CID to produce structurally informative subcomplexes depends on the fate of protein-protein interfaces after transfer into the gas phase. Collapse of low affinity interfaces conjoins the corresponding subunits and favors CID via monomer ejection. Structurally informative subcomplexes are formed only if low affinity interfaces do not undergo a major collapse. However, even in these favorable cases CID is still dominated by monomer ejection, requiring careful analysis of the experimental data for the identification of structurally informative subcomplexes.

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.

Towards the Analysis of High Molecular Weight Proteins and Protein complexes using TIMS-MS

In the present work, we demonstrate the potential and versatility of TIMS for the analysis of proteins, DNA-protein complexes and protein-protein complexes in their native and denatured states. In addition, we show that accurate CCS measurement are possible and in good agreement with previously reported CCS values using other IMS analyzers (<5% difference). The main challenges for the analysis of high mass proteins and protein complexes in the mobility and m/z domain are described. That is, the analysis of high molecular weight systems in their native state may require the use of higher electric fields or a compromise in the TIMS mobility resolution by reducing the bath gas velocity in order to effectively trap at lower electric fields. This is the first report of CCS measurements of high molecular weight biomolecules and biomolecular complexes (~ 150 kDa) using TIMS-MS.

Gas-phase microsolvation of ubiquitin: investigation of crown ether complexation sites using ion mobility-mass spectrometry

The influence of side chain to backbone interactions on the gas-phase structure of ubiquitin and ubiquitin lysine-to-arginine mutants was analysed.

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.

Activation state-selective kinase inhibitor assay based on ion mobility-mass spectrometry

The discovery of activation state dependent kinase inhibitors, which bind specifically to the inactive conformation of the protein, is considered to be a promising pathway to improved cancer treatments. Identifying such inhibitors is challenging, however, because they can have Kd values similar to molecules known to inhibit kinase function by interacting with the active form. Further, while inhibitor induced changes within the kinase tertiary structure are significant, few technologies are able to correctly assign inhibitor binding modes in a high-throughput fashion based exclusively on protein-inhibitor complex formation and changes in local protein structure. We have developed a new assay, using ion mobility-mass spectrometry, capable of both rapidly detecting inhibitor binding and classifying the resultant kinase binding modes. Here, we demonstrate the ability of our approach to classify a broad set of kinase inhibitors, using micrograms of protein, without the need for protein modification or tagging.

The dynamics of water evaporation from partially solvated cytochrome c in the gas phase

The study of evaporation of water from biological macromolecules is important for the understanding of electrospray mass spectrometry experiments. In electrospray ionization (ESI), electrically charged nanoscale droplets are formed from solutions of, for example, proteins. Then evaporation of the solvent leads to dry protein ions that can be analyzed in the mass spectrometer. In this work the dynamics of water evaporation from native cytochrome c covered by a monolayer of water is studied by molecular dynamics (MD) simulations at constant energy. A model of the initial conditions of the process is introduced. The temperature of the protein drops by about 100 K during the 400 picoseconds of the simulations. This sharp drop in temperature causes the water evaporation rate to decrease by about an order of magnitude, leaving the protein with 50% to 90% of the original water molecules, depending on the initial temperature of the simulation. The structural changes of the protein upon desolvation were considered through calculations of the radius of gyration and the root mean square (RMS) of the protein. A variation of 0.4 A in the radius of gyration, together with an RMS value of less than 3 A, indicates only minor changes in the overall shape of the protein structure. The water coordination number of the solvation shell is much smaller than that for bulk water. The mobility of water is high at the beginning of the simulations and drops as the simulation progresses and the temperature decreases. Incomplete desolvation of protein ions was also observed in recent experiments.

Field Asymmetric Waveform Ion Mobility Spectrometry Studies of Proteins: Dipole Alignment in Ion Mobility Spectrometry?

Approaches to separation and characterization of ions based on their mobilities in gases date back to the 1960s. Conventional ion mobility spectrometry (IMS) measures the absolute mobility, and field asymmetric waveform IMS (FAIMS) exploits the difference between mobilities at high and low electric fields. However, in all previous IMS and FAIMS experiments ions experienced an essentially free rotation; thus the separation was based on the orientationally averaged cross-sections Omega(avg) between ions and buffer gas molecules. Virtually all large ions are permanent electric dipoles that will be oriented by a sufficiently strong electric field. Under typical FAIMS conditions this will occur for dipole moments >400 D, found for many macroions including most proteins above approximately 30 kDa. Mobilities of aligned dipoles depend on directional cross-sections Omega(dir) (rather than Omega(avg)), which should have a major effect on FAIMS separation parameters. Here we report the FAIMS behavior of electrospray-ionization-generated ions for 10 proteins up to approximately 70 kDa. Those above 29 kDa exhibit a strong increase of mobility at high field, which is consistent with predicted ion dipole alignment. This effect expands the useful FAIMS separation power by an order of magnitude, allowing separation of up to approximately 10(2) distinct protein conformers and potentially revealing information about Omega(dir) and ion dipole moment that is of utility for structural characterization. Possible approaches to extending dipole alignment to smaller ions are discussed.

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.

A molecular dynamics simulation of SNase and its hydration shell at high temperature and high pressure

Temperature- and pressure-induced unfolding of staphylococcal nuclease (SNase) was studied by Royer, Winter et al. using a variety of experimental techniques (SAXS, FT-IR and fluorescence spectroscopy, DSC, PPC, densimetry). For a more detailed understanding of the underlying mechanistic processes of the different unfolding scenarios, we have carried out a series of molecular dynamics (MD) computer simulations on SNase. We investigated the initial changes of the structure of the protein upon application of pressure (up to 5 kbar) and discuss volumetric and structural differences between the native and pressure pre-denatured state. Additionally, we have obtained the compressibility of the protein and hydration water and compare these data with experimental results. As water plays a crucial role in determining the structure, dynamics and function of proteins, we undertook a detailed analysis of the structure of the interfacial water and the protein-solvent H-bond network as well. Moreover, we report here also MD results on the temperature-induced unfolding of SNase. The time evolution of the protein volume and solvent accessible surface area during thermal unfolding have been investigated, and we present a detailed discussion of the temperature-induced unfolding pathway of SNase in terms of secondary and tertiary structural changes.

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.

Multidimensional separations of ubiquitin conformers in the gas phase: relating ion cross sections to H/D exchange measurements

Investigating gas-phase structures of protein ions can lead to an improved understanding of intramolecular forces that play an important role in protein folding. Both hydrogen/deuterium (H/D) exchange and ion mobility spectrometry provide insight into the structures and stabilities of different gas-phase conformers, but how best to relate the results from these two methods has been hotly debated. Here, high-field asymmetric waveform ion mobility spectrometry (FAIMS) is combined with Fourier-transform ion cyclotron resonance mass spectrometry (FT/ICR MS) and is used to directly relate ubiquitin ion cross sections and H/D exchange extents. Multiple conformers can be identified using both methods. For the 9+ charge state of ubiquitin, two conformers (or unresolved populations of conformers) that have cross sections differing by 10% are resolved by FAIMS, but only one conformer is apparent using H/D exchange at short times. For the 12+ charge state, two conformers (or conformer populations) have cross sections differing by <1%, yet H/D exchange of these conformers differ significantly (6 versus 25 exchanges). These and other results show that ubiquitin ion collisional cross sections and H/D exchange distributions are not strongly correlated and that factors other than surface accessibility appear to play a significant role in determining rates and extents of H/D exchange. Conformers that are not resolved by one method could be resolved by the other, indicating that these two methods are highly complementary and that more conformations can be resolved with this combination of methods than by either method alone.

Using ellipsoids to model charge distributions in gas phase protein complex ion dissociation

Gas phase protein complex dissociation was modelled using ellipsoids on which discrete charges were placed in randomly chosen charge sites. A number of sizes, shapes, orientations, and types of ellipsoids were considered. For each case, charge transfer parameters and electrostatic energies were calculated as a function of the fractional surface area. It was found that to within 10–15%, the charge on product ions after dissociation is distributed according to their fractional surface area. This can imply, for example, that in experiments measuring the dissociation of homodimers, charge asymmetries of greater than 10–15% in the product ions signal that one of the ions has greatly increased its surface area, such as would occur with unfolding. This assumes that proton transfer occurs on a timescale fast enough that the dissociation products adopt the minimum electrostatic energy configuration. Calculation of this energy shows that it is minimized when one or more of the monomers is in an extended conformation.Key words: protein complex dissociation, ellipsoid charge model, charge asymmetry.

Structural details, pathways, and energetics of unfolding apomyoglobin

Protein folding is often difficult to characterize experimentally because of the transience of intermediate states, and the complexity of the protein-solvent system. Atomistic simulations, which could provide more detailed information, have had to employ highly simplified models or high temperatures, to cope with the long time scales of unfolding; direct simulation of folding is even more problematic. We report a fully atomistic simulation of the acid-induced unfolding of apomyoglobin in which the protonation of acidic side-chains to simulate low pH is sufficient to induce unfolding at room temperature with no added biasing forces or other unusual conditions; and the trajectory is validated by comparison to experimental characterization of intermediate states. Novel insights provided by their analysis include: characterization of a dry swollen globule state forming a barrier to initial unfolding or final folding; observation of cooperativity in secondary and tertiary structure formation and its explanation in terms of dielectric environments; and structural details of the intermediate and the completely unfolded states. These insights involve time scales and levels of structural detail that are presently beyond the range of experiment, but come within reach through the simulation methods described here. An implicit solvation model is used to analyze the energetics of protein folding at various pH and ionic strength values, and a reasonable estimate of folding free energy is obtained. Electrostatic interactions are found to disfavor folding.

Charged states of proteins. Reactions of doubly protonated alkyldiamines with NH(3): solvation or deprotonation. Extension of two proton cases to multiply protonated globular proteins observed in the gas phase

The apparent gas-phase basicities (GB(app)'s) of basic sites in multiply protonated molecules, such as proteins, can be approximately predicted. An approach used by Williams and co-workers was to develop an equation for a diprotonated system, NH(3)(CH(2))(7)NH(3)(2+), and then extend it with a summation of pairwise interactions to multiply protonated systems. Experimental determinations of the rates of deprotonation of NH(3)(CH(2))(7)NH(3)(2+) by a variety of bases B, in the present work, showed that GB(app) = GB(NH(3)) = 196 kcal/mol. This result is supported also by determinations of the equilibria: NH(3)(CH(2))(p)NH(3)(2+) + NH(3) = NH(3)(CH(2))(p)NH(3) x NH(3)(2+), for p = 7, 8, 10, 12. The described experimental GB(app) is 14 kcal/mol higher than the value predicted by the equation used by Williams and co-workers but in agreement with an ab initio result by Gronert. Equations based on electrostatics are developed for the two proton and multiproton systems which allow the evaluation of GB(app) of the basic sites on proteins. These are applied for the evaluation of GB(app) of the basic sites and of N(SB), the maximum number of protons that the nondenatured proteins, carbonic anhydrase (CAII), cytochrome c (CYC), and pepsin, can hold. The N(SB) values are compared with the observed charges, Z(obs)'s, when the nondenatured proteins are produced by electrospray and found in agreement with the proposal by de la Mora that Z(obs) is determined by the number of charges provided by the droplet that contains the protein, according to the charge residue model (CRM). The GB(app) values of proteins have many other applications. They can be compared with experimental measurements and are also needed for the understanding of the thermal denaturing of charged proteins and the thermal dissociation of charged protein complexes.

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.

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.