Experiments suggest that simulations may overestimate electrostatic contributions to the mechanical stability of a fibronectin type III domain

Steered molecular dynamics simulations have previously been used to investigate the mechanical properties of the extracellular matrix protein fibronectin. The simulations suggest that the mechanical stability of the tenth type III domain from fibronectin (FNfn10) is largely determined by a number of critical hydrogen bonds in the peripheral strands. Interestingly, the simulations predict that lowering the pH from 7 to approximately 4.7 will increase the mechanical stability of FNfn10 significantly (by approximately 33 %) due to the protonation of a few key acidic residues in the A and B strands. To test this simulation prediction, we used single-molecule atomic force microscopy (AFM) to investigate the mechanical stability of FNfn10 at neutral pH and at lower pH where these key residues have been shown to be protonated. Our AFM experimental results show no difference in the mechanical stability of FNfn10 at these different pH values. These results suggest that some simulations may overestimate the role played by electrostatic interactions in determining the mechanical stability of proteins.     

pH dependence of stability of the 10th human fibronectin type III domain: a computational study

We present detailed computational studies based on electrostatic calculations to evaluate the origins of pKa values and the pH dependence of stability for the 10th type III domain of human fibronectin (FNfn10). One of our goals is to validate the calculation protocols by comparison to experimental data (Koide, A.; Jordan, M. R.; Horner, S.; Batori, V.; Koide, S. Biochemistry 2001, 40, 10326-10333). Another goal is to evaluate the sensitivity of the calculated ionization free energies and apparent pKa values on local structural fluctuations, which do not alter the structural convergence to a particular architecture, by using a complete ensemble of solution NMR structures and the NMR average minimized structure of FNfn10 (Main, A. L.; Harvey, T. S.; Baron, M.; Boyd, J.; Campbell, I. D. Cell 1992, 71, 671-678). Our calculations demonstrate that, at high ionic strength, FNfn10 is more stable at low pH compared to neutral pH, in overall agreement with experimental data. This behavior is attributed to contributions from unfavorable Coulombic interactions in a surface patch for the pairs Asp7-Glu9 and Asp7-Asp23. The unfavorable interactions are decreased at low pH, where the acidic residues become neutral, and are further decreased at high ionic strength because of increased screening by salt ions. Elimination of the unfavorable interactions in the theoretical mutants Asp7Asn (D7N) and Asp7Lys (D7K) produce higher calculated stabilities at neutral pH and any ionic strength compared to the wild-type, in agreement with the experimental data. We also discuss subtleties in the calculated apparent pKa values and ionization free energies, which are not in agreement with the experimental data. This work demonstrates that comparative electrostatic calculations can provide rapid predictions of pH-dependent properties of proteins and can be significant aids in guiding the design of proteins with tailored properties.     

Electrostatic interactions in ubiquitin: stabilization of carboxylates by lysine amino groups

To explore electrostatic interactions in ubiquitin, pK(a) values have been determined by NMR for all 12 carboxyl groups in wild-type ubiquitin and in variants where single lysines have been replaced by neutral residues. Aspartate pK(a) values in ubiquitin range from 3.1 to 3.8 and are generally less than model compound values. Most aspartate pK(a) values are within 0.2 pH unit of those predicted with a simple Tanford-Kirkwood model. Glutamate pK(a) values range from 3.8 to 4.5, close to model compound values and differing by 0.1-0.8 pH unit from calculated values. To determine the role of positive charges in modulating carboxyl pK(a) values, we mutated lysines at positions 11, 29, and 33 to glutamine and threonine. NMR studies with these six single-site mutants reveal significant interactions of Lys 11 and Lys 29 with Glu 34 and Asp 21, respectively: pK(a) values for Glu 34 and Asp 21 increase by approximately 0.5-0.8 pH unit, similar to predicted values, when the lysines are replaced by neutral residues. In contrast, the predicted interaction between Lys 33 and Glu 34 is not observed experimentally. In some instances, substitution of lysine by glutamine and threonine did not lead to the same changes in carboxyl pK(a) values. These may reflect new short-range interactions between the mutated residues and the carboxyl groups. Carboxyl pK(a) shifts > 0.5 pH unit result from mutations at groups that are <5 A from the carboxyl group. No interactions are observed at >10 A.     

