Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Salt bridges destabilize a leucine zipper designed for maximized ion pairing between helices

Interhelical salt bridges are common in leucine zippers and are thought to stabilize the coiled coil conformation. Here we present a detailed thermodynamic investigation of the designed, disulfide-linked leucine zipper AB(SS) whose high-resolution NMR structure shows six interhelical ion pairs between heptad positions g of one helix and e' of the other helix but no ion pairing within single helices. The average pK(a) value of the Glu side chain carboxyl groups of AB(SS) is slightly higher than the pK(a) of a freely accessible Glu in an unfolded peptide [Marti, D. N., Jelesarov, I., and Bosshard, H. R. (2000) Biochemistry 39, 12804-12818]. This indicates that the salt bridges are destabilizing, a prediction we now have confirmed by determining the pH +/- stability profile of AB(SS). Circular dichroism-monitored unfolding by urea and by heating and differential scanning calorimetry show that the coiled coil conformation is approximately 5 kJ/mol more stable when salt bridges are broken by protonation of the carboxyl side chains. Using guanidinium chloride as the denaturant, the increase in the free energy of unfolding on protonation of the carboxyl side chains is larger, approximately 17 kJ/mol. The discrepancy between urea and guanidinium chloride unfolding can be ascribed to the ionic nature of guanidinium chloride, which screens charge-charge interactions. This work demonstrates the difficulty of predicting the energetic contribution of salt bridges from structural data alone even in a case where the ion pairs are seen in high-resolution NMR structures. The reason is that the contribution to stability results from a fine balance between energetically favorable Coulombic attractions and unfavorable desolvation of charges and conformational constraints of the residues involved in ion pairing. The apparent discrepancy between the results presented here and mutational studies indicating stabilization by salt bridges is discussed and resolved. An explanation is proposed for why interhelical salt bridges are frequently found in natural coiled coils despite evidence that they do not directly contribute to stability.     

Interhelical ion pairing in coiled coils: solution structure of a heterodimeric leucine zipper and determination of pKa values of Glu side chains

Residues of opposite charge often populate heptad positions g (heptad i on chain 1) and e' (heptad i + 1 on chain 2) in dimeric coiled coils and may stabilize the dimer by formation of interchain ion pairs. To investigate the contribution to stability of such electrostatic interactions we have designed a disulfide-linked heterodimeric zipper (AB zipper) consisting of the acidic chain Ac-E-VAQLEKE-VAQAEAE-NYQLEQE-VAQLEHE-CG-NH(2) and the basic chain Ac-E-VQALKKR-VQALKAR-NYAAKQK-VQALRHK-CG-NH(2) in which all e and g positions are occupied by either E or K/R to form a maximum of seven interhelical salt bridges. Temperature-induced denaturation experiments monitored by circular dichroism reveal a stable coiled coil conformation below 50 degrees C and in the pH range 1.2-8.0. Stability is highest at pH approximately 4.0 [DeltaG(U) (37 degrees C) = 5.18 +/- 0.51 kcal mol(-)(1)]. The solution structure of the AB zipper at pH 5.65 has been elucidated on the basis of homonuclear (1)H NMR data collected at 800 MHz [heavy atom rmsd's for the ensemble of 50 calculated structures are 0.47 +/- 0.13 A (backbone) and 0.95 +/- 0.16 A (all)]. Both chains of the AB zipper are almost entirely in alpha-helical conformation and form a superhelix with a left-handed twist. Overhauser connectivities reveal close contacts between g position residues (heptad i on chain 1) and residues d/f (heptad i on chain 1), residues a/d (heptad i + 1 on chain 1), and residue a' (heptad i + 1 on chain 2). Residues in position e (heptad i on chain 1) are in contact with residues a/b/d/f (heptad i on chain 1) and residue d' (heptad i on chain 2). These connectivities hint at a relatively defined alignment of the side chains across the helix interface. Partial H-bond formation between the functional groups of residues g and e'(+1) is observed in the calculated structures. NMR pH titration experiments disclose pK(a) values for Glu delta-carboxylate groups: 4.14 +/- 0.02 (E(1)), 4.82 +/- 0.07 (E(6)), 4.52 +/- 0.01 (E(8)), 4.37 +/- 0.03 (E(13)), 4.11 +/- 0.02 (E(15)), 4.41 +/- 0.07 (E(20)), 4.82 +/- 0.03 (E(22)), 4.65 +/- 0.04 (E(27)), 4.63 +/- 0.03 (E(29)), 4.22 +/- 0.02 (E(1)(')). By comparison with pK(a) of Glu in unfolded peptides ( approximately 4. 3 +/- 0.1), our pK(a) data suggest marginal or even unfavorable contribution of charged Glu to the stability of the AB zipper. The electrostatic energy gained from interhelical ion pairs is likely to be surpassed by hydrophobic energy terms upon protonation of Glu, due to increased hydrophobicity of uncharged Glu and, thus, better packing against apolar residues at the chain interface.     

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)     

