Search

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

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.

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.

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.

Pin1 is a phosphorylation-dependent peptidyl-prolyl isomerase that plays a critical role in mediating protein conformational changes involved in signaling processes related to cell cycle control. Pin1 has also been implicated as being neuroprotective in aging-related neurodegenerative disorders including Alzheimer's disease where Pin1 activity is diminished. Notably, recent proteomic analysis of brain samples from patients with mild cognitive impairment revealed that Pin1 is oxidized and also displays reduced activity. Since the Pin1 active site contains a functionally critical cysteine residue (Cys113) with a low predicted pKa, we hypothesized that Cys113 is sensitive to oxidation. Consistent with this hypothesis, we observed that treatment of Pin1 with hydrogen peroxide results in a 32 Da mass increase, likely resulting from the oxidation of Cys113 to sulfinic acid (Cys-SO2H). This modification results in loss of peptidyl-prolyl isomerase activity. Notably, Pin1 with Cys113 substituted by aspartic acid retains activity and is no longer sensitive to oxidation. Structural studies by X-ray crystallography revealed increased electron density surrounding Cys113 following hydrogen peroxide treatment. At lower concentrations of hydrogen peroxide, oxidative inhibition of Pin1 can be partially reversed by treatment with dithiothreitol, suggesting that oxidation could be a reversible modification with a regulatory role. We conclude that the loss of Pin1 activity upon oxidation results from oxidative modification of the Cys113 sulfhydryl to sulfenic (Cys-SOH) or sulfinic acid (Cys-SO2H). Given the involvement of Pin1 in pathological processes related to neurodegenerative diseases and to cancer, these findings could have implications for the prevention or treatment of disease.

© 2015 Published by NRC Research Press. Pin1 is a peptidyl-prolyl isomerase (PPIase) that plays a central role in eukaryotic cell cycle regulation, making this protein an interesting target for cancer therapy. Pin1 exhibits high specificity for substrates where proline is preceded by phosphoserine or phosphothreonine. The protein comprises an N-terminal WW (tryptophan-tryptophan) domain and a C-terminal PPIase domain. The cyclic peptide [CRYPEVEIC] (square brackets are used to denote the cyclic structure) represents a lead compound for a new class of nonphosphorylated Pin1 inhibitors. Unfortunately, it has not been possible thus far to characterize the Pin1-[CRYPEVEIC] complex by X-ray crystallography. Thus, the exact binding mode remains unknown. The current work employs hydrogen/deuterium exchange mass spectrometry for gaining insights into the Pin1-[CRYPEVEIC] interactions. The WW domain shows extensive conformational dynamics, both in the presence and in the absence of ligand. In contrast, profound changes in deuteration kinetics are observed in the PPIase domain after the addition of [CRYPEVEIC]. The secondary structure elements 2,3, and 4 exhibit markedly reduced deuteration, consistent with their postulated involvement in ligand binding. Unexpectedly, [CRYPEVEIC] destabilizes the range of residues 61-86, a segment that comprises basic side chains that normally interact with the substrate phosphate. This destabilization is likely caused by steric clashes with Y3 or E5 of the inhibitor. Ligand-induced destabilization has previously been reported for a few other proteins, but effects of this type are not very common. Our findings suggest that future crystallization trials on Pin1 variants deleted for residues in the 61-86 range might provide a path towards high-resolution X-ray structures of Pin1 bound to cyclic peptide inhibitors.