Analysis of electrostatic interactions in the denatured state ensemble of the N-terminal domain of L9 under native conditions

The pH dependence of protein stability is defined by the difference in the number of protons bound to the folded state and to the denatured state ensemble (DSE) as a function of pH. In many cases, the protonation behavior can be described as the sum of a set of independently titrating residues; in this case, the pH dependence of stability reflects differences in folded and DSE pK(a)'s. pH dependent stability studies have shown that there are energetically important interactions involving charged residues in the DSE of the N-terminal domain of L9 (NTL9), which affect significantly the stability of the protein. The DSE of wild type NTL9 cannot be directly characterized under native conditions because of its high stability. A destabilized double mutant of NTL9, V3AI4A, significantly populates the folded state and the DSE in the absence of denaturant. The two states are in slow exchange on the nuclear magnetic resonance time scale, and diffusion measurements indicate that the DSE is compact. The DSE pK(a)'s of all of the acidic residues were directly determined. The DSE pK(a) of Asp8 and Asp23 are depressed relative to model compounds values. Use of the mutant DSE pK(a)'s together with known native state pK(a)'s leads to a significantly improved agreement between the measured pH dependent stability and that predicted by the Tanford-Wyman linkage relationship. An analysis of the literature suggests that DSE interactions involving charged residues are relatively common and should be considered in discussions of protein stability.     

The determinants of carboxyl pKa values in turkey ovomucoid third domain

A computational methodology for protein pK(a) predictions, based on ab initio quantum mechanical treatment of part of the protein and linear Poisson-Boltzmann equation treatment of the bulk solvent, is presented. The method is used to predict and interpret the pK(a) values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in the serine protease inhibitor turkey ovomucoid third domain. All the predicted pK(a) values are within 0.5 pH units of experiment, with a root-mean-square deviation of 0.31 pH units. We show that the decreased pK(a) values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrogen bonds involving amide protons are shown to be particularly important, and the effect of hydrogen bonding is shown to be nonadditive. Hydrophobic effects are also shown to be important in raising the pK(a). Interactions with charged residues are shown to have relatively little effect on the carboxyl pK(a) values in this protein, in general agreement with experiment.     

Role of charged residues in stabilization of Pyrococcus horikoshii CutA1, which has a denaturation temperature of nearly 150 °C

The CutA1 protein from Pyrococcus horikoshii (PhCutA1), a hyperthermophile, has an unusually high content of charged residues and an unusually high denaturation temperature. To elucidate the role of ion-ion interactions in protein stability, mutant proteins of PhCutA1 in which charged residues were substituted by noncharged residues were comprehensively examined. The denaturation temperatures (T(d)) for 13 of 53 examined mutant proteins were higher than that of the wild-type (148.5 °C at pH 7.0), among which E99Q had the highest T(d) at 154.9 °C. R25A had the largest decrease in T(d) among single mutants at ΔT(d) = -12.4 °C. The average decrease in T(d) of Lys or Arg mutants was greater than that of Glu or Asp mutants, and the average change in T(d) (ΔT(d)) of 21 Glu mutants was negligible, at 0.03 ± 2.05 °C. However, the electrostatic energy (-159.3 kJ·mol(-1)) of PhCutA1 was quite high, compared with that of CutA1 from Escherichia coli (-9.7 kJ·mol(-1)), a mesophile. These results indicate that: (a) many Glu and Asp residues of PhCutA1 should be essential for highly efficient interactions with positively charged residues and for generating high electrostatic energy, although they were forced to be partially repulsive to each other; (b) the changes in stability of mutant proteins with a T(d) value of ~140-150 °C were able to be explained by considering factors important for protein stability and the structural features of mutant sites; and (c) these findings are useful for the design of proteins that are stable at temperatures > 100 °C.     

Insensitivity of perturbed carboxyl pK(a) values in the ovomucoid third domain to charge replacement at a neighboring residue

