Statistical analysis of protein structures suggests that buried ionizable residues in proteins are hydrogen bonded or form salt bridges

It is well known that nonpolar residues are largely buried in the interior of proteins, whereas polar and ionizable residues tend to be more localized on the protein surface where they are solvent exposed. Such a distribution of residues between surface and interior is well understood from a thermodynamic point: nonpolar side chains are excluded from the contact with the solvent water, whereas polar and ionizable groups have favorable interactions with the water and thus are preferred at the protein surface. However, there is an increasing amount of information suggesting that polar and ionizable residues do occur in the protein core, including at positions that have no known functional importance. This is inconsistent with the observations that dehydration of polar and in particular ionizable groups is very energetically unfavorable. To resolve this, we performed a detailed analysis of the distribution of fractional burial of polar and ionizable residues using a large set of ?2600 nonhomologous protein structures. We show that when ionizable residues are fully buried, the vast majority of them form hydrogen bonds and/or salt bridges with other polar/ionizable groups. This observation resolves an apparent contradiction: the energetic penalty of dehydration of polar/ionizable groups is paid off by favorable energy of hydrogen bonding and/or salt bridge formation in the protein interior. Our conclusion agrees well with the previous findings based on the continuum models for electrostatic interactions in proteins. Proteins 2011; © 2011 Wiley‐Liss, Inc.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Stabilization of internal charges in a protein: water penetration or conformational change?

The ionizable amino acid side chains of proteins are usually located at the surface. However, in some proteins an ionizable group is embedded in an apolar internal region. Such buried ionizable groups destabilize the protein and may trigger conformational changes in response to pH variations. Because of the prohibitive energetic cost of transferring a charged group from water to an apolar medium, other stabilizing factors must be invoked, such as ionization-induced water penetration or structural changes. To examine the role of water penetration, we have measured the 17O and 2H magnetic relaxation dispersions (MRD) for the V66E and V66K mutants of staphylococcal nuclease, where glutamic acid and lysine residues are buried in predominantly apolar environments. At neutral pH, where these residues are uncharged, we find no evidence of buried water molecules near the mutation site. This contrasts with a previous cryogenic crystal structure of the V66E mutant, but is consistent with the room-temperature crystal structure reported here. MRD measurements at different pH values show that ionization of Glu-66 or Lys-66 is not accompanied by penetration of long-lived water molecules. On the other hand, the MRD data are consistent with a local conformational change in response to ionization of the internal residues.     

Structural reorganization triggered by charging of Lys residues in the hydrophobic interior of a protein

Structural consequences of ionization of residues buried in the hydrophobic interior of proteins were examined systematically in 25 proteins with internal Lys residues. Crystal structures showed that the ionizable groups are buried. NMR spectroscopy showed that in 2 of 25 cases studied, the ionization of an internal Lys unfolded the protein globally. In five cases, the internal charge triggered localized changes in structure and dynamics, and in three cases, it promoted partial or local unfolding. Remarkably, in 15 proteins, the ionization of the internal Lys had no detectable structural consequences. Highly stable proteins appear to be inherently capable of withstanding the presence of charge in their hydrophobic interior, without the need for specialized structural adaptations. The extent of structural reorganization paralleled loosely with global thermodynamic stability, suggesting that structure-based pK(a) calculations for buried residues could be improved by calculation of thermodynamic stability and by enhanced conformational sampling.     

pH dependence of amide chemical shifts in natively disordered polypeptides detects medium-range interactions with ionizable residues

A growing number of natively disordered proteins undergo a folding/binding process that is essential for their biological function. An interesting question is whether these proteins have incompletely solvated regions that drive the folding/binding process. Although the presence of predominantly hydrophobic buried regions can be easily ascertained by high-sensitivity differential scanning calorimetry analysis, the identification of those residues implicated in the burial requires NMR analysis. We have selected a partially solvated natively disordered fragment of Escherichia coli, thioredoxin, C37 (38-108), for full NMR spectral assignment. The secondary chemical shifts, temperature coefficients, and relaxation rates (R(1) and R(2)) of this fragment indicate the presence of a flexible backbone without a stable hydrogen bond network near neutral pH. (1)H-(15)N heteronuclear single quantum coherence analysis of the pH dependence of amide chemical shifts in fragment C37 within pH 2.0 and 7.0 suggests the presence of interactions between nonionizable residues and the carboxylate groups of four Asp and four Glu residues. The pH midpoints (pH(m)) of the amides in the ionizable residues (Asp or Glu) and, consequently, the shifts in the pH(m) (DeltapH(m)) of these residues with respect to model tetrapeptides, are sequence-dependent; and the nonionizable residues that show pH dependence cluster around the ionizable ones. The same pH dependence has been observed in two fragments: M37 (38-73) and C73 (74-108), ruling out the participation of long-range interactions. Our studies indicate the presence of a 15-residue pH-dependent segment with the highest density of ionizable sites in the disordered ensembles of fragments C37 and M37. The observed correlations between ionizable and nonionizable residues in this segment suggest the organization of the backbone and side chains through local and medium-range interactions up to nine residues apart, in contrast to only a few interactions in fragment C73. These results agree qualitatively with the predominantly hydrophobic buried surface detected only in fragments C37 and M37 by highly sensitive differential scanning calorimetry analysis. This work offers a sensitive and rapid new tool to obtain clues about local and nonlocal interactions between ionizable and nonionizable residues in the growing family of natively disordered small proteins with full NMR assignments.     