Hydrogen/deuterium exchange (HDX) mass spectrometry (MS) has become a key technique for monitoring structural and dynamic aspects of proteins in solution. This approach relies on the fact that exposure of a protein to D(2)O induces rapid amide H → D exchange in disordered regions that lack stable hydrogen-bonding. Tightly folded elements are much more protected from HDX, resulting in slow isotope exchange that is mediated by the structural dynamics ("breathing motions") of the protein. MS-based peptide mapping is a well established technique for measuring the mass shifts of individual protein segments. This tutorial review briefly discusses basic fundamentals of HDX/MS, before highlighting a number of recent developments and applications. Gas phase fragmentation strategies represent a promising alternative to the traditional proteolysis-based approach, but experimentalists have to be aware of scrambling phenomena that can be encountered under certain conditions. Electron-based dissociation methods provide a solution to this problem. We also discuss recent advances that facilitate the applicability of HDX/MS to membrane proteins, and to the characterization of short-lived protein folding intermediates. It is hoped that this review will provide a starting point for novices, as well as a useful reference for practitioners, who require an overview of some recent trends in HDX/MS.

Exposure of the membrane protein bacteriorhodopsin (BR) to SDS induces partial breakdown of the native conformation. The exact structural properties of this SDS state remain a matter of debate, despite its widespread use in BR folding experiments. The current work employs hydroxyl radical (·OH) labeling in conjunction with mass spectrometry (MS)-based peptide mapping for probing the solvent accessibility of individual BR segments in the presence of SDS. Previous work revealed methionine sulfoxide formation to be the dominant oxidative pathway. Those data suggested extensive unfolding of helices A and D in SDS. Unfortunately, the lack of Met residues in helices C and F implies that no direct information on the behavior of the latter two elements could be obtained. Here, we address this problem by employing two variants with additional Met residues, L93M (helix C) and V179M (helix F). The oxidation behavior of the resulting 11 methionines can be grouped into three categories: (1) extensively labeled both in native BR and in SDS (loop residues M32, M68, and M163); (2) protected in the native state but not in SDS (M20, M118); (3) always protected (M56, M60, M93, M145, M179, M209). These data show that a solvent-inaccessible core is retained in SDS. This core consists of partially intact helices B, C, E, F, and G. The termini of these helices are highly dynamic and/or unraveled, particularly on the cytoplasmic side. Overall, this work demonstrates how the use of engineered ·OH labeling sites can provide insights into structural properties of membrane proteins.

This work combines pulsed hydrogen/deuterium exchange (HDX) and top-down mass spectrometry for the structural characterization of short-lived protein folding intermediates. A custom-built flow device with three sequential mixing steps is used for (i) triggering protein folding, (ii) pulsed D(2)O labeling, and (iii) acid quenching. The earliest folding time point that can be studied with this system is 10 ms. The mixing device was coupled online to the electrospray source of a Fourier transform mass spectrometer, where intact protein ions are fragmented by electron capture dissociation (ECD). The viability of this experimental strategy is demonstrated by applying it to the refolding of horse apo-myoglobin (aMb), a reaction known to involve a transient intermediate. Cooling of the mixing device to 0 °C reduces the reaction rate such that the folding process occurs within the experimentally accessible time window. Top-down ECD provides an average spatial resolution of ca. 2 residues, surpassing the resolution typically achieved in traditional proteolytic digestion/HDX studies. Amide back exchange is virtually eliminated by the short (∼1 s) duration of the acid quenching step. The aMb folding intermediate exhibits HDX protection in helices G and H, whereas the remainder of the protein is largely unfolded. Marginal protection is seen for helix A. Overall, the top-down ECD approach used here offers insights into the sequence of events leading from the unfolded state to the native conformation, with envisioned future applications in the areas of protein misfolding and aggregation. The time-resolved experiments reported herein represent an extension of our previous work, where HDX/MS with top-down ECD was employed for monitoring "static" protein structures under equilibrium conditions (Pan et al. J. Am. Chem. Soc. 2009, 131, 12801).

