Correction of both NBD1 energetics and domain interface is required to restore ΔF508 CFTR folding and function.

The folding and misfolding mechanism of multidomain proteins remains poorly understood. Although thermodynamic instability of the first nucleotide-binding domain (NBD1) of ΔF508 CFTR (cystic fibrosis transmembrane conductance regulator) partly accounts for the mutant channel degradation in the endoplasmic reticulum and is considered as a drug target in cystic fibrosis, the link between NBD1 and CFTR misfolding remains unclear. Here, we show that ΔF508 destabilizes NBD1 both thermodynamically and kinetically, but correction of either defect alone is insufficient to restore ΔF508 CFTR biogenesis. Instead, both ΔF508-NBD1 energetic and the NBD1-MSD2 (membrane-spanning domain 2) interface stabilization are required for wild-type-like folding, processing, and transport function, suggesting a synergistic role of NBD1 energetics and topology in CFTR-coupled domain assembly. Identification of distinct structural deficiencies may explain the limited success of ΔF508 CFTR corrector molecules and suggests structure-based combination corrector therapies. These results may serve as a framework for understanding the mechanism of interface mutation in multidomain membrane proteins.

Mechanism of Protein Aggregation Inhibition by Arginine: Blockage of Anionic Side Chains Favors Unproductive Encounter Complexes.

Aggregation refers to the assembly of proteins into nonphysiological higher order structures. While amyloid has been studied extensively, much less is known about amorphous aggregation, a process that interferes with protein expression and storage. Free arginine (Arg+) is a widely used aggregation inhibitor, but its mechanism remains elusive. Focusing on myoglobin (Mb), we recently applied atomistic molecular dynamics (MD) simulations for gaining detailed insights into amorphous aggregation (Ng J. Phys. Chem. B 2021, 125, 13099). Building on that approach, the current work for the first time demonstrates that MD simulations can directly elucidate aggregation inhibition mechanisms. Comparative simulations with and without Arg+ reproduced the experimental finding that Arg+ significantly decreased the Mb aggregation propensity. Our data reveal that, without Arg+, protein-protein encounter complexes readily form salt bridges and hydrophobic contacts, culminating in firmly linked dimeric aggregation nuclei. Arg+ promotes the dissociation of encounter complexes. These "unproductive" encounter complexes are favored because Arg+ binding to D- and E- lowers the tendency of these anionic residues to form interprotein salt bridges. Side chain blockage is mediated largely by the guanidinium group of Arg+, which binds carboxylates through H-bond-reinforced ionic contacts. Our MD data revealed Arg+ self-association into a dynamic quasi-infinite network, but we found no evidence that this self-association is important for protein aggregation inhibition. Instead, aggregation inhibition by Arg+ is similar to that mediated by free guanidinium ions. The computational strategy used here should be suitable for the rational design of aggregation inhibitors with enhanced potency.

Atomistic Details of Peptide Reversed-Phase Liquid Chromatography from Molecular Dynamics Simulations.

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

On the Chemistry of Aqueous Ammonium Acetate Droplets during Native Electrospray Ionization Mass Spectrometry.

Ammonium acetate (NH4Ac) is a widely used solvent additive in native electrospray ionization (ESI) mass spectrometry. NH4Ac can undergo proton transfer to form ammonia and acetic acid (NH4+ + Ac- → NH3 + HAc). The volatility of these products ensures that electrosprayed ions are free of undesired adducts. NH4Ac dissolution in water yields pH 7, providing "physiological" conditions. However, NH4Ac is not a buffer at pH 7 because NH4+ and Ac- are not a conjugate acid/base pair (Konermann, L. J. Am. Soc. Mass Spectrom. 2017, 28, 1827-1835.). In native ESI, it is desirable that analytes experience physiological conditions not only in bulk solution but also while they reside in ESI droplets. Little is known about the internal milieu of NH4Ac-containing ESI droplets. The current work explored the acid/base chemistry of such droplets, starting from a pH 7 analyte solution. We used a two-pronged approach involving evaporation experiments on bulk solutions under ESI-mimicking conditions, as well as molecular dynamics simulations using a newly developed algorithm that allows for proton transfer. Our results reveal that during droplet formation at the tip of the Taylor cone, electrolytically generated protons get neutralized by Ac-, making NH4+ the net charge carriers in the weakly acidic nascent droplets. During the subsequent evaporation, the droplets lose water as well as NH3 and HAc that were generated by proton transfer. NH3 departs more quickly because of its greater volatility, causing the accumulation of HAc. Together with residual Ac-, these HAc molecules form an acetate buffer that stabilizes the average droplet pH at 5.4 ± 0.1, as governed by the Henderson-Hasselbalch equation. The remarkable success of native ESI investigations in the literature implies that this pH drop by ∼1.6 units relative to the initially neutral analyte solution can be tolerated by most biomolecular analytes on the short time scale of the ESI process.