X-ray and thermodynamic studies of staphylococcal nuclease variants I92E and I92K: insights into polarity of the protein interior

We have used crystallography and thermodynamic analysis to study nuclease variants I92E and I92K, in which an ionizable side-chain is placed in the hydrophobic core of nuclease. We find that the energetic cost of burying ionizable groups is rather modest. The X-ray determinations show water molecules solvating the buried glutamic acid under cryo conditions, but not at room temperature. The lysine side-chain does not appear solvated in either case. Guanidine hydrochloride (GnHCl) denaturation of I92E and I92K, done as a function of pH and monitored by tryptophan fluorescence, showed that I92E and I92K are folded in the pH range pH 3.5-9.0 and pH 5.5-9.5, respectively. The stability of the parental protein is independent of pH over a broad range. In contrast, the stabilities of I92E and I92K exhibit a pH dependence, which is quantitatively explained by thermodynamic analysis: the PK(a) value of the buried K92 is 5.6, while that of the buried E92 is 8.65. The free energy difference between burying the uncharged and charged forms of the groups is modest, about 6 kcal/mol. We also found that epsilon(app) for I92K and I92E is in the range approximately 10-12, instead of 2-4 commonly used to represent the protein interior. Side-chains 92E and 92K were uncharged under the conditions of the X-ray experiment. Both are buried completely inside the well-defined hydrophobic core of the variant proteins without forming salt-bridges or hydrogen bonds to other functional groups of the proteins. Under cryo conditions 92E shows a chain of four water molecules, which hydrate one oxygen atom of the carboxyl group of the glutamic acid. Two other water molecules, which are present in the wild-type at all temperatures, are also connected to the water ring observed inside the hydrophobic core. The ready burial of water with an uncharged E92 raises the possibility that solvent excursions into the interior also take place in the wild-type protein, but in a random, dynamic way not detectable by crystallography. Such transient excursions could increase the average polarity, and thus epsilon(app), of the protein interior.     

Very fast empirical prediction and rationalization of protein pKa values

A very fast empirical method is presented for structure-based protein pKa prediction and rationalization. The desolvation effects and intra-protein interactions, which cause variations in pKa values of protein ionizable groups, are empirically related to the positions and chemical nature of the groups proximate to the pKa sites. A computer program is written to automatically predict pKa values based on these empirical relationships within a couple of seconds. Unusual pKa values at buried active sites, which are among the most interesting protein pKa values, are predicted very well with the empirical method. A test on 233 carboxyl, 12 cysteine, 45 histidine, and 24 lysine pKa values in various proteins shows a root-mean-square deviation (RMSD) of 0.89 from experimental values. Removal of the 29 pKa values that are upper or lower limits results in an RMSD = 0.79 for the remaining 285 pKa values.     

Electrostatic contributions to the stability of hyperthermophilic proteins.

Electrostatic contributions to the folding free energy of several hyperthermophilic proteins and their mesophilic homologs are calculated. In all the cases studied, electrostatic interactions are more favorable in the hyperthermophilic proteins. The electrostatic free energy is found not to be correlated with the number of ionizable amino acid residues, ion pairs or ion pair networks in a protein, but rather depends on the location of these groups within the protein structure. Moreover, due to the large free energy cost associated with burying charged groups, buried ion pairs are found to be destabilizing unless they undergo favorable interactions with additional polar groups, including other ion pairs. The latter case involves the formation of stabilizing ion pair networks as is observed in a number of proteins. Ion pairs located on the protein surface also provide stabilizing interactions in a number of cases. Taken together, our results suggest that many hyperthermophilic proteins enhance electrostatic interactions through the optimum placement of charged amino acid residues within the protein structure, although different design strategies are used in different cases. Other physical mechanisms are also likely to contribute, however optimizing electrostatic interactions offers a simple means of enhancing stability without disrupting the core residues characteristic of different protein families.     

Protein Ionizable Groups: pK Values and Their Contribution to Protein Stability and Solubility