A number of carboxyl groups in turkey ovomucoid third domain (OMTKY3) have low pK(a) values. A previous study suggested that neighboring amino groups were primarily responsible for the low carboxyl pK(a) values. However, the expected elevation in pK(a) values for these amino groups was not observed. In the present study, site-directed mutagenesis is used to investigate the origins of perturbed carboxyl pK(a) values in OMTKY3. Electrostatic calculations suggest that Lys 34 has large effects, 0.4-0.6 unit, on Asp 7, Glu 10, and Glu 19 which are 5-11 A away from Lys 34. Two-dimensional (1)H NMR techniques were used to determine pK(a) values of the acidic residues in OMTKY3 mutants in which Lys 34 has been replaced with threonine and glutamine. Surprisingly, the pK(a) values in the mutants are very close to those of the wild-type protein. The insensitivity of the acidic residues to replacement of Lys 34 suggests that long-range electrostatic interactions play less of a role in perturbing carboxyl pK(a) values than originally thought. We hypothesize that hydrogen bonds play a key role in perturbing some of the carboxyl ionization equilibria in OMTKY3.     

Hydrogen bonding markedly reduces the pK of buried carboxyl groups in proteins

The ionizable groups in proteins with the lowest pKs are the carboxyl groups of aspartic acid side-chains. One of the lowest, pK=0.6, is observed for Asp76 in ribonuclease T1. This low pK appeared to result from hydrogen bonds to a water molecule and to the side-chains of Asn9, Tyr11, and Thr91. The results here confirm this by showing that the pK of Asp76 increases to 1.7 in N9A, to 4.0 in Y11F, to 4.2 in T91V, to 4.4 in N9A+Y11F, to 4.9 in N9A+T91V, to 5.9 in Y11F+T91V, and to 6.4 in the triple mutant: N9A+Y11F+T91V. In ribonuclease Sa, the lowest pK=2.4 for Asp33. This pK increases to 3.9 in T56A, which removes the hydrogen bond to Asp33, and to 4.4 in T56V, which removes the hydrogen bond and replaces the -OH group with a -CH(3) group. It is clear that hydrogen bonds are able to markedly lower the pK values of carboxyl groups in proteins. These same hydrogen bonds make large contributions to the conformational stability of the proteins. At pH 7, the stability of D76A ribonuclease T1 is 3.8 kcal mol(-1) less than wild-type, and the stability of D33A ribonuclease Sa is 4.1 kcal mol(-1) less than wild-type. There is a good correlation between the changes in the pK values and the changes in stability. The results suggest that the pK values for these buried carboxyl groups would be greater than 8 in the absence of hydrogen bonds, and that the hydrogen bonds and other interactions of the carboxyl groups contribute over 8 kcal mol(-1) to the stability.     

Determination of pK(a) values of carboxyl groups in the N-terminal domain of rat CD2: anomalous pK(a) of a glutamate on the ligand-binding surface

The ligand-binding surface of the T-lymphocyte glycoprotein CD2 has an unusually high proportion of charged residues, and ionic interactions are thought to play a significant role in defining the ligand specificity and binding affinity of CD2 with the structurally homologous ligands CD48 (in rodents) and CD58 (in humans). The determination of the electrostatic properties of these proteins can therefore contribute to our understanding of structure-activity relationships for these adhesion complexes that underpin T-cell adhesion to antigen-presenting cells. In this study, we investigated the pH titration behavior of the carboxyl groups of the N-terminal domain of rat CD2 (CD2d1) using the chemical shifts of backbone amide nitrogen-15 ((15)N) and proton NMR resonances, and carboxyl carbon-13 ((13)C) signals. The analysis revealed the presence of a glutamate (Glu41) on the binding surface of rat CD2 with an unusually elevated acidity constant (pK(a) = 6.73) for CD2d1 samples at 1.2 mM concentration. pH titration of CD2d1 at low protein concentration (0.1 mM) resulted in a slight decrease of the measured pK(a) of Glu41 to 6.36. The ionization of Glu41 exhibited reciprocal interactions with a second glutamate (Glu29) in a neighboring location, with both residues demonstrating characteristic biphasic titration behavior of the carboxyl (13)C resonances. Measurements at pH 5.5 of the two-bond deuterium isotope shift for the (13)C carboxyl resonances for Glu41 and Glu29 [(2)DeltaC(delta)(O(epsilon)D) = 0.2 and 0.1 ppm, respectively] were consistent with the assignment of the anomalous pK(a) to Glu41, under the strong influence of Glu29. The characterization of single site mutations of CD2d1 residues Glu41 and Glu29 to glutamine confirmed the anomalous pK(a) for Glu41, and indicated that electrostatic interaction with the Glu29 side chain is a significant contributing influence for this behavior in the wild-type protein. The implications of these observations are discussed with respect to recent structural and functional analyses of the interaction of rat CD2 with CD48. In particular, CD2 Glu41 must be a candidate residue to explain the previously reported strong pH dependence of binding of these two proteins in vitro.     

