Proton Transfer Kinetics in Histidine Side Chains Determined by pH-Dependent Multi-Nuclear NMR Relaxation.

Histidine is a key amino-acid residue in proteins with unique properties engendered by its imidazole side chain that can exist in three different states: two different neutral tautomeric forms and a protonated, positively charged one with a pKa value close to physiological pH. Commonly, two or all three states coexist and interchange rapidly, enabling histidine to act as both donor and acceptor of hydrogen bonds, coordinate metal ions, and engage in acid/base catalysis. Understanding the exchange dynamics among the three states is critical for assessing histidine's mechanistic role in catalysis, where the rate of proton exchange and interconversion among tautomers might be rate limiting for turnover. Here, we determine the exchange kinetics of histidine residues with pKa values representative of the accessible range from 5 to 9 by measuring pH-dependent 15N, 13C, and 1H transverse relaxation rate constants for 5 nuclei in each imidazole. Proton exchange between the imidazole and the solvent is mediated by hydronium ions at acidic and neutral pH, whereas hydroxide mediated exchange becomes the dominant mechanism at basic pH. Proton transfer is very fast and reaches the diffusion limit for pKa values near neutral pH. We identify a direct pathway between the two tautomeric forms, likely mediated by a bridging water molecule or, in the case of high pH, hydroxide ion. For histidines with pKa 7, we determine all rate constants (lifetimes) involving protonation over the entire pH range. Our approach should enable critical insights into enzymatic acid/base catalyzed reactions involving histidines in proteins.

Energetics and dynamics of the proton shuttle of carbonic anhydrase II.

Human carbonic anhydrase II catalyzes the reversible reaction of carbon dioxide and water to form bicarbonate and a proton. His64-mediated proton shuttling between the active site and the bulk solvent is rate limiting. Here we investigate the protonation behavior of His64 as well as its structural and dynamic features in a pH dependent way. We derive two pKa values for His64, 6.25 and 7.60, that we were able to assign to its inward and outward conformation. Furthermore, we show that His64 exists in both conformations equally, independent of pH. Both conformations display an equal distribution of their two neutral tautomeric states. The life time of each conformation is short and both states display high flexibility within their orientation. Therefore, His64 is never static, but rather poised to change conformation. These findings support an energetic, dynamic and solution ensemble-based framework for the high enzymatic activity of human carbonic anhydrase II.

Impact of distant peptide substrate residues on enzymatic activity of SlyD.

Peptidyl-prolyl isomerases (PPIases) catalyze intrinsically slow and often rate-limiting isomerization of prolyl-peptide bonds in unfolded or partially folded proteins, thereby speeding up the folding process and preventing misfolding. They often possess binding and chaperone domains in addition to the domain carrying the isomerization activity. Although generally, their substrates display no identity in their amino acid sequence upstream and downstream of the proline with 20 possibilities for each residue, PPIases are efficient enzymes. SlyD is a highly efficient PPIase consisting of an isomerase domain and an additional chaperone domain. The binding of peptide substrates to SlyD and its enzymatic activity depend to some extend on the proline-proximal residues, however, the impact of proline-distant residues has not been investigated so far. Here, we introduce a label-free NMR-based method to measure SlyD activity on different peptide substrates and analysed the data in the context of obtained binding affinities and several co-crystal structures. We show that especially charged and aromatic residues up to eight positions downstream and three positions upstream of the proline and outside the canonical region of similar conformations affect the activity and binding, although they rarely display distinct conformations in our crystal structures. We hypothesize that these positions primarily influence the association reaction. In the absence of the chaperone domain the isomerase activity strongly correlates with substrate affinity, whereas additional factors play a role in its presence. The mutual orientation of isomerase and chaperone domains depends on the presence of substrates in both binding sites, implying allosteric regulation of enzymatic activity.

Ultrastructural evidence for self-replication of Alzheimer-associated Aß42 amyloid along the sides of fibrils.

