The pKa Values of Acidic and Basic Residues Buried at the Same Internal Location in a Protein Are Governed by Different Factors

The pK_a values of internal ionizable groups are usually very different from the normal pK_a values of ionizable groups in water. To examine the molecular determinants of pK_a values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK_a, whereas Asp38 and Glu38 titrate with elevated pK_a values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK_a of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK_a values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK_a values of internal groups in proteins. This study suggests that improvements to computational methods for pK_a calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.     

Thermodynamic principles for the engineering of pH-driven conformational switches and acid insensitive proteins

The general thermodynamic principles behind pH driven conformational transitions of biological macromolecules are well understood. What is less obvious is how they can be used to engineer pH switches in proteins. The acid unfolding of staphylococcal nuclease (SNase) was used to illustrate different factors that can affect pH-driven conformational transitions. Acid unfolding is a structural transition driven by preferential H⁺ binding to the acid unfolded state (U) over the native (N) state of a protein. It is the result of carboxylic groups that titrate with more normal pK_a values in the U state than in the N state. Acid unfolding profiles of proteins reflect a balance between electrostatic and non-electrostatic contributions to stability. Several strategies were used in attempts to turn SNase into an acid insensitive protein: (1) enhancing global stability of the protein with mutagenesis or with osmolytes, (2) use of high salt concentrations to screen Coulomb interactions, (3) stabilizing the N state through specific anion effects, (4) removing Asp or Glu residues that titrate with depressed pK_a values in the N state, and (5) removing basic residues that might have strong repulsive interactions in the N state at low pH. The only effective way to engineer acid resistance in SNase is not through modulation of pK_a values of Asp/Glu but by enhancing the global stability of the protein. Modulation of pH-driven conformational transitions by selective manipulation of the electrostatic component of the switch is an extremely difficult undertaking.     

pH dependence of conformational fluctuations of the protein backbone

Proton binding equilibria (pK_a values) of ionizable groups in proteins are exquisitely sensitive to their microenvironments. Apparent pK_a values measured for individual ionizable residues with NMR spectroscopy are actually population‐weighted averages of the pK_a in different conformational microstates. NMR spectroscopy experiments with staphylococcal nuclease were used to test the hypothesis that pK_a values of surface Glu and Asp residues are affected by pH‐sensitive fluctuations of the backbone between folded and locally unfolded conformations. ¹⁵N spin relaxation studies showed that as the pH decreases from the neutral into the acidic range the amplitudes of backbone fluctuations in the ps‐ns timescale increase near carboxylic residues. Hydrogen exchange experiments suggested that backbone conformational fluctuations promoted by decreasing pH also reflect slower local or sub‐global unfolding near carboxylic groups. This study has implications for structure‐based pK_a calculations: (1) The timescale of the backbone's response to ionization events in proteins can range from ps to ms, and even longer; (2) pH‐sensitive fluctuations of the backbone can be localized to both the segment the ionizable residue is attached to or the one that occludes the ionizable group; (3) Structural perturbations are not necessarily propagated through Coulomb interactions; instead, local fluctuations appear to be coupled through the co‐operativity inherent to elements of secondary structure and to networks of hydrogen bonds. These results are consistent with the idea that local conformational fluctuations and stabilities are important determinants of apparent pK_a values of ionizable residues in proteins. Proteins 2014; 82:3132–3143. © 2014 Wiley Periodicals, Inc.     

Carboxyl pKa Values and Acid Denaturation of BBL

The protein BBL undergoes structural transitions and acid denaturation between pH 1.2 and 8.0. Using NMR spectroscopy, we measured the pK_a values of all the carboxylic residues in this pH range. We employed ¹³C direct-detection two-dimensional IPAP (in-phase antiphase) CACO NMR spectroscopy to monitor the ionization state of different carboxylic groups and demonstrated its advantages over other NMR techniques in measuring pK_a values of carboxylic residues. The two residues Glu161 and Asp162 had significantly lowered pK_a values, showing that these residues are involved in a network of stabilizing electrostatic interactions, as is His166. The other carboxylates had unperturbed values. The pH dependence of the free energy of denaturation was described quantitatively by the ionizations of those three residues of perturbed pK_a, and, using thermodynamic cycles, we could calculate their pK_as in the native and denatured states as well as the equilibrium constants for denaturation of the different protonation states. We also measured ¹³C^α chemical shifts of individual residues as a function of pH. These shifts sense structural transitions rather than ionizations, and they titrated with pH consistent with the change in equilibrium constant for denaturation. Kinetic measurements of the folding of BBL E161Q indicated that, at pH 7, the stabilizing interactions with Glu161 are formed mainly in the transition state. We also found that local interactions still exist in the acid-denatured state of BBL, which attenuate somewhat the flexibility of the acid-denatured state.     