Stabilization of a fibronectin type III domain by the removal of unfavorable electrostatic interactions on the protein surface

It is generally considered that electrostatic interactions on the protein surface, such as ion pairs, contribute little to protein stability, although they may play important roles in conformational specificity. We found that the tenth fibronectin type III domain of human fibronectin (FNfn10) is more stable at acidic pH than neutral pH, with an apparent midpoint of transition near pH 4. Determination of pK(a)'s for all the side chain carboxyl groups of Asp and Glu residues revealed that Asp 23 and Glu 9 have an upshifted pK(a). These residues and Asp 7 form a negatively charged patch on the surface of FNfn10, with Asp 7 centrally located between Asp 23 and Glu 9, suggesting repulsive electrostatic interactions among these residues at neutral pH. Mutant proteins, D7N and D7K, in which Asp 7 was replaced with Asn and Lys, respectively, exhibited a modest but significant increase in stability at neutral pH, compared to the wild type, and they no longer showed pH dependence of stability. The pK(a)'s of Asp 23 and Glu 9 in these mutant proteins shifted closer to their respective unperturbed values, indicating that the unfavorable electrostatic interactions have been reduced in the mutant proteins. Interestingly, the wild-type and mutant proteins were all stabilized to a similar degree by the addition of 1 M sodium chloride at both neutral and acidic pH, suggesting that the repulsive interactions between the carboxyl groups cannot be effectively shielded by 1 M sodium chloride. These results indicate that repulsive interactions between like charges on the protein surface can destabilize a protein, and protein stability can be significantly improved by relieving these interactions.     

Contributions of engineered surface salt bridges to the stability of T4 lysozyme determined by directed mutagenesis.

