Biological insights from hydrogen exchange mass spectrometry

Over the past two decades, hydrogen exchange mass spectrometry (HXMS) has achieved the status of a widespread and routine approach in the structural biology toolbox. The ability of hydrogen exchange to detect a range of protein dynamics coupled with the accessibility of mass spectrometry to mixtures and large complexes at low concentrations result in an unmatched tool for investigating proteins challenging to many other structural techniques. Recent advances in methodology and data analysis are helping HXMS deliver on its potential to uncover the connection between conformation, dynamics and the biological function of proteins and complexes. This review provides a brief overview of the HXMS method and focuses on four recent reports to highlight applications that monitor structure and dynamics of proteins and complexes, track protein folding, and map the thermodynamics and kinetics of protein unfolding at equilibrium. These case studies illustrate typical data, analysis and results for each application and demonstrate a range of biological systems for which the interpretation of HXMS in terms of structure and conformational parameters provides unique insights into function. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.     

Determination of macromolecular folding and structure by synchrotron x-ray radiolysis techniques

Radiolysis of water by synchrotron X-rays generates oxygen-containing radicals that undergo reactions with solvent accessible sites of macromolecules inducing stable covalent modifications or cleavage on millisecond time scales. The extent and site of these reactions are determined by gel electrophoresis and mass spectrometry analysis. These data are used to construct a high-resolution map of solvent accessibility at individual reactive sites. The experiments can be performed in a time-resolved manner to provide kinetic rate constants for dynamic events occurring at individual sites within macromolecules or can provide equilibrium parameters of binding and thermodynamics of folding processes. The application of this synchrotron radiolysis technique to the study of lysozyme protein structure and the equilibrium urea induced unfolding of apomyoglobin are described. The Mg2+-induced folding of Tetrahymena thermophila group I ribozyme shows the capability of the method to study kinetics of folding.     

Protein Dynamics Governed by Interfaces of High Polarity and Low Packing Density

The folding pathway, three-dimensional structure and intrinsic dynamics of proteins are governed by their amino acid sequences. Internal protein surfaces with physicochemical properties appropriate to modulate conformational fluctuations could play important roles in folding and dynamics. We show here that proteins contain buried interfaces of high polarity and low packing density, coined as LIPs: Light Interfaces of high Polarity, whose physicochemical properties make them unstable. The structures of well-characterized equilibrium and kinetic folding intermediates indicate that the LIPs of the corresponding native proteins fold late and are involved in local unfolding events. Importantly, LIPs can be identified using very fast and uncomplicated computational analysis of protein three-dimensional structures, which provides an easy way to delineate the protein segments involved in dynamics. Since LIPs can be retained while the sequences of the interacting segments diverge significantly, proteins could in principle evolve new functional features reusing pre-existing encoded dynamics. Large-scale identification of LIPS may contribute to understanding evolutionary constraints of proteins and the way protein intrinsic dynamics are encoded.     

HDX-analyzer: a novel package for statistical analysis of protein structure dynamics

