Hydrogen exchange in proteins

Considerations in the Analysis of Hydrogen Exchange Mass Spectrometry Data

A major component of a hydrogen exchange mass spectrometry experiment is the analysis of protein and peptide mass spectra to yield information about deuterium incorporation. The processing of data that are produced includes the identification of each peptic peptide to create a master table/array of peptide identity that typically includes sequence, retention time and retention time range, mass range, and undeuterated mass. The amount of deuterium incorporated into each of the peptides in this array must then be determined. Various software platforms have been developed in order to perform this specific type of data analysis. We describe the fundamental parameters to be considered at each step along the way and how data processing, either by an individual or by software, must approach the analysis.

PSX, Protein–Solvent Exchange: software for calculation of deuterium-exchange effects in small-angle neutron scattering measurements from protein coordinates

Recent developments in neutron scattering instrumentation and sample handling have enabled studies of more complex biological samples and measurements at shorter exposure times. The experiments are typically conducted in D2O-based buffers to emphasize or diminish scattering from a particular component or to minimize background noise in the experiment. To extract most information from such experiments it is thus desirable to determine accurate estimates of how and when closely bound hydrogen atoms from the biomolecule exchange with the deuterium in the solvent. This article introduces and documents software, PSX, for exploring the effect of hydrogen–deuterium exchange for proteins solubilized in D2O as well as the underlying bioinformatical models. The software aims to be generally applicable for any atomistic structure of a protein and its surrounding environment, and thus captures effects of both heterogenous exchange rates throughout the protein structure and varying the experimental conditions such as pH and temperature. The paper concludes with examples of applications and estimates of the effect in typical scenarios emerging in small-angle neutron scattering on biological macromolecules in solution. The analysis presented here suggests that the common assumption of 90% exchange is in many cases an overestimate with the rapid sample handling systems currently available, which leads to fitting and calibration issues when analysing the data. Source code for the presented software is available from an online repository in which it is published under version 3 of the GNU publishing licence.

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.

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

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

Interlaboratory Comparison of Hydrogen-Deuterium Exchange Mass Spectrometry Measurements of the Fab Fragment of NISTmAb

Hydrogen-deuterium exchange mass spectrometry (HDX-MS) is an established, powerful tool for investigating protein-ligand interactions, protein folding, and protein dynamics. However, HDX-MS is still an emergent tool for quality control of biopharmaceuticals and for establishing dynamic similarity between a biosimilar and an innovator therapeutic. Because industry will conduct quality control and similarity measurements over a product lifetime and in multiple locations, an understanding of HDX-MS reproducibility is critical. To determine the reproducibility of continuous-labeling, bottom-up HDX-MS measurements, the present interlaboratory comparison project evaluated deuterium uptake data from the Fab fragment of NISTmAb reference material (PDB: 5K8A ) from 15 laboratories. Laboratories reported ∼89 800 centroid measurements for 430 proteolytic peptide sequences of the Fab fragment (∼78 900 centroids), giving ∼100% coverage, and ∼10 900 centroid measurements for 77 peptide sequences of the Fc fragment. Nearly half of peptide sequences are unique to the reporting laboratory, and only two sequences are reported by all laboratories. The majority of the laboratories (87%) exhibited centroid mass laboratory repeatability precisions of ⟨ sLab⟩ ≤ (0.15 ± 0.01) Da (1σx̅). All laboratories achieved ⟨sLab⟩ ≤ 0.4 Da. For immersions of protein at THDX = (3.6 to 25) °C and for D2O exchange times of tHDX = (30 s to 4 h) the reproducibility of back-exchange corrected, deuterium uptake measurements for the 15 laboratories is σreproducibility15 Laboratories( tHDX) = (9.0 ± 0.9) % (1σ). A nine laboratory cohort that immersed samples at THDX = 25 °C exhibited reproducibility of σreproducibility25C cohort( tHDX) = (6.5 ± 0.6) % for back-exchange corrected, deuterium uptake measurements.

Investigating the Conformational Response of the Sortilin Receptor upon Binding Endogenous Peptide- and Protein Ligands by HDX-MS

Sortilin is a multifunctional neuronal receptor involved in sorting of neurotrophic factors and apoptosis signaling. So far, structural characterization of sortilin and its endogenous ligands has been limited to crystallographic studies of sortilin in complex with the neuropeptide neurotensin. Here, we use hydrogen/deuterium exchange mass spectrometry to investigate the conformational response of sortilin to binding biological ligands including the peptides neurotensin and the sortilin propeptide and the proteins progranulin and pro-nerve growth factor-β. The results show that the ligands use two binding sites inside the cavity of the β-propeller of sortilin. However, ligands have distinct differences in their conformational impact on the receptor. Interestingly, the protein ligands induce conformational stabilization in a remote membrane-proximal domain, hinting at an unknown conformational link between the ligand binding region and this membrane-proximal region of sortilin. Our findings improve our structural understanding of sortilin and how it mediates diverse ligand-dependent functions important in neurobiology.

Conformational dynamics of the human serotonin transporter during substrate and drug binding

The serotonin transporter (SERT), a member of the neurotransmitter:sodium symporter family, is responsible for termination of serotonergic signaling by re-uptake of serotonin (5-HT) into the presynaptic neuron. Its key role in synaptic transmission makes it a major drug target, e.g. for the treatment of depression, anxiety and post-traumatic stress. Here, we apply hydrogen-deuterium exchange mass spectrometry to probe the conformational dynamics of human SERT in the absence and presence of known substrates and targeted drugs. Our results reveal significant changes in dynamics in regions TM1, EL3, EL4, and TM12 upon binding co-transported ions (Na+/K+) and ligand-mediated changes in TM1, EL3 and EL4 upon binding 5-HT, the drugs S-citalopram, cocaine and ibogaine. Our results provide a comprehensive direct view of the conformational response of SERT upon binding both biologically relevant substrate/ions and ligands of pharmaceutical interest, thus advancing our understanding of the structure-function relationship in SERT.