Roles of carboxyl groups in the transmembrane insertion of peptides

We have used pHLIP® [pH (low) insertion peptide] to study the roles of carboxyl groups in transmembrane (TM) peptide insertion. pHLIP binds to the surface of a lipid bilayer as a disordered peptide at neutral pH; when the pH is lowered, it inserts across the membrane to form a TM helix. Peptide insertion is reversed when the pH is raised above the characteristic pK_a (6.0). A key event that facilitates membrane insertion is the protonation of aspartic acid (Asp) and/or glutamic acid (Glu) residues, since their negatively charged side chains hinder membrane insertion at neutral pH. In order to gain mechanistic understanding, we studied the membrane insertion and exit of a series of pHLIP variants where the four Asp residues were sequentially mutated to nonacidic residues, including histidine (His). Our results show that the presence of His residues does not prevent the pH-dependent peptide membrane insertion at ∼ pH 4 driven by the protonation of carboxyl groups at the inserting end of the peptide. A further pH drop leads to the protonation of His residues in the TM part of the peptide, which induces peptide exit from the bilayer. We also find that the number of ionizable residues that undergo a change in protonation during membrane insertion correlates with the pH-dependent insertion into the lipid bilayer and exit from the lipid bilayer, and that cooperativity increases with their number. We expect that our understanding will be used to improve the targeting of acidic diseased tissue by pHLIP.     

Predicting protein pKa by environment similarity

A statistical method to predict protein pK_a has been developed by using the 3D structure of a protein and a database of 434 experimental protein pK_a values. Each pK_a in the database is associated with a fingerprint that describes the chemical environment around an ionizable residue. A computational tool, MoKaBio, has been developed to identify automatically ionizable residues in a protein, generate fingerprints that describe the chemical environment around such residues, and predict pK_a from the experimental pK_a values in the database by using a similarity metric. The method, which retrieved the pK_a of 429 of the 434 ionizable sites in the database correctly, was crossvalidated by leave‐one‐out and yielded root mean square error (RMSE) = 0.95, a result that is superior to that obtained by using the Null Model (RMSE 1.07) and other well‐established protein pK_a prediction tools. This novel approach is suitable to rationalize protein pK_a by comparing the region around the ionizable site with similar regions whose ionizable site pK_a is known. The pK_a of residues that have a unique environment not represented in the training set cannot be predicted accurately, however, the method offers the advantage of being trainable to increase its predictive power. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Proton transfer in catalysis and the role of proton shuttles in carbonic anhydrase

The undisputed role of His64 in proton transfer during catalysis by carbonic anhydrases in the α class has raised questions concerning the details of its mechanism. The highly conserved residues Tyr7, Asn62, and Asn67 in the active-site cavity function to fine tune the properties of proton transfer by human carbonic anhydrase II (HCA II). For example, hydrophobic residues at these positions favor an inward orientation of His64 and a low pK_a for its imidazole side chain. It appears that the predominant manner in which this fine tuning is achieved in rate constants for proton transfer is through the difference in pK_a between His64 and the zinc-bound solvent molecule. Other properties of the active-site cavity, such as inward and outward conformers of His64, appear associated with the change in ΔpK_a; however, there is no strong evidence to date that the inward and outward orientations of His64 are in themselves requirements for facile proton transfer in carbonic anhydrase.     

The pKa values of two histidine residues in human haemoglobin, the Bohr effect, and the dipole moments of alpha-helices