Six designed mutants of T4 lysozyme were created in an attempt to create putative salt bridges on the surface of the protein. The first three of the mutants, T115E (Thr 115 to Glu), Q123E, and N144E, were designed to introduce a new charged side chain close to one or more existing charged groups of the opposite sign on the surface of the protein. In each of these cases the putative electrostatic interactions introduced by the mutation include possible salt bridges between residues within consecutive turns of an alpha-helix. Effects of the mutations ranged from no change in stability to a 1.5 degrees C (0.5 kcal/mol) increase in melting temperature. In two cases, secondary (double) mutants were constructed as controls in which the charge partner was removed from the primary mutant structure. These controls proteins indicate that the contributions to stability from each of the engineered salt bridges is very small (about 0.1-0.25 kcal/mol in 0.15 M KCl). The structures of the three primary mutants were determined by X-ray crystallography and shown to be essentially the same as the wild-type structure except at the site of the mutation. Although the introduced charges in the T115E and Q123E structures are within 3-5 A of their intended partner, the introduced side chains and their intended partners were observed to be quite mobile. It has been shown that the salt bridge between His 31 and Asp 70 in T4 lysozyme stabilizes the protein by 3-5 kcal/mol [Anderson, D. E., Becktel, W. J., & Dahlquist, F. W. (1990) Biochemistry 29, 2403-2408]. To test the effectiveness of His...Asp interactions in general, three additional double mutants, K60H/L13D, K83H/A112D, and S90H/Q122D, were created in order to introduce histidine-aspartate charge pairs on the surface of the protein. Each of these mutants destabilizes the protein by 1-3 kcal/mol in 0.15 M KCl at pH values from 2 to 6.5. The X-ray crystallographic structure of the mutant K83H/A112D has been determined and shows that there are backbone conformational changes of 0.3-0.6 A extending over several residues. The introduction of the histidine and aspartate presumably introduces strain into the folded protein that destabilizes this variant. It is concluded that pairs of oppositely charged residues that are on the surface of a protein and have freedom to adopt different conformations do not tend to come together to form structurally localized salt bridges. Rather, such residues tend to remain mobile, interact weakly if at all, and do not contribute significantly to protein stability. It is argued that the entropic cost of localizing a pair of solvent-exposed charged groups on the surface of a protein largely offsets the interaction energy expected from the formation of a defined salt bridge. There are examples of strong salt bridges in proteins, but such interactions require that the folding of the protein provides the requisite driving energy to hold the interacting partners in the correct rigid alignment.     

An electrostatic basis for the stability of thermophilic proteins

Two factors provide key contributions to the stability of thermophilic proteins relative to their mesophilic homologues: electrostatic interactions of charged residues in the folded state and the dielectric response of the folded protein. The dielectric response for proteins in a "thermophilic series" globally modulates the thermal stability of its members, with the calculated dielectric constant for the protein increasing from mesophiles to hyperthermophiles. This variability results from differences in the distribution of charged residues on the surface of the protein, in agreement with structural and genetic observations. Furthermore, the contribution of electrostatic interactions to the stability of the folded state is more favorable for thermophilic proteins than for their mesophilic homologues. This leads to the conclusion that electrostatic interactions play an important role in determining the stability of proteins at high temperatures. The interplay between electrostatic interactions and dielectric response also provides further rationalization for the enhanced stability of thermophilic proteins with respect to cold-denaturation. Taken together, the distribution of charged residues and their fluctuations have been shown to be factors in modulating protein stability over the entire range of biologically relevant temperatures.     

Fluctuations between stabilizing and destabilizing electrostatic contributions of ion pairs in conformers of the c-Myc-Max leucine zipper

In solution proteins often exhibit backbone and side-chain flexibility. Yet electrostatic interactions in proteins are sensitive to motions. Hence, here we study the contribution of ion pairs toward protein stability in a range of conformers which sample the conformational space in solution. Specifically, we focus on the electrostatic contributions of ion pairs to the stability of each of the conformers in the NMR ensemble of the c-Myc-Max leucine zipper and to their average energy minimized structure. We compute the electrostatic contributions of inter- and intra-helical ion pairs and of an ion pair network. We find that the electrostatic contributions vary considerably among the 40 NMR conformers. Each ion pair, and the network, fluctuates between being stabilizing and being destabilizing. This fluctation reflects the variability in the location of the ion pairing residues and in the geometric orientation of these residues, both with respect to each other and with respect to other charged groups in the rest of the protein. Ion pair interactions in the c-Myc-Max leucine zipper in solution depend on the protein conformer which is analyzed. Hence, the overall stabilizing (or destabilizing) contribution of an ion pair is conformer population-dependent. This study indicates that free energy calculations performed using the continuum electrostatics methodology are sensitive to protein conformational details.     

Molecular dynamics simulations of human prion protein: importance of correct treatment of electrostatic interactions.