The structure, stability, solubility, and function of proteins depend on their net charge and on the ionization state of the individual residues. Consequently, biochemists are interested in the pK values of the ionizable groups in proteins and how these pK values depend on their environment. We review what has been learned about pK values of ionizable groups in proteins from experimental studies and discuss the important contributions they make to protein stability and solubility.     

How many ionizable groups can sit on a protein hydrophobic core?

Full or partial burial of ionizable groups in the hydrophobic interior of proteins underlies the large modulation in group properties (modified pK value, high nucleophilicity, enhanced capability of interaction with chemical moieties of the substrate, etc.) linked to biological function. Indeed, the few internal ionizable residues found in proteins are known to play important functional roles in catalysis and, in general, in energy transduction processes. However, ionizable‐group burial is expected to be seriously disruptive and, it is important to note, most functional sites contain not just one, but several ionizable residues. Hence, the adaptations involved in the development of function in proteins (through in vitro engineering or during the course of natural evolution) are not fully understood. Here, we explore experimentally how proteins respond to the accumulation of hydrophobic‐to‐ionizable residue substitutions. For this purpose, we have constructed a combinatorial library targeting a hydrophobic cluster in a consensus‐engineered, stabilized form of a small model protein. Contrary to naïve expectation, half of the variants randomly selected from the library are soluble, folded, and active, despite including up to four mutations. Furthermore, for these variants, the dependence of stability with the number of mutations is not synergistic and catastrophic, but smooth and approximately linear. Clearly, stabilized protein scaffolds may be robust enough to withstand many disruptive hydrophobic‐to‐ionizable residue mutations, even when they are introduced in the same region of the structure. These results should be relevant for protein engineering and may have implications for the understanding of the early evolution of enzymes. Proteins 2012; © 2011 Wiley Periodicals, Inc.     

A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules

Structure and function of macromolecules depend critically on the ionization states of their acidic and basic groups. Most current structure-based theoretical methods that predict pK of ionizable groups in macromolecules include, as one of the key steps, a computation of the partition sum (Boltzmann average) over all possible protonation microstates. As the number of these microstates depends exponentially on the number of ionizable groups present in the molecule, direct computation of the sum is not realistically feasible for many typical proteins that may have tens or even hundreds of ionizable groups. We have tested a simple and robust approximate algorithm for computing these partition sums for macromolecules. The method subdivides the interacting sites into independent clusters, based upon the strength of site-site electrostatic interaction. The resulting partition function is factorizable into computationally manageable components. Two variants of the approach are presented and validated on a representative test set of 602 proteins, by comparing the pK(1/2) values computed by the proposed method with those obtained by the standard Monte Carlo approach used as a reference. With 95% confidence, the relative error introduced by the more accurate of the two methods is less than 0.25 pK units. The algorithms are one to two orders of magnitude faster than the Monte Carlo method, with the typical settings. A graphical representation is introduced that visualizes the clusters of strong site-site interactions in the context of the three-dimensional (3D) structure of the macromolecule, facilitating identification of functionally important clusters of ionizable groups; the approach is exemplified on two proteins, bacteriorhodopsin and myoglobin.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

Multiple-site titration and molecular modeling: Two rapid methods for computing energies and forces for ionizable groups in proteins

Computer models of proteins frequently treat the energies and forces associated with ionizable groups as if they were purely electrostatic. This paper examines the validity of the purely electrostatic approach, and concludes that significant errors in energies can result from the neglect of ionization changes. However, a complete treatment of ionizable groups presents substantial computational obstacles, because of the large number of ionization states which must be examined in systems having multiple interacting titratable groups. In order to address this problem, two novel methods for treating the energetics and forces associated with ionizable groups with a minimum of computer time have been developed. The most rapid method yields approximate energies by computing the free energy of a single highly occupied ionization state. The second method separates ionizable groups into clusters, and treats intracluster interactions exactly, but intercluster interactions approximately. This method yields both accurate energies and fractional charges. Good results are obtained in tests of both methods on proteins having has many as 123 ionizable groups. The more rapid method requires computer times of 0.01 to 0.34 sec, while the more accurate method requires 0.7 to 15 sec. These methods may be fast enough to permit the incorporation of ionization effects in iterative computations, such as energy minimizations and conformational searches. © 1993 Wiley-Liss, Inc.     

Chemical modification of mitochondria. Uncoupler binding by mitochondria in different metabolic states