Studies of abnormal and chemically modified haemoglobins indicate that in 0.1 m-NaCl about 40% of the alkaline Bohr effect of human haemoglobin is contributed by the C-terminal histidine HC3(146)β. In deoxyhaemoglobin, the imidazole of this histidine forms a salt bridge with aspartate FG1(94)β, in oxyhaemoglobin or carbonmonoxyhaemoglobin it accepts a hydrogen bond from its own NH group instead. Kilmartin et al. (1973) showed that in 0.2 m-NaCl + 0.2 m-phosphate this change of ligation lowered the pK_a of the histidine from 8.0 in Hb³ to 7.1 in HbCO, but Russu et al. (1980) claimed that in bis-Tris buffer without added NaCl its pK_a in HbCO dropped no lower than 7.85, and that in this medium the C-terminal histidine made only a negligible contribution to the alkaline Bohr effect.We have compared the histidine resonances of HbCO A with those of three abnormal haemoglobins: HbCO Cowtown (His HC3(146)β → Leu), HbCO Wood (His FG4(97)β → Leu) and HbCO Malmø (His FG4(97)β → Gln). Our results show that the resonance assigned by Russu et al. to His HC3(146)β in fact belongs to His FG4(97)β. Although in Hb the pK_a of His HC3(146)β is 8.05 ± 0.05 independent of ionic strength, in HbCO its pK_a drops sharply with diminishing ionic strength, so that in the buffer employed by Russu et al. it has a pK_a of 6.2 and makes a contribution to the alkaline Bohr effect that is 57% larger than in the phosphate buffer employed by Kilmartin et al. (1973).In HbCO A, His FG4(97)β does not contribute to the Bohr effect, but in HbCO from which His HC3(146)β has been cleaved (HbCO des-His), His FG4(97)β is in equilibrium between two conformations with different pK_a values. This equilibrium varies with ionic strength and pH, and presumably also with degree of ligation of the haem moiety.In HbCO A, His FG4(97)β has a pK_a of 7.8 compared to the pK_a value of about 6.6 characteristic of free histidines at the surface of proteins. This high pK_a is accounted for by its interaction with the negative pole at the C terminus of helices F and FG. It corresponds to a free energy change of the same order as that observed in the interaction of histidines with carboxylate ions and confirms the strongly dipolar character of α-helices, which manifests itself even when they lie on the surface of the protein.     

Ionization constants and thermodynamic parameters of histidine and derivatives

The pK_a values for the proton dissociation of carboxyl, imidazolium, and ammonium groups for histidine and ten of its derivatives were determined electrometrically at seven temperatures in the range 10–40°C. The ΔH and ΔS values were estimated from the temperature dependence of the dissociation constants of histidine and its derivatives. These results and the pK_a values compared in terms of inductive effect suggest an ion-dipole interaction between the protonated amino group and the unprotonated imidazole ring. The charge and the solvation effects of the neighboring groups are the main factors that determine the imidazole group pK_a in histidine and its studied derivatives. The N^τ-H tautomer is favored over the N^π-H by 1.6 kcal/mol, indicating that the inductive substituent effect at position 4 of the imidazole ring is the major component in determining this tautomeric preference.     

Electrostatic contributions to the stability of the GCN4 leucine zipper structure

Ion pairs are ubiquitous in X-ray structures of coiled coils, and mutagenesis of charged residues can result in large stability losses. By contrast, pK_a values determined by NMR in solution often predict only small contributions to stability from charge interactions. To help reconcile these results we used triple-resonance NMR to determine pK_a values for all groups that ionize between pH 1 and 13 in the 33 residue leucine zipper fragment, GCN4p. In addition to the native state we also determined comprehensive pK_a values for two models of the GCN4p denatured state: the protein in 6 M urea, and unfolded peptide fragments of the protein in water. Only residues that form ion pairs in multiple X-ray structures of GCN4p gave large pK_a differences between the native and denatured states. Moreover, electrostatic contributions to stability were not equivalent for oppositely charged partners in ion pairs, suggesting that the interactions between a charge and its environment are as important as those within the ion pair. The pH dependence of protein stability calculated from NMR-derived pK_a values agreed with the stability profile measured from equilibrium urea-unfolding experiments as a function of pH. The stability profile was also reproduced with structure-based continuum electrostatic calculations, although contributions to stability were overestimated at the extremes of pH. We consider potential sources of errors in the calculations, and how pK_a predictions could be improved. Our results show that although hydrophobic packing and hydrogen bonding have dominant roles, electrostatic interactions also make significant contributions to the stability of the coiled coil.     

Assignment of the histidine 12-threonine 45 hydrogen-bonded proton in the NMR spectrum of ribonuclease A in H2O.

Described herein are proton nmr experiments on chemically modified derivatives of ribonuclease A designed to elucidate the origin of an exchangeable resonance, assigned previously to a histidine ring N proton that titrates between 11 to 13 ppm with a pK_a of 6.1 in H₂O solution. Histidines 48 and 105, which are distant from the active site, are eliminated as candidates for this resonance from inhibitor binding studies on the enzyme in acetate–water solutions. This exchangeable resonance titrates with modified pK_a's and constant area over the above pH range in His-119-N₁-carboxymethylated-RNase A and des-(121–124)-RNase A, thus eliminating the imidazole N₃ proton in the His 119-Asp 121 hydrogen bond. In His-12-N₁-carboxymethylated-RNase A, this resonance is also observable, but broadens on raising the pH above 7 and at elevated temperatures above neutrality. It exhibits a pH-independent chemical shift characteristic of the protonated state of histidine. On the basis of these findings, this exchangeable resonance, designated a, is assigned to the imidazole N₁ proton of His 12, which is hydrogen-bonded to the carbonyl oxygen of Thr 45 in the crystal.     