The nucleation of Alzheimer-associated Aß peptide monomers can be catalyzed by preexisting Aß fibrils. This leads to autocatalytic amplification of aggregate mass and underlies self-replication and generation of toxic oligomers associated with several neurodegenerative diseases. However, the nature of the interactions between the monomeric species and the fibrils during this key process, and indeed the ultrastructural localization of the interaction sites have remained elusive. Here we used NMR and optical spectroscopy to identify conditions that enable the capture of transient species during the aggregation and secondary nucleation of the Aß42 peptide. Cryo-electron microscopy (cryo-EM) images show that new aggregates protrude from the entire length of the progenitor fibril. These protrusions are morphologically distinct from the well-ordered fibrils dominating at the end of the aggregation process. The data provide direct evidence that self-replication through secondary nucleation occurs along the sides of fibrils, which become heavily decorated under the current solution conditions (14 µM Aß42, 20 mM sodium phosphate, 200 µM EDTA, pH 6.8).

Monitoring protein unfolding transitions by NMR-spectroscopy.

NMR-spectroscopy has certain unique advantages for recording unfolding transitions of proteins compared e.g. to optical methods. It enables per-residue monitoring and separate detection of the folded and unfolded state as well as possible equilibrium intermediates. This allows a detailed view on the state and cooperativity of folding of the protein of interest and the correct interpretation of subsequent experiments. Here we summarize in detail practical and theoretical aspects of such experiments. Certain pitfalls can be avoided, and meaningful simplification can be made during the analysis. Especially a good understanding of the NMR exchange regime and relaxation properties of the system of interest is beneficial. We show by a global analysis of signals of the folded and unfolded state of GB1 how accurate values of unfolding can be extracted and what limits different NMR detection and unfolding methods. E.g. commonly used exchangeable amides can lead to a systematic under determination of the thermodynamic protein stability. We give several perspectives of how to deal with more complex proteins and how the knowledge about protein stability at residue resolution helps to understand protein properties under crowding conditions, during phase separation and under high pressure.

1H R1ρ relaxation dispersion experiments in aromatic side chains.

Aromatic side chains are attractive probes of protein dynamic, since they are often key residues in enzyme active sites and protein binding sites. Dynamic processes on microsecond to millisecond timescales can be studied by relaxation dispersion experiments that attenuate conformational exchange contributions to the transverse relaxation rate by varying the refocusing frequency of applied radio-frequency fields implemented as either CPMG pulse trains or continuous spin-lock periods. Here we present an aromatic 1H R1ρ relaxation dispersion experiment enabling studies of two to three times faster exchange processes than achievable by existing experiments for aromatic side chains. We show that site-specific isotope labeling schemes generating isolated 1H-13C spin pairs with vicinal 2H-12C moieties are necessary to avoid anomalous relaxation dispersion profiles caused by Hartmann-Hahn matching due to the 3JHH couplings and limited chemical shift differences among 1H spins in phenylalanine, tyrosine and the six-ring moiety of tryptophan. This labeling pattern is sufficient in that remote protons do not cause additional complications. We validated the approach by measuring ring-flip kinetics in the small protein GB1. The determined rate constants, kflip, agree well with previous results from 13C R1ρ relaxation dispersion experiments, and yield 1H chemical shift differences between the two sides of the ring in good agreement with values measured under slow-exchange conditions. The aromatic1H R1ρ relaxation dispersion experiment in combination with the site-selective 1H-13C/2H-12C labeling scheme enable measurement of exchange rates up to kex = 2kflip = 80,000 s-1, and serve as a useful complement to previously developed 13C-based methods.

Slow ring flips in aromatic cluster of GB1 studied by aromatic 13C relaxation dispersion methods.

Ring flips of phenylalanine and tyrosine are a hallmark of protein dynamics. They report on transient breathing motions of proteins. In addition, flip rates also depend on stabilizing interactions in the ground state, like aromatic stacking or cation-π interaction. So far, experimental studies of ring flips have almost exclusively been performed on aromatic rings without stabilizing interactions. Here we investigate ring flip dynamics of Phe and Tyr in the aromatic cluster in GB1. We found that all four residues of the cluster, Y3, F30, Y45 and F52, display slow ring flips. Interestingly, F52, the central residue of the cluster, which makes aromatic contacts with all three others, is flipping significantly faster, while the other rings are flipping with the same rates within margin of error. Determined activation enthalpies and activation volumes of these processes are in the same range of other reported ring flips of single aromatic rings. There is no correlation of the number of aromatic stacking interactions to the activation enthalpy, and no correlation of the ring's extent of burying to the activation volume. Because of these findings, we speculate that F52 is undergoing concerted ring flips with each of the other rings.

Site-selective 1H/2H labeling enables artifact-free 1H CPMG relaxation dispersion experiments in aromatic side chains.