BACKGROUND: HDX mass spectrometry is a powerful platform to probe protein structure dynamics during ligand binding, protein folding, enzyme catalysis, and such. HDX mass spectrometry analysis derives the protein structure dynamics based on the mass increase of a protein of which the backbone protons exchanged with solvent deuterium. Coupled with enzyme digestion and MS/MS analysis, HDX mass spectrometry can be used to study the regional dynamics of protein based on the m/z value or percentage of deuterium incorporation for the digested peptides in the HDX experiments. Various software packages have been developed to analyze HDX mass spectrometry data. Despite the progresses, proper and explicit statistical treatment is still lacking in most of the current HDX mass spectrometry software. In order to address this issue, we have developed the HDXanalyzer for the statistical analysis of HDX mass spectrometry data using R, Python, and RPY2. IMPLEMENTATION AND RESULTS: HDXanalyzer package contains three major modules, the data processing module, the statistical analysis module, and the user interface. RPY2 is employed to enable the connection of these three components, where the data processing module is implemented using Python and the statistical analysis module is implemented with R. RPY2 creates a low-level interface for R and allows the effective integration of statistical module for data processing. The data processing module generates the centroid for the peptides in form of m/z value, and the differences of centroids between the peptides derived from apo and ligand-bound protein allow us to evaluate whether the regions have significant changes in structure dynamics or not. Another option of the software is to calculate the deuterium incorporation rate for the comparison. The two types of statistical analyses are Paired Student's t-test and the linear combination of the intercept for multiple regression and ANCOVA model. The user interface is implemented with wxpython to facilitate the data visualization in graphs and the statistical analysis output presentation. In order to evaluate the software, a previously published xylanase HDX mass spectrometry analysis dataset is processed and presented. The results from the different statistical analysis methods are compared and shown to be similar. The statistical analysis results are overlaid with the three dimensional structure of the protein to highlight the regional structure dynamics changes in the xylanase enzyme. CONCLUSION: Statistical analysis provides crucial evaluation of whether a protein region is significantly protected or unprotected during the HDX mass spectrometry studies. Although there are several other available software programs to process HDX experimental data, HDXanalyzer is the first software program to offer multiple statistical methods to evaluate the changes in protein structure dynamics based on HDX mass spectrometry analysis. Moreover, the statistical analysis can be carried out for both m/z value and deuterium incorporation rate. In addition, the software package can be used for the data generated from a wide range of mass spectrometry instruments.     

Multiscale simulations of protein folding: application to formation of secondary structures

A multiscale simulation method of protein folding is proposed, using atomic representation of protein and solvent, combing genetic algorithms to determine the key protein structures from a global view, with molecular dynamic simulations to reveal the local folding pathways, thus providing an integrated landscape of protein folding. The method is found to be superior to previously investigated global search algorithms or dynamic simulations alone. For secondary structure formation of a selected peptide, RN24, the structures and dynamics produced by this method agree well with corresponding experimental results. Three most populated conformations are observed, including hairpin, β-sheet and α-helix. The energetic barriers separating these three structures are comparable to the kinetic energy of the atoms of the peptide, implying that the transition between these states can be easily triggered by kinetic perturbations, mainly through electrostatic interactions between charged atoms. Transitions between α-helix and β-sheet should jump over at least two energy barriers and may stay in the energetic trap of hairpin. It is proposed that the structure of proteins should be jointly governed by thermodynamic and dynamic factors; free energy is not the exclusive dominant for stability of proteins.     

Does native state topology determine the RNA folding mechanism?

Recent studies in protein folding suggest that native state topology plays a dominant role in determining the folding mechanism, yet an analogous statement has not been made for RNA, most likely due to the strong coupling between the ionic environment and conformational energetics that make RNA folding more complex than protein folding. Applying a distributed computing architecture to sample nearly 5000 complete tRNA folding events using a minimalist, atomistic model, we have characterized the role of native topology in tRNA folding dynamics: the simulated bulk folding behavior predicts well the experimentally observed folding mechanism. In contrast, single-molecule folding events display multiple discrete folding transitions and compose a largely diverse, heterogeneous dynamic ensemble. This both supports an emerging view of heterogeneous folding dynamics at the microscopic level and highlights the need for single-molecule experiments and both single-molecule and bulk simulations in interpreting bulk experimental measurements.     

Pressure as a denaturing agent in studies of single-point mutants of an amyloidogenic protein human cystatin c

Recently, we presented a convenient method combining a deuterium‐hydrogen exchange and electrospray mass spectrometry for studying high‐pressure denaturation of proteins (Stefanowicz et al., Biosci Rep 2009; 30:91–99). Here, we present results of pressure‐induced denaturation studies of an amyloidogenic protein—the wild‐type human cystatin C (hCC) and its single‐point mutants, in which Val57 residue from the hinge region was substituted by Asn, Asp or Pro, respectively. The place of mutation and the substituting residues were chosen mainly on a basis of theoretical calculations. Observation of H/D isotopic exchange proceeding during pressure induced unfolding and subsequent refolding allowed us to detect differences in the proteins stability and folding dynamics. On the basis of the obtained results we can conclude that proline residue at the hinge region makes cystatin C structure more flexible and dynamic, what probably facilitates the dimerization process of this hCC variant. Polar asparagine does not influence stability of hCC conformation significantly, whereas charged aspartic acid in 57 position makes the protein structure slightly more prone to unfolding. Our experiments also point out pressure denaturation as a valuable supplementary method in denaturation studies of mutated proteins. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Prediction of optimal folding routes of proteins that satisfy the principle of lowest entropy loss: dynamic contact maps and optimal control

