The hydrogen exchange core and protein folding

Backbone 1H, 13C and 15N resonance assignment of the ubiquitin specific protease 7 catalytic domain (residues 208-554) in complex with a small molecule ligand

The backbone 1H, 13C and 15N resonance assignment of Ubiquitin Specific Protease 7 catalytic domain (residues 208-554) was performed in its complex with a small molecule ligand and in its apo form as a reference. The amide 1H-15N signal intensities were boosted by an amide hydrogen exchange protocol, where expressed 2H, 13C, 15N-labeled protein was unfolded and re-folded to ensure exchange of amide deuterons to protons. The resonance assignments were used to determine chemical shift perturbations on ligand binding, which are consistent with the binding site observed by crystallography.

Type 2B von Willebrand disease mutations differentially perturb autoinhibition of the A1 domain

Type 2B von Willebrand disease (VWD) is an inherited bleeding disorder in which a subset of point mutations in the von Willebrand factor (VWF) A1 domain and recently identified autoinhibitory module (AIM) cause spontaneous binding to glycoprotein Ibα (GPIbα) on the platelet surface. All reported type 2B VWD mutations share this enhanced binding; however, type 2B VWD manifests as variable bleeding complications and platelet levels in patients, depending on the underlying mutation. Understanding how these mutations localizing to a similar region can result in such disparate patient outcomes is essential for detailing our understanding of VWF regulatory and activation mechanisms. In this study, we produced recombinant glycosylated AIM-A1 fragments bearing type 2B VWD mutations and examined how each mutation affects the A1 domain's thermodynamic stability, conformational dynamics, and biomechanical regulation of the AIM. We found that the A1 domain with mutations associated with severe bleeding occupy a higher affinity state correlating with enhanced flexibility in the secondary GPIbα-binding sites. Conversely, mutation P1266L, associated with normal platelet levels, has similar proportions of high-affinity molecules to wild-type (WT) but shares regions of solvent accessibility with both WT and other type 2B VWD mutations. V1316M exhibited exceptional instability and solvent exposure compared with all variants. Lastly, examination of the mechanical stability of each variant revealed variable AIM unfolding. Together, these studies illustrate that the heterogeneity among type 2B VWD mutations is evident in AIM-A1 fragments.

Rational Design of Protein-Specific Folding Modifiers

Protein-folding can go wrong in vivo and in vitro, with significant consequences for the living organism and the pharmaceutical industry, respectively. Here we propose a design principle for small-peptide-based protein-specific folding modifiers. The principle is based on constructing a "xenonucleus", which is a prefolded peptide that mimics the folding nucleus of a protein. Using stopped-flow kinetics, NMR spectroscopy, Förster resonance energy transfer, single-molecule force measurements, and molecular dynamics simulations, we demonstrate that a xenonucleus can make the refolding of ubiquitin faster by 33 ± 5%, while variants of the same peptide have little or no effect. Our approach provides a novel method for constructing specific, genetically encodable folding catalysts for suitable proteins that have a well-defined contiguous folding nucleus.

Interpreting a black box predictor to gain insights into early folding mechanisms

Protein folding and function are closely connected, but the exact mechanisms by which proteins fold remain elusive. Early folding residues (EFRs) are amino acids within a particular protein that induce the very first stages of the folding process. High-resolution EFR data are only available for few proteins, which has previously enabled the training of a protein sequence-based machine learning 'black box' predictor (EFoldMine). Such a black box approach does not allow a direct extraction of the 'early folding rules' embedded in the protein sequence, whilst such interpretation is essential to improve our understanding of how the folding process works. We here apply and investigate a novel 'grey box' approach to the prediction of EFRs from protein sequence to gain mechanistic residue-level insights into the sequence determinants of EFRs in proteins. We interpret the rule set for three datasets, a default set comprised of natural proteins, a scrambled set comprised of the scrambled default set sequences, and a set of de novo designed proteins. Finally, we relate these data to the secondary structure adopted in the folded protein and provide all information online via http://xefoldmine.bio2byte.be/, as a resource to help understand and steer early protein folding.

Functional and Regulatory Roles of Fold-Switching Proteins

Fold-switching proteins respond to cellular stimuli by remodeling their secondary structures and changing their functions. Whereas several previous reviews have focused on various structural, physical-chemical, and evolutionary aspects of this newly emerging class of proteins, this minireview focuses on how fold switching modulates protein function and regulates biological processes. It first compares and contrasts fold switchers with other known types of proteins. Second, it presents examples of how various proteins can change their functions through fold switching. Third, it demonstrates that fold switchers can regulate biological processes by discussing two proteins, RfaH and KaiB, whose dramatic secondary structure remodeling events directly affect gene expression and a circadian clock, respectively. Finally, this minireview discusses how the field of protein fold switching might advance.

Paramagnetic relaxation enhancement-assisted structural characterization of a partially disordered conformation of ubiquitin

Nuclear magnetic resonance (NMR) is a powerful tool to study three-dimensional structures as well as protein conformational fluctuations in solution, but it is compromised by increases in peak widths and missing signals. We previously reported that ubiquitin has two folded conformations, N₁ and N₂ and plus another folded conformation, I, in which some amide group signals of residues 33-41 almost disappeared above 3 kbar at pH 4.5 and 273 K. Thus, well-converged structural models could not be obtained for this region owing to the absence of distance restraints. Here, we reexamine the problem using the ubiquitin Q41N variant as a model for this locally disordered conformation, I. We demonstrate that the variant shows pressure-induced loss of backbone amide group signals at residues 28, 33, 36, and 39-41 like the wild-type, with a similar but smaller effect on CαH and CβH signals. In order to characterize this I structure, we measured paramagnetic relaxation enhancement (PRE) under high pressure to obtain distance restraints, and calculated the structure assisted by Bayesian inference. We conclude that the more disordered I conformation observed at pH 4.0, 278 K, and 2.5 kbar largely retained the N₂ conformation, although the amide groups at residues 33-41 have more heterogeneous conformations and more contact with water, which differ from the N₁ and N₂ states. The PRE-assisted strategy has the potential to improve structural characterization of proteins that lack NMR signals, especially for relatively more open and hydrated protein conformations.

StructureDistiller: Structural relevance scoring identifies the most informative entries of a contact map