Aromatic side chains are often key residues in enzyme active sites and protein binding sites, making them attractive probes of protein dynamics on the millisecond timescale. Such dynamic processes can be studied by aromatic 13C or 1H CPMG relaxation dispersion experiments. Aromatic 1H CPMG relaxation dispersion experiments in phenylalanine, tyrosine and the six-ring moiety of tryptophan, however, are affected by 3J 1H-1H couplings which are causing anomalous relaxation dispersion profiles. Here we show that this problem can be addressed by site-selective 1H/2H labeling of the aromatic side chains and that artifact-free relaxation dispersion profiles can be acquired. The method has been further validated by measuring folding-unfolding kinetics of the small protein GB1. The determined rate constants and populations agree well with previous results from 13C CPMG relaxation dispersion experiments. Furthermore, the CPMG-derived chemical shift differences between the folded and unfolded states are in excellent agreement with those obtained directly from the spectra. In summary, site-selective 1H/2H labeling enables artifact-free aromatic 1H CPMG relaxation dispersion experiments in phenylalanine and the six-ring moiety of tryptophan, thereby extending the available methods for studying millisecond dynamics in aromatic protein side chains.

Experimental pKa Value Determination of All Ionizable Groups of a Hyperstable Protein.

Electrostatic interactions significantly contribute to the stability and function of proteins. The stabilizing or destabilizing effect of local charge is reflected in the perturbation of the pKa value of an ionizable group from the intrinsic pKa value. Herein, the charge network of a hyperstable dimeric protein (ribbon-helix-helix (rhh) protein from plasmid pRN1 from Sulfolobus islandicus) is studied through experimental determination of the pKa values of all ionizable groups. Transitions were monitored by multiple NMR signals per ionizable group between pH 0 and 12.5, prior to a global analysis, which accounted for the effects of neighboring residues. It is found that for several residues involved in salt bridges (four Asp and one Lys) the pKa values are shifted in favor of the charged state. Furthermore, the pKa values of residues C40 and Y47, both located in the hydrophobic dimer interface, are shifted beyond 13.7. The necessary energy for such a shift is about two-thirds of the total stability of the protein, which confirms the importance of the hydrophobic core to the overall stability of the rhh protein.

Conserved S/T Residues of the Human Chaperone DNAJB6 Are Required for Effective Inhibition of Aß42 Amyloid Fibril Formation.

The human molecular chaperone DNAJB6, an oligomeric protein with a conserved S/T-rich region, is an efficient suppressor of amyloid fibril formation by highly aggregation-prone peptides such as the Aß and polyQ peptides associated with Alzheimer's and Huntington's disease, respectively. We previously showed that DNAJB6 can inhibit the processes through which amyloid fibrils are formed via strong interactions with aggregated forms of Aß42 that become sequestered. Here we report that the concentration-dependent capability of DNAJB6 to suppress fibril formation in thioflavin T fluorescence assays decreases progressively with an increasing number of S/T substitutions, with an almost complete loss of suppression when 18 S/T residues are substituted. The kinetics of primary nucleation in particular are affected. No detectable changes in the structure are caused by the substitutions. Also, the level of binding of DNAJB6 to Aß42 decreases with the S/T substitutions, as determined by surface plasmon resonance and microscale thermophoresis. The aggregation process monitored using nuclear magnetic resonance spectroscopy showed that DNAJB6, in contrast to a mutational variant with 18 S/T residues substituted, can keep monomeric Aß42 soluble for an extended time. The inhibition of the primary nucleation is likely to depend on hydroxyl groups in side chains of the S/T residues, and hydrogen bonding with Aß42 is one plausible molecular mechanism, although other possibilities cannot be excluded. The loss of the ability to suppress fibril formation upon S/T to A substitution was previously observed also for polyQ peptides, suggesting that the S/T residues in the DNAJB6-like chaperones have a general ability to inhibit amyloid fibril formation by different aggregation-prone peptides.

Conformational exchange of aromatic side chains by 1H CPMG relaxation dispersion.