Molecular dynamics simulations have been used to investigate the dynamical and structural behavior of a homology model of human prion protein HuPrP(90-230) generated from the NMR structure of the Syrian hamster prion protein ShPrP(90-231) and of ShPrP(<90-231) itself. These PrPs have a large number of charged residues on the protein surface. At the simulation pH 7, HuPrP(90-230) has a net charge of -1 eu from 15 positively and 14 negatively charged residues. Simulations for both PrPs, using the AMBER94 force field in a periodic box model with explicit water molecules, showed high sensitivity to the correct treatment of the electrostatic interactions. Highly unstable behavior of the structured region of the PrPs (127-230) was found using the truncation method, and stable trajectories could be achieved only by including all the long-range electrostatic interactions using the particle mesh Ewald (PME) method. The instability using the truncation method could not be reduced by adding sodium and chloride ions nor by replacing some of the sodium ions with calcium ions. The PME simulations showed, in accordance with NMR experiments with ShPrP and mouse PrP, a flexibly disordered N-terminal part, PrP(90-126), and a structured C-terminal part, PrP(127-230), which includes three alpha-helices and a short antiparallel beta-strand. The simulations showed some tendency for the highly conserved hydrophobic segment PrP(112-131) to adopt an alpha-helical conformation and for helix C to split at residues 212-213, a known disease-associated mutation site (Q212P). Three highly occupied salt bridges could be identified (E146/D144<-->R208, R164<-->D178, and R156<-->E196) which appear to be important for the stability of PrP by linking the stable main structured core (helices B and C) with the more flexible structured part (helix A and strands A and B). Two of these salt bridges involve disease-associated mutations (R208H and D178N). Decreased PrP stability shown by protein unfolding experiments on mutants of these residues and guanidinium chloride or temperature-induced unfolding studies indicating reduced stability at low pH are consistent with stabilization by salt bridges. The fact that electrostatic interactions, in general, and salt bridges, in particular, appear to play an important role in PrP stability has implications for PrP structure and stability at different pHs it may encounter physiologically during normal or abnormal recycling from the pH neutral membrane surface into endosomes or lysomes (acidic pHs) or in NMR experiments (5.2 for ShPrP and 4.5 for mouse PrP).     

Attractive interhelical electrostatic interactions in the proline- and acidic-rich region (PAR) leucine zipper subfamily preclude heterodimerization with other basic leucine zipper subfamilies

Basic region-leucine zipper (B-ZIP) proteins homo- or heterodimerize to bind sequence-specific double-stranded DNA. We present circular dichroism (CD) thermal denaturation data on vitellogenin promoter-binding protein (VBP), a member of the PAR subfamily of B-ZIP proteins that also includes thyroid embryonic factor, hepatocyte leukemia factor, and albumin site D-binding protein. VBP does not heterodimerize with B-ZIP domains from C/EBP alpha, JUND, or FOS. We describe a dominant negative protein, A-VBP, that contains the VBP leucine zipper and an acidic amphipathic protein sequence that replaces the basic region critical for DNA binding. The acidic extension forms a coiled coil structure with the VBP basic region in the VBP.A-VBP heterodimer. This new alpha-helical structure extends the leucine zipper N-terminally, stabilizing the complex by 2.0 kcal/mol. A-VBP abolishes DNA binding of VBP in an equimolar competition assay, but does not affect DNA binding even at 100-fold excess of CREB, C/EBP alpha, or FOS/JUND. Likewise, proteins containing the acidic extension appended to seven other leucine zippers do not inhibit VBP DNA binding. We show that conserved g <--> e' or i, i' +5 salt bridges are sufficient to confer specificity to VBP by mutating the C/EBPalpha leucine zipper to contain the g <--> e' salt bridges that characterize VBP. A-VBP heterodimerizes with this mutant C/EBP, preventing it from binding to DNA. These conserved g <--> e' electrostatic interactions define the specificity of the PAR subfamily of B-ZIP proteins and preclude interaction with other B-ZIP subfamilies.     

Relationship between ion pair geometries and electrostatic strengths in proteins

