Search

Differential scanning calorimetry (DSC) is a standard tool for probing the resilience of monoclonal antibodies (mAbs) and other protein therapeutics against thermal degradation. Unfortunately, DSC usually only provides insights into global unfolding, although sequential steps are sometimes discernible for multidomain proteins. Temperature-dependent hydrogen/deuterium exchange (HDX) mass spectrometry (MS) has the potential to probe heat-induced events at a much greater level of detail. We recently proposed a strategy to deconvolute temperature-dependent HDX data into contributions from local dynamics, global unfolding/refolding, as well as chemical labeling. However, that strategy was validated only for a small protein (Tajoddin, N. N.; Konermann, L.

Excess protons in water exhibit unique transport properties because they can rapidly hop along H-bonded water wires. Considerable progress has been made in unraveling this Grotthuss diffusion mechanism using quantum mechanical-based computational techniques. Unfortunately, high computational cost tends to restrict those techniques to small systems and short times. Molecular dynamics (MD) simulations can be applied to much larger systems and longer time windows. However, standard MD methods do not permit the dissociation/formation of covalent bonds, such that Grotthuss diffusion cannot be captured. Here, we bridge this gap by combining atomistic MD simulations (using Gromacs and TIP4P/2005 water) with proton hopping. Excess protons are modeled as hydronium ions that undergo H

Peptide separations by reversed-phase liquid chromatography (RPLC) are an integral part of bottom-up proteomics. These separations typically employ C18 columns with water/acetonitrile gradient elution in the presence of formic acid. Despite the widespread use of such workflows, the exact nature of peptide interactions with the stationary and mobile phases are poorly understood. Here we employ microsecond molecular dynamics (MD) simulations to uncover details of peptide RPLC. We examined two tryptic peptides, a hydrophobic and a hydrophilic species, in a slit pore lined with C18 chains that were grafted onto SiO2 support. Our simulations explored peptide trapping, followed by desorption and elution. Trapping in an aqueous mobile phase was initiated by C18 contacts with Lys butyl moieties. This was followed by extensive anchoring of nonpolar side chains (Leu/Ile/Val) in the C18 layer. Exposure to water/acetonitrile triggered peptide desorption in a stepwise fashion; charged sites close to the termini were the first to lift off, followed by the other residues. During water/acetonitrile elution, both peptides preferentially resided close to the pore center. The hydrophilic peptide exhibited no contacts with the stationary phase under these conditions. In contrast, the hydrophobic species underwent multiple transient Leu/Ile/Val binding interactions with C18 chains. These nonpolar interactions represent the foundation of differential peptide retention, in agreement with the experimental elution behavior of the two peptides. Extensive peptide/formate ion pairing was observed in water/acetonitrile, particularly at N-terminal sites. Overall, this work uncovers an unprecedented level of RPLC molecular details, paving the way for MD simulations as a future tool for improving retention prediction algorithms, and for the design of novel column materials.

Events taking place during electrospray ionization (ESI) can trigger the self-assembly of various nanoclusters. These products are often dominated by magic number clusters (MNCs) that have highly symmetrical structures. The literature rationalizes the dominance of MNCs by noting their high stability. However, this argument is not necessarily adequate because thermodynamics cannot predict the outcome of kinetically controlled reactions. Thus, the mechanisms responsible for MNC dominance remain poorly understood. Molecular dynamics (MD) simulations can provide atomistic insights into self-assembly reactions, but even this approach has thus far failed to provide pertinent answers. The current work overcomes this limitation. We focused on salt clusters formed from aqueous NaCl solutions during ESI. The corresponding mass spectra are dominated by the Na

The transfer of peptide ions from solution into the gas phase by electrospray ionization (ESI) is an integral component of mass spectrometry (MS)-based proteomics. The mechanisms whereby gaseous peptide ions are released from charged ESI nanodroplets remain unclear. This is in contrast to intact protein ESI, which has been the focus of detailed investigations using molecular dynamics (MD) simulations and other methods. Under acidic liquid chromatography/MS conditions, many peptides carry a solution charge of 3+ or 2+. Because of this pre-existing charge and their relatively small size, prevailing views suggest that peptides follow the ion evaporation mechanism (IEM). The IEM entails analyte ejection from ESI droplets, driven by electrostatic repulsion between the analyte and droplet. Surprisingly, recent peptide MD investigations reported a different behavior, that is, the release of peptide ions

N/A

Oxidation of key residues in cytochrome c (cyt c) by chloramine T (CT) converts the protein from an electron transporter to a peroxidase. This peroxidase-activated state represents an important model system for exploring the early steps of apoptosis. CT-induced transformations include oxidation of the distal heme ligand Met80 (MetO, +16 Da) and carbonylation (LysCHO, -1 Da) in the range of Lys53/55/72/73. Remarkably, the 15 remaining Lys residues in cyt c are not susceptible to carbonylation. The cause of this unusual selectivity is unknown. Here we applied top-down mass spectrometry (MS) to examine whether CT-induced oxidation is catalyzed by heme. To this end, we compared the behavior of cyt c with (holo-cyt c) and without heme (apoSS-cyt c). CT caused MetO formation at Met80 for both holo- and apoSS-cyt c, implying that this transformation can proceed independently of heme. The aldehyde-specific label Girard's reagent T (GRT) reacted with oxidized holo-cyt c, consistent with the presence of several LysCHO. In contrast, oxidized apo-cyt c did not react with GRT, revealing that LysCHO forms only in the presence of heme. The heme dependence of LysCHO formation was further confirmed using microperoxidase-11 (MP11). CT exposure of apoSS-cyt c in the presence of MP11 caused extensive nonselective LysCHO formation. Our results imply that the selectivity of LysCHO formation at Lys53/55/72/73 in holo-cyt c is caused by the spatial proximity of these sites to the reactive (distal) heme face. Overall, this work highlights the utility of top-down MS for unravelling complex oxidative modifications.

H/D exchange (HDX) mass spectrometry (MS) is a widely used technique for interrogating protein structure and dynamics. Backbone HDX is mediated by opening/closing (unfolding/refolding) fluctuations. In traditional HDX-MS, proteins are incubated in D2O as a function of time at constant temperature (T). There is an urgent need to complement this traditional approach with experiments that probe proteins in a T-dependent fashion, e.g., for assessing the stability of therapeutic antibodies. A key problem with such studies is the absence of strategies for interpreting HDX-MS data in the context of T-dependent protein dynamics. Specifically, it has not been possible thus far to separate T-induced changes of the chemical labeling step (kch) from thermally enhanced protein fluctuations. Focusing on myoglobin, the current work solves this problem by dissecting T-dependent HDX-MS profiles into contributions from kch(T), as well as local and global protein dynamics. Experimental profiles started off with surprisingly shallow slopes that seemed to defy the quasi-exponential kch(T) dependence. Just below the melting temperature (Tm) the profiles showed a sharp increase. Our analysis revealed that local dynamics dominate at low T, while global events become prevalent closer to Tm. About half of the backbone NH sites exhibited a canonical scenario, where local opening/closing was associated with positive ΔH and ΔS. Many of the remaining sites had negative ΔH and ΔS, thereby accounting for the shallowness of the experimental HDX-MS profiles at low T. In summary, this work provides practitioners with the tools to analyze proteins over a wide temperature range, paving the way toward T-dependent high-throughput screening applications by HDX-MS.

Covalent modifications by reactive oxygen species can modulate the function and stability of proteins. Thermal unfolding experiments in solution are a standard tool for probing oxidation-induced stability changes. Complementary to such solution investigations, the stability of electrosprayed protein ions can be assessed in the gas phase by collision-induced unfolding (CIU) and ion-mobility spectrometry. A question that remains to be explored is whether oxidation-induced stability alterations in solution are mirrored by the CIU behavior of gaseous protein ions. Here, we address this question using chloramine-T-oxidized cytochrome c (CT-cyt c) as a model system. CT-cyt c comprises various proteoforms that have undergone MetO formation (+16 Da) and Lys carbonylation (LysCH2-NH2 → LysCHO, -1 Da). We found that CT-cyt c in solution was destabilized, with a ∼5 °C reduced melting temperature compared to unmodified controls. Surprisingly, CIU experiments revealed the opposite trend, i.e., a stabilization of CT-cyt c in the gas phase. To pinpoint the source of this effect, we performed proteoform-resolved CIU on CT-cyt c fractions that had been separated by cation exchange chromatography. In this way, it was possible to identify MetO formation at residue 80 as the key modification responsible for stabilization in the gas phase. Possibly, this effect is caused by newly formed contacts of the sulfoxide with aromatic residues in the protein core. Overall, our results demonstrate that oxidative modifications can affect protein stability in solution and in the gas phase very differently.

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 La3+ 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 La3+ 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 La3+ showed a very different behavior. The trivalent nature of La3+ 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 La3+ 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.

Proteins that encounter unfavorable solvent conditions are prone to aggregation, a phenomenon that remains poorly understood. This work focuses on myoglobin (Mb) as a model protein. Upon heating, Mb produces amorphous aggregates. Thermal unfolding experiments at low concentration (where aggregation is negligible), along with centrifugation assays, imply that Mb aggregation proceeds via globally unfolded conformers. This contrasts studies on other proteins that emphasized the role of partially folded structures as aggregate precursors. Molecular dynamics (MD) simulations were performed to gain insights into the mechanism by which heat-unfolded Mb molecules associate with one another. A prerequisite for these simulations was the development of a method for generating monomeric starting structures. Periodic boundary condition artifacts necessitated the implementation of a partially immobilized water layer lining the walls of the simulation box. Aggregation simulations were performed at 370 K to track the assembly of monomeric Mb into pentameric species. Binding events were preceded by multiple unsuccessful encounters. Even after association, protein-protein contacts remained in flux. Binding was mediated by hydrophobic contacts, along with salt bridges that involved hydrophobically embedded Lys residues. Overall, this work illustrates that atomistic MD simulations are well suited for garnering insights into protein aggregation mechanisms.

Nuclear factor erythroid 2-related factor 2 (Nrf2) is a transcription regulator that plays a pivotal role in coordinating the cellular response to oxidative stress. Through interactions with other proteins, such as Kelch-like ECH-associated protein 1 (Keap1), CREB-binding protein (CBP), and retinoid X receptor alpha (RXRα), Nrf2 mediates the transcription of cytoprotective genes critical for removing toxicants and preventing DNA damage, thereby playing a significant role in chemoprevention. Dysregulation of Nrf2 is linked to tumorigenesis and chemoresistance, making Nrf2 a promising target for anticancer therapeutics. However, despite the physiological importance of Nrf2, the molecular details of this protein and its interactions with most of its targets remain unknown, hindering the rational design of Nrf2-targeted therapeutics. With this in mind, we used a combined bioinformatics and experimental approach to characterize the structure of full-length Nrf2 and its interaction with Keap1. Our results show that Nrf2 is partially disordered, with transiently structured elements in its Neh2, Neh7, and Neh1 domains. Moreover, interaction with the Kelch domain of Keap1 leads to protection of the binding motifs in the Neh2 domain of Nrf2, while the rest of the protein remains highly dynamic. This work represents the first detailed structural characterization of full-length Nrf2 and provides valuable insights into the molecular basis of Nrf2 activity modulation in oxidative stress response.

Proteins that are unfolded in solution produce higher charge states during electrospray ionization (ESI) than their natively folded counterparts. Protein charge states can be further increased by the addition of supercharging agents (SCAs) such as sulfolane. The mechanism whereby these supercharged [M + zH] z+ ions are formed under unfolded conditions remains unclear. Here we employed a combination of mass spectrometry (MS), ion mobility spectrometry (IMS), and molecular dynamics (MD) simulations for probing the ESI mechanism under denatured supercharging conditions. ESI of acid-unfolded apo-myoglobin (aMb) in the presence of sulfolane produced charge states around 27+, all the way to fully protonated (33+) aMb. MD simulations of aMb 27+ to 33+ in Rayleigh-charged water/sulfolane droplets culminated in electrostatically driven protein expulsion, consistent with the chain ejection model (CEM). The electrostatically stretched conformations predicted by these simulations were in agreement with IMS experiments. The CEM involves partitioning of mobile H+ between the droplet and the departing protein. Our results imply that supercharging of unfolded proteins is caused by residual sulfolane that stabilizes protonated sites on the protruding chains, thereby promoting H+ retention on the protein. The stabilization of charged sites is due to charge-dipole interactions mediated by the large dipole moment and the low volatility of sulfolane. Support for this mechanism comes from the experimental observation of sulfolane adducts on the most highly charged ions, a phenomenon previously noted by Venter ( J. Am. Soc. Mass Spectrom. 2012, 23, 489-497). The "CEM supercharging model" proposed here for unfolded proteins is distinct from the charge trapping mechanism believed to be operative during native ESI supercharging.

Experiments and molecular dynamics (MD) simulations in the literature indicate that gaseous proteins generated by electrospray ionization (ESI) can retain native-like structures. However, the exact properties of these ions remain to be explored. Focusing on ubiquitin and lysozyme, we examined several pertinent questions. (1) We applied solvent MD runs to test whether the X-ray structures of both proteins are affected by crystal packing. Main and side-chain orientations were retained in solution, providing a justification for the hitherto unscrutinized approach of relying on crystal data for "solution" versus gas-phase comparisons. (2) Most earlier gas-phase protein MD investigations employed short (ns) simulation windows. By extending this time frame to 1 μs, we were able to observe rare unfolding/folding transitions in ubiquitin. These predicted fluctuations were consistent with a semi-unfolded subpopulation detected by ion mobility spectrometry (IMS). (3) Most earlier modeling studies did not account for the high H+ mobility in gaseous proteins. For the first time, we compared static and mobile H+ simulations, focusing on both positively and negatively charged ions. The MD runs revealed a strong preference for retention of a solution-like backbone fold, whereas titratable/polar side chains collapsed onto the protein surface. This side-chain collapse was caused by a multitude of intramolecular salt bridges, H-bonds, and charge-dipole interactions. Our results generalize the findings of Steinberg et al. ( ChemBioChem, 2008, 9, 2417-2423) who had first proposed the occurrence of such side-chain contacts on the basis of short-term simulations with static H+. (4) Calculated collision cross sections of the MD conformers were in close agreement with IMS experiments. Overall, this study supports the view that solution-like protein structures can be retained because of kinetic trapping on the time scale of typical ESI-IMS experiments.

Electrospray ionization (ESI) is an essential technique for transferring proteins from solution into the gas phase for mass spectrometry and ion mobility spectrometry. The mechanisms whereby [M + zH]z+ protein ions are released from charged nanodroplets during ESI have been controversial for many years. Here we discuss recent computational and experimental studies that have shed light on many of the mysteries in this area. Four types of protein ESI experiments can be distinguished, each of which appears to be associated with a specific mechanism. (i) Native ESI proceeds according to the charged residue model (CRM) that entails droplet evaporation to dryness, generating compact protein ions in low charge states. (ii) Native ESI supercharging is also a CRM process, but the dried-out proteins accumulate additional charge because supercharging agents such as sulfolane interfere with the ejection of small ions (Na+, NH4+, etc.) from the shrinking droplets. (iii) Denaturing ESI follows the chain ejection model (CEM), where protein ions are gradually expelled from the droplet surface. H+ equilibration between the droplets and the protruding chains culminates in highly charged gaseous proteins, analogous to the collision-induced dissociation of multi-protein complexes. (iv) Denatured ESI supercharging also generates protein ions via the CEM. Supercharging agents stabilize protonated sites on the protein tail via charge-dipole interactions, causing the chain to acquire additional charge. There will likely be scenarios that fall outside of these four models, but it appears that the framework outlined here covers most of the experimentally relevant conditions.

The peroxidase activity of cytochrome c (cyt c) plays a key role during apoptosis. Peroxidase catalysis requires a vacant Fe coordination site, i.e., cyt c must undergo an activation process involving structural changes that rupture the native Met80-Fe contact. A common strategy for dissociating this bond is the conversion of Met80 to sulfoxide (MetO). It is widely believed that this MetO formation in itself is sufficient for cyt c activation. This notion originates from studies on chloramine-T-treated cyt c (CT-cyt c) which represents a standard model for the peroxidase activated state. CT-cyt c is considered to be a "clean" species that has undergone selective MetO formation, without any other modifications. Using optical, chromatographic, and mass spectrometry techniques, the current work demonstrates that CT-induced activation of cyt c is more complicated than previously thought. MetO formation alone results in only marginal peroxidase activity, because dissociation of the Met80-Fe bond triggers alternative ligation scenarios where Lys residues interfere with access to the heme. We found that CT causes not only MetO formation, but also carbonylation of several Lys residues. Carbonylation is associated with -1 Da mass shifts that have gone undetected in the CT-cyt c literature. Proteoforms possessing both MetO and Lys carbonylation exhibit almost fourfold higher peroxidase activity than those with MetO alone. Carbonylation abrogates the capability of Lys to coordinate the heme, thereby freeing up the distal site as required for an active peroxidase. Previous studies on CT-cyt c may have inadvertently examined carbonylated proteoforms, potentially misattributing effects of carbonylation to solely MetO formation.

Previous structural studies of osteoprotegerin (OPG), a crucial negative regulator of bone remodeling and osteoclastogenesis, were mostly limited to the N-terminal ligand-binding domains. It is now known that the three C-terminal domains of OPG also play essential roles in its function by mediating OPG dimerization, OPG–heparan sulfate (HS) interactions, and formation of the OPG–HS–receptor activator of nuclear factor κB ligand (RANKL) ternary complex. Employing hydrogen–deuterium exchange MS methods, here we investigated the structure of full-length OPG in complex with HS or RANKL in solution. Our data revealed two noteworthy aspects of the OPG structure. First, we found that the interconnection between the N- and C-terminal domains is much more rigid than previously thought, possibly because of hydrophobic interactions between the fourth cysteine-rich domain and the first death domain. Second, we observed that two hydrophobic clusters located in two separate C-terminal domains directly contribute to OPG dimerization, likely by forming a hydrophobic dimerization interface. Aided by site-directed mutagenesis, we further demonstrated that an intact dimerization interface is essential for the biological activity of OPG. Our study represents an important step toward deciphering the structure–function relationship of the full-length OPG protein.

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 [NanI(n-1)]+ clusters. The addition of sulfolane resulted in abundant doubly charged [NanI(n-2)Sulfolanes]2+ 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 [NanIi]z+ and [NanIiSulfolanes]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.

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 [(αβ)2, (αβ)3, and (αβ)4] 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.

Oxidation of key residues in cytochrome

The mechanisms whereby protein ions are released into the gas phase from charged droplets during electrospray ionization (ESI) continue to be controversial. Several pathways have been proposed. For native ESI the charged residue model (CRM) is favored; it entails the liberation of proteins via solvent evaporation to dryness. Unfolded proteins likely follow the chain ejection model (CEM), which involves the gradual expulsion of stretched-out chains from the droplet. According to the ion evaporation model (IEM) ions undergo electrostatically driven desorption from the droplet surface. The IEM is well supported for small precharged species such as Na+. However, it is unclear whether proteins can show IEM behavior as well. We examined this question using molecular dynamics (MD) simulations, mass spectrometry (MS), and ion mobility spectrometry (IMS) in positive ion mode. Ubiquitin was chosen as the model protein because of its structural stability which allows the protein charge in solution to be controlled via pH adjustment without changing the protein conformation. MD simulations on small ESI droplets (3 nm radius) showed CRM behavior regardless of the protein charge in solution. Surprisingly, many MD runs on larger droplets (5.5 nm radius) culminated in IEM ejection of ubiquitin, as long as the protein carried a sufficiently large positive solution charge. MD simulations predicted that nonspecific salt adducts are less prevalent for IEM-generated protein ions than for CRM products. This prediction was confirmed experimentally. Also, collision cross sections of MD structures were in good agreement with IMS data. Overall, this work reveals that the CRM, CEM, and IEM all represent viable pathways for generating gaseous protein ions during ESI. The IEM is favored for proteins that are tightly folded and highly charged in solution and for droplets in a suitable size regime.

Covalent modifications by reactive oxygen species can modulate the function and stability of proteins. Thermal unfolding experiments in solution are a standard tool for probing oxidation-induced stability changes. Complementary to such solution investigations, the stability of electrosprayed protein ions can be assessed in the gas phase by collision-induced unfolding (CIU) and ion-mobility spectrometry. A question that remains to be explored is whether oxidation-induced stability alterations in solution are mirrored by the CIU behavior of gaseous protein ions. Here, we address this question using chloramine-T-oxidized cytochrome

The properties of electrosprayed protein ions continue to be enigmatic, owing to the absence of high-resolution structure determination methods in the gas phase. There is considerable evidence that under properly optimized conditions these ions preserve solution-like conformations and interactions. However, it is unlikely that these solution-like conformers represent the "intrinsic" structural preferences of gaseous proteins. In an effort to uncover what such intrinsically preferred conformers might look like, we performed molecular dynamics (MD) simulations of gaseous ubiquitin. Our work was inspired by recent gas phase experiments, where highly extended 13+ ubiquitin ions were transformed to compact 3+ species by proton stripping (Laszlo, K. J.; Munger, E. B.; Bush, M. F. J. Am. Chem. Soc. 2016, 138, 9581-9588). Our simulations covered several microseconds and used a mobile-proton algorithm to account for the fact that a H+ in gaseous proteins can migrate between different titratable sites. Proton stripping caused folding of ubiquitin into heterogeneous "inside-out" structures. The hydrophilic core of these conformers was stabilized by charge-charge and polar interactions, while hydrophobic residues were located on the protein surface. Collision cross sections of these MD structures were in good agreement with experimental results. The inside-out structures generated during gas phase folding are in striking contrast to the solution behavior which is dominated by the hydrophobic effect, i.e., the tendency to bury hydrophobic side chains in the core (instead of exposing them to the surface). We do not dispute that native-like proteins can be transferred into the gas phase as kinetically trapped species. However, those metastable conformers do not represent the intrinsic structural preferences of gaseous proteins. Our work for the first time provides detailed insights into the properties of intrinsically preferred gas phase conformers, and we unequivocally find them to have inside-out architectures.

Recent advances in molecular dynamics (MD) simulations have made it possible to examine the behavior of large charged droplets that contain analytes such as proteins or polymers, thereby providing insights into electrospray ionization (ESI) mechanisms. In the present study, we use this approach to investigate the release of polylactide (PLA) ions from water/acetonitrile ESI droplets. We found that cationized gaseous PLA ions can be formed via various competing pathways. Some MD runs showed extrusion and subsequent separation of polymer chains from the droplet, as envisioned by the chain ejection model (CEM). On other occasions the PLA chains remained inside the droplets and were released after solvent evaporation to dryness, consistent with the charge residue model (CRM). Following their release from ESI droplets, the nascent gaseous PLA ions were subjected to structural relaxation for several μs in vacuo. The MD conformations generated in this way for various PLA charge states compared favorably to experimental results obtained by ion mobility spectrometry-mass spectrometry (IMS-MS). The structures of all PLA ions evolved during relaxation in the gas phase. However, some macroion species retained features that resembled their nascent structures. For this subset of ions, the IMS-MS response appears to be strongly correlated with the ESI release mechanism (CEM vs. CRM). The former favored extended structures, whereas the latter preferentially generated compact conformers.

F1-ATPase uses ATP hydrolysis to drive rotation of the γ subunit. The γ C-terminal helix constitutes the rotor tip that is seated in an apical bearing formed by α3β3. It remains uncertain to what extent the γ conformation during rotation differs from that seen in rigid crystal structures. Existing models assume that the entire γ subunit participates in every rotation. Here we interrogated E. coli F1-ATPase by hydrogen–deuterium exchange (HDX) mass spectrometry. Rotation of γ caused greatly enhanced deuteration in the γ C-terminal helix. The HDX kinetics implied that most F1 complexes operate with an intact rotor at any given time, but that the rotor tip is prone to occasional unfolding. A molecular dynamics (MD) strategy was developed to model the off-axis forces acting on γ. MD runs showed stalling of the rotor tip and unfolding of the γ C-terminal helix. MD-predicted H-bond opening events coincided with experimental HDX patterns. Our data suggest that in vitro operation of F1-ATPase is associated with s...

The ion evaporation model (IEM) and the charged residue model (CRM) represent cornerstones of any discussion related to the mechanism of electrospray ionization (ESI). Molecular dynamics (MD) simulations have confirmed that small ions such as Na+ are ejected from the surface of aqueous ESI droplets (IEM), while folded proteins in native ESI are released by water evaporation to dryness (CRM). ESI of unfolded proteins yields [M + zH] z+ ions that are much more highly charged than their folded counterparts. A chain ejection model (CEM) has been proposed to account for the protein ESI behavior under such non-native conditions (Konermann, L., et al. Anal. Chem. 2013, 85, 2-9). The CEM envisions that unfolded proteins are driven to the droplet surface by hydrophobic and electrostatic factors, followed by gradual ejection via intermediates where droplets carry extended protein tails. Thus far, it has not been possible to support the CEM through MD simulations using realistic protein models and atomistic force fields. Such endeavors require much larger droplets than in previous MD studies. Also, the incorporation of CEM-related H+ migration is difficult. This work overcomes these challenges in MD simulations on unfolded apo-myoglobin (aMb) in droplets with a 5.5 nm radius (∼22500 water molecules). We focused on solutions at pH ∼4 where the aMb solution charge coincides with the charge on some of the electrosprayed ions (22+ to 27+), such that H+ migration could be neglected. Na+ ions were added to ensure a droplet charge close to the Rayleigh limit. We found that 16 of 17 MD runs on various protonation patterns produced [M + zH] z+ ions via chain ejection. The predicted stretched-out aMb conformations were consistent with experimental collision cross sections. These results support the view that unfolded proteins follow the CEM. Overall, the IEM/CRM/CEM triad can account for a wide range of ESI scenarios involving various types of analytes.

Recent molecular dynamics (MD) simulations from various laboratories have advanced the general understanding of electrospray ionization (ESI)-related processes. Unfortunately, computational cost has limited most of those previous endeavors to ESI droplets with radii of ∼3 nm or less, which represent the low end of the size distribution in the ESI plume. The current work extends this range by conducting simulations on aqueous ESI droplets with radii of 5.5 nm (∼23 000 water molecules). Considering that computational cost increases with r

"Native" electrospray ionization (ESI) mass spectrometry (MS) aims to transfer proteins from solution into the gas phase while maintaining solution-like structures and interactions. The ability to control the charge states of protein ions produced in these experiments is of considerable importance. Supercharging agents (SCAs) such as sulfolane greatly elevate charge states without significantly affecting the protein structure in bulk aqueous solution. The origin of native ESI supercharging remains contentious. According to one model, SCAs trigger unfolding within ESI droplets. In contrast, the "charge trapping model" envisions that SCAs impede the ejection of charge carriers (e.g., NH4+ or Na+) from the droplet. We addressed this controversy experimentally and computationally by employing 18C6 crown ether as a mechanistic probe in native ESI-MS experiments on holo-myoglobin. Remarkably, 18C6 suppressed the supercharging capability of sulfolane. Molecular dynamics (MD) simulations reproduced the experimental charge states. The MD data revealed that 18C6 altered the location of charge carriers in the ESI droplets. Without 18C6, sulfolane covered the droplets in an ionophobic layer that impeded charge carrier access to the surface. In contrast, 18C6 complexation caused charge carrier enrichment in this surface layer, thereby promoting charge ejection. For late droplets, all the water had left and the protein was encapsulated in sulfolane; charge ejection at this stage continued only in the presence of 18C6. As a result, evaporation to dryness of charge-depleted water/sulfolane/18C6 droplets produced low protein charge states, whereas charge-abundant water/sulfolane droplets generated high charge states. Our data support the view that native ESI supercharging is caused by charge trapping. Unfolding within the droplet may play an ancillary role under some conditions, but for the cases examined here, protein structural changes are not a causative factor for supercharging. Our conclusions are bolstered by dendrimer supercharging experiments.

In addition to serving as respiratory electron shuttle, ferri-cytochrome c (cyt c) acts as a peroxidase; i.e., it catalyzes the oxidation of organic substrates by H2O2. This peroxidase function plays a key role during apoptosis. Typical peroxidases have a five-coordinate heme with a vacant distal coordination site that permits the iron center to interact with H2O2. In contrast, native cyt c is six-coordinate, as the distal coordination site is occupied by Met80. It thus seems counterintuitive that native cyt c would exhibit peroxidase activity. The current work scrutinizes the origin of this structure–function mismatch. Cyt c-catalyzed peroxidase reactions show an initial lag phase that is consistent with the in situ conversion of a precatalyst to an active peroxidase. Using mass spectrometry, we demonstrate the occurrence of cyt c self-oxidation in the presence of H2O2. The newly generated oxidized proteoforms are shown to possess significantly enhanced peroxidase activity. H2O2-induced modifications com...

Allosteric proteins possess dynamically coupled residues for the propagation of input signals to distant target binding sites. The input signals usually correspond to "effector is present" or "effector is not present". Many aspects of allosteric regulation remain incompletely understood. This work focused on S100A11, a dimeric EF-hand protein with two hydrophobic target binding sites. An annexin peptide (Ax) served as the target. Target binding is allosterically controlled by Ca2+ over a distance of ∼26 Å. Ca2+ promotes formation of a [Ca4 S100 Ax2] complex, where the Ax peptides are accommodated between helices III/IV and III'/IV'. Without Ca2+ these binding sites are closed, precluding interactions with Ax. The allosteric mechanism was probed by microsecond MD simulations in explicit water, complemented by hydrogen exchange mass spectrometry (HDX/MS). Consistent with experimental data, MD runs in the absence of Ca2+ and Ax culminated in target binding site closure. In simulations on [Ca4 S100] the target binding sites remained open. These results capture the essence of allosteric control, revealing how Ca2+ prevents binding site closure. Both HDX/MS and MD data showed that the metalation sites become more dynamic after Ca2+ loss. However, these enhanced dynamics do not represent the primary trigger of the allosteric cascade. Instead, a labile salt bridge acts as an incessantly active "agitator" that destabilizes the packing of adjacent residues, causing a domino chain of events that culminates in target binding site closure. This agitator represents the starting point of the allosteric signal propagation pathway. Ca2+ binding rigidifies elements along this pathway, thereby blocking signal transmission. This blocking mechanism does not conform to the commonly held view that allosteric communication pathways generally originate at the sites where effectors interact with the protein.

It is believed that electrosprayed proteins and protein complexes can retain solution-like conformations in the gas phase. However, the lack of high-resolution structure determination methods for gaseous protein ions implies that their properties remain poorly understood. Many practitioners tackle this difficulty by complementing mass spectrometry-based experiments with molecular dynamics (MD) simulations. It is a potential problem that the standard MD force fields used for this purpose (such as OPLS-AA/L and CHARMM) were optimized for solution conditions. The question whether these force fields produce meaningful gas-phase data has received surprisingly little attention. Standard force fields are overpolarized to account for an aqueous environment, i.e., atomic charges and intramolecular dipole moments are ∼20% larger than predicted by gas-phase ab initio methods. Here, we examined the implications of this overpolarization by conducting a series of MD simulations on electrosprayed proteins. Force fields were modified via a charge scaling factor (CSF), while ensuring that the net protein charge remained unchanged. CSF = 0.8 should roughly eliminate water-associated overpolarization. Gas-phase CHARMM simulations on myoglobin with CSF = 0.8 and with unmodified parameters (CSF = 1) yielded similar results, preserving a compact structure that was consistent with ion mobility experiments. Major structural changes caused by weakened charge-dipole and dipole-dipole contacts occurred only when lowering CSF to physically unreasonable values (0.5 and 0.1). Similar results were obtained in mobile-proton OPLS-AA/L simulations on the collision-induced dissociation of transthyretin. Our data support the view that gas-phase MD simulations with standard (solution) force fields are suitable for modeling gaseous protein ions in a semiquantitative manner. Although this is welcome news for the mass spectrometry community, it is hoped that dedicated gas-phase MD force fields will become available in the near future.

The E3 ligase parkin ubiquitinates outer mitochondrial membrane\ud proteins during oxidative stress and is linked to early-onset\ud Parkinson’s disease. Parkin is autoinhibited but is activated by the\ud kinase PINK1 that phosphorylates ubiquitin leading to parkin\ud recruitment, and stimulates phosphorylation of parkin’s N-terminal\ud ubiquitin-like (pUbl) domain. How these events alter the\ud structure of parkin to allow recruitment of an E2~Ub conjugate\ud and enhanced ubiquitination is an unresolved question. We\ud present a model of an E2~Ub conjugate bound to the phosphoubiquitin-loaded\ud C-terminus of parkin, derived from NMR chemical\ud shift perturbation experiments. We show the UbcH7~Ub conjugate\ud binds in the open state whereby conjugated ubiquitin binds to the\ud RING1/IBR interface. Further, NMR and mass spectrometry experiments\ud indicate the RING0/RING2 interface is re-modelled,\ud remote from the E2 binding site, and this alters the reactivity of\ud the RING2(Rcat) catalytic cysteine, needed for ubiquitin transfer.\ud Our experiments provide evidence that parkin phosphorylation\ud and E2~Ub recruitment act synergistically to enhance a weak\ud interaction of the pUbl domain with the RING0 domain and rearrange\ud the location of the RING2(Rcat) domain to drive parkin\ud activity.

Infusion of NaCl solutions into an electrospray ionization (ESI) source produces [Na(n+1)Cl n ]+ and other gaseous clusters. The n = 4, 13, 22 magic number species have cuboid ground state structures and exhibit elevated abundance in ESI mass spectra. Relatively few details are known regarding the mechanisms whereby these clusters undergo collision-induced dissociation (CID). The current study examines to what extent molecular dynamics (MD) simulations can be used to garner insights into the sequence of events taking place during CID. Experiments on singly charged clusters reveal that the loss of small neutrals is the dominant fragmentation pathway. MD simulations indicate that the clusters undergo extensive structural fluctuations prior to decomposition. Consistent with the experimentally observed behavior, most of the simulated dissociation events culminate in ejection of small neutrals ([NaCl] i , with i = 1, 2, 3). The MD data reveal that the prevalence of these dissociation channels is linked to the presence of short-lived intermediates where a relatively compact core structure carries a small [NaCl] i protrusion. The latter can separate from the parent cluster via cleavage of a single Na-Cl contact. Fragmentation events of this type are kinetically favored over other dissociation channels that would require the quasi-simultaneous rupture of multiple electrostatic contacts. The CID behavior of NaCl cluster ions bears interesting analogies to that of collisionally activated protein complexes. Overall, it appears that MD simulations represent a valuable tool for deciphering the dissociation of noncovalently bound systems in the gas phase.

Much remains to be learned about the way in which bound metal ions modulate the response of electrosprayed proteins and protein complexes to collisional excitation. Nonspecific metal adducts can affect the extent of collision-induced unfolding (CIU) and collision-induced dissociation (CID). Here, we examine how Na(+) and Ca(2+) adducts alter the CIU response of monomeric proteins under native electrospray conditions. Both of these metals are commonly encountered in biological samples. Measured collision cross sections are largely independent of metal adduction as long as in-source excitation is minimized. In contrast, under CIU conditions, the metal-adducted proteins are markedly more compact than their metal-free counterparts. This phenomenon is particularly pronounced for Ca(2+) binding, but Na(+) adducts have significant effects as well. Molecular dynamics simulations reproduce the experimentally observed trends. The simulations show that structural expansion of the collisionally unfolded proteins is limited by multidentate metal contacts that restrict the conformational freedom of the polypeptide chains. Multidentate interactions with carboxylates and other electron-rich moieties are to be anticipated for divalent metals such as Ca(2+). It is surprising that Na(+) also engages in multidentate ligation. Electrostatic mapping reveals that the propensity of both Na(+) and Ca(2+) to interact with multiple electron-rich groups is caused by ineffective charge shielding during ion pairing. Despite their compactness, the CIU structures of metalated proteins do not retain native-like elements. Instead, CIU generates inside-out conformations where previously surface-exposed hydrophilic side chains get buried along with most of the metal ions. Our findings caution that the observation of compact conformers after collisional excitation does not imply the survival of solution-like structural features. We also discuss possible implications of adduct-mediated effects for CIU fingerprinting studies.

The protein Chibby (Cby) is an antagonist of the Wnt signaling pathway, where it inhibits the binding between the transcriptional coactivator β-catenin and the Tcf/Lef transcription factors. The 126 residue Cby is partially disordered; its N-terminal half is unstructured while its C-terminal half comprises a coiled-coil domain. Previous structural analyses of Cby using NMR spectroscopy suffered from severe line broadening for residues within the protein's C-terminal half, hindering detailed characterization of the coiled-coil domain. Here, we use hydrogen/deuterium exchange-mass spectrometry (HDX-MS) to examine Cby's C-terminal half. Results reveal that Cby is divided into three structural elements: a disordered N-terminal half, a coiled-coil domain, and a C-terminal unstructured extension consisting of the last ∼ 25 residues (which we term C-terminal extension). A series of truncation constructs were designed to assess the roles of individual structural elements in protein stability and Cby binding to TC-1, a positive regulator of the Wnt signaling pathway. CD and NMR data show that Cby maintains coiled-coil structure upon deletion of either disordered region. NMR and ITC binding experiments between Cby and TC-1 illustrate that the interaction is retained upon deletion of either Cby's N-terminal half or its C-terminal extension. Intriguingly, Cby's C-terminal half alone binds to TC-1 with significantly greater affinity compared to full-length Cby, implying that target binding of the coiled-coil domain is affected by the flanking disordered regions.

Electrospray ionization (ESI) allows the production of intact gas-phase ions from proteins in solution. Nondenaturing solvent conditions usually culminate in low ESI charge states. However, many mass spectrometric applications benefit from protein ions that are more highly charged. One way to boost protein charge is the addition of supercharging agents (SCAs) such as sulfolane or m-nitrobenzyl alcohol (m-NBA) to the aqueous solution. The supercharging mechanism remains controversial. We use molecular dynamics (MD) simulations to examine how SCAs affect the behavior of ESI nanodroplets. Simulations were conducted on myoglobin in water, water/sulfolane, and water/m-NBA. Na(+) ions served as surrogate charge carriers instead of H(+). We focus on conditions where the protein initially adopts its native conformation. MD-generated charge states show remarkable agreement with experimental data. Droplet shrinkage is accompanied by Na(+) ejection, consistent with the ion evaporation model (IEM). The droplets segregate into an outer SCA shell and an aqueous core. This core harbors protein and Na(+). Unfavorable SCA solvation restricts Na(+) access to the droplet surface, thereby impeding IEM ejection. Rapid water loss causes SCA enrichment, ultimately forcing all remaining Na(+) to bind the protein. IEM ejection is no longer feasible after this point, such that the protein becomes supercharged by Na(+) trapping. SCA-free droplets produce lower charge states because the aqueous environment ensures a higher IEM efficiency. For all scenarios examined here, proteins are released via solvent evaporation to dryness, as envisioned by the charged residue model. Our data provide the first atomistic view of the supercharging mechanism.

F

The assembly of proteins into dimers and oligomers is a necessary step for the proper function of transcription factors, muscle proteins, and proteases. In uncontrolled states, oligomerization can also contribute to illnesses such as Alzheimer's disease. The S100 protein family is a group of dimeric proteins that have important roles in enzyme regulation, cell membrane repair, and cell growth. Most S100 proteins have been examined in their homodimeric state, yet some of these important proteins are found in similar tissues implying that heterodimeric molecules can also be formed from the combination of two different S100 members. In this work, we have established co-expression methods in order to identify and quantify the distribution of homo- and heterodimers for four specific pairs of S100 proteins in their calcium-free states. The split GFP trap methodology was used in combination with other GFP variants to simultaneously quantify homo- and heterodimeric S100 proteins in vitro and in living cells. For the specific S100 proteins examined, NMR, mass spectrometry, and GFP trap experiments consistently show that S100A1:S100B, S100A1:S100P, and S100A11:S100B heterodimers are the predominant species formed compared to their corresponding homodimers. We expect the tools developed here will help establish the roles of S100 heterodimeric proteins and identify how heterodimerization might alter the specificity for S100 protein action in cells.

The ability to transfer intact proteins and protein complexes into the gas phase by electrospray ionization (ESI) has opened up numerous mass spectrometry (MS)-based avenues for exploring biomolecular structure and function. However, many details regarding the ESI process and the properties of gaseous analyte ions are difficult to decipher when relying solely on experimental data. Molecular dynamics (MD) simulations can provide additional insights into the behavior of ESI droplets and protein ions. This review is geared primarily towards experimentalists who wish to adopt MD simulations as a complementary research tool. We touch on basic points such as force fields, the choice of a proper water model, GPU-acceleration, possible artifacts, as well as shortcomings of current MD models. Following this technical overview, we highlight selected applications. Simulations on aqueous droplets confirm that "native" ESI culminates in protein ion release via the charged residue model. MD-generated charge states and collision cross sections match experimental data. Gaseous protein ions produced by native ESI retain much of their solution structure. Moving beyond classical fixed-charge algorithms, we discuss a simple strategy that captures the mobile nature of H+ within gaseous biomolecules. These mobile proton simulations confirm the high propensity of gaseous proteins to form salt bridges, as well as the occurrence of charge migration during collision-induced unfolding and dissociation. It is hoped that this review will promote the use of MD simulations in ESI-related research. We also hope to encourage the development of improved algorithms for charged droplets and gaseous biomolecular ions.

It is believed that enzyme catalysis is facilitated by conformational dynamics of the protein scaffold that surrounds the active site, yet the exact nature of catalytically relevant protein motions remains largely unknown. Hydrogen/deuterium exchange (HDX) mass spectrometry (MS) reports on backbone H-bond fluctuations. HDX/MS thus represents a promising avenue for probing the relationship between enzyme dynamics and catalysis. A seemingly straightforward strategy for such studies involves comparative measurements during substrate turnover and in the resting state. We examined the feasibility of this approach using rabbit muscle pyruvate kinase (rM1-PK) which catalyzes the conversion of phosphoenolpyruvate and Mg-ADP to pyruvate and Mg-ATP. HDX/MS revealed that catalytically active rM1-PK undergoes significant rigidification in the active site. This finding is counterintuitive, considering the purported correlation between dynamics and catalysis. Interestingly, virtually the same rigidification was seen upon exposing rM1-PK to substrates or products in the absence of turnover. These data imply that the active site dynamics during turnover are dominated by protein-ligand binding interactions. These interactions stabilize H-bonds in the vicinity of the active site, thereby masking subtle dynamic features that might be uniquely associated with catalysis. Our data uncover an inherent problem with side-by-side turnover/resting state measurements, i.e., the difficulty to design a suitable reference state against which the working enzyme can be compared. Comparative HDX/MS experiments on enzyme dynamics should therefore be interpreted with caution.

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

Gas/water interfaces (such as air bubbles or foam) are detrimental to the stability of proteins, often causing aggregation. This represents a potential problem for industrial processes, for example, the production and handling of protein drugs. Proteins possess surfactant-like properties, resulting in a high affinity for gas/water interfaces. The tendency of previously buried nonpolar residues to maximize contact with the gas phase can cause significant structural distortion. Most earlier studies in this area employed spectroscopic tools that could only provide limited information. Here we use hydrogen/deuterium exchange (HDX) mass spectrometry (MS) for probing the conformational dynamics of the model protein myoglobin (Mb) in the presence of N2 bubbles. HDX/MS relies on the principle that unfolded and/or highly dynamic regions undergo faster deuteration than tightly folded segments. In bubble-free solution Mb displays EX2 behavior, reflecting the occurrence of short-lived excursions to partially unfolded conformers. A dramatically different behavior is seen in the presence of N2 bubbles; EX2 dynamics still take place, but in addition the protein shows EX1 behavior. The latter results from interconversion of the native state with conformers that are globally unfolded and long-lived. These unfolded species likely correspond to Mb that is adsorbed to the surface of gas bubbles. N2 sparging also induces aggregation. To explain the observed behavior we propose a simple model, that is, “semi-unfolded” ↔ “native” ↔ “globally unfolded” → “aggregated”. This model quantitatively reproduces the experimentally observed kinetics. To the best of our knowledge, the current study marks the first exploration of surface denaturation phenomena by HDX/MS.

In addition to serving as respiratory electron shuttle, ferri-cytochrome c (cyt c) acts as a peroxidase; i.e., it catalyzes the oxidation of organic substrates by H

Molecular dynamics (MD) simulations have become a key tool for examining the properties of electrosprayed protein ions. Traditional force fields employ static charges on titratable sites, whereas in reality, protons are highly mobile in gas-phase proteins. Earlier studies tackled this problem by adjusting charge patterns during MD runs. Within those algorithms, proton redistribution was subject to energy minimization, taking into account electrostatic and proton affinity contributions. However, those earlier approaches described (de)protonated moieties as point charges, neglecting charge solvation, which is highly prevalent in the gas phase. Here, we describe a mobile proton algorithm that considers the electrostatic contributions from all atoms, such that charge solvation is explicitly included. MD runs were broken down into 50 ps fixed-charge segments. After each segment, the electrostatics was reanalyzed and protons were redistributed. Challenges associated with computational cost were overcome by devising a streamlined method for electrostatic calculations. Avidin (a 504-residue protein complex) maintained a nativelike fold over 200 ns. Proton transfer and side chain rearrangements produced extensive salt bridge networks at the protein surface. The mobile proton technique introduced here should pave the way toward future studies on protein folding, unfolding, collapse, and subunit dissociation in the gas phase.

Native ESI-MS involves the transfer of intact proteins and biomolecular complexes from solution into the gas phase. One potential pitfall is the occurrence of pH-induced changes that can affect the analyte while it is still surrounded by solvent. Most native ESI-MS studies employ neutral aqueous ammonium acetate solutions. It is a widely perpetuated misconception that ammonium acetate buffers the analyte solution at neutral pH. By definition, a buffer consists of a weak acid and its conjugate weak base. The buffering range covers the weak acid pKa ± 1 pH unit. NH4 + and CH3-COO− are not a conjugate acid/base pair, which means that they do not constitute a buffer at pH 7. Dissolution of ammonium acetate salt in water results in pH 7, but this pH is highly labile. Ammonium acetate does provide buffering around pH 4.75 (the pKa of acetic acid) and around pH 9.25 (the pKa of ammonium). This implies that neutral ammonium acetate solutions electrosprayed in positive ion mode will likely undergo acidification down to pH 4.75 ± 1 in the ESI plume. Ammonium acetate nonetheless remains a useful additive for native ESI-MS. It is a volatile electrolyte that can mimic the solvation properties experienced by proteins under physiological conditions. Also, a drop from pH 7 to around pH 4.75 is less dramatic than the acidification that would take place in pure water. It is hoped that the habit of referring to pH 7 solutions as ammonium acetate “buffer” will disappear from the literature. Ammonium acetate “solution” should be used instead.

Royal jelly (RJ) triggers the development of female honeybee larvae into queens. This effect has been attributed to the presence of major royal jelly protein 1 (MRJP1) in RJ. MRJP1 isolated from royal jelly is tightly associated with apisimin, a 54-residue α-helical peptide that promotes the noncovalent assembly of MRJP1 into multimers. No high-resolution structural data are available for these complexes, and their binding stoichiometry remains uncertain. We examined MRJP1/apisimin using a range of biophysical techniques. We also investigated the behavior of deglycosylated samples, as well as samples with reduced apisimin content. Our mass spectrometry (MS) data demonstrate that the native complexes predominantly exist in a (MRJP14 apisimin4) stoichiometry. Hydrogen/deuterium exchange MS reveals that MRJP1 within these complexes is extensively disordered in the range of residues 20-265. Marginally stable secondary structure (likely antiparallel β-sheet) exists around residues 266-432. These weakly structured regions interchange with conformers that are extensively unfolded, giving rise to bimodal (EX1) isotope distributions. We propose that the native complexes have a "dimer of dimers" quaternary structure in which MRJP1 chains are bridged by apisimin. Specifically, our data suggest that apisimin acts as a linker that forms hydrophobic contacts involving the MRJP1 segment 316VLFFGLV322. Deglycosylation produces large soluble aggregates, highlighting the role of glycans as aggregation inhibitors. Samples with reduced apisimin content form dimeric complexes with a (MRJP12 apisimin1) stoichiometry. The information uncovered in this work will help pave the way toward a better understanding of the unique physiological role played by MRJP1 during queen differentiation.

F0F1 ATP synthase harnesses a transmembrane electrochemical gradient for the production of ATP. When operated in reverse, this multiprotein complex catalyzes ATP hydrolysis. In bacteria, the ε subunit is involved in regulating this ATPase activity. Also, ε is essential for coupling ATP hydrolysis (or synthesis) to proton translocation. The ε subunit consists of a β sandwich and two C-terminal helices, α1 and α2. The protein can switch from a compact fold to an alternate conformation where α1 and α2 are separated, resulting in an extended structure. ε from the thermophile Bacillus PS3 (Tε) binds ATP with high affinity such that this protein may function as an intracellular ATP level sensor. ATP binding to isolated Tε triggers a major conformational transition. Earlier data were interpreted in terms of an ATP + Tεextended → ATP·Tεcompact transition that may mimic aspects of the regulatory switching within F0F1 (Yagi et al. (2007) Proc. Natl. Acad. Sci. U.S.A., 104, 11233–11238). In this work, we employ complementary biophysical techniques for examining the ATP-induced conformational switching of isolated Tε. CD spectroscopy confirmed the occurrence of a large-scale conformational transition upon ATP binding, consistent with the formation of stable helical structure. Hydrogen/deuterium exchange (HDX) mass spectrometry revealed that this transition is accompanied by a pronounced stabilization in the vicinity of the ATP-binding pocket. Surprisingly, dramatic stabilization is also seen in the β8−β9 region, which is remote from the site of ATP interaction. Analytical ultracentrifugation uncovered a previously unrecognized feature of Tε: a high propensity to undergo dimerization in the presence of ATP. Comparison with existing crystallography data strongly suggests that the unexpected β8−β9 HDX protection is due to newly formed protein–protein contacts. Hence, ATP binding to isolated Tε proceeds according to 2ATP + 2Tεextended → (ATP·Tεcompact)2. Implications of this dimerization propensity for the possible role of Tε as an antibiotic target are discussed.

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

Hydrogen/deuterium exchange (HDX) mass spectrometry (MS) is a widely used technique for probing protein structure and dynamics. Exposure to D2O induces the deuteration of backbone N-H groups via a process that involves transient excursions to partially unfolded protein conformers. The resulting mass shifts can be probed by MS, usually in combination with proteolytic digestion and/or electron-based fragmentation. Studies on protein-ligand complexes represent a particularly important HDX/MS application. The prevailing view is that ligand binding should reduce deuteration rates, and it is often expected that this reduction will be most pronounced in the vicinity of the interaction site. Many protein-ligand systems do indeed behave in a fashion that is consistent with this paradigm. In this review we point out that the opposite effect may be encountered as well. Also, mixed scenarios are possible where ligand binding induces elevated HDX rates in some protein regions, whereas rates in other segments are reduced. We present a framework that links ligand-induced changes in HDX kinetics to alterations in the occupancy of excited protein conformers. Spontaneous ligand binding will always lower the free energy of the ground state. In contrast, the corresponding free energy shifts of excited states are largely unpredictable, giving rise to a range of possible HDX responses. "Type 1" scenarios, characterized by a reduction of HDX rates are just as feasible as "Type 2" behavior where deuteration is accelerated. Even "Type 0" phenomena may be encountered, where HDX rates are unaffected by the presence of ligand. Type 0/1/2 scenarios can coexist in the same protein (these terms are not to be confused with the EX1/EX2 expressions which refer to a different aspect of protein HDX). Allosteric effects and ligand-induced protein-protein contacts can affect the outcome of protein-ligand binding studies as well. In summary, comparative HDX measurements conducted in the presence and in the absence ligand provide a detailed fingerprint of biomolecular interactions. However, protein-ligand interactions can elicit a wide range of responses, and the interpretation of binding site mapping experiments may not always be straightforward.