Abstract
High-resolution energetic information about protein conformational ensembles is essential for understanding protein function, yet remains challenging to obtain. Here we present PIGEON-FEATHER, a method for calculating ensemble free energies of opening (âGop) at single-amino-acid or near-single-amino-acid resolution for proteins of all sizes from hydrogen exchangeâmass spectrometry (HXâMS) data. PIGEON-FEATHER disambiguates and reconstructs all experimentally measured HXâMS isotopic mass envelopes using a Bayesian Monte Carlo sampling approach. We applied PIGEON-FEATHER to reveal how Escherichia coli and human dihydrofolate reductases (ecDHFR and hDHFR) have evolved distinct ensembles. We show how two competitive inhibitors bind these orthologs differently, solving the longstanding mystery of why both therapeutic molecules inhibit ecDHFR but only one inhibits hDHFR. Extending PIGEON-FEATHER to a large proteinâDNA complex, we mapped ligand-induced ensemble reweighting in the E.âcoli lac repressor to describe the functional switching mechanism crucial for transcriptional regulation.

This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$32.99 /Â 30Â days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout






Similar content being viewed by others
Data availability
Data supporting the main findings are included in the article and its Supplementary Information. All data needed to reproduce this work are available from the PRIDE database84 (dataset identifier PXD057539). Source data are provided with this paper.
Code availability
The source code for the software is available from GitHub (https://github.com/glasgowlab/PIGEON-FEATHER), along with the synthetic datasets used for benchmarking, an example HXâMS dataset and a tutorial.
References
Boehr, D. D., Nussinov, R. & Wright, P. E. The role of dynamic conformational ensembles in biomolecular recognition. Nat. Chem. Biol. 5, 789â796 (2009).
Cooper, A. Thermodynamic fluctuations in protein molecules. Proc. Natl Acad. Sci. USA 73, 2740â2741 (1976).
Hilser, V. J. An ensemble view of allostery. Science 327, 653â654 (2010).
Balbach, J. et al. Following protein folding in real time using NMR spectroscopy. Nat. Struct. Biol. 2, 865â870 (1995).
Roder, H., Elöve, G. A. & Englander, S. W. Structural characterization of folding intermediates in cytochrome c by H-exchange labelling and proton NMR. Nature 335, 700â704 (1988).
Shaw, D. E. et al. Anton, a special-purpose machine for molecular dynamics simulation. Commun. ACM 51, 91â97 (2008).
Kish, M. et al. Online fully automated system for hydrogen/deuterium-exchange mass spectrometry with millisecond time resolution. Anal. Chem. 95, 5000â5008 (2023).
Espada, A. et al. A decoupled automation platform for hydrogen/deuterium exchange mass spectrometry experiments. J. Am. Soc. Mass. Spectrom. 30, 2580â2583 (2019).
Krishna, R. et al. Generalized biomolecular modeling and design with RoseTTAFold All-Atom. Science 384, eadl2528 (2024).
Abramson, J. et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature 630, 493â500 (2024).
Konermann, L., Pan, J. & Liu, Y.-H. Hydrogen exchange mass spectrometry for studying protein structure and dynamics. Chem. Soc. Rev. 40, 1224â1234 (2011).
Bai, Y., Sosnick, T. R., Mayne, L. & Englander, S. W. Protein folding intermediates: native-state hydrogen exchange. Science 269, 192â197 (1995).
Skinner, J. J., Lim, W. K., Bédard, S., Black, B. E. & Englander, S. W. Protein hydrogen exchange: testing current models. Protein Sci. 21, 987â995 (2012).
Bai, Y., Milne, J. S., Mayne, L. & Englander, S. W. Primary structure effects on peptide group hydrogen exchange. Proteins 17, 75â86 (1993).
Masson, G. R. et al. Recommendations for performing, interpreting and reporting hydrogen deuterium exchange mass spectrometry (HDXâMS) experiments. Nat. Methods 16, 595â602 (2019).
Weis, D. D. Recommendations for the propagation of uncertainty in hydrogen exchangeâmass spectrometric measurements. J. Am. Soc. Mass. Spectrom. 32, 1610â1617 (2021).
Hilser, V. J. & Freire, E. Structure-based calculation of the equilibrium folding pathway of proteins. correlation with hydrogen exchange protection factors. J. Mol. Biol. 262, 756â772 (1996).
Liu, T. et al. Quantitative assessment of protein structural models by comparison of H/D exchange MS data with exchange behavior accurately predicted by DXCOREX. J. Am. Soc. Mass. Spectrom. 23, 43â56 (2012).
Eggertson, M. J., Fadgen, K., Engen, J. R. & Wales, T. E. Considerations in the analysis of hydrogen exchange mass spectrometry data. In Mass Spectrometry Data Analysis in Proteomics (ed. Matthiesen, R.)(Springer, 2020).
Kan, Z., Ye, X., Skinner, J. J., Mayne, L. & Englander, S. W. ExMS2: an integrated solution for hydrogenâdeuterium exchange mass spectrometry data analysis. Anal. Chem. 91, 7474â7481 (2019).
Guttman, M., Weis, D. D., Engen, J. R. & Lee, K. K. Analysis of overlapped and noisy hydrogen/deuterium exchange mass spectra. J. Am. Soc. Mass. Spectrom. 24, 1906â1912 (2013).
Saltzberg, D. J. et al. A residue-resolved Bayesian approach to quantitative interpretation of hydrogenâdeuterium exchange from mass spectrometry: application to characterizing proteinâligand interactions. J. Phys. Chem. B 121, 3493â3501 (2017).
Smit, J. H. et al. Probing universal protein dynamics using hydrogenâdeuterium exchange mass spectrometry-derived residue-level Gibbs free energy. Anal. Chem. 93, 12840â12847 (2021).
Stofella, M., Skinner, S. P., Sobott, F., Houwing-Duistermaat, J. & Paci, E. High-resolution hydrogenâdeuterium protection factors from sparse mass spectrometry data validated by nuclear magnetic resonance measurements. J. Am. Soc. Mass. Spectrom. 33, 813â822 (2022).
Kan, Z.-Y., Walters, B. T., Mayne, L. & Englander, S. W. Protein hydrogen exchange at residue resolution by proteolytic fragmentation mass spectrometry analysis. Proc. Natl Acad. Sci. USA 110, 16438â16443 (2013).
Zhang, Z., Zhang, A. & Xiao, G. Improved protein hydrogen/deuterium exchange mass spectrometry platform with fully automated data processing. Anal. Chem. 84, 4942â4949 (2012).
Ferrari, Ã. J. R. et al. Large-scale discovery, analysis, and design of protein energy landscapes. Preprint at bioRxiv https://doi.org/10.1101/2025.03.20.644235 (2025).
Huennekens, F. M. The methotrexate story: a paradigm for development of cancer chemotherapeutic agents. Adv. Enzym. Regul. 34, 397â419 (1994).
Oyeyemi, O. A. et al. Comparative hydrogenâdeuterium exchange for a mesophilic vs thermophilic dihydrofolate reductase at 25â°C: identification of a single active site region with enhanced flexibility in the mesophilic protein. Biochemistry 50, 8251â8260 (2011).
Yamamoto, T. Mass spectrometry on segment-specific hydrogen exchange of dihydrofolate reductase. J. Biochem. 135, 17â24 (2004).
Wales, T. E. et al. Resolving chaperone-assisted protein folding on the ribosome at the peptide level. Nat. Struct. Mol. Biol. 31, 1888â1897 (2024).
Pancsa, R., Varadi, M., Tompa, P. & Vranken, W. F. Start2Fold: a database of hydrogen/deuterium exchange data on protein folding and stability. Nucleic Acids Res. 44, D429âD434 (2016).
Jones, B. E. & Matthews, C. R. Early intermediates in the folding of dihydrofolate reductase from escherichia coli detected by hydrogen exchange and NMR. Protein Sci. 4, 167â177 (1995).
Schnell, J. R., Dyson, H. J. & Wright, P. E. Structure, dynamics, and catalytic function of dihydrofolate reductase. Annu. Rev. Biophys. 33, 119â140 (2004).
Wan, Q. et al. Toward resolving the catalytic mechanism of dihydrofolate reductase using neutron and ultrahigh-resolution X-ray crystallography. Proc. Natl Acad. Sci. USA 111, 18225â18230 (2014).
Vendruscolo, M., Paci, E., Dobson, C. M. & Karplus, M. Rare fluctuations of native proteins sampled by equilibrium hydrogen exchange. J. Am. Chem. Soc. 125, 15686â15687 (2003).
Alford, R. F. et al. The Rosetta All-Atom energy function for macromolecular modeling and design. J. Chem. Theory Comput. 13, 3031â3048 (2017).
Nguyen, T. T., Marzolf, D. R., Seffernick, J. T., Heinze, S. & Lindert, S. Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR. Structure 30, 313â320 (2022).
Hollien, J. & Marqusee, S. Structural distribution of stability in a thermophilic enzyme. Proc. Natl Acad. Sci. USA 96, 13674â13678 (1999).
Liu, C. T. et al. Functional significance of evolving protein sequence in dihydrofolate reductase from bacteria to humans. Proc. Natl Acad. Sci. USA 110, 10159â10164 (2013).
Bhabha, G. et al. Divergent evolution of protein conformational dynamics in dihydrofolate reductase. Nat. Struct. Mol. Biol. 20, 1243â1249 (2013).
Bystroff, C. & Kraut, J. Crystal structure of unliganded Escherichia coli dihydrofolate reductase. Ligand-induced conformational changes and cooperativity in binding. Biochemistry 30, 2227â2239 (1991).
Bhabha, G. et al. A dynamic knockout reveals that conformational fluctuations influence the chemical step of enzyme catalysis. Science 332, 234â238 (2011).
Cameron, C. E. & Benkovic, S. J. Evidence for a functional role of the dynamics of glycine-121 of Escherichia coli dihydrofolate reductase obtained from kinetic analysis of a site-directed mutant. Biochemistry 36, 15792â15800 (1997).
Miller, G. P. & Benkovic, S. J. Strength of an interloop hydrogen bond determines the kinetic pathway in catalysis by Escherichia coli dihydrofolate reductase. Biochemistry 37, 6336â6342 (1998).
Appleman, J. R. et al. Atypical transient state kinetics of recombinant human dihydrofolate reductase produced by hysteretic behavior. Comparison with dihydrofolate reductases from other sources. J. Biol. Chem. 264, 2625â2633 (1989).
Blakley, R. L. Eukaryotic dihydrofolate reductase. Adv. Enzymol. Relat. Areas Mol. Biol. 70, 23â102 (1995).
Carroll, M. J. et al. Evidence for dynamics in proteins as a mechanism for ligand dissociation. Nat. Chem. Biol. 8, 246â252 (2012).
Sawaya, M. R. & Kraut, J. Loop and subdomain movements in the mechanism of Escherichia coli dihydrofolate reductase: crystallographic evidence. Biochemistry 36, 586â603 (1997).
Manna, M. S. et al. A trimethoprim derivative impedes antibiotic resistance evolution. Nat. Commun. 12, 2949 (2021).
Mauldin, R. V., Carroll, M. J. & Lee, A. L. Dynamic dysfunction in dihydrofolate reductase results from antifolate drug binding: modulation of dynamics within a structural state. Structure 17, 386â394 (2009).
Chen, J., Dima, R. I. & Thirumalai, D. Allosteric communication in dihydrofolate reductase: signaling network and pathways for closed to occluded transition and back. J. Mol. Biol. 374, 250â266 (2007).
Krucinska, J. et al. Structure-guided functional studies of plasmid-encoded dihydrofolate reductases reveal a common mechanism of trimethoprim resistance in Gram-negative pathogens. Commun. Biol. 5, 1â14 (2022).
Rodriguez, D. C. P., Weber, K. C., Sundberg, B. & Glasgow, A. MAGPIE: an interactive tool for visualizing and analyzing proteinâligand interactions. Protein Sci. 33, e5027 (2024).
Dion-Schultz, A. & Howell, E. E. Effects of insertions and deletions in a beta-bulge region of Escherichia coli dihydrofolate reductase. Protein Eng. Des. Sel. 10, 263â272 (1997).
Bullerjahn, A. M. & Freisheim, J. H. Site-directed deletion mutants of a carboxyl-terminal region of human dihydrofolate reductase. J. Biol. Chem. 267, 864â870 (1992).
Ratnam, M., Tan, X., Prendergast, N. J., Smith, P. L. & Freisheim, J. H. Ligand-induced structural constraints in human dihydrofolate reductase revealed by peptide-specific antibodies. Biochemistry 27, 4800â4804 (1988).
Howell, E. E., Booth, C., Farnum, M., Kraut, J. & Warren, M. S. A second-site mutation at phenylalanine-137 that increases catalytic efficiency in the mutant aspartate-27 .fwdarw. serine Escherichia coli dihydrofolate reductase. Biochemistry 29, 8561â8569 (1990).
Brown, K. A., Howell, E. E. & Kraut, J. Long-range structural effects in a second-site revertant of a mutant dihydrofolate reductase. Proc. Natl Acad. Sci. USA 90, 11753â11756 (1993).
Dion, A., Linn, C. E., Bradrick, T. D., Georghiou, S. & Howell, E. E. How do mutations at phenylalanine-153 and isoleucine-155 partially suppress the effects of the aspartate-27 .fwdarw. serine mutation in Escherichia coli dihydrofolate reductase?. Biochemistry 32, 3479â3487 (1993).
Watson, M., Liu, J.-W. & Ollis, D. Directed evolution of trimethoprim resistance in Escherichia coli. FEBS J. 274, 2661â2671 (2007).
Toprak, E. Evolutionary paths to strong antibiotic resistance under dynamically sustained drug stress. Nat. Genet. 44, 101â105 (2011).
Daber, R., Stayrook, S., Rosenberg, A. & Lewis, M. Structural analysis of lac repressor bound to allosteric effectors. J. Mol. Biol. 370, 609â619 (2007).
Glasgow, A. et al. Ligand-specific changes in conformational flexibility mediate long-range allostery in the lac repressor. Nat. Commun. 14, 1179 (2023).
Markiewicz, P., Kleina, L. G., Cruz, C., Ehret, S. & Miller, J. H. Genetic studies of the lac repressor. XIV. Analysis of 4000 altered Escherichia coli lac repressors reveals essential and non-essential residues, as well as âspacersâ which do not require a specific sequence. J. Mol. Biol. 240, 421â433 (1994).
Suckow, J. et al. Genetic studies of the lac repressor XV: 4000 single amino acid substitutions and analysis of the resulting phenotypes on the basis of the protein structure. J. Mol. Biol. 261, 509â523 (1996).
Romanuka, J. et al. Genetic switching by the lac repressor is based on two-state MonodâWymanâChangeux allostery. Proc. Natl Acad. Sci. USA 120, e2311240120 (2023).
Kalodimos, C. G., Folkers, G. E., Boelens, R. & Kaptein, R. Strong DNA binding by covalently linked dimeric Lac headpiece: evidence for the crucial role of the hinge helices. Proc. Natl Acad. Sci. USA 98, 6039â6044 (2001).
Zhong, E. D., Bepler, T., Berger, B. & Davis, J. H. CryoDRGN: reconstruction of heterogeneous cryo-EM structures using neural networks. Nat. Methods 18, 176â185 (2021).
Lin, Z. et al. Evolutionary-scale prediction of atomic-level protein structure with a language model. Science 379, 1123â1130 (2023).
van Kempen, M. et al. Fast and accurate protein structure search with Foldseek. Nat. Biotechnol. 42, 243â246 (2024).
Halabi, N., Rivoire, O., Leibler, S. & Ranganathan, R. Protein sectors: evolutionary units of three-dimensional structure. Cell 138, 774â786 (2009).
Steinegger, M. & Söding, J. MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat. Biotechnol. 35, 1026â1028 (2017).
Leaver-Fay, A. et al. ROSETTA3: an object-oriented software suite for the simulation and design of macromolecules. Methods Enzymol. 487, 545â574 (2011).
Jumper, J. et al. Highly accurate protein structure prediction with AlphaFold. Nature 596, 583â589 (2021).
Bell, C. E. & Lewis, M. A closer view of the conformation of the lac repressor bound to operator. Nat. Struct. Biol. 7, 209â214 (2000).
Chambers, M. C. et al. A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 30, 918â920 (2012).
Goloborodko, A. A., Levitsky, L. I., Ivanov, M. V. & Gorshkov, M. V. Pyteomicsâa Python framework for exploratory data analysis and rapid software prototyping in proteomics. J. Am. Soc. Mass. Spectrom. 24, 301â304 (2013).
Levitsky, L. I., Klein, J. A., Ivanov, M. V. & Gorshkov, M. V. Pyteomics 4.0: five years of development of a Python proteomics framework. J. Proteome Res. 18, 709â714 (2019).
Skinner, J. J., Lim, W. K., Bédard, S., Black, B. E. & Englander, S. W. Protein dynamics viewed by hydrogen exchange. Protein Sci. 21, 996â1005 (2012).
Smit, J. Jhsmit/HDXrate: HDXrate version 0.2.2. Zenodo https://doi.org/10.5281/zenodo.10022160 (2023).
Nguyen, D., Mayne, L., Phillips, M. C. & Walter Englander, S. Reference parameters for protein hydrogen exchange rates. J. Am. Soc. Mass. Spectrom. 29, 1936â1939 (2018).
Röst, H. L., Schmitt, U., Aebersold, R. & Malmström, L. pyOpenMS: a Python-based interface to the OpenMS mass-spectrometry algorithm library. Proteomics 14, 74â77 (2014).
Perez-Riverol, Y. et al. The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 50, D543âD552 (2022).
Acknowledgements
We acknowledge A.G. Lab members for requesting, adding, discussing, testing and troubleshooting various features of the software. PIGEON-FEATHER was tested by D. Swingle at the Advanced Science Research Center (ASRC) at the City University of New York. We thank G. Rocklinâs group for providing the EEHEE_rd4_087 sample and NMR dataset and for useful discussions about these data. We thank A. Piserchio, K. Gardner and R. Abzalimov from the ASRC, as well as S. Costello, S. Shoemaker and N. Latorraca, for feedback on the method. We also thank R. Abzalimov in his role as the MS facility manager at the ASRC, where the HXâMS experiments were performed. We are grateful to S. Pratihar, H. Al-Hashimi and A. Palmer at Columbia University Medical Center for valuable discussions and S. Costello, H. Al-Hashimi, R. Abzalimov, R. Azad, J. Glasgow, T. Kortemme and S. Lomvardas for critical feedback on this manuscript. This work was supported by National Institutes of Health grants R00GM135529 and R21EB035208 to A.G. and a National Science Foundation Graduate Research Fellowship to S.K.M.
Author information
Authors and Affiliations
Contributions
C.L., M.L.W. and A.G. developed the computational methods. M.L.W., A.R. and S.K.M. collected the HXâMS data. All authors tested the methods and analyzed and interpreted the data. C.L., M.L.W. and A.G. wrote the manuscript with input from all authors. A.G. conceptualized and supervised the study.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Chemical Biology thanks the anonymous reviewer(s) for their contribution to the peer review of this work.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Peptide disambiguation with PIGEON.
a, HXâMS peptides are frequently degenerate in mass. b, The number of degenerate peptides in E. coli DHFR (ecDHFR) given measurement uncertainty (Îm/z ~0.001 within experiments, ~0.01 between experiments). c, ecDHFR peptides are duplicatively matched to the same peak. d, Disambiguation of peptides in panel c using tandem MS fragment information. The natural isotopic distribution corresponding to a single MS compound is shown, as extracted from the marked portion of the base peak intensity chromatogram. MS2 ions from this compound are compared to theoretical ions for the matched peptides, and the matches scored based on the degree of fragment ion support. e,(i), MS2 data are pooled for all replicates and matched to theoretical peptides for each protein in the dataset. e(ii), Systematic mass errors are corrected by fitting a trendline to the highest scoring matches and re-thresholding around the trendline. All but the best peak for each peptide are discarded. e(iii), Each peptide for which there exists an identical match and which is not adequately supported by MS fragments is dropped. e(iv), Each charge state for which there exists an identical match (as in iii.) but which has fragment support is either kept (âKEEPâ), or discarded (âDROPâ). These peptides co-elute and both contribute to the signal. e(v), The final cleaned dataset contains peptides with scores distributed over the full range. f, Initial score and m/z error distribution for a representative pooled dataset (as in E, i). g, A duplicate match where the fragment ions support one peptide but not the other. h, Impact of each PIGEON step (systematic error correction and disambiguation) and peptide selection mode for four ecDHFR HXâMS datasets (Supplementary Tables 1 and 8). The use of multiple proteases increases peptide coverage: fungal protease (FP), pepsin, nepenthesin (NP), and alanyl aminopeptidase (AP). The m/z and score distribution for one resulting peptide pool (red box) are shown.
Extended Data Fig. 3 FEATHER benchmark against HXâNMR.
a, Absolute ÎGop for EEHEE_rd4_087 projected onto the AlphaFold3-predicted structure. b, Overlay of FEATHER-derived exchange rates and HXâNMR data, formatted and calculated as in Fig. 2e. HXâMS data were pooled for FEATHER analysis from three technical replicate timecourses (see Methods and Supplementary Table 8). c, Comparison of ÎGop for 16 protons from both HXâMS-PIGEON-FEATHER and HXâNMR. R and r.m.s.e are calculated using only site-resolvable residues. Error bars reflect FEATHER bootstrapped uncertainties as in panel b. Points are colored according to peptide coverage using the scale in b.
Supplementary information
Supplementary Information
Review of HX theory, Recommendations for PIGEON-FEATHER usage and Supplementary Figs. 1â31, Tables 1â17, Methods and References.
Supplementary Data
Experimental deuterium uptake plots and model fitting plots.
Source data
Source Data Fig. 2
FEATHER-derived âlog(kex) values versus true values for a representative simulated dataset.
Source Data Fig. 3
FEATHER-derived ÎGop and âlog(kex) values for the apo states of ecDHFR and hDHFR. HXâMS datasets include peptide centroid deuterium uptake and isotopic mass envelopes. Reconstruction errors for both isotopic mass envelopes and centroids are also provided.
Source Data Fig. 4
FEATHER-derived ÎGop and âlog(kex) values for the apo, TMP and MTX states of ecDHFR.
Source Data Fig. 6
FEATHER-derived ÎGop and âlog(kex) values for the apo, DNA and IPTG states of LacI. HXâMS datasets include peptide centroid deuterium uptake and isotopic mass envelopes of LacI.
Source Data Extended Data Fig. 1
Peptide coverage and matching results using PIGEON. Includes statistics on degenerate peptides, illustrated with the ecDHFR dataset as an example.
Source Data Extended Data Fig. 3
FEATHER-derived ÎGop and âlog(kex) values for EEHEE_rd4_087.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Lu, C., Wells, M.L., Reckers, A. et al. Site-resolved energetic information from HXâMS experiments. Nat Chem Biol (2025). https://doi.org/10.1038/s41589-025-02049-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41589-025-02049-1