Global Implications of Local Unfolding Phenomena, Probed by Cysteine Reactivity in Human Frataxin

Local events that affect specific regions of proteins are of utmost relevance for stability and function. The aim of this study is to quantitatively assess the importance of locally-focused dynamics by means of a simple chemical modification procedure. Taking human Frataxin as a working model, we investigated local fluctuations of the C-terminal region (the last 16 residues of the protein) by means of three L → C replacement mutants: L98C, L200C and L203C. The conformation and thermodynamic stability of each variant was assessed. All the variants exhibited native features and high stabilities: 9.1 (wild type), 8.1 (L198C), 7.0 (L200C) and 10.0 kcal mol^-1 (L203C). In addition, kinetic rates of Cys chemical modification by DTNB and DTDPy were measured, conformational dynamics data were extracted and free energy for the local unfolding of the C-terminal region was estimated. The analysis of these results indicates that the conformation of the C-terminal region fluctuates with partial independence from global unfolding events. Additionally, numerical fittings of the kinetic model of the process suggest that the local transition occurs in the seconds to minutes timescale. In fact, standard free energy differences for local unfolding were found to be significantly lower than those of the global unfolding reaction, showing that chemical modification results may not be explained in terms of the global unfolding reaction alone. These results provide unequivocal experimental evidence of local phenomena with global effects and contribute to understanding how global and local stability are linked to protein dynamics.

High-Resolution HDX-MS of Cytochrome c Using Pepsin/Fungal Protease Type XIII Mixed Bed Column

A pepsin/FPXIII (protease from Aspergillus saitoi, type XIII) mixed bed column significantly improved the resolution of bottom-up hydrogen/deuterium exchange mass spectrometry (HDX-MS) data compared with a pepsin-only column. The HDX-MS method using the mixed bed column determined 65 amide hydrogen exchange rates out of one hundred cytochrome c backbone amide hydrogens. Different cleavage specificities of the two enzymes generated 138 unique high-quality peptic fragments, which allows fine sub-localization of deuterium. The exchange rates determined in this method are consistent within the current study as well as with the previous HDX-NMR study. High-resolution HDX-MS data can determine the exchange rate of each residue not the deuterium buildup curve of a peptic fragment. The exchange rates provide more precise and quantitative measurements of protein dynamics in a more reproducible manner. Graphical Abstract ᅟ.

Role of domain interactions in the aggregation of full-length immunoglobulin light chains

Amyloid light-chain (LC) amyloidosis is a protein misfolding disease in which the aggregation of an overexpressed antibody LC from a clonal plasma cell leads to organ toxicity and patient death if left untreated. While the overall dimeric architecture of LC molecules is established, with each LC composed of variable (V_L) and constant (C_L) domains, the relative contributions of LC domain-domain interfaces and intrinsic domain stabilities to protection against LC aggregation are not well understood. To address these topics we have engineered a number of domain-destabilized LC mutants and used solution NMR spectroscopy to characterize their structural properties and intrinsic stabilities. Moreover, we used fluorescence spectroscopy to assay their aggregation propensities. Our results point to the importance of both dimerization strength and intrinsic monomer stability in stabilizing V_L domains against aggregation. Notably, in all cases considered V_L domains aggregate at least 10-fold faster than full-length LCs, establishing the important protective role of C_L domains. A strong protective coupling is found between V_L-V_L and C_L-C_L dimer interfaces, with destabilization of one interface adversely affecting the stability of the other. Fibril formation is observed when either the V_L or C_L domain in the full-length protein is severely destabilized (i.e., where domain unfolding free energies are less than 2 kcal/mol). The important role of C_L domains in preventing aggregation highlights the potential of the C_L-C_L interface as a target for the development of drugs to stabilize the dimeric LC structure and hence prevent LC amyloidosis.

Hydrogen/deuterium exchange in mass spectrometry

The isotopic exchange approach is in use since the first observation of such reactions in 1933 by Lewis. This approach allows the investigation of the pathways of chemical and biochemical reactions, determination of structure, composition, and conformation of molecules. Mass spectrometry has now become one of the most important analytical tools for the monitoring of the isotopic exchange reactions. Investigation of conformational dynamics of proteins, quantitative measurements, obtaining chemical, and structural information about individual compounds of the complex natural mixtures are mainly based on the use of isotope exchange in combination with high resolution mass spectrometry. The most important reaction is the Hydrogen/Deuterium exchange, which is mainly performed in the solution. Recently we have developed the approach allowing performing of the Hydrogen/Deuterium reaction on-line directly in the ionization source under atmospheric pressure. Such approach simplifies the sample preparation and can accelerate the exchange reaction so that certain hydrogens that are considered as non-labile will also participate in the exchange. The use of in-ionization source H/D exchange in modern mass spectrometry for structural elucidation of molecules serves as the basic theme in this review. We will focus on the mechanisms of the isotopic exchange reactions and on the application of in-ESI, in-APCI, and in-APPI source Hydrogen/Deuterium exchange for the investigation of petroleum, natural organic matter, oligosaccharides, and proteins including protein-protein complexes. The simple scenario for adaptation of H/D exchange reactions into mass spectrometric method is also highlighted along with a couple of examples collected from previous studies.