Aromatic side chains are attractive probes of protein dynamics on the millisecond time scale, because they are often key residues in enzyme active sites and protein binding sites. Further they allow to study specific processes, like histidine tautomerization and ring flips. Till now such processes have been studied by aromatic 13C CPMG relaxation dispersion experiments. Here we investigate the possibility of aromatic 1H CPMG relaxation dispersion experiments as a complementary method. Artifact-free dispersions are possible on uniformly 1H and 13C labeled samples for histidine δ2 and ε1, as well as for tryptophan δ1. The method has been validated by measuring fast folding-unfolding kinetics of the small protein CspB under native conditions. The determined rate constants and populations agree well with previous results from 13C CPMG relaxation dispersion experiments. The CPMG-derived chemical shift differences between the folded and unfolded states are in good agreement with those obtained directly from the spectra. In contrast, the 1H relaxation dispersion profiles in phenylalanine, tyrosine and the six-ring moiety of tryptophan, display anomalous behavior caused by 3J 1H-1H couplings and, if present, strong 13C-13C couplings. Therefore they require site-selective 1H/2H and, in case of strong couplings, 13C/12C labeling. In summary, aromatic 1H CPMG relaxation dispersion experiments work on certain positions (His δ2, His ε1 and Trp δ1) in uniformly labeled samples, while other positions require site-selective isotope labeling.

Ganglioside lipids accelerate α-synuclein amyloid formation.

The deposition of α-synuclein fibrils is one hallmark of Parkinson's disease. Here, we investigate how ganglioside lipids, present in high amounts in neurons and exosomes, influence the aggregation kinetics of α-synuclein. Gangliosides, as well as, other anionic lipid species with small or large headgroups were found to induce conformational changes of α-synuclein monomers and catalyse their aggregation at mildly acidic conditions. Although the extent of this catalytic effect was slightly higher for gangliosides, the results imply that charge interactions are more important than headgroup chemistry in triggering aggregation. In support of this idea, uncharged lipids with large headgroups were not found to induce any conformational change and only weakly catalyse aggregation. Intriguingly, aggregation was also triggered by free ganglioside headgroups, while these caused no conformational change of α-synuclein monomers. Our data reveal that partially folded α-synuclein helical intermediates are not required species in triggering of α-synuclein aggregation.

Structure-Based Insights into the Dynamics and Function of Two-Domain SlpA from Escherichia coli.

SlpA (SlyD-like protein A) comprises two domains, a FK506 binding domain (FKBP fold) of moderate prolyl cis/trans-isomerase activity and an inserted in flap (IF) domain that hosts its chaperone activity. Here we present the nuclear magnetic resonance (NMR) solution structure of apo Escherichia coli SlpA determined by NMR that mirrors the structural properties seen for various SlyD homologues. Crucial structural differences in side-chain orientation arise for F37, which points directly into the hydrophobic core of the active site. It forms a prominent aromatic stacking with F15, one of the key residues for PPIase activity, thus giving a possible explanation for the inherently low PPIase activity of SlpA. The IF domain reveals the highest stability within the FKBP-IF protein family, most likely arising from an aromatic cluster formed by four phenylalanine residues. Both the thermodynamic stability and the PPIase and chaperone activity let us speculate that SlpA is a backup system for homologous bacterial systems under unfavorable conditions.

Dynamics of Aromatic Side Chains in the Active Site of FKBP12.

FKBP12, a small human enzyme, aids protein folding by catalyzing cis-trans isomerization of peptidyl-prolyl bonds, and is involved in cell signaling pathways, calcium regulation, and the immune response. The underlying molecular mechanisms are not fully understood, but it is well-known that aromatic residues in the active site and neighboring loops are important for substrate binding and catalysis. Here we report micro- to millisecond exchange dynamics of aromatic side chains in the active site region of ligand-free FKBP12, involving a minor state population of 0.5% and an exchange rate of 3600 s-1, similar to previous results for the backbone and methyl-bearing side chains. The exchange process involves tautomerization of H87. In the major state H87 is highly flexible and occupies the common HNε2 tautomer, while in the minor state it occupies the rare HNδ1 tautomer, which typically requires stabilization by specific interactions, such as hydrogen bonds. This finding suggests that the exchange process is coupled to a rearrangement of the hydrogen bond network around H87. Upon addition of the active-site inhibitor FK506 the exchange of all aromatic residues is quenched, with exception of H87. The H87 resonances are broadened beyond detection, suggesting that interconversion between tautomers prevail in the FK506-bound state. While key active-site residues undergo conformational exchange in the apo state, the exchange rate is considerably faster than the catalytic turnover, as determined herein by Michaelis-Menten type analysis of NMR line shapes and chemical shifts. We discuss alternative interpretations of this observation in terms of FKBP12 function.