Synergistic recruitment of UbcH7~Ub and phosphorylated Ubl domain triggers parkin activation.

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

Formation of Gaseous Peptide Ions from Electrospray Droplets: Competition between the Ion Evaporation Mechanism and Charged Residue Mechanism.

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 via droplet evaporation to dryness which represents the hallmark of the charged residue mechanism (CRM). Here, we resolved this conundrum by performing MD simulations on a common model peptide (bradykinin) in Rayleigh-charged aqueous droplets. The primary focus was on pH 2 conditions (bradykinin solution charge = 3+), but we also verified that our MD strategy captured pH-dependent charge state shifts seen in ESI-MS experiments. In agreement with earlier simulations, we found that droplets with initial radii of 1.5-3 nm predominantly release peptide ions via the CRM. In contrast, somewhat larger radii (4-5 nm) favor IEM behavior. It appears that these are the first MD data to unequivocally demonstrate the viability of peptide IEM events. Electrostatic arguments can account for the observed droplet size dependence. In summary, both CRM and IEM can be operative in peptide ESI-MS. The prevalence of one over the other mechanism depends on the droplet size distribution in the ESI plume.

Mechanism of Magic Number NaCl Cluster Formation from Electrosprayed Water Nanodroplets.

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 Na14Cl13+ MNC. Simulations of ESI droplets showed nonspecific association of Na+ and Cl-, culminating in gaseous clusters via solvent evaporation to dryness (charged residue mechanism). These nascent clusters did not show any preference for MNCs. In mass spectrometry experiments, analyte ions undergo in-source activation prior to detection. We emulated in-source activation by heating nascent clusters in our MD runs. Heating triggered structural fluctuations and dissociation events, generating MNC-dominated product distributions. Why are MNCs preferred after in-source activation? Thermally excited clusters frequently adopt structures consisting of a preformed MNC and a stringlike protrusion that contains the surplus ions. Facile separation of these protrusions releases the MNC (Clusterhot → MNC-protrusion → MNC + protrusion). This work marks the first time that MD simulations were able to capture cluster self-assembly with subsequent "molecular pruning", generating MNC-dominated product distributions that agree with experiments.

Structural Dynamics of a Thermally Stressed Monoclonal Antibody Characterized by Temperature-Dependent H/D Exchange Mass Spectrometry.

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. Anal. Chem. 2020, 92, 10058). The current work explores the applicability of this HDX framework to the NIST reference mAb (NISTmAb), a large multidomain protein that is prone to aggregation and has three melting points. Using global fitting, we were able to model HDX profiles across the NISTmAb sequence between zero and 95 °C, and for time points between 15 s and 20 min. We uncovered the enthalpic and entropic contributions of local fluctuations that govern the conformational dynamics at low temperatures. The CH2 and CH3 domains were found to be increasingly affected by global unfolding/refolding in the vicinity of their melting points, although the transiently unfolded protein displayed significant residual protection. Global dynamics were not involved in the deuteration of the Fab domains (which have the highest melting point). Instead, global Fab unfolding was followed immediately by irreversible aggregation. Our results reveal that the thermodynamic HDX-MS strategy applied in this work is well suited for probing spatially resolved dynamics of thermally stressed large proteins such as mAbs, complementing data obtained by DSC.

Recommendations for performing, interpreting and reporting hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments.

Hydrogen deuterium exchange mass spectrometry (HDX-MS) is a powerful biophysical technique being increasingly applied to a wide variety of problems. As the HDX-MS community continues to grow, adoption of best practices in data collection, analysis, presentation and interpretation will greatly enhance the accessibility of this technique to nonspecialists. Here we provide recommendations arising from community discussions emerging out of the first International Conference on Hydrogen-Exchange Mass Spectrometry (IC-HDX; 2017). It is meant to represent both a consensus viewpoint and an opportunity to stimulate further additions and refinements as the field advances.