Protein folding and structure prediction are two sides of the same coin. Contact maps and the related techniques of constraint-based structure reconstruction can be considered as unifying aspects of both processes. We present the Structural Relevance (SR) score which quantifies the information content of individual contacts and residues in the context of the whole native structure. The physical process of protein folding is commonly characterized with spatial and temporal resolution: some residues are Early Folding while others are Highly Stable with respect to unfolding events. We employ the proposed SR score to demonstrate that folding initiation and structure stabilization are subprocesses realized by distinct sets of residues. The example of cytochrome c is used to demonstrate how StructureDistiller identifies the most important contacts needed for correct protein folding. This shows that entries of a contact map are not equally relevant for structural integrity. The proposed StructureDistiller algorithm identifies contacts with the highest information content; these entries convey unique constraints not captured by other contacts. Identification of the most informative contacts effectively doubles resilience toward contacts which are not observed in the native contact map. Furthermore, this knowledge increases reconstruction fidelity on sparse contact maps significantly by 0.4 Å.

How internal cavities destabilize a protein

Although many proteins possess a distinct folded structure lying at a minimum in a funneled free energy landscape, thermal energy causes any protein to continuously access lowly populated excited states. The existence of excited states is an integral part of biological function. Although transitions into the excited states may lead to protein misfolding and aggregation, little structural information is currently available for them. Here, we show how NMR spectroscopy, coupled with pressure perturbation, brings these elusive species to light. As pressure acts to favor states with lower partial molar volume, NMR follows the ensuing change in the equilibrium spectroscopically, with residue-specific resolution. For T4 lysozyme L99A, relaxation dispersion NMR was used to follow the increase in population of a previously identified "invisible" folded state with pressure, as this is driven by the reduction in cavity volume by the flipping-in of a surface aromatic group. Furthermore, multiple partly disordered excited states were detected at equilibrium using pressure-dependent H/D exchange NMR spectroscopy. Here, unfolding reduced partial molar volume by the removal of empty internal cavities and packing imperfections through subglobal and global unfolding. A close correspondence was found for the distinct pressure sensitivities of various parts of the protein and the amount of internal cavity volume that was lost in each unfolding event. The free energies and populations of excited states allowed us to determine the energetic penalty of empty internal protein cavities to be 36 cal⋅Å^-3.

Application of an interpretable classification model on Early Folding Residues during protein folding

Background

Machine learning strategies are prominent tools for data analysis. Especially in life sciences, they have become increasingly important to handle the growing datasets collected by the scientific community. Meanwhile, algorithms improve in performance, but also gain complexity, and tend to neglect interpretability and comprehensiveness of the resulting models.

Results

Generalized Matrix Learning Vector Quantization (GMLVQ) is a supervised, prototype-based machine learning method and provides comprehensive visualization capabilities not present in other classifiers which allow for a fine-grained interpretation of the data. In contrast to commonly used machine learning strategies, GMLVQ is well-suited for imbalanced classification problems which are frequent in life sciences. We present a Weka plug-in implementing GMLVQ. The feasibility of GMLVQ is demonstrated on a dataset of Early Folding Residues (EFR) that have been shown to initiate and guide the protein folding process. Using 27 features, an area under the receiver operating characteristic of 76.6% was achieved which is comparable to other state-of-the-art classifiers. The obtained model is accessible at https://biosciences.hs-mittweida.de/efpred/.

Conclusions

The application on EFR prediction demonstrates how an easy interpretation of classification models can promote the comprehension of biological mechanisms. The results shed light on the special features of EFR which were reported as most influential for the classification: EFR are embedded in ordered secondary structure elements and they participate in networks of hydrophobic residues. Visualization capabilities of GMLVQ are presented as we demonstrate how to interpret the results.

Electronic supplementary material

The online version of this article (10.1186/s13040-018-0188-2) contains supplementary material, which is available to authorized users.

The Singular NMR Fingerprint of a Polyproline II Helical Bundle

Polyproline II (PPII) helices play vital roles in biochemical recognition events and structures like collagen and form part of the conformational landscapes of intrinsically disordered proteins (IDPs). Nevertheless, this structure is generally hard to detect and quantify. Here, we report the first thorough NMR characterization of a PPII helical bundle protein, the Hypogastrura harveyi "snow flea" antifreeze protein (sfAFP). J-couplings and nuclear Overhauser enhancement spectroscopy confirm a natively folded structure consisting of six PPII helices. NMR spectral analyses reveal quite distinct Hα2 versus Hα3 chemical shifts for 28 Gly residues as well as ¹³Cα, ¹⁵N, and ¹HN conformational chemical shifts (Δδ) unique to PPII helical bundles. The ¹⁵N Δδ and ¹HN Δδ values and small negative ¹HN temperature coefficients evince hydrogen-bond formation. ¹H-¹⁵N relaxation measurements reveal that the backbone structure is generally highly rigid on ps-ns time scales. NMR relaxation parameters and biophysical characterization reveal that sfAFP is chiefly a dimer. For it, a structural model featuring the packing of long, flat hydrophobic faces at the dimer interface is advanced. The conformational stability, measured by amide H/D exchange to be 6.24 ± 0.2 kcal·mol^-1, is elevated. These are extraordinary findings considering the great entropic cost of fixing Gly residues and, together with the remarkable upfield chemical shifts of 28 Gly ¹Hα, evidence significant stabilizing contributions from CαHα ||| O═C hydrogen bonds. These stabilizing interactions are corroborated by density functional theory calculations and natural bonding orbital analysis. The singular conformational chemical shifts, J-couplings, high hNOE ratios, small negative temperature coefficients, and slowed H/D exchange constitute a unique set of fingerprints to identify PPII helical bundles, which may be formed by hundreds of Gly-rich motifs detected in sequence databases. These results should aid the quantification of PPII helices in IDPs, the development of improved antifreeze proteins, and the incorporation of PPII helices into novel designed proteins.

Characterizing the relation of functional and Early Folding Residues in protein structures using the example of aminoacyl-tRNA synthetases

Proteins are chains of amino acids which adopt a three-dimensional structure and are then able to catalyze chemical reactions or propagate signals in organisms. Without external influence, many proteins fold into their native structure, and a small number of Early Folding Residues (EFR) have previously been shown to initiate the formation of secondary structure elements and guide their respective assembly. Using the two diverse superfamilies of aminoacyl-tRNA synthetases (aaRS), it is shown that the position of EFR is preserved over the course of evolution even when the corresponding sequence conservation is small. Folding initiation sites are positioned in the center of secondary structure elements, independent of aaRS class. In class I, the predicted position of EFR resembles an ancient structural packing motif present in many seemingly unrelated proteins. Furthermore, it is shown that EFR and functionally relevant residues in aaRS are almost entirely disjoint sets of residues. The Start2Fold database is used to investigate whether this separation of EFR and functional residues can be observed for other proteins. EFR are found to constitute crucial connectors of protein regions which are distant at sequence level. Especially, these residues exhibit a high number of non-covalent residue-residue contacts such as hydrogen bonds and hydrophobic interactions. This tendency also manifests as energetically stable local regions, as substantiated by a knowledge-based potential. Despite profound differences regarding how EFR and functional residues are embedded in protein structures, a strict separation of structurally and functionally relevant residues cannot be observed for a more general collection of proteins.