Evidence for a Partially Stalled γ Rotor in F1-ATPase from Hydrogen-Deuterium Exchange Experiments and Molecular Dynamics Simulations

F₁-ATPase uses ATP hydrolysis to drive rotation of the γ subunit. The γ C-terminal helix constitutes the rotor tip that is seated in an apical bearing formed by α₃β₃. It remains uncertain to what extent the γ conformation during rotation differs from that seen in rigid crystal structures. Existing models assume that the entire γ subunit participates in every rotation. Here we interrogated E. coli F₁-ATPase by hydrogen-deuterium exchange (HDX) mass spectrometry. Rotation of γ caused greatly enhanced deuteration in the γ C-terminal helix. The HDX kinetics implied that most F₁ complexes operate with an intact rotor at any given time, but that the rotor tip is prone to occasional unfolding. A molecular dynamics (MD) strategy was developed to model the off-axis forces acting on γ. MD runs showed stalling of the rotor tip and unfolding of the γ C-terminal helix. MD-predicted H-bond opening events coincided with experimental HDX patterns. Our data suggest that in vitro operation of F₁-ATPase is associated with significant rotational resistance in the apical bearing. These conditions cause the γ C-terminal helix to get "stuck" (and unfold) sporadically while the remainder of γ continues to rotate. This scenario contrasts the traditional "greasy bearing" model that envisions smooth rotation of the γ C-terminal helix. The fragility of the apical rotor tip in F₁-ATPase is attributed to the absence of a c₁₀ ring that stabilizes the rotation axis in intact F_oF₁. Overall, the MD/HDX strategy introduced here appears well suited for interrogating the inner workings of molecular motors.

Interpreting Hydrogen-Deuterium Exchange Events in Proteins Using Atomistic Simulations: Case Studies on Regulators of G-Protein Signaling Proteins

Hydrogen-deuterium exchange (HDX) experiments are widely used in studies of protein dynamics. To predict the propensity of amide hydrogens for exchange with deuterium, several models have been reported in which computations of amide-hydrogen protection factors are carried out using molecular dynamics (MD) simulations. Given significant variation in the criteria used in different models, the robustness and broader applicability of these models to other proteins, especially homologous proteins showing distinct amide-exchange patterns, remains unknown. The sensitivity of the predictions when MD simulations are conducted with different force-fields is yet to tested and quantified. Using MD simulations and experimental HDX data on three homologous signaling proteins, we report detailed studies quantifying the performance of seven previously reported models (M1-M7) of two general types: empirical and fractional-population models. We find that the empirical models show inconsistent predictions but predictions of the fractional population models are robust. Contrary to previously reported work, we find that the solvent-accessible surface area of amide hydrogens is a useful metric when combined with a new metric defining the distances of amide hydrogens from the first polar atoms in proteins. On the basis of this, we report two new models, one empirical (M8) and one population-based (M9). We find strong protection of amide hydrogens from solvent exchange both within the stable helical motifs and also in the interhelical loops. We further observe that the exchange-competent states of amide hydrogens occur on the sub 100 ps time-scale via localized fluctuations, and such states among amides of a given protein do not appear to show any cooperativity or allosteric coupling.

Modulating long-range energetics via helix stabilization: A case study using T4 lysozyme

Cooperative protein folding requires distant regions of a protein to interact and provide mutual stabilization. The mechanism of this long-distance coupling remains poorly understood. Here, we use T4 lysozyme (T4L*) as a model to investigate long-range communications across two subdomains of a globular protein. T4L* is composed of two structurally distinct subdomains, although it behaves in a two-state manner at equilibrium. The subdomains of T4L* are connected via two topological connections: the N-terminal helix that is structurally part of the C-terminal subdomain (the A-helix) and a long helix that spans both subdomains (the C-helix). To understand the role that the C-helix plays in cooperative folding, we analyzed a circularly permuted version of T4L* (CP13*), whose subdomains are connected only by the C-helix. We demonstrate that when isolated as individual fragments, both subdomains of CP13* can fold autonomously into marginally stable conformations. The energetics of the N-terminal subdomain depend on the formation of a salt bridge known to be important for stability in the full-length protein. We show that the energetic contribution of the salt bridge to the stability of the N-terminal fragment increases when the C-helix is stabilized, such as occurs upon folding of the C-terminal subdomain. These results suggest a model where long-range energetic coupling is mediated by helix stabilization and not specific tertiary interactions.

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.

Site-Specific Hydrogen Exchange in a Membrane Environment Analyzed by Infrared Spectroscopy

Hydrogen exchange is a powerful method to examine macromolecules. In membrane proteins, exchange can distinguish between solvent-accessible and -inaccessible residues due to shielding by the hydrophobic environment of the lipid bilayer. Herein, rather than examining which residues undergo hydrogen exchange, we employ a protocol that enables the full deuteration of all polar hydrogens in a membrane protein. We then measure the impact of hydrogen exchange on the shift of the amide I vibrational mode of individually labeled sites. The results enable us to correlate polarity with vibrational shifts, thereby providing a powerful tool to examine specific locations within a membrane protein in its native membrane environment.

Effects of Disease-Causing Mutations on the Conformation of Human Apolipoprotein A-I in Model Lipoproteins