Electrostatic analysis of the hepatitis C virus NS3 helicase reveals both active and allosteric site locations

Multi-conformation continuum electrostatics (MCCE) was used to analyze various structures of the NS3 RNA helicase from the hepatitis C virus in order to determine the ionization state of amino acid side chains and their pK_as. In MCCE analyses of HCV helicase structures that lacked ligands, several active site residues were identified to have perturbed pK_as in both the nucleic acid binding site and in the distant ATP-binding site, which regulates helicase movement. In all HCV helicase structures, Glu493 was unusually basic and His369 was abnormally acidic. Both these residues are part of the HCV helicase nucleic acid binding site, and their roles were analyzed by examining the pH profiles of site-directed mutants. Data support the accuracy of MCCE predicted pK_a values, and reveal that Glu493 is critical for low pH enzyme activation. Several key residues, which were previously shown to be involved in helicase-catalyzed ATP hydrolysis, were also identified to have perturbed pK_as including Lys210 in the Walker-A motif and the DExD/H-box motif residues Asp290 and His293. When DNA was present in the structure, the calculated pK_as shifted for both Lys210 and Asp290, demonstrating how DNA binding might lead to electrostatic changes that stimulate ATP hydrolysis.     

Structural origins of high apparent dielectric constants experienced by ionizable groups in the hydrophobic core of a protein

The side chains of Lys66, Asp66, and Glu66 in staphylococcal nuclease are fully buried and surrounded mainly by hydrophobic matter, except for internal water molecules associated with carboxylic oxygen atoms. These ionizable side chains titrate with pK_a values of 5.7, 8.8, and 8.9, respectively. To reproduce these pK_a values with continuum electrostatics calculations, we treated the protein with high dielectric constants. We have examined the structural origins of these high apparent dielectric constants by using NMR spectroscopy to characterize the structural response to the ionization of these internal side chains. Substitution of Val66 with Lys66 and Asp66 led to increased conformational fluctuations of the microenvironments surrounding these groups, even under pH conditions where Lys66 and Asp66 are neutral. When Lys66, Asp66, and Glu66 are charged, the proteins remain almost fully folded, but resonances for a few backbone amides adjacent to the internal ionizable residues are broadened. This suggests that the ionization of the internal groups promotes a local increase in dynamics on the intermediate timescale, consistent with either partial unfolding or increased backbone fluctuations of helix 1 near residue 66, or, less likely, with increased fluctuations of the charged side chains at position 66. These experiments confirm that the high apparent dielectric constants reported by internal Lys66, Asp66, and Glu66 reflect localized changes in conformational fluctuations without incurring detectable global structural reorganization. To improve structure-based pK_a calculations in proteins, we will need to learn how to treat this coupling between ionization of internal groups and local changes in conformational fluctuations explicitly.     

Modulation of an Active-Site Cysteine pKa Allows PDI to Act as a Catalyst of both Disulfide Bond Formation and Isomerization

Protein disulfide isomerase (PDI) plays a central role in disulfide bond formation in the endoplasmic reticulum. It is implicated both in disulfide bond formation and in disulfide bond reduction and isomerization. To be an efficient catalyst of all three reactions requires complex mechanisms. These include mechanisms to modulate the pK_a values of the active-site cysteines of PDI. Here, we examined the role of arginine 120 in modulating the pK_a values of these cysteines. We find that arginine 120 plays a significant role in modulating the pK_a of the C-terminal active-site cysteine in the a domain of PDI and plays a role in determining the reactivity of the N-terminal active-site cysteine but not via direct modulation of its pK_a. Mutation of arginine 120 and the corresponding residue, arginine 461, in the a′ domain severely reduces the ability of PDI to catalyze disulfide bond formation and reduction but enhances the ability to catalyze disulfide bond isomerization due to the formation of more stable PDI-substrate mixed disulfides. These results suggest that the modulation of pK_a of the C-terminal active cysteine by the movement of the side chain of these arginine residues into the active-site locales has evolved to allow PDI to efficiently catalyze both oxidation and isomerization reactions.     