Water-Protein Interactions Coupled with Protein Conformational Transition

Conformational fluctuations of proteins are crucially important for their functions. However, changes in the location and dynamics of hydrated water in many proteins accompanied by the conformational transition have not been fully understood. Here, we used phase-modulated clean chemical exchange NMR approach to investigate pressure-induced changes in water-to-amide proton exchange occurring at sub-second time scale. With the transition of ubiquitin from its native conformation (N1) to an alternative conformation (N2) at 250 MPa, proton exchange rates of residues 32-35, 40-41, and 71, which are located at the C-terminal side of the protein, were significantly increased. These observations can be explained by the destabilization of the hydrogen bonds in the backbone and partial exposure of those amide groups to solvent in N2. We conclude that phase-modulated clean chemical exchange NMR approach coupled with pressure perturbation will be a useful tool for investigations of more open and hydrated protein structures.

Extant fold-switching proteins are widespread

It is commonly thought that each globular protein has a single 3D structure, or fold, that fosters its function. In contrast, recent studies have identified several fold-switching proteins whose secondary structures can be remodeled in response to cellular stimuli. Although thought to be rare, we found 96 literature-validated fold-switching proteins by exhaustively searching the database of protein structures [Protein Data Bank (PDB)]. Characterizing these proteins led us to hypothesize that their abundance may be underrepresented in the PDB. Thus, we developed a computational method that identifies fold-switching proteins and used it to estimate that 0.5–4% of PDB proteins switch folds. These results suggest that proteins switch folds with significant frequency, which has implications for cell biology, genomics, and human health.

A central tenet of biology is that globular proteins have a unique 3D structure under physiological conditions. Recent work has challenged this notion by demonstrating that some proteins switch folds, a process that involves remodeling of secondary structure in response to a few mutations (evolved fold switchers) or cellular stimuli (extant fold switchers). To date, extant fold switchers have been viewed as rare byproducts of evolution, but their frequency has been neither quantified nor estimated. By systematically and exhaustively searching the Protein Data Bank (PDB), we found ∼100 extant fold-switching proteins. Furthermore, we gathered multiple lines of evidence suggesting that these proteins are widespread in nature. Based on these lines of evidence, we hypothesized that the frequency of extant fold-switching proteins may be underrepresented by the structures in the PDB. Thus, we sought to identify other putative extant fold switchers with only one solved conformation. To do this, we identified two characteristic features of our ∼100 extant fold-switching proteins, incorrect secondary structure predictions and likely independent folding cooperativity, and searched the PDB for other proteins with similar features. Reassuringly, this method identified dozens of other proteins in the literature with indication of a structural change but only one solved conformation in the PDB. Thus, we used it to estimate that 0.5–4% of PDB proteins switch folds. These results demonstrate that extant fold-switching proteins are likely more common than the PDB reflects, which has implications for cell biology, genomics, and human health.

A discontinuous autoinhibitory module masks the A1 domain of von Willebrand factor

Essentials The mechanism for the auto-inhibition of von Willebrand factor (VWF) remains unclear. Hydrogen exchange of two VWF A1 fragments with disparate activities was measured and compared. Discontinuous residues flanking A1 form a structural module that blocks A1 binding to the platelet. Our results suggest a potentially unified model of VWF activation. Click to hear an ISTH Academy presentation on the domain architecture of VWF and activation by elongational flow by Dr Springer SUMMARY: Background How von Willebrand factor (VWF) senses and responds to shear flow remains unclear. In the absence of shear flow, VWF or its fragments can be induced to bind spontaneously to platelet GPIbα. Objectives To elucidate the auto-inhibition mechanism of VWF. Methods Hydrogen-deuterium exchange (HDX) of two recombinant VWF fragments expressed from baby hamster kidney cells were measured and compared. Results The shortA1 protein contains VWF residues 1261-1472 and binds GPIbα with a significantly higher affinity than the longA1 protein that contains VWF residues 1238-1472. Both proteins contain the VWF A1 domain (residues 1272-1458). Many residues in longA1, particularly those in the N- and C-terminal sequences flanking the A1 domain, and in helix α1, loops α1β2 and β3α2, demonstrated markedly reduced HDX compared with their counterparts in shortA1. The HDX-protected region in longA1 overlaps with the GPIbα-binding interface and is clustered with type 2B von Willebrand disease (VWD) mutations. Additional comparison with the HDX of denatured longA1 and ristocetin-bound longA1 indicates the N- and C-terminal sequences flanking the A1 domain form cooperatively an integrated autoinhibitory module (AIM) that interacts with the HDX-protected region. Binding of ristocetin to the C-terminal part of the AIM desorbs the AIM from A1 and enables longA1 binding to GPIbα. Conclusion The discontinuous AIM binds the A1 domain and prevents it from binding to GPIbα, which has significant implications for the pathogenesis of type 2B VWD and the shear-induced activation of VWF activity.

Exploring the Sequence-based Prediction of Folding Initiation Sites in Proteins

Protein folding is a complex process that can lead to disease when it fails. Especially poorly understood are the very early stages of protein folding, which are likely defined by intrinsic local interactions between amino acids close to each other in the protein sequence. We here present EFoldMine, a method that predicts, from the primary amino acid sequence of a protein, which amino acids are likely involved in early folding events. The method is based on early folding data from hydrogen deuterium exchange (HDX) data from NMR pulsed labelling experiments, and uses backbone and sidechain dynamics as well as secondary structure propensities as features. The EFoldMine predictions give insights into the folding process, as illustrated by a qualitative comparison with independent experimental observations. Furthermore, on a quantitative proteome scale, the predicted early folding residues tend to become the residues that interact the most in the folded structure, and they are often residues that display evolutionary covariation. The connection of the EFoldMine predictions with both folding pathway data and the folded protein structure suggests that the initial statistical behavior of the protein chain with respect to local structure formation has a lasting effect on its subsequent states.