This review discusses various mass spectrometry (MS)-based approaches for exploring structural aspects of proteins in solution. Electrospray ionization (ESI)-MS, in particular, has found fascinating applications in this area. For example, when used in conjunction with solution-phase hydrogen/deuterium exchange (HDX), ESI-MS is a highly sensitive tool for probing conformational dynamics. The main focus of this article is a technique that is complementary to HDX, that is, the covalent labeling of proteins by hydroxyl radicals. The reactivity of individual amino acid side chains with *OH is strongly affected by their degree of solvent exposure. Thus, analysis of the oxidative labeling pattern by peptide mapping and tandem mass spectrometry provides detailed structural information. A convenient method for *OH production is the photolysis of H(2)O(2) by a pulsed UV laser, resulting in oxidative labeling on the microsecond time scale. Selected examples demonstrate the use of this technique for structural studies on membrane proteins, and the combination with rapid mixing devices for characterizing the properties of short-lived protein (un)folding intermediates.

Water nanodroplets charged with excess protons play a central role during electrospray ionization (ESI). In the current study molecular dynamics (MD) simulations were used for gaining insights into the nanodroplet behavior based on classical mechanics. The SPC/E water model was modified to permit the inclusion of protons as highly mobile point charges at minimum computational cost. A spherical trapping potential was assigned to every SPC/E oxygen, thereby allowing the formation of protonated water molecules. Within a tightly packed nanodroplet the individual potential wells merge to form a three-dimensional energy landscape that facilitates rapid proton hopping between water molecules. This approach requires short-range modifications to the standard Coulomb potential for modeling electrostatic proton-water interactions. Simulations on nanodroplets consisting of 1248 water molecules and 10 protons (radius, ca. 21 A) result in a proton diffusion coefficient that is in agreement with the value measured in bulk solution. Radial proton distributions extracted from 1 ns MD runs exhibit a large peak around 14 A, in addition to substantial population density closer to the droplet center. Similar radial distributions were found for nanodroplets charged with Na+ ions. This behavior is dramatically different from that expected on the basis of continuum electrostatic theory, which predicts that excess charge should be confined to a thin layer on the droplet surface. One important contributor to this effect seems to be the ordering of water molecules at the liquid/vacuum interface. This ordering results in an electrical double layer, generating a potential gradient that tends to pull positive charge carriers (such as protons, but also others such as Na+ ions) toward the droplet interior. This deviation from the widely assumed surface charge paradigm could have implications for the mechanism by which protonated analyte ions are formed during ESI.

The S100 proteins comprise a group of EF-hand proteins that undergo a calcium-induced conformational change allowing them to interact with other proteins and produce a biological response. A unique feature of these proteins is the fact that they can form both homo- and heterodimers independent of calcium binding. The reported dissociation constants for several S100 proteins span a very large range, from 1-4 microM to

The structural characterization of short-lived intermediates provides insights into the mechanisms of protein folding and unfolding. Using holo-myoglobin as a model system, this work reports the application of oxidative pulse labeling for experiments of this kind. Protein unfolding is triggered by a pH jump from 6.5 to 3.2 in 150 mM NaCl. Subsequent (.-)OH exposure at various time points using laser photolysis of H2O2 leads to covalent modifications of solvent-exposed side chains within approximately 1 mus (Hambly, D. M.; Gross, M. L. J. Am. Soc. Mass Spectrom. 2005, 16, 2057-2063). Most of these modifications appear as 16 Da adducts in the mass spectrum of the intact protein. The overall extent of labeling increases with time, reflecting the exposure of reactive side chains that had previously been buried. Unfolding and disruption of heme-protein interactions go to completion within approximately 10 s. Spatially resolved information is obtained by monitoring the signal intensities of unmodified tryptic peptides. After 50 ms, many regions have lost most of their protection, whereas structure is retained in the B, E, F, and G helices. The BEF core remains partially folded, even after 500 ms, at which point helix G is fully unprotected. The observation of an "early" (BEFG) and a "late" (BEF) intermediate is in accord with optical stopped-flow measurements. Formation of these transient species is attributed to the persistence of heme-protein interactions during the early stages of the reaction. Overall, this work demonstrates the feasibility of laser-induced oxidative labeling as a tool for characterizing the structure of short-lived protein conformers. The combination of this approach with ultrarapid mixing or photochemical triggering should allow folding experiments in the submillisecond range.