pH-dependent random coil 1H, 13C,and 15N chemical shifts of the ionizable amino acids: a guide for protein pK a measurements

The pK _a values and charge states of ionizable residues in polypeptides and proteins are frequently determined via NMR-monitored pH titrations. To aid the interpretation of the resulting titration data, we have measured the pH-dependent chemical shifts of nearly all the ¹H, ¹³C, and ¹⁵N nuclei in the seven common ionizable amino acids (X = Asp, Glu, His, Cys, Tyr, Lys, and Arg) within the context of a blocked tripeptide, acetyl-Gly-X-Gly-amide. Alanine amide and N-acetyl alanine were used as models of the N- and C-termini, respectively. Together, this study provides an essentially complete set of pH-dependent intra-residue and nearest-neighbor reference chemical shifts to help guide protein pK _a measurements. These data should also facilitate pH-dependent corrections in algorithms used to predict the chemical shifts of random coil polypeptides. In parallel, deuterium isotope shifts for the side chain ¹⁵N nuclei of His, Lys, and Arg in their positively-charged and neutral states were also measured. Along with previously published results for Asp, Glu, Cys, and Tyr, these deuterium isotope shifts can provide complementary experimental evidence for defining the ionization states of protein residues.     

Molecular Asymmetry in Pigeonpea Urease: pH Inactivation Studies

Pigeonpea (Cajanus cajan) urease was inactivated by incubating it in buffer of low pHs i.e., 4.8 and 4.5. The pattern of inactivation at both pHs was found to be biphasic, in which half of the activity was destroyed more rapidly than the remaining half. This distribution of active site into two categories is suggestive of site-site heterogeneity, or more specifically, the half-site reactivity of the enzyme moiety. Our pH studies on the rate of reaction showed the presence of two ionizable groups of pK _a values 6.2 ± 0.1 and 8.8 ± 0.1, respectively (Srivastava PK & Kayastha AM, J Mol Catal B: Enz, 16 (2001) 81-89). The later group corresponds to the pK _a value of cysteine group. Here we correlate the loss of urease activity by low pH treatment is due to the effect on essential thiol residues.     

Origin of the pKa shift of the catalytic lysine in acetoacetate decarboxylase

The pK_a value of Lys115, the catalytic residue in acetoacetate decarboxylate, was calculated using atomic coordinates of the X-ray crystal structure with consideration of the protonation states of all titratable sites in the protein. The calculated pK_a value of Lys115 (pK_a(Lys115)) was unusually low (≈6) in agreement with the experimentally measured value. Although charged residues impact pK_a(Lys115) considerably in the native protein, the significant pK_a(Lys115) downshift in the protein with respect to aqueous solution was mainly due to loss of the solvation energy in the catalytic active site relative to bulk water.     

Influence of p<Emphasis Type="Italic">K</Emphasis><Subscript>a</Subscript> Shifts on the Calculated Dipole Moments of Proteins

The protein dipole moment is a low-resolution parameter that characterizes the second-order charge organization of a biomolecule. Theoretical approaches to calculate protein dipole moments rely on pK _a values, which are either computed individually for each ionizable residue or obtained from model compounds. The influence of pK _a shifts are evaluated first by comparing calculated and measured dipole moments of β-lactoglobulin. Second, calculations are made on a dataset of 66 proteins from the Protein Data Bank, and average differences are determined between dipole moments calculated with model pK _as, pK _as derived using a Poisson–Boltzmann approach, and empirically-calculated pK _as. Dipole moment predictions that neglect pK _a shifts are consistently larger than predictions in which they are included. The importance of pK _a shifts are observed to vary with protein size, internal permittivity, and solution pH.     

Chemically accurate protein structures: Validation of protein NMR structures by comparison of measured and predicted pK a values

A new method is presented for evaluating the quality of protein structures obtained by NMR. This method exploits the dependence between measurable chemical properties of a protein, namely pK _a values of acidic residues, and protein structure. The accurate and fast empirical computational method employed by the PROPKA program () allows the user to test the ability of a given structure to reproduce known pK _a values, which in turn can be used as a criterion for the selection of more accurate structures. We demonstrate the feasibility of this novel idea for a series of proteins for which both␣NMR and X-ray structures, as well as pK _a values of all ionizable residues, have been determined. For the 17 NMR ensembles used in this study, this criterion is shown effective in the elimination of a large number of NMR structure ensemble members.