Are acidic and basic groups in buried proteins predicted to be ionized?

Ionizable residues play essential roles in proteins, modulating protein stability, fold and function. Asp, Glu, Arg, and Lys make up about a quarter of the residues in an average protein. Multi-conformation continuum electrostatic (MCCE) calculations were used to predict the ionization states of all acidic and basic residues in 490 proteins. Of all 36,192 ionizable residues, 93.5% were predicted to be ionized. Thirty-five percent have lost 4.08 kcal/mol solvation energy (DeltaDeltaG(rxn)) sufficient to shift a pK(a) by three pH units in the absence of other interactions and 17% have DeltaDeltaG(rxn) sufficient to shift pK(a) by five pH units. Overall 85% of these buried residues (DeltaDeltaG(rxn)>5DeltapK units) are ionized, including 92% of the Arg, 86% of the Asp, 77% of the Glu, and 75% of the Lys. Ion-pair interactions stabilize the ionization of both acids and bases. The backbone dipoles stabilize anions more than cations. The interactions with polar side-chains are also different for acids and bases. Asn and Gln stabilize all charges, Ser and Thr stabilize only acids while Tyr rarely stabilize Lys. Thus, hydroxyls are better hydrogen bond donors than acceptors. Buried ionized residues are more likely to be conserved than those on the surface. There are 3.95 residues buried per 100 residues in an average protein.     

Hydrogen bonding is the prime determinant of carboxyl pKa values at the N-termini of alpha-helices

Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.     

Stabilization of apoflavodoxin by replacing hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or Gln

Knowledge of protein stability principles provides a means to increase protein stability in a rational way. Here we explore the feasibility of stabilizing proteins by replacing solvent-exposed hydrogen-bonded charged Asp or Glu residues by the neutral isosteric Asn or GLN: The rationale behind this is a previous observation that, in some cases, neutral hydrogen bonds may be more stable that charged ones. We identified, in the apoflavodoxin from Anabaena PCC 7119, three surface-exposed aspartate or glutamate residues involved in hydrogen bonding with a single partner and we mutated them to asparagine or glutamine, respectively. The effect of the mutations on apoflavodoxin stability was measured by both urea and temperature denaturation. We observed that the three mutant proteins are more stable than wild-type (on average 0.43 kcal/mol from urea denaturation and 2.8 degrees C from a two-state analysis of fluorescence thermal unfolding data). At high ionic strength, where potential electrostatic repulsions in the acidic apoflavodoxin should be masked, the three mutants are similarly more stable (on average 0.46 kcal/mol). To rule out further that the stabilization observed is due to removal of electrostatic repulsions in apoflavodoxin upon mutation, we analysed three control mutants and showed that, when the charged residue mutated to a neutral one is not hydrogen bonded, there is no general stabilizing effect. Replacing hydrogen-bonded charged Asp or Glu residues by Asn or Gln, respectively, could be a straightforward strategy to increase protein stability.     

Cumulative site-directed charge-change replacements in bacteriophage T4 lysozyme suggest that long-range electrostatic interactions contribute little to protein stability