Atomistic Insights into the Formation of Nonspecific Protein Complexes during Electrospray Ionization.

Native electrospray ionization (ESI)-mass spectrometry (MS) is widely used for the detection and characterization of multi-protein complexes. A well-known problem with this approach is the possible occurrence of nonspecific protein clustering in the ESI plume. This effect can distort the results of binding affinity measurements, and it can even generate gas-phase complexes from proteins that are strictly monomeric in bulk solution. By combining experiments and molecular dynamics (MD) simulations, the current work for the first time provides detailed insights into the ESI clustering of proteins. Using ubiquitin as a model system, we demonstrate how the entrapment of more than one protein molecule in an ESI droplet can generate nonspecific clusters (e.g., dimers or trimers) via solvent evaporation to dryness. These events are in line with earlier proposals, according to which protein clustering is associated with the charged residue model (CRM). MD simulations on cytochrome c (which carries a large intrinsic positive charge) confirmed the viability of this CRM avenue. In addition, the cytochrome c data uncovered an alternative mechanism where protein-protein contacts were formed early within ESI droplets, followed by cluster ejection from the droplet surface. This second pathway is consistent with the ion evaporation model (IEM). The observation of these IEM events for large protein clusters is unexpected because the IEM has been thought to be associated primarily with low-molecular-weight analytes. In all cases, our MD simulations produced protein clusters that were stabilized by intermolecular salt bridges. The MD-generated charge states agreed with experiments. Overall, this work reveals that ESI-induced protein clustering does not follow a tightly orchestrated pathway but can proceed along different avenues.

Hydrogen/Deuterium Exchange Measurements May Provide an Incomplete View of Protein Dynamics: a Case Study on Cytochrome c.

Many aspects of protein function rely on conformational fluctuations. Hydrogen/deuterium exchange (HDX) mass spectrometry (MS) provides a window into these dynamics. Despite the widespread use of HDX-MS, it remains unclear whether this technique provides a truly comprehensive view of protein dynamics. HDX is mediated by H-bond-opening/closing events, implying that HDX methods provide an H-bond-centric view. This raises the question if there could be fluctuations that leave the H-bond network unaffected, thereby rendering them undetectable by HDX-MS. We explore this issue in experiments on cytochrome c (cyt c). Compared to the Fe(II) protein, Fe(III) cyt c shows enhanced deuteration on both the distal and proximal sides of the heme. Previous studies have attributed the enhanced dynamics of Fe(III) cyt c to the facile and reversible rupture of the distal M80-Fe(III) bond. Using molecular dynamics (MD) simulations, we conducted a detailed analysis of various cyt c conformers. Our MD data confirm that rupture of the M80-Fe(III) contact triggers major reorientation of the distal Ω loop. Surprisingly, this event takes place with only miniscule H-bonding alterations. In other words, the distal loop dynamics are almost "HDX-silent". Moreover, distal loop movements cannot account for enhanced dynamics on the opposite (proximal) side of the heme. Instead, enhanced deuteration of Fe(III) cyt c is attributed to sparsely populated conformers where both the distal (M80) and proximal (H18) coordination bonds have been ruptured, along with opening of numerous H-bonds on both sides of the heme. We conclude that there can be major structural fluctuations that are only weakly coupled to changes in H-bonding, making them virtually impossible to track by HDX-MS. In such cases, HDX-MS may provide an incomplete view of protein dynamics.

Nrf2, the Major Regulator of the Cellular Oxidative Stress Response, is Partially Disordered.

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.

Delineating Heme-Mediated versus Direct Protein Oxidation in Peroxidase-Activated Cytochrome c by Top-Down Mass Spectrometry.

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.

Formation of Gaseous Proteins via the Ion Evaporation Model (IEM) in Electrospray Mass Spectrometry.

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.

Analysis of Temperature-Dependent H/D Exchange Mass Spectrometry Experiments.

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.

Effects of electrospray mechanisms and structural relaxation on polylactide ion conformations in the gas phase: insights from ion mobility spectrometry and molecular dynamics simulations.

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.

Dimerization interface of osteoprotegerin revealed by hydrogen-deuterium exchange mass spectrometry.

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.

Mechanism of Electrospray Supercharging for Unfolded Proteins: Solvent-Mediated Stabilization of Protonated Sites During Chain Ejection.

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.

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

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.

How to run molecular dynamics simulations on electrospray droplets and gas phase proteins: Basic guidelines and selected applications.

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.