An optimization model is introduced in which proteins try to evade high energy regions of the folding landscape, and prefer low entropy loss routes during folding. We make use of the framework of optimal control whose convenient solution provides practical and useful insight into the sequence of events during folding. We assume that the native state is available. As the protein folds, it makes different set of contacts at different folding steps. The dynamic contact map is constructed from these contacts. The topology of the dynamic contact map changes during the course of folding and this information is utilized in the dynamic optimization model. The solution is obtained using the optimal control theory. We show that the optimal solution can be cast into the form of a Gaussian Network that governs the optimal folding dynamics. Simulation results on three examples (CI2, Sso7d and Villin) show that folding starts by the formation of local clusters. Non-local clusters generally require the formation of several local clusters. Non-local clusters form cooperatively and not sequentially. We also observe that the optimal controller prefers "zipping" or small loop closure steps during folding. The folding routes predicted by the proposed method bear strong resemblance to the results in the literature.     

Conformational transitions in the membrane scaffold protein of phospholipid bilayer nanodiscs

Phospholipid bilayer nanodiscs are model membrane systems that provide an environment where membrane proteins are highly stable and monodisperse without the use of detergents or liposomes. Nanodiscs consist of a discoidal phospholipid bilayer encircled by two copies of an amphipathic alpha helical membrane scaffold protein, which is modeled from apolipoprotein A-1. Hydrogen exchange mass spectrometry was used to probe the structure and dynamics of the scaffold protein in the presence and absence of lipid. On nanodisc self-assembly, the entire scaffold protein gained significant protection from exchange, consistent with a large, protein-wide, structural rearrangement. This protection was short-lived and the scaffold protein was highly deuterated within 2 h. Several regions of the scaffold protein, in both the lipid-free and lipid-associated states, displayed EX1 unfolding kinetics. The rapid deuteration of the scaffold protein and the presence of correlated unfolding events both indicate that nanodiscs are dynamic rather than rigid bodies in solution. This work provides a catalog of the expected scaffold protein peptic peptides in a nanodisc-hydrogen exchange mass spectrometry experiment and their deuterium uptake signatures, data that can be used as a benchmark to verify correct assembly and nanodisc structure. Such reference data will be useful control data for all hydrogen exchange mass spectrometry experiments involving nanodiscs in which transmembrane or lipid-associated proteins are the primary molecule(s) of interest.     

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.     

Structural proteomics of macromolecular assemblies using oxidative footprinting and mass spectrometry

Understanding the composition, structure and dynamics of macromolecules and their assemblies is at the forefront of biological science today. Hydroxyl-radical-mediated protein footprinting using mass spectrometry can define macromolecular structure, macromolecular assembly and conformational changes of macromolecules in solution based on measurements of reactivity of amino acid side-chain groups with covalent-modification reagents. Subsequent to oxidation by reactive oxygen species, proteins are digested by specific proteases to generate peptides for analysis by mass spectrometry. Accurate measurements of side-chain reactivity are achieved using quantitative liquid-chromatography-coupled mass spectrometry, whereas the side-chain sites within the macromolecular probes are identified using tandem mass spectrometry. In addition, the use of footprinting data in conjunction with computational modeling approaches is a powerful new method for testing and refining structural models of macromolecules and their complexes.     

New insights into the mechanisms of protein misfolding and aggregation in amyloidogenic diseases derived from pressure studies