Early Folding Events, Local Interactions, and Conservation of Protein Backbone Rigidity

Protein folding is in its early stages largely determined by the protein sequence and complex local interactions between amino acids, resulting in lower energy conformations that provide the context for further folding into the native state. We compiled a comprehensive data set of early folding residues based on pulsed labeling hydrogen deuterium exchange experiments. These early folding residues have corresponding higher backbone rigidity as predicted by DynaMine from sequence, an effect also present when accounting for the secondary structures in the folded protein. We then show that the amino acids involved in early folding events are not more conserved than others, but rather, early folding fragments and the secondary structure elements they are part of show a clear trend toward conserving a rigid backbone. We therefore propose that backbone rigidity is a fundamental physical feature conserved by proteins that can provide important insights into their folding mechanisms and stability.

Fold and flexibility: what can proteins' mechanical properties tell us about their folding nucleus?

The determination of a protein's folding nucleus, i.e. a set of native contacts playing an important role during its folding process, remains an elusive yet essential problem in biochemistry. In this work, we investigate the mechanical properties of 70 protein structures belonging to 14 protein families presenting various folds using coarse-grain Brownian dynamics simulations. The resulting rigidity profiles combined with multiple sequence alignments show that a limited set of rigid residues, which we call the consensus nucleus, occupy conserved positions along the protein sequence. These residues' side chains form a tight interaction network within the protein's core, thus making our consensus nuclei potential folding nuclei. A review of experimental and theoretical literature shows that most (above 80%) of these residues were indeed identified as folding nucleus member in earlier studies.

Aberrant increase of NMR signal in hydrogen exchange experiments. Observation and explanation

Hydrogen exchange (HX) NMR spectroscopy is widely used for monitoring structure, stability and dynamics of proteins at the level of individual residues. The stochastic replacement of protons by deuterons typically leads to an exponential decrease of the NMR signals. However, an unusual signal increase was observed in HX of several amides for T4 lysozyme L99A. This effect can be attributed to peak sharpening as a result of reduced dipolar relaxation from proximal amide protons that experience more rapid hydrogen/deuterium (H/D) exchange. The behavior was specifically observed at the termini of secondary structure elements, where large differences in protection against H/D exchange are observed. This effect is expected to be more widespread in NMR HX studies, and is important for the accurate determination of protection factors.

Adaptive local learning in sampling based motion planning for protein folding

BACKGROUND

Simulating protein folding motions is an important problem in computational biology. Motion planning algorithms, such as Probabilistic Roadmap Methods, have been successful in modeling the folding landscape. Probabilistic Roadmap Methods and variants contain several phases (i.e., sampling, connection, and path extraction). Most of the time is spent in the connection phase and selecting which variant to employ is a difficult task. Global machine learning has been applied to the connection phase but is inefficient in situations with varying topology, such as those typical of folding landscapes.

RESULTS

We develop a local learning algorithm that exploits the past performance of methods within the neighborhood of the current connection attempts as a basis for learning. It is sensitive not only to different types of landscapes but also to differing regions in the landscape itself, removing the need to explicitly partition the landscape. We perform experiments on 23 proteins of varying secondary structure makeup with 52-114 residues. We compare the success rate when using our methods and other methods. We demonstrate a clear need for learning (i.e., only learning methods were able to validate against all available experimental data) and show that local learning is superior to global learning producing, in many cases, significantly higher quality results than the other methods.

CONCLUSIONS

We present an algorithm that uses local learning to select appropriate connection methods in the context of roadmap construction for protein folding. Our method removes the burden of deciding which method to use, leverages the strengths of the individual input methods, and it is extendable to include other future connection methods.

Autoinhibition of ETV6 DNA Binding Is Established by the Stability of Its Inhibitory Helix

The ETS transcriptional repressor ETV6 (or TEL) is autoinhibited by an α-helix that sterically blocks its DNA-binding ETS domain. The inhibitory helix is marginally stable and unfolds when ETV6 binds to either specific or non-specific DNA. Using NMR spectroscopy, we show that folding of the inhibitory helix requires a buried charge-dipole interaction with helix H1 of the ETS domain. This interaction also contributes directly to autoinhibition by precluding a highly conserved dipole-enhanced hydrogen bond between the phosphodiester backbone of bound DNA and the N terminus of helix H1. To probe further the thermodynamic basis of autoinhibition, ETV6 variants were generated with amino acid substitutions introduced along the solvent exposed surface of the inhibitory helix. These changes were designed to increase the intrinsic helical propensity of the inhibitory helix without perturbing its packing interactions with the ETS domain. NMR-monitored amide hydrogen exchange measurements confirmed that the stability of the folded inhibitory helix increases progressively with added helix-promoting substitutions. This also results in progressively reinforced autoinhibition and decreased DNA-binding affinity. Surprisingly, locking the inhibitory helix onto the ETS domain by a disulfide bridge severely impairs, but does not abolish DNA binding. Weak interactions still occur via an interface displaced from the canonical ETS domain DNA-binding surface. Collectively, these studies establish a direct thermodynamic linkage between inhibitory helix stability and ETV6 autoinhibition, and demonstrate that helix unfolding does not strictly precede DNA binding. Modulating inhibitory helix stability provides a potential route for the in vivo regulation of ETV6 activity.

Characterization of folding cores in the cyclophilin A-cyclosporin A complex

Determining the folding core of a protein yields information about its folding process and dynamics. The experimental procedures for identifying the amino acids that make up the folding core include hydrogen-deuterium exchange and Φ-value analysis and can be expensive and time consuming. Because of this, there is a desire to improve upon existing methods for determining protein folding cores theoretically. We have obtained HDX data for the complex of cyclophilin A with the immunosuppressant cyclosporin A. We compare these data, as well as literature values for uncomplexed cyclophilin A, to theoretical predictions using a combination of rigidity analysis and coarse-grained simulations of protein motion. We find that in this case, the most specific prediction of folding cores comes from a combined approach that models the rigidity of the protein using the first software suite and the dynamics of the protein using the froda tool.

Clusters of isoleucine, leucine, and valine side chains define cores of stability in high-energy states of globular proteins: Sequence determinants of structure and stability

