QUDeX-MS: hydrogen/deuterium exchange calculation for mass spectra with resolved isotopic fine structure

Computational Tools for Hydrogen-Deuterium Exchange Mass Spectrometry Data Analysis

Hydrogen-deuterium exchange (HDX) has become a pivotal method for investigating the structural and dynamic properties of proteins. The versatility and sensitivity of mass spectrometry (MS) made the technique the ideal companion for HDX, and today HDX-MS is addressing a growing number of applications in both academic research and industrial settings. The prolific generation of experimental data has spurred the concurrent development of numerous computational tools, designed to automate parts of the workflow while employing different strategies to achieve common objectives. Various computational methods are available to perform automated peptide searches and identification; different statistical tests have been implemented to quantify differences in the exchange pattern between two or more experimental conditions; alternative strategies have been developed to deconvolve and analyze peptides showing multimodal behavior; and different algorithms have been proposed to computationally increase the resolution of HDX-MS data, with the ultimate aim to provide information at the level of the single residue. This review delves into a comprehensive examination of the merits and drawbacks associated with the diverse strategies implemented by software tools for the analysis of HDX-MS data.

Advances in Hydrogen/Deuterium Exchange Mass Spectrometry and the Pursuit of Challenging Biological Systems

Solution-phase hydrogen/deuterium exchange (HDX) coupled to mass spectrometry (MS) is a widespread tool for structural analysis across academia and the biopharmaceutical industry. By monitoring the exchangeability of backbone amide protons, HDX-MS can reveal information about higher-order structure and dynamics throughout a protein, can track protein folding pathways, map interaction sites, and assess conformational states of protein samples. The combination of the versatility of the hydrogen/deuterium exchange reaction with the sensitivity of mass spectrometry has enabled the study of extremely challenging protein systems, some of which cannot be suitably studied using other techniques. Improvements over the past three decades have continually increased throughput, robustness, and expanded the limits of what is feasible for HDX-MS investigations. To provide an overview for researchers seeking to utilize and derive the most from HDX-MS for protein structural analysis, we summarize the fundamental principles, basic methodology, strengths and weaknesses, and the established applications of HDX-MS while highlighting new developments and applications.

Structure-function analyses of the G729R 2-oxoadipate dehydrogenase genetic variant associated with a disorder of l-lysine metabolism

2-Oxoadipate dehydrogenase (E1a, also known as DHTKD1, dehydrogenase E1, and transketolase domain-containing protein 1) is a thiamin diphosphate-dependent enzyme and part of the 2-oxoadipate dehydrogenase complex (OADHc) in l-lysine catabolism. Genetic findings have linked mutations in the DHTKD1 gene to several metabolic disorders. These include α-aminoadipic and α-ketoadipic aciduria (AMOXAD), a rare disorder of l-lysine, l-hydroxylysine, and l-tryptophan catabolism, associated with clinical presentations such as developmental delay, mild-to-severe intellectual disability, ataxia, epilepsy, and behavioral disorders that cannot currently be managed by available treatments. A heterozygous missense mutation, c.2185G→A (p.G729R), in DHTKD1 has been identified in most AMOXAD cases. Here, we report that the G729R E1a variant when assembled into OADHc in vitro displays a 50-fold decrease in catalytic efficiency for NADH production and a significantly reduced rate of glutaryl-CoA production by dihydrolipoamide succinyl-transferase (E2o). However, the G729R E1a substitution did not affect any of the three side-reactions associated solely with G729R E1a, prompting us to determine the structure-function effects of this mutation. A multipronged systematic analysis of the reaction rates in the OADHc pathway, supplemented with results from chemical cross-linking and hydrogen-deuterium exchange MS, revealed that the c.2185G→A DHTKD1 mutation affects E1a-E2o assembly, leading to impaired channeling of OADHc intermediates. Cross-linking between the C-terminal region of both E1a and G729R E1a with the E2o lipoyl and core domains suggested that correct positioning of the C-terminal E1a region is essential for the intermediate channeling. These findings may inform the development of interventions to counter the effects of pathogenic DHTKD1 mutations.

Poisson Model To Generate Isotope Distribution for Biomolecules

We introduce a simplified computational algorithm for computing isotope distributions (relative abundances and masses) of biomolecules. The algorithm is based on Poisson approximation to binomial and multinomial distributions. It leads to a small number of arithmetic operations to compute isotope distributions of molecules. The approach uses three embedded loops to compute the isotope distributions, as compared with the eight embedded loops in exact calculations. The speed improvement is about 3-fold compared to the fast Fourier transformation-based isotope calculations, often termed as ultrafast isotope calculation. The approach naturally incorporates the determination of the masses of each molecular isotopomer. It is applicable to high mass accuracy and resolution mass spectrometry data. The application to tryptic peptides in a UniProt protein database revealed that the mass accuracy of the computed isotopomers is better than 1 ppm. Even better mass accuracy (below 1 ppm) is achievable when the method is paired with the exact calculations, which we term a hybrid approach. The algorithms have been implemented in a freely available C/C++ code.

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.

Parsing disease-relevant protein modifications from epiphenomena: perspective on the structural basis of SOD1-mediated ALS

Conformational change and modification of proteins are involved in many cellular functions. However, they can also have adverse effects that are implicated in numerous diseases. How structural change promotes disease is generally not well-understood. This perspective illustrates how mass spectrometry (MS), followed by toxicological and epidemiological validation, can discover disease-relevant structural changes and therapeutic strategies. We (with our collaborators) set out to characterize the structural and toxic consequences of disease-associated mutations and post-translational modifications (PTMs) of the cytosolic antioxidant protein Cu/Zn-superoxide dismutase (SOD1). Previous genetic studies discovered >180 different mutations in the SOD1 gene that caused familial (inherited) amyotrophic lateral sclerosis (fALS). Using hydrogen-deuterium exchange with mass spectrometry, we determined that diverse disease-associated SOD1 mutations cause a common structural defect - perturbation of the SOD1 electrostatic loop. X-ray crystallographic studies had demonstrated that this leads to protein aggregation through a specific interaction between the electrostatic loop and an exposed beta-barrel edge strand. Using epidemiology methods, we then determined that decreased SOD1 stability and increased protein aggregation are powerful risk factors for fALS progression, with a combined hazard ratio > 300 (for comparison, a lifetime of smoking is associated with a hazard ratio of ~15 for lung cancer). The resulting structural model of fALS etiology supported the hypothesis that some sporadic ALS (sALS, ~80% of ALS is not associated with a gene defect) could be caused by post-translational protein modification of wild-type SOD1. We developed immunocapture antibodies and high sensitivity top-down MS methods and characterized PTMs of wild-type SOD1 using human tissue samples. Using global hydrogen-deuterium exchange, X-ray crystallography and neurotoxicology, we then characterized toxic and protective subsets of SOD1 PTMs. To cap this perspective, we present proof-of-concept that post-translational modification can cause disease. We show that numerous mutations (N➔D; Q➔E), which result in the same chemical structure as the PTM deamidation, cause multiple diseases. Copyright © 2017 John Wiley & Sons, Ltd.