Hydrostatic pressure is a robust tool for studying the thermodynamics of protein folding and protein interactions, as well as the dynamics and structure of folding intermediates. One of the main innovations obtained from using high pressure is the stabilization of folding intermediates such as molten-globule conformations, thus providing a unique opportunity for characterizing their structure and dynamics. Equally important is the prospect of understanding protein misfolding diseases by using pressure to populate partially folded intermediates at the junction between productive and off-pathway folding, which may give rise to misfolded proteins, aggregates, and amyloids. High hydrostatic pressure (HHP) has also been used to dissociate nonamyloid aggregates and inclusion bodies. In many proteins, the competition between correct folding and misfolding can lead to formation of insoluble aggregates, an important problem for the biotechnology industry and for human pathologies such as amyloidosis, Alzheimer's, Parkinson's, prion's, and tumor diseases. The diversity of diseases that result from protein misfolding has made this theme an important research focus for pharmaceutical and biotechnology companies. The use of high-pressure promises to contribute to the identification of the mechanisms behind these defects and creation of therapies against these diseases.     

Snapshots of protein folding problem: implications of folding and misfolding studies

Deciphering the code that determines the three-dimensional structure of proteins and the ability to predict the final folded form of a protein is still elusive to molecular biophysists. In the case of several proteins a similar tertiary structure is not accompanied by any significant sequence similarity. The question now remains whether a code beyond the genetic code that describes the arrangement of the amino acid within a three dimensional protein structure. The available data undoubtedly demonstrates that the redundancy of this code must be tremendous. Several techniques such as nuclear magnetic resonance spectroscopy and laser detection techniques, coupled with fast initiation of the folding reaction, can now probe the folding events in milliseconds or even faster and provide highly relevant information. The thermodynamic analysis of the folding process and of kinetic intermediates opens whole new avenue of understanding. Breaking the protein folding code would enable scientists to look at a gene whose function is unknown and predict the three-dimensional structure of the protein it encodes. This would give them a very good idea of what the gene does. In this review we hope to bring together the information available about protein folding with particular emphasis on folding intermediate(s). Additionally, the practical consequences of the solution of the protein folding problem in medicine and biotechnology are also discussed.     

Protein folding mechanisms studied by pulsed oxidative labeling and mass spectrometry

Deciphering the mechanisms of protein folding remains a considerable challenge. In this review we discuss the application of pulsed oxidative labeling for tracking protein structural changes in a time-resolved fashion. Exposure to a microsecond OH pulse at selected time points during folding induces the oxidation of solvent-accessible side chains, whereas buried residues are protected. Oxidative modifications can be detected by mass spectrometry. Folding is associated with dramatic accessibility changes, and therefore this method can provide detailed mechanistic insights. Solvent accessibility patterns are complementary to H/D exchange investigations, which report on the extent of hydrogen bonding. This review highlights the application of pulsed OH labeling to soluble proteins as well as membrane proteins.     

Autonomous subdomains in protein folding.

Proteolytic dissection of native trp repressor and horse heart cytochrome c has been used to infer some of the steps in the folding pathways of the intact proteins. For both proteins, small fragments are capable of undergoing spontaneous noncovalent association to form subdomains with native-like secondary and/or tertiary structural features, suggesting that dissection/reassembly may be a general method to gain insight into the structures of folding intermediates. The importance of this approach is its simplicity and potential applicability to studying the folding pathways of a wide range of proteins. The proteases report on the structure and dynamics of the native state, circumventing the need for prior knowledge of the structures of folding intermediates. The observation that small fragments of proteins can associated noncovalently suggests that protein folding can be viewed as an intramolecular "recognition" process. The results imply that substantial information about protein structure and folding is encoded at the level of subdomains, and that chain connectivity has only a minor role in determining the fold.     

Mapping the protein-protein interface between a toxin and its cognate antitoxin from the bacterial pathogen Streptococcus pyogenes