Measurements of protection against exchange of main chain amide hydrogens (NH) with solvent hydrogens in globular proteins have provided remarkable insights into the structures of rare high-energy states that populate their folding free-energy surfaces. Lacking, however, has been a unifying theory that rationalizes these high-energy states in terms of the structures and sequences of their resident proteins. The Branched Aliphatic Side Chain (BASiC) hypothesis has been developed to explain the observed patterns of protection in a pair of TIM barrel proteins. This hypothesis supposes that the side chains of isoleucine, leucine, and valine (ILV) residues often form large hydrophobic clusters that very effectively impede the penetration of water to their underlying hydrogen bond networks and, thereby, enhance the protection against solvent exchange. The linkage between the secondary and tertiary structures enables these ILV clusters to serve as cores of stability in high-energy partially folded states. Statistically significant correlations between the locations of large ILV clusters in native conformations and strong protection against exchange for a variety of motifs reported in the literature support the generality of the BASiC hypothesis. The results also illustrate the necessity to elaborate this simple hypothesis to account for the roles of adjacent hydrocarbon moieties in defining stability cores of partially folded states along folding reaction coordinates.

Start2Fold: a database of hydrogen/deuterium exchange data on protein folding and stability

Proteins fulfil a wide range of tasks in cells; understanding how they fold into complex three-dimensional (3D) structures and how these structures remain stable while retaining sufficient dynamics for functionality is essential for the interpretation of overall protein behaviour. Since the 1950's, solvent exchange-based methods have been the most powerful experimental means to obtain information on the folding and stability of proteins. Considerable expertise and care were required to obtain the resulting datasets, which, despite their importance and intrinsic value, have never been collected, curated and classified. Start2Fold is an openly accessible database (http://start2fold.eu) of carefully curated hydrogen/deuterium exchange (HDX) data extracted from the literature that is open for new submissions from the community. The database entries contain (i) information on the proteins investigated and the underlying experimental procedures and (ii) the classification of the residues based on their exchange protection levels, also allowing for the instant visualization of the relevant residue groups on the 3D structures of the corresponding proteins. By providing a clear hierarchical framework for the easy sharing, comparison and (re-)interpretation of HDX data, Start2Fold intends to promote a better understanding of how the protein sequence encodes folding and structure as well as the development of new computational methods predicting protein folding and stability.

Cavity as a source of conformational fluctuation and high-energy state: high-pressure NMR study of a cavity-enlarged mutant of T4 lysozyme

Although the structure, function, conformational dynamics, and controlled thermodynamics of proteins are manifested by their corresponding amino acid sequences, the natural rules for molecular design and their corresponding interplay remain obscure. In this study, we focused on the role of internal cavities of proteins in conformational dynamics. We investigated the pressure-induced responses from the cavity-enlarged L99A mutant of T4 lysozyme, using high-pressure NMR spectroscopy. The signal intensities of the methyl groups in the (1)H/(13)C heteronuclear single quantum correlation spectra, particularly those around the enlarged cavity, decreased with the increasing pressure, and disappeared at 200 MPa, without the appearance of new resonances, thus indicating the presence of heterogeneous conformations around the cavity within the ground state ensemble. Above 200 MPa, the signal intensities of >20 methyl groups gradually decreased with the increasing pressure, without the appearance of new resonances. Interestingly, these residues closely matched those sensing a large conformational change between the ground- and high-energy states, at atmospheric pressure. (13)C and (1)H NMR line-shape simulations showed that the pressure-induced loss in the peak intensity could be explained by the increase in the high-energy state population. In this high-energy state, the aromatic side chain of F114 gets flipped into the enlarged cavity. The accommodation of the phenylalanine ring into the efficiently packed cavity may decrease the partial molar volume of the high-energy state, relative to the ground state. We suggest that the enlarged cavity is involved in the conformational transition to high-energy states and in the volume fluctuation of the ground state.

How amide hydrogens exchange in native proteins

Amide hydrogen exchange (HX) is widely used in protein biophysics even though our ignorance about the HX mechanism makes data interpretation imprecise. Notably, the open exchange-competent conformational state has not been identified. Based on analysis of an ultralong molecular dynamics trajectory of the protein BPTI, we propose that the open (O) states for amides that exchange by subglobal fluctuations are locally distorted conformations with two water molecules directly coordinated to the N-H group. The HX protection factors computed from the relative O-state populations agree well with experiment. The O states of different amides show little or no temporal correlation, even if adjacent residues unfold cooperatively. The mean residence time of the O state is ∼100 ps for all examined amides, so the large variation in measured HX rate must be attributed to the opening frequency. A few amides gain solvent access via tunnels or pores penetrated by water chains including native internal water molecules, but most amides access solvent by more local structural distortions. In either case, we argue that an overcoordinated N-H group is necessary for efficient proton transfer by Grotthuss-type structural diffusion.

High-Pressure NMR Spectroscopy Reveals Functional Sub-states of Ubiquitin and Ubiquitin-Like Proteins

High-pressure nuclear magnetic resonance (NMR) spectroscopy has revealed that ubiquitin has at least two high-energy states--an alternatively folded state N2 and a locally disordered state I--between the basic folded state N1 and totally unfolded U state. The high-energy states are conserved among ubiquitin-like post-translational modifiers, ubiquitin, NEDD8, and SUMO-2, showing the E1-E2-E3 cascade reaction. It is quite intriguing that structurally similar high-energy states are evolutionally conserved in the ubiquitin-like modifiers, and the thermodynamic stabilities vary among the proteins. To investigate atomic details of the high-energy states, a Q41N mutant of ubiquitin was created as a structural model of N2, which is 71% populated even at atmospheric pressure. The convergent structure of the "pure" N2 state was obtained by nuclear Overhauser effect (NOE)-based structural analysis of the Q41N mutant at 2.5 kbar, where the N2 state is 97% populated. The N2 state of ubiquitin is closely similar to the conformation of the protein bound to the ubiquitin-activating enzyme E1. The recognition of E1 by ubiquitin is best explained by conformational selection rather than by induced-fit motion.

Proline substitution of dimer interface β-strand residues as a strategy for the design of functional monomeric proteins