At low uncoupler concentrations the binding of carbonyl-cyanide-m-chlorophenyl-hydrazone to mitochondria was found to depend sensitively on the metabolic state of mitochondria. The binding data are consistent with the assumption that at low concentrations and pH 7.4 the uncoupler is bound mainly in anionic form to the inner mitochondrial membrane and that upon energization the inner membrane undergoes conformation change, exposes buried ionizable groups and hence acquires a negative net membrane charge. Deenergization of the inner membrane by a small amount of uncoupler removes the negative net membrane charge and consequently increases the apparent binding constants. Based upon the present results on uncoupler binding and previous observations on the physiological properties of alkylating uncouplers, a possible molecular mechanism involving electron carriers and coupling factors is suggested for coupling electron transport to phosphorylation.     

Radioiodination of proteins by reductive alkylation

The use of the aliphatic aldehyde, para-hydroxyphenylacetaldehyde as the reactive moiety in the radioiodination of proteins by reductive alkylation is described. The para-hydroxyphenyl group is radiolabeled with 125I, reacted through its aliphatic aldehyde group with primary amino groups on proteins to form a reversible Schiff base linkage which can then be stabilized with the mild reducing agent NaCNBH3. The introduction of the methylene group between the benzene ring and the aldehyde group increases its reactivity with protein amino groups permitting efficient labeling at low aldehyde concentrations. Using this method, radioiodinated proteins with high specific activity can be produced. The reductive alkylation procedure is advantageous in that the labeling conditions are mild, the reaction is specific for lysyl residues, and the modification of the epsilon-ammonium group of lysine results in ionizable secondary amino groups avoiding major changes in protein charge.     

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.     

Completely buried, non-ion-paired glutamic acid contributes favorably to the conformational stability of pyrrolidone carboxyl peptidases from hyperthermophiles

Pyrrolidone carboxyl peptidases (PCPs) from hyperthermophiles have a structurally conserved and completely buried Glu192 in the hydrophobic core; in contrast, the corresponding residue in the mesophile protein is a hydrophobic residue, Ile. Does the buried ionizable residue contribute to stabilization or destabilization of hyperthermophile PCPs? To elucidate the role of the buried glutamic acid in stabilizing PCP from hyperthermophiles, we constructed five Glu192 mutants of PCP-0SH (C142S/C188S, Cys-free double mutant of PCP) from Pyrococcus furiosus and examined their thermal and pH-induced unfolding and crystal structures and compared them with those of PCP-0SH. The stabilities of apolar (E192A/I/V) and polar (E192D/Q) mutants were less than PCP-0SH at acidic pH values. In the alkaline region, the mutant proteins, except for E192D, were more stable than PCP-0SH. The thermal stability data and theoretical calculations indicated an apparent pKa value > or = 7.3 for Glu192. Present results confirmed that the protonated Glu192 in PCP-0SH forms strong hydrogen bonds with the carbonyl oxygen and peptide nitrogen of Pro168. New intermolecular hydrogen bonds in the E --> A/D mutants were formed by a water molecule introduced into the cavity created around position 192, whereas the hydrogen bonds disappeared in the E --> I/V mutants. Structure-based empirical stability of mutant proteins was in good agreement with the experimental results. The results indicated that (1) completely buried Glu192 contributes to the stabilization of PCP-0SH because of the formation of strong intramolecular hydrogen bonds and (2) the hydrogen bonds by the nonionized and buried Glu can contribute more than the burial of hydrophobic groups to the conformational stability of proteins.     

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.     

Improved pKa calculations through flexibility based sampling of a water-dominated interaction scheme

Ionizable groups play critical roles in biological processes. Computation of pK(a)s is complicated by model approximations and multiple conformations. Calculated and experimental pK(a)s are compared for relatively inflexible active-site side chains, to develop an empirical model for hydration entropy changes upon charge burial. The modification is found to be generally small, but large for cysteine, consistent with small molecule ionization data and with partial charge distributions in ionized and neutral forms. The hydration model predicts significant entropic contributions for ionizable residue burial, demonstrated for components in the pyruvate dehydrogenase complex. Conformational relaxation in a pH-titration is estimated with a mean-field assessment of maximal side chain solvent accessibility. All ionizable residues interact within a low protein dielectric finite difference (FD) scheme, and more flexible groups also access water-mediated Debye-Hückel (DH) interactions. The DH method tends to match overall pH-dependent stability, while FD can be more accurate for active-site groups. Tolerance for side chain rotamer packing is varied, defining access to DH interactions, and the best fit with experimental pK(a)s obtained. The new (FD/DH) method provides a fast computational framework for making the distinction between buried and solvent-accessible groups that has been qualitatively apparent from previous work, and pK(a) calculations are significantly improved for a mixed set of ionizable residues. Its effectiveness is also demonstrated with computation of the pH-dependence of electrostatic energy, recovering favorable contributions to folded state stability and, in relation to structural genomics, with substantial improvement (reduction of false positives) in active-site identification by electrostatic strain.