Bacteriophage T4 lysozyme is a basic molecule with an isoelectric point above 9.0, and an excess of nine positive charges at neutral pH. It might be expected that it would be energetically costly to bring these out-of-balance charges from the extended, unfolded, form of the protein into the compact folded state. To determine the contribution of such long-range electrostatic interactions to the stability of the protein, five positively charged surface residues, Lys16, Arg119, Lys135, Lys147 and Arg154, were individually replaced with glutamic acid. Eight selected double, triple and quadruple mutants were also constructed so as to sequentially reduce the out-of-balance formal charge on the molecule from +9 to +1 units. Each of the five single variant proteins was crystallized and high-resolution X-ray analysis confirmed that each mutant structure was, in general, very similar to the wild-type. In the case of R154E, however, the Arg154 to Glu replacement caused a rearrangement in which Asp127 replaced Glu128 as the capping residue of a nearby alpha-helix. The thermal stabilities of all 13 variant proteins were found to be fairly similar, ranging from 0.5 kcal/mol more stable than wild-type to 1.7 kcal/mol less stable than wild-type. In the case of the five single charge-change variants, for which the structures were determined, the changes in stability can be rationalized in terms of changes in local interactions at the site of the replacement. There is no evidence that the reduction in the out-of-balance charge on the molecule increases the stability of the folded relative to the unfolded form, either at pH 2.8 or at pH 5.3. This indicates that long-range electrostatic interactions between the substituted amino acid residues and other charged groups on the surface of the molecule are weak or non-existent. Furthermore, the relative stabilities of the multiple charge replacement mutant proteins were found to be almost exactly equal to the sums of the relative stabilities of the constituent single mutant proteins. This also clearly indicates that the electrostatic interactions between the replaced charges are negligibly small. The activities of the charge-change mutant lysozymes, as measured by the rate of hydrolysis of cell wall suspensions, are essentially equal to that of the wild-type lysozyme, but on a lysoplate assay the mutant enzymes appear to have higher activity.(ABSTRACT TRUNCATED AT 400 WORDS)     

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.     

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.     

Dissecting the electrostatic interactions and pH-dependent activity of a family 11 glycosidase

Previous studies of the low molecular mass family 11 xylanase from Bacillus circulans show that the ionization state of the nucleophile (Glu78, pK(a) 4.6) and the acid/base catalyst (Glu172, pK(a) 6.7) gives rise to its pH-dependent activity profile. Inspection of the crystal structure of BCX reveals that Glu78 and Glu172 are in very similar environments and are surrounded by several chemically equivalent and highly conserved active site residues. Hence, there are no obvious reasons why their apparent pK(a) values are different. To address this question, a mutagenic approach was implemented to determine what features establish the pK(a) values (measured directly by (13)C NMR and indirectly by pH-dependent activity profiles) of these two catalytic carboxylic acids. Analysis of several BCX variants indicates that the ionized form of Glu78 is preferentially stabilized over that of Glu172 in part by stronger hydrogen bonds contributed by two well-ordered residues, namely, Tyr69 and Gln127. In addition, theoretical pK(a) calculations show that Glu78 has a lower pK(a) value than Glu172 due to a smaller desolvation energy and more favorable background interactions with permanent partial charges and ionizable groups within the protein. The pK(a) value of Glu172 is in turn elevated due to electrostatic repulsion from the negatively charged glutamate at position 78. The results also indicate that all of the conserved active site residues act concertedly in establishing the pK(a) values of Glu78 and Glu172, with no particular residue being singly more important than any of the others. In general, residues that contribute positive charges and hydrogen bonds serve to lower the pK(a) values of Glu78 and Glu172. The degree to which a hydrogen bond lowers a pK(a) value is largely dependent on the length of the hydrogen bond (shorter bonds lower pK(a) values more) and the chemical nature of the donor (COOH > OH > CONH(2)). In contrast, neighboring carboxyl groups can either lower or raise the pK(a) values of the catalytic glutamic acids depending upon the electrostatic linkage of the ionization constants of the residues involved in the interaction. While the pH optimum of BCX can be shifted from -1.1 to +0.6 pH units by mutating neighboring residues within the active site, activity is usually compromised due to the loss of important ground and/or transition state interactions. These results suggest that the pH optima of an enzyme might be best engineered by making strategic amino acid substitutions, at positions outside of the "core" active site, that electrostatically influence catalytic residues without perturbing their immediate structural environment.     

Intrinsically disordered inhibitor of glutamine synthetase is a functional protein with random‐coil‐like pKa values