Proteins that exist in monomer-dimer equilibrium can be found in all organisms ranging from bacteria to humans; this facilitates fine-tuning of activities from signaling to catalysis. However, studying the structural basis of monomer function that naturally exists in monomer-dimer equilibrium is challenging, and most studies to date on designing monomers have focused on disrupting packing or electrostatic interactions that stabilize the dimer interface. In this study, we show that disrupting backbone H-bonding interactions by substituting dimer interface β-strand residues with proline (Pro) results in fully folded and functional monomers, by exploiting proline's unique feature, the lack of a backbone amide proton. In interleukin-8, we substituted Pro for each of the three residues that form H-bonds across the dimer interface β-strands. We characterized the structures, dynamics, stability, dimerization state, and activity using NMR, molecular dynamics simulations, fluorescence, and functional assays. Our studies show that a single Pro substitution at the middle of the dimer interface β-strand is sufficient to generate a fully functional monomer. Interestingly, double Pro substitutions, compared to single Pro substitution, resulted in higher stability without compromising native monomer fold or function. We propose that Pro substitution of interface β-strand residues is a viable strategy for generating functional monomers of dimeric, and potentially tetrameric and higher-order oligomeric proteins.

Solution structure of CXCL5--a novel chemokine and adipokine implicated in inflammation and obesity

The chemokine CXCL5 is selectively expressed in highly specialized cells such as epithelial type II cells in the lung and white adipose tissue macrophages in muscle, where it mediates diverse functions from combating microbial infections by regulating neutrophil trafficking to promoting obesity by inhibiting insulin signaling. Currently very little is known regarding the structural basis of how CXCL5 mediates its novel functions. Towards this missing knowledge, we have solved the solution structure of the CXCL5 dimer by NMR spectroscopy. CXCL5 is a member of a subset of seven CXCR2-activating chemokines (CAC) that are characterized by the highly conserved ELR motif in the N-terminal tail. The structure shows that CXCL5 adopts the typical chemokine fold, but also reveals several distinct differences in the 30 s loop and N-terminal residues; not surprisingly, crosstalk between N-terminal and 30 s loop residues have been implicated as a major determinant of receptor activity. CAC function also involves binding to highly sulfated glycosaminoglycans (GAG), and the CXCL5 structure reveals a distinct distribution of positively charged residues, suggesting that differences in GAG interactions also influence function. The availability of the structure should now facilitate the design of experiments to better understand the molecular basis of various CXCL5 functions, and also serve as a template for the design of inhibitors for use in a clinical setting.

Epitope mapping of inhibitory antibodies targeting the C2 domain of coagulation factor VIII by hydrogen-deuterium exchange mass spectrometry

BACKGROUND

The development of anti-factor VIII antibodies (inhibitors) is a significant complication in the management of patients with hemophilia A, leading to significant increases in morbidity and treatment cost. Using a panel of mAbs against different epitopes on FVIII, we have recently shown that epitope specificity, inhibitor kinetics and time to maximum inhibition are more important than inhibitor titer in predicting responses to FVIII and the combination of FVIII and recombinant FVIIa. In particular, a subset of high-titer inhibitors responded to high-dose FVIII, which would not be predicted on the basis of their inhibitor titer alone. Thus, the ability to quickly map the epitope spectrum of patient plasma with a clinically feasible assay may fundamentally change how clinicians approach the treatment of high-titer inhibitor patients.

OBJECTIVES

To map the epitopes of anti-FVIII mAbs, three of which are classic inhibitors and one of which is a non-classic inhibitor, by the use of hydrogen-deuterium exchange coupled with mass spectrometry (HDX-MS).

METHODS

The binding epitopes of four mAbs targeting the FVIII C2 domain were mapped with HDX-MS.

RESULTS

The epitopes determined with HDX-MS are consistent with those obtained earlier through structural characterization and antibody competition assays. In addition, classic and non-classic inhibitor epitopes could be distinguished by the use of a limited subset of C2 domain-derived peptic fragments.

CONCLUSION

Our results demonstrate the effectiveness and robustness of the HDX-MS method for epitope mapping, and suggest a potential role of rapid mapping of FVIII inhibitor epitopes in facilitating individualized treatment of inhibitor patients.

Autoinhibition of ETV6 (TEL) DNA binding: appended helices sterically block the ETS domain

ETV6 (or TEL), a transcriptional repressor belonging to the ETS family, is frequently involved in chromosomal translocations linked with human cancers. It displays a DNA-binding mode distinct from other ETS proteins due to the presence of a self-associating PNT domain. In this study, we used NMR spectroscopy to dissect the structural and dynamic bases for the autoinhibition of ETV6 DNA binding by sequences C-terminal to its ETS domain. The C-terminal inhibitory domain (CID) contains two helices, H4 and H5, which sterically block the DNA-binding interface of the ETS domain. Importantly, these appended helices are only marginally stable as revealed by amide hydrogen exchange and (15)N relaxation measurements. The CID is thus poised to undergo a facile conformational change as required for DNA binding. The CID also dampens millisecond timescale motions of the ETS domain hypothesized to be critical for the recognition of specific ETS target sequences. This work illustrates the use of appended sequences on conserved structural domains to generate biological diversity and complements previous studies of the allosteric mechanism of ETS1 autoinhibition to reveal both common and divergent features underlying the regulation of DNA binding by ETS transcription factors.

A delicate interplay of structure, dynamics, and thermodynamics for function: a high pressure NMR study of outer surface protein A

Outer surface protein A (OspA) is a crucial protein in the infection of Borrelia burgdorferi causing Lyme disease. We studied conformational fluctuations of OspA with high-pressure (15)N/(1)H two-dimensional NMR along with high-pressure fluorescence spectroscopy. We found evidence within folded, native OspA for rapid local fluctuations of the polypeptide backbone in the nonglobular single layer β-sheet connecting the N- and C-terminal domains with τ << ms, which may give the two domains certain independence in mobility and thermodynamic stability. Furthermore, we found that folded, native OspA is in equilibrium (τ >> ms) with a minor conformer I, which is almost fully disordered and hydrated for the entire C-terminal part of the polypeptide chain from β8 to the C-terminus. Conformer I is characterized with ΔG(0) = 32 ± 9 kJ/mol and ΔV(0) = -140 ± 40 mL/mol, populating only ∼0.001% at 40°C at 0.1 MPa, pH 5.9. Because in the folded conformer the receptor binding epitope of OspA is buried in the C-terminal domain, its transition into conformer I under in vivo conditions may be critical for the infection of B. burgdorferi. The formation and stability of the peculiar conformer I are apparently supported by a large packing defect or cavity located in the C-terminal domain.

Electrostatics of hydrogen exchange for analyzing protein flexibility

Electrostatic interactions at the protein-aqueous interface modulate the reactivity of solvent-exposed backbone amides by a factor of at least a billion fold. The brief (∼10 ps) lifetime of the peptide anion formed during the hydroxide-catalyzed exchange reaction helps enable the experimental rates to be robustly predictable by continuum dielectric methods. Since this ability to predict the structural dependence of exchange reactivity also applies to the protein amide hydrogens that are only rarely exposed to the bulk solvent phase, electrostatic analysis of the experimental exchange rates provides an effective assessment of whether a given model ensemble is consistent with the properly weighted Boltzmann conformational distribution of the protein native state.