Plasma high-density lipoproteins (HDLs) are protein-lipid nanoparticles that transport lipids and protect against atherosclerosis. Human apolipoprotein A-I (apoA-I) is the principal HDL protein whose mutations can cause either aberrant lipid metabolism or amyloid disease. Hydrogen-deuterium exchange (HDX) mass spectrometry (MS) was used to study the apoA-I conformation in model discoidal lipoproteins similar in size to large plasma HDL. We examined how point mutations associated with hereditary amyloidosis (F71Y and L170P) or atherosclerosis (L159R) influence the local apoA-I conformation in model lipoproteins. Unlike other apoA-I forms, the large particles showed minimal conformational heterogeneity, suggesting a fully extended protein conformation. Mutation-induced structural perturbations in lipid-bound protein were attenuated compared to the free protein and indicated close coupling between the two belt-forming apoA-I molecules. These perturbations propagated to distant lipoprotein sites, either increasing or decreasing their protection. This HDX MS study of large model HDL, compared with previous studies of smaller particles, ascertained that apoA-I's central region helps accommodate the protein conformation to lipoproteins of various sizes. This study also reveals that the effects of mutations on lipoprotein conformational dynamics are much weaker than those in a lipid-free protein. Interestingly, the mutation-induced perturbations propagate to distant sites nearly 10 nm away and alter their protection in ways that cannot be predicted from the lipoprotein structure and stability. We propose that long-range mutational effects are mediated by both protein and lipid and can influence lipoprotein functionality.

Defining Hsp33's Redox-regulated Chaperone Activity and Mapping Conformational Changes on Hsp33 Using Hydrogen-deuterium Exchange Mass Spectrometry

Living organisms regularly need to cope with fluctuating environments during their life cycle, including changes in temperature, pH, the accumulation of reactive oxygen species, and more. These fluctuations can lead to a widespread protein unfolding, aggregation, and cell death. Therefore, cells have evolved a dynamic and stress-specific network of molecular chaperones, which maintain a "healthy" proteome during stress conditions. ATP-independent chaperones constitute one major class of molecular chaperones, which serve as first-line defense molecules, protecting against protein aggregation in a stress-dependent manner. One feature these chaperones have in common is their ability to utilize structural plasticity for their stress-specific activation, recognition, and release of the misfolded client. In this paper, we focus on the functional and structural analysis of one such intrinsically disordered chaperone, the bacterial redox-regulated Hsp33, which protects proteins against aggregation during oxidative stress. Here, we present a toolbox of diverse techniques for studying redox-regulated chaperone activity, as well as for mapping conformational changes of the chaperone, underlying its activity. Specifically, we describe a workflow which includes the preparation of fully reduced and fully oxidized proteins, followed by an analysis of the chaperone anti-aggregation activity in vitro using light-scattering, focusing on the degree of the anti-aggregation activity and its kinetics. To overcome frequent outliers accumulated during aggregation assays, we describe the usage of Kfits, a novel graphical tool which allows easy processing of kinetic measurements. This tool can be easily applied to other types of kinetic measurements for removing outliers and fitting kinetic parameters. To correlate the function with the protein structure, we describe the setup and workflow of a structural mass spectrometry technique, hydrogen-deuterium exchange mass spectrometry, that allows the mapping of conformational changes on the chaperone and substrate during different stages of Hsp33 activity. The same methodology can be applied to other protein-protein and protein-ligand interactions.

Dynamics of Scabin toxin. A proposal for the binding mode of the DNA substrate

Scabin is a mono-ADP-ribosyltransferase enzyme and is a putative virulence factor produced by the plant pathogen, Streptomyces scabies. Previously, crystal structures of Scabin were solved in the presence and absence of substrate analogues and inhibitors. Herein, experimental (hydrogen-deuterium exchange), simulated (molecular dynamics), and theoretical (Gaussian Network Modeling) approaches were systematically applied to study the dynamics of apo-Scabin in the context of a Scabin·NAD⁺·DNA model. MD simulations revealed that the apo-Scabin solution conformation correlates well with the X-ray crystal structure, beyond the conformation of the exposed, mobile regions. In turn, the MD fluctuations correspond with the crystallographic B-factors, with the fluctuations derived from a Gaussian network model, and with the experimental H/D exchange rates. An Essential Dynamics Analysis identified the dynamic aspects of the toxin as a crab-claw-like mechanism of two topological domains, along with coupled deformations of exposed motifs. The “crab-claw” movement resembles the motion of C3-like toxins and emerges as a property of the central β scaffold of catalytic single domain toxins. The exposure and high mobility of the cis side motifs in the Scabin β-core suggest involvement in DNA substrate binding. A ternary Scabin·NAD⁺·DNA model was produced via an independent docking methodology, where the intermolecular interactions correspond to the region of high mobility identified by dynamics analyses and agree with binding and kinetic data reported for wild-type and Scabin variants. Based on data for the Pierisin-like toxin group, the sequence motif R_β1–R_La–N_Lc–STT_β2–W_PN–W_ARTT–(QxE)_ARTT emerges as a catalytic signature involved in the enzymatic activity of these DNA-acting toxins. However, these results also show that Scabin possesses a unique DNA-binding motif within the Pierisin-like toxin group.

Hydrogen-deuterium exchange reveals long-range dynamical allostery in soybean lipoxygenase