The sequential action of glutamine synthetase (GS) and glutamate synthase (GOGAT) in cyanobacteria allows the incorporation of ammonium into carbon skeletons. In the cyanobacterium Synechocystis sp. PCC 6803, the activity of GS is modulated by the interaction with proteins, which include a 65‐residue‐long intrinsically disordered protein (IDP), the inactivating factor IF7. This interaction is regulated by the presence of charged residues in both IF7 and GS. To understand how charged amino acids can affect the binding of an IDP with its target and to provide clues on electrostatic interactions in disordered states of proteins, we measured the pK_a values of all IF7 acidic groups (Glu32, Glu36, Glu38, Asp40, Asp58, and Ser65, the backbone C‐terminus) at 100 mM NaCl concentration, by using NMR spectroscopy. We also obtained solution structures of IF7 through molecular dynamics simulation, validated them on the basis of previous experiments, and used them to obtain theoretical estimates of the pK_a values. Titration values for the two Asp and three Glu residues of IF7 were similar to those reported for random‐coil models, suggesting the lack of electrostatic interactions around these residues. Furthermore, our results suggest the presence of helical structure at the N‐terminus of the protein and of conformational changes at acidic pH values. The overall experimental and in silico findings suggest that local interactions and conformational equilibria do not play a role in determining the electrostatic features of the acidic residues of IF7.     

Inverse electrostatic effect: electrostatic repulsion in the unfolded state stabilizes a leucine zipper

The pH-dependent stability of a protein is strongly affected by electrostatic interactions between ionizable residues in the folded as well as unfolded state. Here we characterize the individual contributions of charged Glu and His residues to stability and determine the NMR structure of the designed, heterodimeric leucine zipper AB consisting of an acidic A chain and a basic B chain. Thermodynamic parameters are compared with those of the homologous leucine zipper AB(SS) in which the A and B chains are disulfide-linked. NMR structures of AB based on (1)H NMR data collected at 600 MHz converge, and formation of the same six interchain salt bridges found previously in disulfide-linked AB(SS) [Marti, D. N., and Bosshard, H. R. (2003) J. Mol. Biol. 330, 621-637] is indicated. While the structures of AB and AB(SS) are very similar, their pH-dependent relative stabilities are strikingly different. The stability of AB peaks at pH approximately 4.5 and is higher at pH 8 than at pH 2. In contrast, AB(SS) is most stable at acidic pH where no interhelical salt bridges are formed. The different energetic contributions of charged Glu and His residues to stability of the two coiled coil structures were evaluated from pK(a) shifts induced by folding. The six charged Glu residues involved in salt bridges stabilize leucine zipper AB by 4.5 kJ/mol yet destabilize disulfide-linked AB(SS) by -1.1 kJ/mol. Two non-ion-paired Glu charges destabilize AB by only -1.8 kJ/mol but AB(SS) by -5.6 kJ/mol. The higher relative stability of AB at neutral pH is not caused by more favorable electrostatic interactions in the folded leucine zipper. It is due mainly to unfavorable electrostatic interactions in the unfolded A and B chains and may therefore be called an inverse electrostatic effect. This study illustrates the importance of residual interactions in the unfolded state and how the energetics of the unfolded state affect the stability of the folded protein.     

Rapid Calculation of Protein pK(a) Values Using Rosetta

We developed a Rosetta-based Monte Carlo method to calculate the pK(a) values of protein residues that commonly exhibit variable protonation states (Asp, Glu, Lys, His, and Tyr). We tested the technique by calculating pK(a) values for 264 residues from 34 proteins. The standard Rosetta score function, which is independent of any environmental conditions, failed to capture pK(a) shifts. After incorporating a Coulomb electrostatic potential and optimizing the solvation reference energies for pK(a) calculations, we employed a method that allowed side-chain flexibility and achieved a root mean-square deviation (RMSD) of 0.83 from experimental values (0.68 after discounting 11 predictions with an error over 2 pH units). Additional degrees of side-chain conformational freedom for the proximal residues facilitated the capture of charge-charge interactions in a few cases, resulting in an overall RMSD of 0.85 pH units. The addition of backbone flexibility increased the overall RMSD to 0.93 pH units but improved relative pK(a) predictions for proximal catalytic residues. The method also captures large pK(a) shifts of lysine and some glutamate point mutations in staphylococcal nuclease. Thus, a simple and fast method based on the Rosetta score function and limited conformational sampling produces pK(a) values that will be useful when rapid estimation is essential, such as in docking, design, and folding.     

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.