The electrostatic free energy contribution of an ion pair in a protein depends on two factors, geometrical orientation of the side-chain charged groups with respect to each other and the structural context of the ion pair in the protein. Conformers in NMR ensembles enable studies of the relationship between geometry and electrostatic strengths of ion pairs, because the protein structural contexts are highly similar across different conformers. We have studied this relationship using a dataset of 22 unique ion pairs in 14 NMR conformer ensembles for 11 nonhomologous proteins. In different NMR conformers, the ion pairs are classified as salt bridges, nitrogen-oxygen (N-O) bridges and longer-range ion pairs on the basis of geometrical criteria. In salt bridges, centroids of the side-chain charged groups and at least a pair of side-chain nitrogen and oxygen atoms of the ion-pairing residues are within a 4 A distance. In N-O bridges, at least a pair of the side-chain nitrogen and oxygen atoms of the ion-pairing residues are within 4 A distance, but the distance between the side-chain charged group centroids is greater than 4 A. In the longer-range ion pairs, the side-chain charged group centroids as well as the side-chain nitrogen and oxygen atoms are more than 4 A apart. Continuum electrostatic calculations indicate that most of the ion pairs have stabilizing electrostatic contributions when their side-chain charged group centroids are within 5 A distance. Hence, most (approximately 92%) of the salt bridges and a majority (68%) of the N-O bridges are stabilizing. Most (approximately 89%) of the destabilizing ion pairs are the longer-range ion pairs. In the NMR conformer ensembles, the electrostatic interaction between side-chain charged groups of the ion-pairing residues is the strongest for salt bridges, considerably weaker for N-O bridges, and the weakest for longer-range ion pairs. These results suggest empirical rules for stabilizing electrostatic interactions in proteins.     

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.     

Protein destabilization by electrostatic repulsions in the two-stranded alpha-helical coiled-coil/leucine zipper.

The destabilizing effect of electrostatic repulsions on protein stability has been studied by using synthetic two-stranded alpha-helical coiled-coils as a model system. The native coiled-coil consists of two identical 35-residue polypeptide chains with a heptad repeat QgVaGbAcLdQeKf and a Cys residue at position 2 to allow formation of an interchain disulfide bridge. This peptide, designed to contain no intrahelical or interhelical electrostatic interactions, forms a stable coiled-coil structure at 20 degrees C in benign medium (50 mM KCl, 25 mM PO4, pH 7) with a [urea]1/2 value of 6.1 M. Four mutant coiled-coils were designed to contain one or two Glu substitutions for Gln per polypeptide chain. The resulting coiled-coils contained potential i to i' + 5 Glu-Glu interchain repulsions (denoted as peptide E2(15,20)), i to i' + 2 Glu-Glu interchain repulsions (denoted E2(20,22)), or no interchain ionic interactions (denoted E2(13,22) and E1(20)). The stabilities of the coiled-coils were determined by measuring the ellipticities at 222 nm as a function of urea or guanidine hydrochloride concentration at 20 degrees C in the presence and absence of an interchain disulfide bridge. At pH 7, in the presence of urea, the stabilities of E2(13,22) and E2(20,22) were identical suggesting that the potential i to i' + 2 interchain Glu-Glu repulsion in the E2(20,22) coiled-coil does not occur. In contrast, the mutant E2(15,20) is substantially less stable than E2(13,22) or E2(15,20) by 0.9 kcal/mol due to the presence of two i to i' + 5 interchain Glu-Glu repulsions, which destabilize the coiled-coil by 0.45 kcal/mol each. At pH 3 the coiled-coils were found to increase in stability as the number of Glu substitutions were increased. This, combined with reversed-phase HPLC results at pH 7 and pH 2, supports the conclusion that the protonated Glu side chains present at low pH are significantly more hydrophobic than Gln side chains which are in turn more hydrophobic than the ionized Glu side chains present at neutral pH. The protonated Glu residues increase the hydrophobicity of the coiled-coil interface leading to higher coiled-coil stability. The guanidine hydrochloride results at pH 7 show similar stabilities between the native and mutant coiled-coils indicating that guanidine hydrochloride masks electrostatic repulsions due to its ionic nature and that Glu and Gln in the e and g positions of the heptad repeat have very similar effects on coiled-coil stability in the presence of GdnHCl.     

