Method to estimate amide base-catalyzed back exchange rates in H2O/D2O mixtures

Intrinsic Hydrogen Deuterium Exchange Rates in H2O/D2O Mixtures | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Intrinsic Hydrogen Deuterium Exchange Rates in H 2 O/D 2 O Mixtures View ORCID Profile Antonio Grimaldi , View ORCID Profile Michele Stofella , View ORCID Profile Emanuele Paci doi: https://doi.org/10.1101/2025.09.11.674334 Antonio Grimaldi 1 Department of Physics and Astronomy, University of Bologna Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Antonio Grimaldi For correspondence: a.grimaldi{at}unibo.it e.paci{at}unibo.it Michele Stofella 1 Department of Physics and Astronomy, University of Bologna Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Michele Stofella Emanuele Paci 1 Department of Physics and Astronomy, University of Bologna Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Emanuele Paci For correspondence: a.grimaldi{at}unibo.it e.paci{at}unibo.it Abstract Full Text Info/History Metrics Data/Code Preview PDF Abstract Hydrogen–deuterium exchange (HDX) measurements are widely used to probe protein structural dynamics. Quantitative interpretation of HDX data relies on the concept of an intrinsic exchange rate, which is well characterized in isotopically pure H 2 O or D 2 O but does not explicitly account for the back exchange that necessarily occurs in H 2 O/D 2 O mixtures: in this case, both the approach-to-equilibrium rate and the equilibrium deuterium enrichment of amides depend nontrivially on solvent composition and acidity. A practical method is presented to predict intrinsic forward and reverse amide exchange rates in H 2 O/D 2 O mixtures. The approach combines known second-order reference rates measured in pure solvents with established empirical descriptions of H 2 O/D 2 O mixtures. The resulting framework yields explicit expressions for forward and back exchange rates as functions of solvent composition and acidity, and correctly recovers the known limits in pure H 2 O and pure D 2 O. The model predicts composition-dependent kinetic isotope effects and an equilibrium amide fractionation factor of ϕ = 1.20 for unstructured peptides in base-catalyzed conditions, in close agreement with the experimental value 1.22 reported for poly-D,L-alanine. By providing a physically motivated description of exchange in mixed solvents, this method offers a practical starting point for quantitatively correcting back exchange in HDX–MS and HDX–NMR experiments. Introduction Biophysical techniques measuring amide hydrogen-deuterium exchange (HDX) are broadly used to fingerprint structural dynamics of proteins. 1 – 4 The Linderstrøm-Lang model describes HDX of protein backbone amides according to the reaction 5 where the amide can assume either a closed (NH cl ) or open (NH op ) conformation, the former being exchange-incompetent. An amide in the open state exchanges with a pseudo-first-order rate k int , called the intrinsic (chemical) exchange rate, which can be estimated as a function of sequence, temperature and pH (or pD) in pure H 2 O (or D 2 O). 6 – 8 The protection factor P = k cl /k op is the equilibrium constant for site closing 3 that one aims to determine to characterize structural dynamics. 9 – 14 In EX2 conditions ( k cl ≫ k op , k int ), the protection factor is related to the observed exchange rate k obs as P = k int /k obs . HDX is typically detected using nuclear magnetic resonance (NMR) spectroscopy or mass spectrometry (MS). NMR experiments and the labeling step in HDX-MS are commonly conducted in buffers containing 80-95% D 2 O. 15 While mixed-solvent equilibrium and kinetic isotope effects are well established, 16 , 17 they are often treated approximately or neglected in practical analyses, despite a body of work that used HDX in H 2 O/D 2 O mixtures to quantify amide fractionation and probe hydrogen-bonded structure. 18 – 21 Mixed H 2 O/D 2 O solvents are also deliberately used in NMR protocols that quantify amide–solvent exchange, including CLEANEX-PM 22 and SOLEXSY. 23 In H 2 O/D 2 O mixtures, exchange is intrinsically bidirectional and cannot be described by a single rate constant. n this work, we make this explicit by decomposing intrinsic exchange into forward and reverse components, enabling a quantitative treatment of kinetic and equilibrium isotope effects in mixed solvents that is not accessible within standard intrinsic-rate models. This work presents a practical method to predict forward and reverse amide exchange rates in H 2 O/D 2 O mixtures as functions of the degree of deuteration and acidity. The approach combines published reference parameters measured in pure solvents with an operational acidity scale for mixtures, and a probabilistic treatment for reprotonation inserting H or D. The resulting expressions recover pure-solvent limits and provide explicit mixture-dependent rates. Methods Proton transfer theory The chemistry of HDX reactions is described by proton transfer theory: 3 , 24 A proton donor XH collides with a proton acceptor Y, forming a short-lived encounter complex (bracketed in Eq. 2 ) in which proton redistribution is faster than dissociation. Here, XH and Y are generic acid/base labels and may be neutral or charged depending on the specific mechanism, the corresponding conjugate-acid pairs being XH/X and YH/Y. Complex formation occurs with a second-order (diffusion-limited) rate constant k d estimated as 10 10 M −1 s −1 . 25 The complex is regarded as highly dynamical but always weakly populated. Focusing on the forward process XH + Y →X + HY, the proton transfer rate constant k t can be expressed as 3 where Δp K a = p K a (XH) − p K a (YH) is the difference between the acidity of the donor XH and that of the conjugate acid YH of the acceptor. The acidity constant of a generic acid HA is defined as with K a (HA) being the acid dissociation constant. Intrinsic exchange rates In 100% D 2 O, amide HDX involves removal of a proton (H + ) from the amide group and the transfer of a deuteron (D + ) from bulk solvent to the amide. The reaction can be acid-, base-, or water-catalyzed, and the pseudo-first-order HDX intrinsic exchange rate of Eq. 1 can be written as where k acid , k base and k water are the corresponding second-order rate constants. Temperature and pH dependence of HDX second-order rate constants were measured by Englander’s group for poly-D,L-alanine (PDLA) and 3-alanine (3-Ala). 6 – 8 Below, we dissect the base-catalyzed mechanism of HDX, which dominates by orders of magnitude at near-neutral conditions. 3 Extension to other mechanisms (acid- and water-catalyzed) is reported in the Appendix. For base-catalyzed HDX in pure D 2 O, the kinetic bottleneck is the first step, i.e. abstraction of the amide proton by OD − . The amide is rapidly deuterated by a solvent D 2 O molecule, and a new OD − ion is produced. In this case, the intrinsic rate can be estimated as B λ and B ρ are tabulated factors that depend on amino acid type and position relative to the amide group ( λ for left, ρ for right, see Table A1 ). 6 These account for steric and inductive effects from neighboring residues, and were demonstrated to be simply additive. 26 k B,ref ( T ) is the reference rate (i.e. the exchange rate for PDLA or 3-Ala) at temperature T , which exhibits an Arrhenius-type dependence: where E B = 17 kcal mol −1 , R is the gas constant, and k B,ref ( T ref ) is the second-order rate constant measured. Values for k B,ref (20°C) for H and D in H 2 O and H in D 2 O are given in Table 1 . 6 – 8 View this table: View inline View popup Download powerpoint Table 1: Reference rates (PDLA) for the base-catalyzed exchange of amides in disordered peptides, at 20°C. 8 Reference rates are specified by two pedices, the first indicating the amide-bound isotope being removed, the second indicating the isotope in the catalytic base abstracting the proton, cfr Forward and back exchange rates Section. Rates for 3-Ala are obtained multiplying the reference rates by a factor 1.35. 8 H 2 O/D 2 O mixtures This section introduces an operational acidity scale for mixed H 2 O/D 2 O solvents and then analyzes the influence of solvent composition on hydroxide availability and reprotonation probabilities that dictate intrinsic exchange rates. Let L denote either protium (H) or deuterium (D) isotope. An isotopic exchange reaction in a solvent with D 2 O mole fraction x (H 2 O mole fraction 1 − x ) involving a solute XL is Here, XH and XD are the solute protiated and deuterated species, while HOL and DOL are solvent molecules containing at least one H or D atom, respectively. The equilibrium constant ϕ XL of the reaction (8) is called the fractionation factor of XL which quantifies the enrichment in D of solute XL with respect to the solvent. 16 The concentrations [H + ] and [OH − ] in pure H 2 O follow from the pH and the ionic product, , through Where . Analogous relations hold in pure D 2 O defining [D + ] and [OD − ] in terms of pD, pOD and . In H 2 O/D 2 O mixtures, multiple isotopologue ionization equilibria coexist (involving H 2 O, HDO, D 2 O and the corresponding ions). Rather than assigning separate equilibrium constants to each elementary process, it is convenient to work with an operational acidity scale based on the total concentration of [L + ] = [H + ] + [D + ], introducing the pL as used in studies of mixed H 2 O/D 2 O solvents by glass electrodes. 27 Similarly one can define the pOL via where [OL − ] = [OH − ] + [OD − ]. At 25°C, the ionic product of D 2 O is lower than that of H 2 O . 28 , 29 In mixed solvents, the corresponding effective ionic product pK w ( x ) = pL( x ) +pOL( x ) varies nonlinearly with the deuterium mole fraction x , and can be described empirically as 30 To relate the mixture pL( x ) to the glass electrode pH-meter reading, pH ∗ , the following relation can be used: 31 For D 2 O at 25°C and pH ∗ = 7, Eq. 13 yields pD = pL(1) ≃ 7.43, i.e., the traditional “+0.4” correction widely used by HDX practitioners. The temperature dependence of has been empirically determined as 32 At given composition x , the p K ( x ) can be computed via Eq. 12 . Once the mixture pH ∗ is measured, the pL( x ) can be computed using Eq. 13 and the pOL( x ) is given by subtraction of the two values as stated above. Reported values for the fractionation factor ϕ OL (0.47 31 or 0.43 33 ) can be used to compute [OH − ]( x ) and [OD − ]( x ), cfr Eq. A4 . Examples are shown in the upper and middle panels of Figure 1 . Download figure Open in new tab Figure 1: Acidity, concentration of ions and intrinsic rates in mixtures, expressed as a function of the deuterium fraction in the mixture, in three different scenarios: constant pH read (left), constant pL (center), constant pOL (right). Upper panels: mixture pH read, pL and pOL. Middle panels: concentrations of the ions relevant for base-catalyzed HDX: [OH − ], [OD − ], and [OL − ] = [OH − ] + [OD − ]. Lower panels: intrinsic forward k forw , back k back and mixture intrinsic k int ( x ) = k forw ( x ) + k back ( x ) exchange rates of PDLA. Magnitude and direction of kinetic isotope effects depend on the quantity fixed while exploring different H 2 O/D 2 O compositions. Data for 25°C. Forward and back exchange rates The reaction considered for HDX of an open amide in a H 2 O/D 2 O mixture is where and are the base-catalyzed reactions rate constants. Reference rates in pure solvents ( Table 1 ) are denoted k HH , k DH , k HD , k DD respectively. Here, the first index indicates the isotope initially bound to the amide (NH or ND) and the second indicates the isotope of the catalytic base (OH − or OD − ) performing proton abstraction. These correspond to the base-catalyzed rate removal of the amide-bound proton, followed by rapid reprotonation from the bulk solvent. The unmeasured rate k DD was estimated as that is, assuming that reaction rates are governed entirely by the zero-point energies of the participating species. 34 The rates of reaction (15) are obtained as a sum of terms, cfr Eq. 5 , each having the same structure as in Eq. 6 : The coefficients of the products ( B λ × B ρ ) are tabulated, cfr Appendix. The scheme (15) can be simplified to a pseudo-first order reaction where the pseudo-first order forward and back exchange rates are weighted by the probability of encountering a reactive solvent molecule: that is, upon removal of the proton, in the forward reaction reprotonation occurs by DOL with probability x = [DOL] / [L 2 O], and by HOL with probability 1 − x = [HOL] / [L 2 O] in the reverse. The solution to the kinetics in Eq. 18 , with the constraint d ([NH] + [ND]) / d t = 0 and calling D ( t ) the normalized concentration of deuterated amides over time, is where D 0 is the initial condition, is the fraction of deuterated amides at equilibrium, and the intrinsic rate in the mixture is defined as Results and discussion The rates k forw and k back are generally different and depend on the fraction x of D 2 O in the solvent. The stationary state ( Eq. 21 ), as well as the approach-to-equilibrium rate ( Eq. 22 ), are functions of x through k forw ( x ) and k back ( x ), cfr Eq. 19 . In other words, kinetic and equilibrium isotope effects are present. Kinetic isotope effects Kinetic isotope effects cause the reaction rate k int ( x ) to vary as a function of the solvent composition. These are due not only to the different second order rates of Table 1 , but also to the solvent composition x , which influences the mixture ion product p K w ( x ) and acidity pL( x ), cfr Eqs. 12 and 13 . By construction, the measured reference rates are recovered in pure H 2 O and D 2 O: by Eqs. 17 and 19 , for x = 0 (pure H 2 O), while for x = 1 (pure D 2 O), The question “how does k int ( x ) vary with x ?” has no unique answer: it depends on the acidity measure held constant while varying x . At fixed temperature, the ionic product of the mixture p K w ( x ) is assumed to depend on x only, cfr Eq. 12 . The observed exchange rate depends on the acidity of the solution through the concentration of ions OH − and OD − , cfr Eq. 17 . However, different definitions of acidity are possible in a mixture, e.g. based on pH ∗ , pL or pOL. The behavior of k int ( x ) in these scenarios is illustrated in the lower panels of Figure 1 and commented below. These results extend the observations made for pure solvents 7 to mixtures. In the first column of Figure 1 , mixtures of different compositions that yield same pH ∗ are considered. Starting from neutral pure H 2 O at 25°C (pH = 7), for x > 0 one finds ( Eq. 12 ). For varying composition, the solution remains neutral, thus both pL and pOL increase. A higher pOL implies fewer available OL − ions for base-catalyzed HDX, hence k int ( x ) decreases for increasing x . As a result, the intrinsic rate of H → D in pure D 2 O is about 2-fold slower than D → H in pure H 2 O. In the second example, the effective acidity of the mixture, pL, was fixed to 7.43 (because pL(1) = 7.43 at pH ∗ = 7 and T = 25°C). Here, decreasing x causes p K w ( x ) to decrease. Because pL is fixed, pOL decreases too, and a more substantial kinetic isotope effect is observed. D → H in H 2 O results about 5-fold faster than H → D in D 2 O on a pL scale. The third scenario evaluates the effect of x at fixed pOL. In this case, k int ( x ) increases with x . This is because the concentration of catalysts is fixed hence the HDX rates depends on the rates of Table 1 . Since k HH > k DH and k HD > k DD , that is, it is easier to extract H than D from the amide group, the HDX rate is higher for increasing D 2 O content. In this case, the extent of the isotope effect (moving from pure D 2 O to H 2 O) is smaller than the other cases (1.5 times faster in D 2 O). As an example, Figure 2 displays the approach-to-equilibrium rates k int ( x ) ( Eq. 22 ) for two short sequences at fixed pH ∗ , which is the experimentally controlled parameter, and varying D 2 O content in the mixture. Kinetic isotope effects lead to a decrease of the reaction rate constant when the D 2 O content is increased, as discussed above. Download figure Open in new tab Figure 2: Intrinsic exchange rates k int ( x ) for short polypeptide sequences ‘ AAAAAAA ’ (left) and ‘ PEPTIDE ’ (right), calculated for varying D 2 O fraction in the mixture x, T = 25°C, pH ∗ = 7, and using PDLA reference rates. Equilibrium isotope effects An equilibrium isotope effect is described by the fractionation factor ( Eq. 9 ), which for reaction (15) is simply . Within the assumptions of the base-catalyzed model developed above, the fractionation factor reduces to with k HH and k DH from Table 1 . This value is unaffected by neighboring side chains, because the local sequence multiplicative factors ( B λ × B ρ ) cancel in the ratio, cfr Appendix. A corollary of Eq. 23 is that, for random coil peptides in conditions where base catalysis dominates, the equilibrium deuterium enrichment of amides will exceed the bulk D 2 O content, i.e. D eq ( x ) > x . This is shown in Figure 3 . A physical basis of this difference is given by the difference in the p K a between amide deuterated (ND) and protiated (NH) forms: p K a (ND) > p K a (NH), thus D tends to be bound more tightly to the amide, and the amide is enriched with D at equilibrium. 7 The value ϕ = 1.20, as well as the predicted deuterium enrichment, is insensitive to temperature and pH (provided that the base-catalyzed mechanism is the dominant one). Download figure Open in new tab Figure 3: Equilibrium isotope effects. At equilibrium, the amide is enriched in D with respect to solvent ( D eq ( x ) > x , as ϕ = 1.20 > 1). The theoretical value ϕ = 1.20 ( Eq. 23 ) agrees closely with the measured one for PDLA ( ϕ PDLA = 1.22). 19 The larger fractionation factor reported for poly-D,L-lysine (PDLK) within the same study 19 ( ϕ PDLK = 1.46) may be attributed to a residual helical propensity of PDLA. 35 It is thus possible that the equilibrium fractionation expected for unstructured peptides is underestimated, since structured amides generally display lower fractionation (usually, ϕ ≲ 1) than unstructured ones. 18 , 19 Residual uncertainties remain also because acid- and water-catalyzed exchange, solvent–side-chain interactions, and buffer composition can all shift the apparent fractionation. Notably, the equilibrium fractionation factor emerges directly from the kinetic framework without introducing additional adjustable parameters or sequence-specific assumptions. This indicates that the observed fractionation is an inherent consequence of intrinsic exchange chemistry in mixed solvents rather than an independent thermodynamic input. Conclusions This work provides a practical method to estimate intrinsic exchange rates in H 2 O/D 2 O mixtures as a function of solvent composition and acidity, recovering known limits for pure solvents. 6 – 8 The model yields closed expressions for forward and back exchange rates, that for base-catalyzed exchange incorporate mixture-dependent catalyst availability and probabilistic insertion of H or D. For completeness, the corresponding acid- and water-catalyzed contributions are described in the Appendix. Kinetic and equilibrium isotope effects emerge as natural consequences of the chemistry of the mixture and are not determined by composition alone: both the magnitude and the direction of the kinetic isotope effect depend on how acidity is specified or experimentally controlled (e.g. fixed pH ∗ versus fixed operational acidity). Under base-catalyzed conditions, the framework predicts equilibrium amide enrichment in D, in good agreement with measurements performed on PDLA ( ϕ ≃ 1.20). 19 Mismatch with PDLK ( ϕ ≃ 1.46) 19 suggests that residual structure or other unmodeled factors can affect fractionation and should be assessed by systematic benchmarking against random-coil standards. Beyond providing explicit intrinsic exchange rates in mixtures, this framework emphasizes that HDX is bidirectional and that should be described by both forward and back exchange rates to capture kinetic and equilibrium isotope effects. The same framework can be applied to HDX-NMR data, where mixture-dependent kinetics and final enrichment must be accounted for to obtain unbiased estimates of local stability (in terms of protection factors). 36 This is also relevant to HDX-MS workflows, where back exchange can occur both during labeling (due to incomplete buffer deuteration) and in subsequent handling steps (e.g. quenching, usually performed at intermediate D 2 O content, low temperature and pH). 15 A consistent treatment of exchange in mixture is the first step towards physically-grounded back exchange corrections. In this case, the chemistry of the mixture defines baseline effects, that must be integrated with the environment (e.g. salts, co-solvents, denaturants, and whatever additive that can alter exchange kinetics). Code availability A Python 37 implementation of the method presented, allowing to compute k forw and k back for arbitrary sequences as a function of solvent temperature, glass electrode pH read, and composition, is available at https://github.com/pacilab/hdx-rates-mixtures . Acknowledgement AG was supported by the Italian Ministry of University and Research (MUR) under the National Recovery and Resilience Plan (PNRR), through a PhD scholarship at the University of Bologna (DM 118/2023, CUP J33C23002540002). Appendix Multiplicative factors for intrinsic exchange rates View this table: View inline View popup Download powerpoint Table A1: Multiplicative factors for intrinsic exchange rates. 6 – 8 Coefficients A λ , A ρ , B λ , B ρ are given in Table A1 in both H 2 O and D 2 O for all residues except Asp ( D ), Glu ( E ) and His ( H ). 8 In these cases, the coefficients A ζ ,L , B ζ ,L ( X ) of residue X , where ζ is λ or ρ , and L is H (for D → H exchange) or D (H → D), are given by where and parameters C L ( X ), E a,L ( X ), T ref given in Table A2 . View this table: View inline View popup Download powerpoint Table A2: Parameters for the computation of B λ and B ρ coefficients for Asp, Glu, His. 8 Calculation of ϕ From definitions of Eqs. A1 and A2 , one has that ( B λ × B ρ ) HH = ( B λ × B ρ ) DH and ( B λ × B ρ ) DD = ( B λ × B ρ ) HD . Calling ξ = ( B λ × B ρ ) HD / ( B λ × B ρ ) HH ( ξ = 1 if Asp, Glu or His are not neighboring to the amide group), the rates and of Eq. 17 are written [OH − ]( x ) and [OD − ]( x ) can be written in terms of [OL − ]( x ) using the definition of fractional abundance of OD − in OL − : 16 Hence Eq. A3 becomes where the subscript “HH” of ( B λ × B ρ ) has been dropped to ease the notation. Recalling Eq. 16 , and defining ψ = k HD /k HH = k DD /k DH , from which one has Extension to other mechanisms In H 2 O/D 2 O mixtures (D 2 O mole fraction x ), the intrinsic exchange rates k forw ( x ), k back ( x ) are given as a sum of contributions (of acid-, base-, and water-catalyzed exchange): In the main text, the derivation focused on the base-catalyzed contribution, which dominates near-neutral conditions and therefore governs most HDX measurements. For completeness, this appendix outlines how the other two pathways can be incorporated into the same framework as presented in details in the main text for the case of the base-catalyzed pathway. View this table: View inline View popup Download powerpoint Table A3: Reference rates for acid- and water-catalyzed exchange of at 20°C. 8 Water-catalyzed exchange In the water-catalyzed pathway, a solvent water molecule acts as a base, abstracting the amide H or D. This contribution is essentially pH-independent and usually negligible compared to the other terms, except in the region of the minimum. Sequence effects are addressed by the coefficients of Table A1 , as k W = k W,ref ( B λ × B ρ ), with reference rates given in Table A3 . Literature tabulates the rates of H → D exchange in D 2 O and the rate of D → H in H 2 O . In a mixed solvent with deuterium mole fraction x , it is assumed that the rate of abstraction depends only on the isotope being removed, while the isotope incorporated upon reprotonation is determined solely by the solvent composition. In this case, the forward and back water-catalyzed intrinsic exchange rates are Acid-catalyzed exchange In the acid-catalyzed pathway, exchange proceeds through protonation of the amide by solvated hydronium species before loss of the bound isotope. In this case, the isotopic identity of the proton-donating species, i.e. D + or H + , determines whether the event contributes to H → D or D → H exchange. Acid-catalyzed second-order rates can be calculated from reference values ( Table A3 ) and accounting for the sequence dependence with the corresponding multiplicative factors A λ , A ρ from Table A1 , e.g. . Since the key isotopic selectivity lies in the isotope donated to the amide, it is assumed that the intrinsic rate remains the same as in pure solvent, but the relevant concentrations [H + ]( x ) and [D + ]( x ) should be used: The relative abundance of [D + ] is given, analogously to Eq. A4 , as where ℓ = 0.69 16 is the fractionation parameter for L 3 O + and [L + ]( x ) = [H + ]( x ) + [D + ]( x ), cfr main text. Footnotes Introductory and Method sections revised to avoid confusion between back exchange naturally occurring in a mixture, from back exchange occurring after labeling in HDX-MS workflows. Includes acid- and water-catalysed exchange (Appendix). https://github.com/pacilab/hdx-rates-mixtures References (1). ↵ Masson , G. R. ; Burke , J. E. ; Ahn , N. G. ; Anand , G. S. ; Borchers , C. ; Brier , S. ; Bou-Assaf , G. M. ; Engen , J. R. ; Englander , S. W. ; Faber , J. ; others Recommendations for performing, interpreting and reporting hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments . Nature methods 2019 , 16 , 595 – 602 . OpenUrl PubMed (2). Narang , D. ; Lento , C. ; J. Wilson , D. HDX-MS: an analytical tool to capture protein motion in action . Biomedicines 2020 , 8 , 224 . OpenUrl PubMed (3). ↵ Hamuro , Y. Tutorial: Chemistry of Hydrogen/Deuterium Exchange Mass Spectrometry . Journal of the American Society for Mass Spectrometry 2020 , 32 , 133 – 151 . OpenUrl PubMed (4). ↵ James , E. I. ; Murphree , T. A. ; Vorauer , C. ; Engen , J. R. ; Guttman , M. Advances in hydrogen/deuterium exchange mass spectrometry and the pursuit of challenging biological systems . Chemical reviews 2021 , 122 , 7562 – 7623 . OpenUrl CrossRef PubMed (5). ↵ Hvidt , A. ; Nielsen , S. O. Hydrogen exchange in proteins . Advances in protein chemistry 1966 , 21 , 287 – 386 . OpenUrl CrossRef PubMed (6). ↵ Bai , Y. ; Milne , J. S. ; Mayne , L. ; Englander , S. W. Primary structure effects on peptide group hydrogen exchange . Proteins: Structure, Function, and Bioinformatics 1993 , 17 , 75 – 86 . OpenUrl (7). ↵ Connelly , G. P. ; Bai , Y. ; Jeng , M.-F. ; Englander , S. W. Isotope effects in peptide group hydrogen exchange . Proteins: Structure, Function, and Bioinformatics 1993 , 17 , 87 – 92 . OpenUrl (8). ↵ Nguyen , D. ; Mayne , L. ; Phillips , M. C. ; Walter Englander , S. Reference parameters for protein hydrogen exchange rates . Journal of the American Society for Mass Spectrometry 2018 , 29 , 1936 – 1939 . OpenUrl CrossRef PubMed (9). ↵ Hilser , V. J. ; Freire , E. Structure-based calculation of the equilibrium folding pathway of proteins. Correlation with hydrogen exchange protection factors . Journal of molecular biology 1996 , 262 , 756 – 772 . OpenUrl CrossRef PubMed Web of Science (10). ↵ Vendruscolo , M. ; Paci , E. ; Dobson , C. M. ; Karplus , M. Rare fluctuations of native proteins sampled by equilibrium hydrogen exchange . Journal of the American Chemical Society 2003 , 125 , 15686 – 15687 . OpenUrl CrossRef PubMed Web of Science (11). Best , R. B. ; Vendruscolo , M. Structural interpretation of hydrogen exchange protection factors in proteins: characterization of the native state fluctuations of CI2 . Structure 2006 , 14 , 97 – 106 . OpenUrl CrossRef PubMed (12). Bradshaw , R. T. ; Marinelli , F. ; Faraldo-Gómez , J. D. ; Forrest , L. R. Interpretation of HDX data by maximum-entropy reweighting of simulated structural ensembles . Biophysical Journal 2020 , 118 , 1649 – 1664 . OpenUrl CrossRef PubMed (13). 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 2022 , 30 , 313 – 320 . OpenUrl CrossRef (14). ↵ Devaurs , D. ; Antunes , D. A. ; Borysik , A. J. Computational modeling of molecular structures guided by hydrogen-exchange data . Journal of the American Society for Mass Spectrometry 2022 , 33 , 215 – 237 . OpenUrl CrossRef PubMed (15). ↵ Stofella , M. ; Grimaldi , A. ; Smit , J. H. ; Claesen , J. ; Paci , E. ; Sobott , F. Computational Tools for Hydrogen–Deuterium Exchange Mass Spectrometry Data Analysis . Chemical Reviews 2024 , 124 , 12242 – 12263 . OpenUrl CrossRef PubMed (16). ↵ Gold , V. Advances in Physical Organic Chemistry ; Elsevier , 1969 ; Vol. 7 ; pp 259 – 331 . OpenUrl (17). ↵ Schowen , K. B. ; Schowen , R. L. [29] Solvent isotope effects on enzyme systems . Methods in enzymology 1982 , 87 , 551 – 606 . OpenUrl CrossRef PubMed Web of Science (18). ↵ Loh , S. N. ; Markley , J. L. Hydrogen bonding in proteins as studied by amide hydrogen D/H fractionation factors: application to staphylococcal nuclease . Biochemistry 1994 , 33 , 1029 – 1036 . OpenUrl CrossRef PubMed (19). ↵ Bowers , P. M. ; Klevit , R. E. Hydrogen bonding and equilibrium isotope enrichment in histidine-containing proteins . Nature Structural Biology 1996 , 3 , 522 – 531 . OpenUrl CrossRef PubMed (20). LiWang , A. C. ; Bax , A. Equilibrium protium/deuterium fractionation of backbone amides in U13C/15N labeled human ubiquitin by triple resonance NMR . Journal of the American Chemical Society 1996 , 118 , 12864 – 12865 . OpenUrl (21). ↵ Shi , Z. ; Olson , C. A. ; Kallenbach , N. R. ; Sosnick , T. R. D/H amide isotope effect in model α-helical peptides . Journal of the American Chemical Society 2002 , 124 , 13994 – 13995 . OpenUrl PubMed (22). ↵ Hwang , T.-L. ; Van Zijl , P. C. ; Mori , S. Accurate quantitation of water-amide proton exchange rates using the phase-modulated CLEAN chemical EXchange (CLEANEXPM) approach with a Fast-HSQC (FHSQC) detection scheme . Journal of biomolecular NMR 1998 , 11 , 221 – 226 . OpenUrl CrossRef PubMed Web of Science (23). ↵ Chevelkov , V. ; Xue , Y. ; Krishna Rao , D. ; Forman-Kay , J. D. ; Skrynnikov , N. R. 15NH/D-SOLEXSY experiment for accurate measurement of amide solvent exchange rates: application to denatured drkN SH3 . Journal of Biomolecular NMR 2010 , 46 , 227 – 244 . OpenUrl CrossRef PubMed Web of Science (24). ↵ Eigen , M. Proton transfer, acid-base catalysis, and enzymatic hydrolysis. Part I: elementary processes . Angewandte Chemie International Edition in English 1964 , 3 , 1 – 19 . OpenUrl CrossRef Web of Science (25). ↵ Zhou , G.-Q. ; Zhong , W.-Z. Diffusion-Controlled Reactions of Enzymes: A Comparison between Chou’s Model and Alberty-Hammes-Eigen’s Model . European journal of biochemistry 1982 , 128 , 383 – 387 . OpenUrl CrossRef PubMed (26). ↵ Molday , R. ; Englander , S. ; Kallen , R. Primary structure effects on peptide group hydrogen exchange . Biochemistry 1972 , 11 , 150 – 158 . OpenUrl CrossRef PubMed Web of Science (27). ↵ Mikkelsen , K. ; Nielsen , S. O. Acidity measurements with the glass electrode in H2O-D2O mixtures . The Journal of Physical Chemistry 1960 , 64 , 632 – 637 . OpenUrl CrossRef (28). ↵ Abel , E. ; Bratu , E. ; Redlich , O. Die elektrolytische Dissoziation des schweren Wassers . Zeitschrift für Physikalische Chemie 1935 , 173 , 353 – 364 . OpenUrl (29). ↵ Shoesmith , D. W. ; Lee , W. The ionization constant of heavy water (D2O) in the temperature range 298 to 523 K . Canadian Journal of Chemistry 1976 , 54 , 3553 – 3558 . OpenUrl CrossRef (30). ↵ Gold , V. ; Lowe , B. Measurement of solvent isotope effects with the glass electrode. Part The ionic product of D2O and D2O–H2O mixtures . Journal of the Chemical Society A: Inorganic, Physical, Theoretical 1967 , 936 – 943 . (31). ↵ Salomaa , P. ; Schaleger , L. L. ; Long , F. Solvent deuterium isotope effects on acid-base equilibria . Journal of the American Chemical Society 1964 , 86 , 1 – 7 . OpenUrl CrossRef (32). ↵ Sweeton , F. H. ; Mesmer , R. ; Baes , C. Acidity measurements at elevated temperatures. VII. Dissociation of water . Journal of Solution Chemistry 1974 , 3 , 191 – 214 . OpenUrl (33). ↵ Gold , V. ; Grist , S. Deuterium solvent isotope effects on reactions involving the aqueous hydroxide ion . Journal of the Chemical Society, Perkin Transactions 1972 , 2 , 89 – 95 . OpenUrl (34). ↵ Kiefer , P. M. ; Hynes , J. T. Kinetic isotope effects for adiabatic proton transfer reactions in a polar environment . The Journal of Physical Chemistry A 2003 , 107 , 9022 – 9039 . OpenUrl (35). ↵ Stofella , M. ; Seetaloo , N. ; St John , A. N. ; Paci , E. ; Phillips , J. J. ; Sobott , F. Recalibrating Protection Factors Using Millisecond Hydrogen/Deuterium Exchange Mass Spectrometry . Analytical Chemistry 2025 , 97 , 2648 – 2657 . OpenUrl (36). ↵ Grimaldi , A. ; Stofella , M. ; Hobbs , B. ; Karamanos , T. K. ; Paci , E. Amide Hydrogen Deuterium Exchange in Isotopically Mixed Waters . arXiv preprint arxiv: 2510.24860 2025 , (37). ↵ Van Rossum , G. ; Drake Jr , F. L. Python reference manual ; Centrum voor Wiskunde en Informatica Amsterdam , 1995 . View the discussion thread. Back to top Previous Next Posted January 22, 2026. Download PDF Data/Code Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Intrinsic Hydrogen Deuterium Exchange Rates in H2O/D2O Mixtures Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Intrinsic Hydrogen Deuterium Exchange Rates in H 2 O/D 2 O Mixtures Antonio Grimaldi , Michele Stofella , Emanuele Paci bioRxiv 2025.09.11.674334; doi: https://doi.org/10.1101/2025.09.11.674334 Share This Article: Copy Citation Tools Intrinsic Hydrogen Deuterium Exchange Rates in H 2 O/D 2 O Mixtures Antonio Grimaldi , Michele Stofella , Emanuele Paci bioRxiv 2025.09.11.674334; doi: https://doi.org/10.1101/2025.09.11.674334 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Biophysics Biochemistry Subject Areas All Articles Animal Behavior and Cognition (7620) Biochemistry (17642) Bioengineering (13865) Bioinformatics (41863) Biophysics (21410) Cancer Biology (18548) Cell Biology (25437) Clinical Trials (138) Developmental Biology (13359) Ecology (19865) Epidemiology (2067) Evolutionary Biology (24288) Genetics (15587) Genomics (22468) Immunology (17704) Microbiology (40301) Molecular Biology (17142) Neuroscience (88448) Paleontology (666) Pathology (2825) Pharmacology and Toxicology (4815) Physiology (7634) Plant Biology (15109) Scientific Communication and Education (2042) Synthetic Biology (4285) Systems Biology (9812) Zoology (2268)