Investigating homology between proteins using energetic profiles

Accumulated experimental observations demonstrate that protein stability is often preserved upon conservative point mutation. In contrast, less is known about the effects of large sequence or structure changes on the stability of a particular fold. Almost completely unknown is the degree to which stability of different regions of a protein is generally preserved throughout evolution. In this work, these questions are addressed through thermodynamic analysis of a large representative sample of protein fold space based on remote, yet accepted, homology. More than 3,000 proteins were computationally analyzed using the structural-thermodynamic algorithm COREX/BEST. Estimated position-specific stability (i.e., local Gibbs free energy of folding) and its component enthalpy and entropy were quantitatively compared between all proteins in the sample according to all-vs.-all pairwise structural alignment. It was discovered that the local stabilities of homologous pairs were significantly more correlated than those of non-homologous pairs, indicating that local stability was indeed generally conserved throughout evolution. However, the position-specific enthalpy and entropy underlying stability were less correlated, suggesting that the overall regional stability of a protein was more important than the thermodynamic mechanism utilized to achieve that stability. Finally, two different types of statistically exceptional evolutionary structure-thermodynamic relationships were noted. First, many homologous proteins contained regions of similar thermodynamics despite localized structure change, suggesting a thermodynamic mechanism enabling evolutionary fold change. Second, some homologous proteins with extremely similar structures nonetheless exhibited different local stabilities, a phenomenon previously observed experimentally in this laboratory. These two observations, in conjunction with the principal conclusion that homologous proteins generally conserved local stability, may provide guidance for a future thermodynamically informed classification of protein homology.

Protein structure and function are fundamentally determined by thermodynamics. However, for technical as well as historical reasons, current evolutionary classification schemes and bioinformatics tools do not fully utilize thermodynamic information to describe or analyze proteins. In this work, we address this deficiency by computationally estimating the position-specific thermodynamic quantities of stability (ΔG), enthalpy (ΔH), and entropy (TΔS) for a large and diverse representative sample of protein structures. The sample was drawn from an expertly curated database, such that accepted evolutionary relationships existed for all protein pairs. Importantly, trivial relationships between pairs highly similar in amino acid sequence were explicitly excluded. We found that all position-specific thermodynamic quantities ΔG, ΔH, and TΔS were more similar between proteins that were evolutionarily related (i.e., homologous), and were less similar between proteins that were not evolutionarily related (i.e., non-homologous), with stability being particularly similar between homologous proteins. However, interesting statistically significant exceptions to these trends were observed, exceptions that could indicate novel processes of functional adaptation or evolutionary fold change, mediated by thermodynamics, for the proteins involved. Taken together, these results expand our understanding of the role of thermodynamics in protein evolution and suggest an organizational framework for a future thermodynamically-informed classification of protein homology.

Investigating solution-phase protein structure and dynamics by hydrogen exchange mass spectrometry

By taking advantage of labeling methods such as hydrogen exchange (HX), many details about protein conformation, dynamics, and interactions can be revealed by mass spectrometry. In this unit, hydrogen exchange theory is discussed as it applies to HX-MS protocols, the practice of HX-MS including data analysis and interpretation is explained in detail, and recent advancements in technology which greatly increase the depth of information gained from the technique are highlighted.

Accessibility changes within diphtheria toxin T domain when in the functional molten globule state, as determined using hydrogen/deuterium exchange measurements

The translocation domain (T domain) of diphtheria toxin adopts a partially folded state, the so-called molten globule state, to become functional at acidic pH. We compared, using hydrogen/deuterium exchange experiments associated with MS, the structures of the T domain in its soluble folded state at neutral pH and in its functional molten globule state at acidic pH. In the native state, the alpha-helices TH5 and TH8 are identified as the core of the domain. Based on the high-resolution structure of the T domain, we propose that TH8 is highly protected because it is buried within the native structure. According to the same structure, TH5 is partly accessible at the surface of the T domain. We propose that its high protection is caused by the formation of dimers. Within the molten globule state, high protection is still observed within the helical hairpin TH8-TH9, which is responsible for the insertion of the T domain into the membrane. In the absence of the lipid bilayer, this hydrophobic part of the domain self-assembles, leading to the formation of oligomers. Overall, hydrogen/deuterium-exchange measurements allow the analysis of interaction contacts within small oligomers made of partially folded proteins. Such information, together with crystal structure data, are particularly valuable for using to analyze the self-assembly of proteins.

Pro-interleukin (IL)-1beta shares a core region of stability as compared with mature IL-1beta while maintaining a distinctly different configurational landscape: a comparative hydrogen/deuterium exchange mass spectrometry study

Interleukin-1beta (IL-1beta) is a master cytokine involved in initiating the innate immune response in vertebrates (Dinarello, C. A. (1994) FASEB J. 8, 1314-1325). It is first synthesized as an inactive 269-residue precursor (pro-interleukin-1beta or pro-IL-1beta). Pro-IL-1beta requires processing by caspase-1 to generate the active, mature 153-residue cytokine. In this study, we combined hydrogen/deuterium exchange mass spectrometry, circular dichroism spectroscopy, and enzymatic digestion comparative studies to investigate the configurational landscape of pro-IL-1beta and the role the N terminus plays in modulating the landscape. We find that the N terminus keeps pro-IL-1beta in a protease-labile state while maintaining a core region of stability in the C-terminal region, the eventual mature protein. In mature IL-1beta, this highly protected region maps back to the area protected earliest in the NMR studies characterizing an on-route kinetic refolding intermediate. This protected region also encompasses two important functional loops that participate in the IL-1beta/receptor binding interface required for biological activity. We propose that the purpose of the N-terminal precursor region in pro-IL-1beta is to suppress the function of the eventual mature region while keeping a structurally and also functionally important core region primed for the final folding into the native, active state of the mature protein. The presence of the self-inhibiting precursor region provides yet another layer of regulation in the life cycle of this important cytokine.

Predicting protein folding cores by empirical potential functions