Protein stabilization by salt bridges: concepts, experimental approaches and clarification of some misunderstandings

Salt bridges in proteins are bonds between oppositely charged residues that are sufficiently close to each other to experience electrostatic attraction. They contribute to protein structure and to the specificity of interaction of proteins with other biomolecules, but in doing so they need not necessarily increase a protein's free energy of unfolding. The net electrostatic free energy of a salt bridge can be partitioned into three components: charge-charge interactions, interactions of charges with permanent dipoles, and desolvation of charges. Energetically favorable Coulombic charge-charge interaction is opposed by often unfavorable desolvation of interacting charges. As a consequence, salt bridges may destabilize the structure of the folded protein. There are two ways to estimate the free energy contribution of salt bridges by experiment: the pK(a) approach and the mutation approach. In the pK(a) approach, the contribution of charges to the free energy of unfolding of a protein is obtained from the change of pK(a) of ionizable groups caused by altered electrostatic interactions upon folding of the protein. The pK(a) approach provides the relative free energy gained or lost when ionizable groups are being charged. In the mutation approach, the coupling free energy between interacting charges is obtained from a double mutant cycle. The coupling free energy is an indirect and approximate measure of the free energy of charge-charge interaction. Neither the pK(a) approach nor the mutation approach can provide the net free energy of a salt bridge. Currently, this is obtained only by computational methods which, however, are often prone to large uncertainties due to simplifying assumptions and insufficient structural information on which calculations are based. This state of affairs makes the precise thermodynamic quantification of salt bridge energies very difficult. This review is focused on concepts and on the assessment of experimental methods and does not cover the vast literature.     

Contribution of surface salt bridges to protein stability: guidelines for protein engineering

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.     

Formyl-coenzyme A (CoA):oxalate CoA-transferase from the acidophile Acetobacter aceti has a distinctive electrostatic surface and inherent acid stability

Bacterial formyl-CoA:oxalate CoA-transferase (FCOCT) and oxalyl-CoA decarboxylase work in tandem to perform a proton-consuming decarboxylation that has been suggested to have a role in generalized acid resistance. FCOCT is the product of uctB in the acidophilic acetic acid bacterium Acetobacter aceti. As expected for an acid-resistance factor, UctB remains folded at the low pH values encountered in the A. aceti cytoplasm. A comparison of crystal structures of FCOCTs and related proteins revealed few features in UctB that would distinguish it from nonacidophilic proteins and thereby account for its acid stability properties, other than a strikingly featureless electrostatic surface. The apparently neutral surface is a result of a "speckled" charge decoration, in which charged surface residues are surrounded by compensating charges but do not form salt bridges. A quantitative comparison among orthologs identified a pattern of residue substitution in UctB that may be a consequence of selection for protein stability by constant exposure to acetic acid. We suggest that this surface charge pattern, which is a distinctive feature of A. aceti proteins, creates a stabilizing electrostatic network without stiffening the protein or compromising protein-solvent interactions.     

Single molecule force spectroscopy reveals that electrostatic interactions affect the mechanical stability of proteins

It is well known that electrostatic interactions play important roles in determining the thermodynamic stability of proteins. However, the investigation into the role of electrostatic interactions in mechanical unfolding of proteins has just begun. Here we used single molecule atomic force microscopy techniques to directly evaluate the effect of electrostatic interactions on the mechanical stability of a small protein GB1. We engineered a bi-histidine motif into the force-bearing region of GB1. By varying the pH, histidine residues can switch between protonated and deprotonated states, leading to the change of the electrostatic interactions between the two histidine residues. We found that the mechanical unfolding force of the engineered protein decreased by ∼34% (from 115 pN to 76 pN) on changing the pH from 8.5 to 3, due to the increased electrostatic repulsion between the two positively charged histidines at acidic pH. Our results demonstrated that electrostatic interactions can significantly affect the mechanical stability of elastomeric proteins, and modulating the electrostatic interactions of key charged residues can become a promising method for regulating the mechanical stability of elastomeric proteins.     