In lipoxygenases, the topologically conserved C-terminal domain catalyzes the oxidation of polyunsaturated fatty acids, generating an assortment of biologically relevant signaling mediators. Plant and animal lipoxygenases also contain a 100-150-amino acid N-terminal C2-like domain that has been implicated in interactions with isolated fatty acids and at the phospholipid bilayer. These interactions may lead to increased substrate availability and contribute to the regulation of active-site catalysis. Because of a lack of structural information, a molecular understanding of this lipid-protein interaction remains unresolved. Herein, we employed hydrogen-deuterium exchange MS (HDXMS) to spatially resolve changes in protein conformation upon interaction of soybean lipoxygenase with a fatty acid surrogate, oleyl sulfate (OS), previously shown to act at a site separate from the substrate-binding site. Specific, OS-induced conformational changes are detected both at the N-terminal domain and within the substrate portal nearly 30 Å away. Combining previously measured kinetic properties in the presence of OS with its impact on the Kd for linoleic acid substrate binding, we conclude that OS binding brings about an increase in rate constants for both the ingress and egress of substrate. We discuss the role of OS-induced changes in protein flexibility in the context of changes in the mechanism of substrate acquisition.

Computational methods and challenges in hydrogen/deuterium exchange mass spectrometry

Hydrogen/Deuterium exchange (HDX) has been applied, since the 1930s, as an analytical tool to study the structure and dynamics of (small) biomolecules. The popularity of using HDX to study proteins increased drastically in the last two decades due to the successful combination with mass spectrometry (MS). Together with this growth in popularity, several technological advances have been made, such as improved quenching and fragmentation. As a consequence of these experimental improvements and the increased use of protein-HDXMS, large amounts of complex data are generated, which require appropriate analysis. Computational analysis of HDXMS requires several steps. A typical workflow for proteins consists of identification of (non-)deuterated peptides or fragments of the protein under study (local analysis), or identification of the deuterated protein as a whole (global analysis); determination of the deuteration level; estimation of the protection extent or exchange rates of the labile backbone amide hydrogen atoms; and a statistically sound interpretation of the estimated protection extent or exchange rates. Several algorithms, specifically designed for HDX analysis, have been proposed. They range from procedures that focus on one specific step in the analysis of HDX data to complete HDX workflow analysis tools. In this review, we provide an overview of the computational methods and discuss outstanding challenges. © 2016 Wiley Periodicals, Inc. Mass Spec Rev 36:649-667, 2017.

Probing Transient Release of Membrane-Sequestered Tyrosine-Based Signaling Motif by Solution NMR Spectroscopy

Membrane sequestration of tyrosine-based signaling motifs of antigen receptors effectively restricts the signaling activities in resting lymphocytes. However, low level of basal signaling in resting cells is required for lymphocyte survival and antigen responsiveness, of which the molecular mechanism remains obscure. Here we probe the transient release of the cytoplasmic domain of the membrane-bound IgG heavy chain (mIgG-tail) by hydrogen exchange NMR spectroscopy, illustrating a dynamic molecular basis for its basal signaling activity. To solve the severe resonance overlap problem in the 2D spectra of mIgG-tail, a non-uniformly sampled pseudo-4D hydrogen exchange NMR experiment has been exploited to quantitatively measure site-specific hydrogen exchange rates. Our solution NMR study reveals transient solvent exposure of the ITT signaling motif that can be further enhanced by calcium ion, and provides insight into the mechanism of lymphocyte basal signaling.

Calcium-Mediated Control of S100 Proteins: Allosteric Communication via an Agitator/Signal Blocking Mechanism

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

NMR Analysis of Amide Hydrogen Exchange Rates in a Pentapeptide-Repeat Protein from A. thaliana

At2g44920 from Arabidopsis thaliana is a pentapeptide-repeat protein (PRP) composed of 25 repeats capped by N- and C-terminal α-helices. PRP structures are dominated by four-sided right-handed β-helices typically consisting of mixtures of type II and type IV β-turns. PRPs adopt repeated five-residue (Rfr) folds with an Rfr consensus sequence (STAV)(D/N)(L/F)(S/T/R)(X). Unlike other PRPs, At2g44920 consists exclusively of type II β-turns. At2g44920 is predicted to be located in the thylakoid lumen although its biochemical function remains unknown. Given its unusual structure, we investigated the biophysical properties of At2g44920 as a representative of the β-helix family to determine if it had exceptional global stability, backbone dynamics, or amide hydrogen exchange rates. Circular dichroism measurements yielded a melting point of 62.8°C, indicating unexceptional global thermal stability. Nuclear spin relaxation measurements indicated that the Rfr-fold core was rigid with order parameters ranging from 0.7 to 0.9. At2g44920 exhibited a striking range of amide hydrogen exchange rates spanning 10 orders of magnitude, with lifetimes ranging from minutes to several months. A weak correlation was found among hydrogen exchange rates, hydrogen bonding energies, and amino acid solvent-accessible areas. Analysis of contributions from fast (approximately picosecond to nanosecond) backbone dynamics to amide hydrogen exchange rates revealed that the average order parameter of amides undergoing fast exchange was significantly smaller compared to those undergoing slow exchange. Importantly, the activation energies for amide hydrogen exchange were found to be generally higher for the slowest exchanging amides in the central Rfr coil and decreased toward the terminal coils. This could be explained by assuming that the concerted motions of two preceding or following coils required for hydrogen bond disruption and amide hydrogen exchange have a higher activation energy compared to that required for displacement of a single coil to facilitate amide hydrogen exchange in either the terminal or penultimate coils.

Solid-State NMR Provides Evidence for Small-Amplitude Slow Domain Motions in a Multispanning Transmembrane α-Helical Protein