Protein--protein interactions are ubiquitous and essential for most biological processes. Although new proteomic technologies have generated large catalogs of interacting proteins, considerably less is known about these interactions at the molecular level, information that would aid in predicting protein interactions, designing therapeutics to alter these interactions, and understanding the effects of disease-producing mutations. Here we describe mapping the interacting surfaces of the bacterial toxin SPN (Streptococcus pyogenes NAD(+) hydrolase) in complex with its antitoxin IFS (immunity factor for SPN) by using hydrogen-deuterium amide exchange and electrospray ionization mass spectrometry. This approach affords data in a relatively short time for small amounts of protein, typically 5-7 pmol per analysis. The results show a good correspondence with a recently determined crystal structure of the IFS--SPN complex but additionally provide strong evidence for a folding transition of the IFS protein that accompanies its binding to SPN. The outcome shows that mass-based chemical footprinting of protein interaction surfaces can provide information about protein dynamics that is not easily obtained by other methods and can potentially be applied to large, multiprotein complexes that are out of range for most solution-based methods of biophysical analysis.     

Study of early events in the protein folding of villin headpiece using molecular dynamics simulation

Protein folding is scientifically and computationally challenging problem. The early phases of protein folding are interesting due to various events like nascent secondary structure formation, hydrophobic collapse leading to formation of non-native or meta-stable conformations. These events occur within a very short time span of 100ns as compared to total folding time of few microseconds. It is highly difficult to observe these events experimentally due to very short lifetime. Molecular dynamics simulation technique can efficiently probe the detailed atomic level understanding about these events. In the present paper, all atom molecular dynamics simulation trajectory of nearly 200ns was carried out for fully solvated villin headpiece with PME treatment using AMBER 7 package. Initial hydrophobic collapse along with secondary structure formation resulted into formation of partially stable non-native conformations. The formation of secondary structural elements and hydrophobic collapse takes place simultaneously in the folding process.     

Conservation of folding pathways in evolutionarily distant globin sequences

To test the hypothesis that the folding pathways of evolutionarily related proteins with similar three-dimensional structures but widely different sequences should be similar, the folding pathway of apoleghemoglobin has been characterized using stopped-flow circular dichroism, heteronuclear NMR pulse labeling techniques and mass spectrometry. The pathway of folding was found to differ significantly from that of a protein of the same family, apomyoglobin, although both proteins appear to fold through helical burst phase intermediates. For leghemoglobin, the burst phase intermediate exhibits stable helical structure in the G and H helices, together with a small region in the center of the E helix. The A and B helices are not stabilized until later stages of the folding process. The structure of the burst phase folding intermediate thus differs from that of apomyoglobin, in which stable helical structure is formed in the A, B, G and H helix regions.     

Identification and characterization of key substructures involved in the early folding events of a (beta/alpha)8-barrel protein as studied by experimental and computational methods

A number of studies have examined the structural properties of late folding intermediates of (beta/alpha)8-barrel proteins involved in tryptophan biosynthesis, whereas there is little information available about the early folding events of these proteins. To identify the contiguous polypeptide segments important to the folding of the (beta/alpha)8-barrel protein Escherichia coli N-(5'-phosphoribosyl)anthranilate isomerase, we structurally characterized fragments and circularly permuted forms of the protein. We also simulated thermal unfolding of the protein using molecular dynamics. Our fragmentation experiments demonstrate that the isolated (beta/alpha)(1-4)beta5 fragment is almost as stable as the full-length protein. The far and near-UV CD spectra of this fragment are indicative of native-like secondary and tertiary structures. Structural analysis of the circularly permutated proteins shows that if the protein is cleaved within the two N-terminal betaalpha modules, the amount of secondary structure is unaffected, whereas, when cleaved within the central (beta/alpha)(3-4)beta5 segment, the protein simply cannot fold. An ensemble of the denatured structures produced by thermal unfolding simulations contains a persistent local structure comprised of beta3, beta4 and beta5. The presence of this three-stranded beta-barrel suggests that it may be an important early-stage folding intermediate. Interactions found in (beta/alpha)(3-4)beta5 may be essential for the early events of ePRAI folding if they provide a nucleation site that directs folding.