Site-selective 13C labeling of histidine and tryptophan using ribose.

Experimental studies on protein dynamics at atomic resolution by NMR-spectroscopy in solution require isolated 1H-X spin pairs. This is the default scenario in standard 1H-15N backbone experiments. Side chain dynamic experiments, which allow to study specific local processes like proton-transfer, or tautomerization, require isolated 1H-13C sites which must be produced by site-selective 13C labeling. In the most general way this is achieved by using site-selectively 13C-enriched glucose as the carbon source in bacterial expression systems. Here we systematically investigate the use of site-selectively 13C-enriched ribose as a suitable precursor for 13C labeled histidines and tryptophans. The 13C incorporation in nearly all sites of all 20 amino acids was quantified and compared to glucose based labeling. In general the ribose approach results in more selective labeling. 1-13C ribose exclusively labels His δ2 and Trp δ1 in aromatic side chains and helps to resolve possible overlap problems. The incorporation yield is however only 37% in total and 72% compared to yields of 2-13C glucose. A combined approach of 1-13C ribose and 2-13C glucose maximizes 13C incorporation to 75% in total and 150% compared to 2-13C glucose only. Further histidine positions ß, α and CO become significantly labeled at around 50% in total by 3-, 4- or 5-13C ribose. Interestingly backbone CO of Gly, Ala, Cys, Ser, Val, Phe and Tyr are labeled at 40-50% in total with 3-13C ribose, compared to 5% and below for 1-13C and 2-13C glucose. Using ribose instead of glucose as a source for site-selective 13C labeling enables a very selective labeling of certain positions and thereby expanding the toolbox for customized isotope labeling of amino-acids.

Site-selective 13C labeling of proteins using erythrose.

NMR-spectroscopy enables unique experimental studies on protein dynamics at atomic resolution. In order to obtain a full atom view on protein dynamics, and to study specific local processes like ring-flips, proton-transfer, or tautomerization, one has to perform studies on amino-acid side chains. A key requirement for these studies is site-selective labeling with 13C and/or 1H, which is achieved in the most general way by using site-selectively 13C-enriched glucose (1- and 2-13C) as the carbon source in bacterial expression systems. Using this strategy, multiple sites in side chains, including aromatics, become site-selectively labeled and suitable for relaxation studies. Here we systematically investigate the use of site-selectively 13C-enriched erythrose (1-, 2-, 3- and 4-13C) as a suitable precursor for 13C labeled aromatic side chains. We quantify 13C incorporation in nearly all sites in all 20 amino acids and compare the results to glucose based labeling. In general the erythrose approach results in more selective labeling. While there is only a minor gain for phenylalanine and tyrosine side-chains, the 13C incorporation level for tryptophan is at least doubled. Additionally, the Phe ζ and Trp η2 positions become labeled. In the aliphatic side chains, labeling using erythrose yields isolated 13C labels for certain positions, like Ile ß and His ß, making these sites suitable for dynamics studies. Using erythrose instead of glucose as a source for site-selective 13C labeling enables unique or superior labeling for certain positions and is thereby expanding the toolbox for customized isotope labeling of amino-acid side-chains.

Computational design of a leucine-rich repeat protein with a predefined geometry.

Structure-based protein design offers a possibility of optimizing the overall shape of engineered binding scaffolds to match their targets better. We developed a computational approach for the structure-based design of repeat proteins that allows for adjustment of geometrical features like length, curvature, and helical twist. By combining sequence optimization of existing repeats and de novo design of capping structures, we designed leucine-rich repeats (LRRs) from the ribonuclease inhibitor (RI) family that assemble into structures with a predefined geometry. The repeat proteins were built from self-compatible LRRs that are designed to interact to form highly curved and planar assemblies. We validated the geometrical design approach by engineering a ring structure constructed from 10 self-compatible repeats. Protein design can also be used to increase our structural understanding of repeat proteins. We use our design constructs to demonstrate that buried Cys play a central role for stability and folding cooperativity in RI-type LRR proteins. The computational procedure presented here may be used to develop repeat proteins with various geometrical shapes for applications where greater control of the interface geometry is desired.

Molecular insights into substrate recognition and catalytic mechanism of the chaperone and FKBP peptidyl-prolyl isomerase SlyD.