Proteins are dynamic entities and populate ensembles of conformations. Transitions between states within a conformational ensemble occur over a broad spectrum of amplitude and time scales, and are often related to biological function. Whereas solid-state NMR (SSNMR) spectroscopy has recently been used to characterize conformational ensembles of proteins in the microcrystalline states, its applications to membrane proteins remain limited. Here we use SSNMR to study conformational dynamics of a seven-helical transmembrane (TM) protein, Anabaena Sensory Rhodopsin (ASR) reconstituted in lipids. We report on site-specific measurements of the 15N longitudinal R1 and rotating frame R1ρ relaxation rates at two fields of 600 and 800 MHz and at two temperatures of 7 and 30 °C. Quantitative analysis of the R1 and R1ρ values and of their field and temperature dependencies provides evidence of motions on at least two time scales. We modeled these motions as fast local motions and slower collective motions of TM helices and of structured loops, and used the simple model-free and extended model-free analyses to fit the data and estimate the amplitudes, time scales and activation energies. Faster picosecond (tens to hundreds of picoseconds) local motions occur throughout the protein and are dominant in the middle portions of the TM helices. In contrast, the amplitudes of the slower collective motions occurring on the nanosecond (tens to hundreds of nanoseconds) time scales, are smaller in the central parts of helices, but increase toward their cytoplasmic sides as well as in the interhelical loops. ASR interacts with a soluble transducer protein on its cytoplasmic surface, and its binding affinity is modulated by light. The larger amplitude of motions on the cytoplasmic side of the TM helices correlates with the ability of ASR to undergo large conformational changes in the process of binding/unbinding the transducer.

Soft interactions and volume exclusion by polymeric crowders can stabilize or destabilize transient structure in disordered proteins depending on polymer concentration

The effects of macromolecular crowding on the transient structure of intrinsically disordered proteins is not well-understood. Crowding by biological molecules inside cells could modulate transient structure and alter IDP function. Volume exclusion theory and observations of structured proteins suggest that IDP transient structure would be stabilized by macromolecular crowding. Amide hydrogen exchange (HX) of IDPs in highly concentrated polymer solutions would provide valuable insights into IDP transient structure under crowded conditions. Here, we have used mass spectrometry to measure HX by a transiently helical random coil domain of the activator of thyroid and retinoid receptor (ACTR) in solutions containing 300 g L^-1 and 400 g L^-1 of Ficoll, a synthetic polysaccharide, using a recently-developed strong cation exchange-based cleanup method [Rusinga, et al., Anal Chem 2017;89:1275-1282]. Transiently helical regions of ACTR exchanged faster in 300 g L^-1 Ficoll than in dilute buffer. In contrast, one transient helix exchanged more slowly in 400 g L^-1 Ficoll. Nonspecific interactions destabilize ACTR helicity in 300 g L^-1 Ficoll because ACTR engages with the Ficoll polymer mesh. In contrast, 400 g L^-1 Ficoll is a semi-dilute solution where ACTR cannot engage the Ficoll mesh. At this higher concentration, volume exclusion stabilizes ACTR helicity because ACTR is compacted in interstitial spaces between Ficoll molecules. Our results suggest that the interplay between nonspecific interactions and volume exclusion in different cellular compartments could modulate IDP function by altering the stability of IDP transient structures. Proteins 2017; 85:1468-1479. © 2017 Wiley Periodicals, Inc.

Role of water in the formation of macromolecular structures

This review shows that water in biological systems is not just a passive liquid solvent but also a partner in the formation of the structure of proteins, nucleic acids and their complexes, thereby contributing to the stability and flexibility required for their proper function. Reciprocally, biological macromolecules affect the state of the water contacting them, so that it is only partly in the normal liquid state, being somewhat ordered when bound to macromolecules. While the compaction of globular proteins results from the reluctance of their hydrophobic groups to interact with water, the collagen superhelix is maintained by water forming a hydroxyproline-controlled frame around this coiled-coil macromolecule. As for DNA, its stability and rigidity are linked to water fixed by AT pairs in the minor groove: this leads to the enthalpic contribution of AT pairs exceeding that of GC pairs, but this is overbalanced by their greater entropy contribution, with the result that AT pairs melt at lower temperatures than GCs. Loss of this water drives transcription factor binding to the minor groove.

Functional Studies on Membrane Proteins by Means of H/D Exchange in Infrared: Structural Changes in Na+ NQR from V. cholerae in the Presence of Lipids

H/D exchange kinetics at the level of the amide proton in the mid infrared (1700-1500 cm^-1) make it possible to study the conformational flexibility of membrane proteins, independent of size or the presence of detergent or lipids. Slow, medium, and fast exchanging domains are distinguished, which reveal a different accessibility to the solvent. Whereas amide hydrogens undergo rapid exchange with solvent in an open structure, hydrogens experience much slower exchange when involved in H-bonded structures or when sterically inaccessible to the solvent. Here, we describe the protocol that was used to study the effect of phospholipids on the overall structure of the Na⁺ NQR from V. cholerae, a sodium pumping membrane protein.

Automated Strong Cation-Exchange Cleanup To Remove Macromolecular Crowding Agents for Protein Hydrogen Exchange Mass Spectrometry