Contribution of salt bridges near the surface of a protein to the conformational stability

Salt bridges play important roles in the conformational stability of proteins. However, the effect of a surface salt bridge on the stability remains controversial even today; some reports have shown little contribution of a surface salt bridge to stability, whereas others have shown a favorable contribution. In this study, to elucidate the net contribution of a surface salt bridge to the conformational stability of a protein, systematic mutant human lysozymes, containing one Glu to Gln (E7Q) and five Asp to Asn mutations (D18N, D49N, D67N, D102N, and D120N) at residues where a salt bridge is formed near the surface in the wild-type structure, were examined. The thermodynamic parameters for denaturation between pH 2.0 and 4.8 were determined by use of a differential scanning calorimeter, and the crystal structures were analyzed by X-ray crystallography. The denaturation Gibbs energy (DeltaG) of all mutant proteins was lower than that of the wild-type protein at pH 4, whereas there was little difference between them near pH 2. This is caused by the fact that the Glu and Asp residues are ionized at pH 4 but protonated at pH 2, indicating a favorable contribution of salt bridges to the wild-type structure at pH 4. Each contribution was not equivalent, but we found that the contributions correlate with the solvent inaccessibility of the salt bridges; the salt bridge contribution was small when 100% accessible, while it was about 9 kJ/mol if 100% inaccessible. This conclusion indicates how to reconcile a number of conflicting reports about role of surface salt bridges in protein stability. Furthermore, the effect of salts on surface salt bridges was also examined. In the presence of 0.2 M KCl, the stability at pH 4 decreased, and the differences in stability between the wild-type and mutant proteins were smaller than those in the absence of salts, indicating the compensation to the contribution of salt bridges with salts. Salt bridges with more than 50% accessibility did not contribute to the stability in the presence of 0.2 M KCl.     

Dynamic properties of pH-dependent structural organization of the amyloidogenic ��-protein (1�C40)

The structural organization of the amyloidogenic β-protein containing 40 amino acid residues (Aβ40) was studied by the high temperature molecular dynamics simulations in the acidic (pH ∼ 3) and basic (pH ∼ 8) pH regions. The obtained data suggest that the central Ala21-Gly29 segment of Aβ40 can adopt folded and partially unfolded structures. At the basic pH, this segment forms folded structures stabilized by electrostatic interactions and hydrogen bonds. At the acidic pH, it forms partially unfolded structures. Two other segments flanking to the central segment exhibit the propensity to adopt unstable interconverting α-helical, 3₁₀-helical and turn-like structures. One of these segments is comprised of the Ala30-Val36 residues at both of the considered pHs. The second segment is comprised of the Glu11-Phe20 at the basic pH and of the Glu11-Val24 residues at the acidic pHs. The revealed pH-dependent structuration of the Aβ40 allowed us to suggest a possible scenario for initial Aβ aggregation. According to this scenario, the occurrence of the partially unfolded states of the Ala21-Gly29 segment plays main role in the Aβ oligomerization process.Key words: amyloid-β protein, Alzheimer disease, oligomerization, fibril, electrostatic interactions, molecular dynamics simulations     

Removal of surface charge-charge interactions from ubiquitin leaves the protein folded and very stable

The contribution of solvent-exposed charged residues to protein stability was evaluated using ubiquitin as a model protein. We combined site-directed mutagenesis and specific chemical modifications to first replace all Arg residues with Lys, followed by carbomylation of Lys-amino groups. Under the conditions in which all carboxylic groups are protonated (at pH 2), the chemically modified protein is folded and very stable (DeltaG = 18 kJ/mol). These results indicate that surface charge-charge interactions are not an essential fundamental force for protein folding and stability.