Theoretical and in vitro experiments suggest that protein folding cores form early in the process of folding, and that proteins may have evolved to optimize both folding speed and native-state stability. In our previous work (Chen et al., Structure, 14 (2006) 1401), we developed a set of empirical potential functions and used them to analyze interaction energies among secondary-structure elements in two beta-sandwich proteins. Our work on this group of proteins demonstrated that the predicted folding core also harbors residues that form native-like interactions early in the folding reaction. In the current work, we have tested our empirical potential functions on structurally-different proteins for which the folding cores have been revealed by protein hydrogen-deuterium exchange experiments. Using a set of 29 unrelated proteins, which have been extensively studied in the literature, we demonstrate that the average prediction result from our method is significantly better than predictions based on other computational methods. Our study is an important step towards the ultimate goal of understanding the correlation between folding cores and native structures.

Dynameomics: large-scale assessment of native protein flexibility

Structure is only the first step in understanding the interactions and functions of proteins. In this paper, we explore the flexibility of proteins across a broad database of over 250 solvated protein molecular dynamics simulations in water for an aggregate simulation time of approximately 6 micros. These simulations are from our Dynameomics project, and these proteins represent approximately 75% of all known protein structures. We employ principal component analysis of the atomic coordinates over time to determine the primary axis and magnitude of the flexibility of each atom in a simulation. This technique gives us both a database of flexibility for many protein fold families and a compact visual representation of a particular protein's native-state conformational space, neither of which are available using experimental methods alone. These tools allow us to better understand the nature of protein motion and to describe its relationship to other structural and dynamical characteristics. In addition to reporting general properties of protein flexibility and detailing many dynamic motifs, we characterize the relationship between protein native-state flexibility and early events in thermal unfolding and show that flexibility predicts how a protein will begin to unfold. We provide evidence that fold families have conserved flexibility patterns, and family members who deviate from the conserved patterns have very low sequence identity. Finally, we examine novel aspects of highly inflexible loops that are as important to structural integrity as conventional secondary structure. These loops, which are difficult if not impossible to locate without dynamic data, may constitute new structural motifs.

Residue-wise conformational stability of DLC8 dimer from native-state hydrogen exchange

Dynein light chain (DLC8) is the smallest subunit of the dynein motor complex, which is known to act as a cargo adaptor in intracellular trafficking. The protein exists as a pure dimer at physiological pH and a completely folded monomer below pH 4. Here, we have determined the energy landscape of the dimeric protein using a combination of optical techniques and native-state hydrogen exchange of amide groups, the former giving the global features and the latter yielding the residue level details. The data indicated the presence of intermediates along the equilibrium unfolding transition. The hydrogen exchange data suggested that the molecule has differential stability in its various segments. We deduce from the free energy data that the antiparallel beta-sheets (beta4 and beta5) that form the hydrophobic core of the protein and the alpha2 helix, all of which are highly protected with regard to hydrogen exchange, contribute significantly to the initial step of the protein folding mechanism. Denaturant-dependent hydrogen exchange indicated further that some amides exchange via local fluctuations, whereas there are others which exchange via global unfolding events. Implications of these to cargo adaptability of the dimer are discussed.

Predicting protein folding pathways at the mesoscopic level based on native interactions between secondary structure elements

Background

Since experimental determination of protein folding pathways remains difficult, computational techniques are often used to simulate protein folding. Most current techniques to predict protein folding pathways are computationally intensive and are suitable only for small proteins.

Results

By assuming that the native structure of a protein is known and representing each intermediate conformation as a collection of fully folded structures in which each of them contains a set of interacting secondary structure elements, we show that it is possible to significantly reduce the conformation space while still being able to predict the most energetically favorable folding pathway of large proteins with hundreds of residues at the mesoscopic level, including the pig muscle phosphoglycerate kinase with 416 residues. The model is detailed enough to distinguish between different folding pathways of structurally very similar proteins, including the streptococcal protein G and the peptostreptococcal protein L. The model is also able to recognize the differences between the folding pathways of protein G and its two structurally similar variants NuG1 and NuG2, which are even harder to distinguish. We show that this strategy can produce accurate predictions on many other proteins with experimentally determined intermediate folding states.

Conclusion

Our technique is efficient enough to predict folding pathways for both large and small proteins at the mesoscopic level. Such a strategy is often the only feasible choice for large proteins. A software program implementing this strategy (SSFold) is available at .

A billion-fold range in acidity for the solvent-exposed amides of Pyrococcus furiosus rubredoxin

The exchange rates of the static solvent-accessible amide hydrogens of Pyrococcus furiosus rubredoxin range from near the diffusion-limited rate to a billion-fold slower for the non-hydrogen-bonded Val 38 (eubacterial numbering). Hydrogen exchange directly monitors the kinetic acidity of the peptide nitrogen. Electrostatic solvation free energies were calculated by Poisson-Boltzmann methods for the individual peptide anions that form during the hydroxide-catalyzed exchange reaction to examine how well the predicted thermodynamic acidities match the experimentally determined kinetic acidities. With the exception of the Ile 12 amide, the differential exchange rate constant for each solvent-exposed amide proton that is not hydrogen bonded to a backbone carbonyl can be predicted within a factor of 6 (10 (0.78)) root-mean-square deviation (rmsd) using the CHARMM22 electrostatic parameter set and an internal dielectric value of 3. Under equivalent conditions, the PARSE parameter set yields a larger rmsd value of 1.28 pH units, while the AMBER parm99 parameter set resulted in a considerably poorer correlation. Either increasing the internal dielectric value to 4 or reducing it to a value of 2 significantly degrades the quality of the prediction. Assigning the excess charge of the peptide anion equally between the peptide nitrogen and the carbonyl oxygen also reduces the correlation to the experimental data. These continuum electrostatic calculations were further analyzed to characterize the specific structural elements that appear to be responsible for the wide range of peptide acidities observed for these solvent-exposed amides. The striking heterogeneity in the potential at sites along the protein-solvent interface should prove germane to the ongoing challenge of quantifying the contribution that electrostatic interactions make to the catalytic acceleration achieved by enzymes.

Folding and assembly of large macromolecular complexes monitored by hydrogen-deuterium exchange and mass spectrometry

Recent advances in protein mass spectrometry (MS) have enabled determinations of hydrogen deuterium exchange (HDX) in large macromolecular complexes. HDX-MS became a valuable tool to follow protein folding, assembly and aggregation. The methodology has a wide range of applications in biotechnology ranging from quality control for over-expressed proteins and their complexes to screening of potential ligands and inhibitors. This review provides an introduction to protein folding and assembly followed by the principles of HDX and MS detection, and concludes with selected examples of applications that might be of interest to the biotechnology community.