Measuring amide hydrogen exchange (HX) of intrinsically disordered proteins (IDPs) in solutions containing high concentrations of macromolecular crowding agents would give new insights into the structure and dynamics of these proteins under crowded conditions. High concentrations of artificial crowders, required to simulate cellular crowding, introduce overwhelming interferences to mass spectrometry (MS) analysis. We have developed a fully automated, dual-stage online cleanup that uses strong cation-exchange (SCX) followed by reversed-phase desalting to remove Ficoll, a synthetic polymer, for HX-MS analysis of proteins under crowded conditions. We tested the efficiency of our method by measuring the HX-MS signal intensities of myoglobin peptides from crowded samples containing 300 g L^-1 Ficoll and from uncrowded samples. Although there was loss of abundance relative to uncrowded myoglobin analyzed using conventional HX-MS, 97% coverage of the myoglobin sequence was still obtained. Control HX-MS experiments using unstructured peptides labeled at pD 4.0 under crowded and uncrowded conditions confirmed that Ficoll does not alter chemical exchange and that the same extent of HX is achieved in uncrowded solutions as in solutions containing 300 g L^-1 of predeuterated Ficoll. We validated our method by measuring HX of CBP, the intrinsically disordered nuclear coactivator binding domain of CREB binding protein (UniProt CBP_MOUSE P45481 ), residues 2059-2117, at pD 6.5 under crowded and uncrowded conditions. Ficoll induced both protection and deprotection from HX in different regions of CBP, with the greatest deprotection occurring at the edges of helices. These results are consistent with previous observation of IDPs under the influence of synthetic polymers.

A Residue-Resolved Bayesian Approach to Quantitative Interpretation of Hydrogen-Deuterium Exchange from Mass Spectrometry: Application to Characterizing Protein-Ligand Interactions

Characterization of interactions between proteins and other molecules is crucial for understanding the mechanisms of action of biological systems and, thus, drug discovery. An increasingly useful approach to mapping these interactions is measurement of hydrogen/deuterium exchange (HDX) using mass spectrometry (HDX-MS), which measures the time-resolved deuterium incorporation of peptides obtained by enzymatic digestion of the protein. Comparison of exchange rates between apo- and ligand-bound conditions results in a mapping of the differential HDX (ΔHDX) of the ligand. Residue-level analysis of these data, however, must account for experimental error, sparseness, and ambiguity due to overlapping peptides. Here, we propose a Bayesian method consisting of a forward model, noise model, prior probabilities, and a Monte Carlo sampling scheme. This method exploits a residue-resolved exponential rate model of HDX-MS data obtained from all peptides simultaneously, and explicitly models experimental error. The result is the best possible estimate of ΔHDX magnitude and significance for each residue given the data. We demonstrate the method by revealing richer structural interpretation of ΔHDX data on two nuclear receptors: vitamin D-receptor (VDR) and retinoic acid receptor gamma (RORγ). The method is implemented in HDX Workbench and as a standalone module of the open source Integrative Modeling Platform.

Biophysical Evidence for Intrinsic Disorder in the C-terminal Tails of the Epidermal Growth Factor Receptor (EGFR) and HER3 Receptor Tyrosine Kinases

The epidermal growth factor receptor (EGFR)/ErbB family of receptor tyrosine kinases includes oncogenes important in the progression of breast and other cancers, and they are targets for many drug development strategies. Each member of the ErbB family possesses a unique, structurally uncharacterized C-terminal tail that plays an important role in autophosphorylation and signal propagation. To determine whether these C-terminal tails are intrinsically disordered regions, we conducted a battery of biophysical experiments on the EGFR and HER3 tails. Using hydrogen/deuterium exchange mass spectrometry, we measured the conformational dynamics of intracellular half constructs and compared the tails with the ordered kinase domains. The C-terminal tails demonstrate more rapid deuterium exchange behavior when compared with the kinase domains. Next, we expressed and purified EGFR and HER3 tail-only constructs. Results from circular dichroism spectroscopy, size exclusion chromatography with multiangle light scattering, dynamic light scattering, analytical ultracentrifugation, and small angle X-ray scattering each provide evidence that the EGFR and HER3 C-terminal tails are intrinsically disordered with extended, non-globular structure in solution. The intrinsic disorder and extended conformation of these tails may be important for their function by increasing the capture radius and reducing the thermodynamic barriers for binding of downstream signaling proteins.

Protein dynamics and function from solution state NMR spectroscopy

It is well-established that dynamics are central to protein function; their importance is implicitly acknowledged in the principles of the Monod, Wyman and Changeux model of binding cooperativity, which was originally proposed in 1965. Nowadays the concept of protein dynamics is formulated in terms of the energy landscape theory, which can be used to understand protein folding and conformational changes in proteins. Because protein dynamics are so important, a key to understanding protein function at the molecular level is to design experiments that allow their quantitative analysis. Nuclear magnetic resonance (NMR) spectroscopy is uniquely suited for this purpose because major advances in theory, hardware, and experimental methods have made it possible to characterize protein dynamics at an unprecedented level of detail. Unique features of NMR include the ability to quantify dynamics (i) under equilibrium conditions without external perturbations, (ii) using many probes simultaneously, and (iii) over large time intervals. Here we review NMR techniques for quantifying protein dynamics on fast (ps-ns), slow (μs-ms), and very slow (s-min) time scales. These techniques are discussed with reference to some major discoveries in protein science that have been made possible by NMR spectroscopy.

Charge-mediated Fab-Fc interactions in an IgG1 antibody induce reversible self-association, cluster formation, and elevated viscosity