BACKGROUND: Peptidyl-prolyl isomerases (PPIases) catalyze cis/trans isomerization of peptidyl-prolyl bonds, which is often rate-limiting for protein folding. SlyD is a two-domain enzyme containing both a PPIase FK506-binding protein (FKBP) domain and an insert-in-flap (IF) chaperone domain. To date, the interactions of these domains with unfolded proteins have remained rather obscure, with structural information on binding to the FKBP domain being limited to complexes involving various inhibitor compounds or a chemically modified tetrapeptide. RESULTS: We have characterized the binding of 15-residue-long unmodified peptides to SlyD from Thermus thermophilus (TtSlyD) in terms of binding thermodynamics and enzyme kinetics through the use of isothermal titration calorimetry, nuclear magnetic resonance spectroscopy, and site-directed mutagenesis. We show that the affinities and enzymatic activity of TtSlyD towards these peptides are much higher than for the chemically modified tetrapeptides that are typically used for activity measurements on FKBPs. In addition, we present a series of crystal structures of TtSlyD with the inhibitor FK506 bound to the FKBP domain, and with 15-residue-long peptides bound to either one or both domains, which reveals that substrates bind in a highly adaptable fashion to the IF domain through ß-strand augmentation, and can bind to the FKBP domain as both types VIa1 and VIb-like cis-proline ß-turns. Our results furthermore provide important clues to the catalytic mechanism and support the notion of inter-domain cross talk. CONCLUSIONS: We found that 15-residue-long unmodified peptides can serve as better substrate mimics for the IF and FKBP domains than chemically modified tetrapeptides. We furthermore show how such peptides are recognized by each of these domains in TtSlyD, and propose a novel general model for the catalytic mechanism of FKBPs that involves C-terminal rotation around the peptidyl-prolyl bond mediated by stabilization of the twisted transition state in the hydrophobic binding site.

A Microfluidic Platform for Real-Time Detection and Quantification of Protein-Ligand Interactions.

The key steps in cellular signaling and regulatory pathways rely on reversible noncovalent protein-ligand binding, yet the equilibrium parameters for such events remain challenging to characterize and quantify in solution. Here, we demonstrate a microfluidic platform for the detection of protein-ligand interactions with an assay time on the second timescale and without the requirement for immobilization or the presence of a highly viscous matrix. Using this approach, we obtain absolute values for the electrophoretic mobilities characterizing solvated proteins and demonstrate quantitative comparison of results obtained under different solution conditions. We apply this strategy to characterize the interaction between calmodulin and creatine kinase, which we identify as a novel calmodulin target. Moreover, we explore the differential calcium ion dependence of calmodulin ligand-binding affinities, a system at the focal point of calcium-mediated cellular signaling pathways. We further explore the effect of calmodulin on creatine kinase activity and show that it is increased by the interaction between the two proteins. These findings demonstrate the potential of quantitative microfluidic techniques to characterize binding equilibria between biomolecules under native solution conditions.

Site-specific protonation kinetics of acidic side chains in proteins determined by pH-dependent carboxyl (13)C NMR relaxation.

Proton-transfer dynamics plays a critical role in many biochemical processes, such as proton pumping across membranes and enzyme catalysis. The large majority of enzymes utilize acid-base catalysis and proton-transfer mechanisms, where the rates of proton transfer can be rate limiting for the overall reaction. However, measurement of proton-exchange kinetics for individual side-chain carboxyl groups in proteins has been achieved in only a handful of cases, which typically have involved comparative analysis of mutant proteins in the context of reaction network modeling. Here we describe an approach to determine site-specific protonation and deprotonation rate constants (kon and koff, respectively) of carboxyl side chains, based on (13)C NMR relaxation measurements as a function of pH. We validated the method using an extensively studied model system, the B1 domain of protein G, for which we measured rate constants koff in the range (0.1-3) × 10(6) s(-1) and kon in the range (0.6-300) × 10(9) M(-1) s(-1), which correspond to acid-base equilibrium dissociation constants (Ka) in excellent agreement with previous results determined by chemical shift titrations. Our results further reveal a linear free-energy relationship between log kon and pKa, which provides information on the free-energy landscape of the protonation reaction, showing that the variability among residues in these parameters arises primarily from the extent of charge stabilization of the deprotonated state by the protein environment. We find that side-chain carboxyls with extreme values of koff or kon are involved in hydrogen bonding, thus providing a mechanistic explanation for the observed stabilization of the protonated or deprotonated state.