Concentration-dependent reversible self-association (RSA) of monoclonal antibodies (mAbs) poses a challenge to their pharmaceutical development as viable candidates for subcutaneous delivery. While the role of the antigen-binding fragment (Fab) in initiating RSA is well-established, little evidence supports the involvement of the crystallizable fragment (Fc). In this report, a variety of biophysical tools, including hydrogen exchange mass spectrometry, are used to elucidate the protein interface of such non-covalent protein-protein interactions. Using dynamic and static light scattering combined with viscosity measurements, we find that an IgG1 mAb (mAb-J) undergoes RSA primarily through electrostatic interactions and forms a monomer-dimer-tetramer equilibrium. We provide the first direct experimental mapping of the interface formed between the Fab and Fc domains of an antibody at high protein concentrations. Charge distribution heterogeneity between the positively charged interface spanning complementarity-determining regions CDR3H and CDR2L in the Fab and a negatively charged region in C_H3/Fc domain mediates the RSA of mAb-J. When arginine and NaCl are added, they disrupt RSA of mAb-J and decrease the solution viscosity. Fab-Fc domain interactions between mAb monomers may promote the formation of large transient antibody complexes that ultimately cause increases in solution viscosity. Our findings illustrate how limited specific arrangements of amino-acid residues can cause mAbs to undergo RSA at high protein concentrations and how conserved regions in the Fc portion of the antibody can also play an important role in initiating weak and transient protein-protein interactions.

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.

Protein Folding-How and Why: By Hydrogen Exchange, Fragment Separation, and Mass Spectrometry

Advanced hydrogen exchange (HX) methodology can now determine the structure of protein folding intermediates and their progression in folding pathways. Key developments over time include the HX pulse labeling method with nuclear magnetic resonance analysis, the fragment separation method, the addition to it of mass spectrometric (MS) analysis, and recent improvements in the HX MS technique and data analysis. Also, the discovery of protein foldons and their role supplies an essential interpretive link. Recent work using HX pulse labeling with MS analysis finds that a number of proteins fold by stepping through a reproducible sequence of native-like intermediates in an ordered pathway. The stepwise nature of the pathway is dictated by the cooperative foldon unit construction of the protein. The pathway order is determined by a sequential stabilization principle; prior native-like structure guides the formation of adjacent native-like structure. This view does not match the funneled energy landscape paradigm of a very large number of folding tracks, which was framed before foldons were known and is more appropriate for the unguided residue-level search to surmount an initial kinetic barrier rather than for the overall unfolded-state to native-state folding pathway.

Evolution of the concept of conformational dynamics of enzyme functions over half of a century: A personal view

To most physicists, it was always evident that conformational fluctuation is an inherent property of all molecules. Its existence in proteins was mentioned first by Linderström-Lang and Schellman in 1959 based on their hydrogen-deuterium exchange experiments. The "induced fit" mechanism to explain ligand-induced conformational changes was suggested by Koshland in 1958. Straub combined these two concepts in his "fluctuation fit" theory in 1964. The era of protein X-ray crystallography imposed a static view of protein structures. With proteins becoming accessible to NMR analysis, conformational dynamics could be mapped, and a new wave of dynamic interpretation of enzymatic catalysis and molecular recognition appeared. Energy landscapes, energy funnels, conformational selection, conformational distribution shifts are now frequent terms in interpreting biomolecular recognition and enzymatic catalysis. All these interpretations are based on the concept that evolution uses the conformational fluctuations of enzymes to develop efficient and dynamic catalytic machines. In a resurrection of the original "fluctuation fit" concept, it is generally recognized now that spatial and temporal events of catalysis are equally important to describe its mechanism. This special issue, dedicated to the memory of Henryk Eisenberg, prompted us to look back at the last 50 years of development of a concept that-like other important concepts-appeared, evolved and became accepted during the period covered by the scientific lifespan of Henryk.

Load-dependent destabilization of the γ-rotor shaft in FOF1 ATP synthase revealed by hydrogen/deuterium-exchange mass spectrometry

FoF1 is a membrane-bound molecular motor that uses proton-motive force (PMF) to drive the synthesis of ATP from ADP and Pi. Reverse operation generates PMF via ATP hydrolysis. Catalysis in either direction involves rotation of the γε shaft that connects the α3β3 head and the membrane-anchored cn ring. X-ray crystallography and other techniques have provided insights into the structure and function of FoF1 subcomplexes. However, interrogating the conformational dynamics of intact membrane-bound FoF1 during rotational catalysis has proven to be difficult. Here, we use hydrogen/deuterium exchange mass spectrometry to probe the inner workings of FoF1 in its natural membrane-bound state. A pronounced destabilization of the γ C-terminal helix during hydrolysis-driven rotation was observed. This behavior is attributed to torsional stress in γ, arising from γ⋅⋅⋅α3β3 interactions that cause resistance during γ rotation within the apical bearing. Intriguingly, we find that destabilization of γ occurs only when FoF1 operates against a PMF-induced torque; the effect disappears when PMF is eliminated by an uncoupler. This behavior resembles the properties of automotive engines, where bearings inflict greater forces on the crankshaft when operated under load than during idling.

Infrared spectroscopic studies on reaction induced conformational changes in the NADH ubiquinone oxidoreductase (complex I)

Redox-dependent conformational changes are currently discussed to be a crucial part of the reaction mechanism of the respiratory complex I. Specialized difference Fourier transform infrared techniques allow the detection of side-chain movements and minute secondary structure changes. For complex I, (1)H/(2)H exchange kinetics of the amide modes revealed a better accessibility of the backbone in the presence of NADH and quinone. Interestingly, the presence of phospholipids, that is crucial for the catalytic activity of the isolated enzyme complex, changes the overall conformation. When comparing complex I samples from different species, very similar electrochemically induced FTIR difference spectra and very similar rearrangements are reported. Finally, the information obtained with variants and from Zn(2+) inhibited samples for the conformational reorganization of complex I upon electron transfer are discussed in this review. This article is part of a Special Issue entitled Respiratory complex I, edited by Volker Zickermann and Ulrich Brandt.