Electrostatic contributions to T4 lysozyme stability: solvent-exposed charges versus semi-buried salt bridges

Effects of Charge-to-Alanine Substitutions on the Stability of Ribosomal Protein L30e from Thermococcus celer

The ability to rationally engineer a protein with altered stability depends upon the detailed understanding of the role of noncovalent interactions in defining thermodynamic properties of proteins. In this paper, we used T. celer L30e as a model to address the question of the role of charge-charge interactions in defining the stability of this protein. A total of 26 single-site charge-to-alanine variants of this protein were generated, and the stability of these proteins was determined using thermal- and denaturant-induced unfolding. It was found that, although L30e is isolated from a thermophilic organism and is highly thermostable, some of the substitutions lead to a further increase in the transition temperature. Analysis of the effects of high ionic strength on the stabilities of L30e variants shows that the long-range charge-charge interactions are as important as the short-range (salt bridge) interactions. The changes in stabilities of the T. celer L30e protein variants were compared with the changes in the energy of charge-charge interactions calculated using different computational models. It was found that there is a good qualitative agreement between experimental and calculated data: for 70-80% (19-21 of 26, confidence p < 0.003) of the variants, computational models predict correctly the sign of the stability changes. In particular, computational models identify correctly those charged amino acid residue substitutions of which led to enhancement in thermostability. Thus, optimization of the charge-charge interactions might be a useful approach for the rational increase in protein stability.

How to improve nature: study of the electrostatic properties of the surface of alpha-lactalbumin

It was recently shown that alpha-lactalbumin interacts with histones and simple models of histone proteins such as positively charged polyamino acids, suggesting that some fundamental aspects of the protein surface electrostatics may come into play. In the present work, the energies of charge-charge interaction in apo- and Ca(2+)-loaded alpha-lactalbumin were calculated using a Tanford-Kirkwood algorithm with either solvent accessibility correction or using a finite difference Poisson-Boltzmann method. The analysis revealed two major regions of alpha-lactalbumin that possessed highly unfavorable electrostatic potentials: (a) the Ca(2+)-binding loop and its neighboring residues and (b) the N-terminal region of the protein. Several individual mutants were prepared to neutralize specific individual surface acidic amino acids at both the N-terminus and Ca(2+)-binding loop of bovine alpha-lactalbumin. These mutants were characterized by intrinsic fluorescence, differential scanning microcalorimetry and circular dichroism. The structural and thermodynamic data agree in every case with the theoretical predictions, confirming that the N-terminal region is very sensitive to changes in charge. For example, desMet D14N mutation destabilizes protein and decreases its calcium affinity. On the other hand, desMet E1M and desMet D37N substitutions increase the thermal stability and calcium affinity. The Met E1Q is characterized by a marked increase in protein stability, whereas desMet E7Q and desMet E11L display a slight increase in calcium affinity and thermal stability. Examination of the unfavorable energy contributed by Glu1 and the energetically favorable consequences of neutralizing this residue suggests that nature may have made an error with bovine alpha-lactalbumin from the viewpoint of stabilizing structure and conformation.

Electrostatic Recognition and Induced Fit in the κ-PVIIA Toxin Binding to Shaker Potassium Channel

Brownian dynamics (BD) and molecular dynamics (MD) simulations and electrostatic calculations were performed to study the binding process of kappa-PVIIA to the Shaker potassium channel and the structure of the resulting complex. BD simulations, guided by electrostatic interactions, led to an initial alignment between the toxin and the channel protein. MD simulations were then carried out to allow for rearrangements from this initial structure. After approximately 4 ns, a critical "induced fit" process was observed to last for approximately 2 ns. In this process, the interface was reorganized, and side chains were moved so that favorable atomic contacts were formed or strengthened, while unfavorable contacts were eliminated. The final complex structure was stabilized through electrostatic interactions with the positively charged side chain of Lys7 of kappa-PVIIA deeply inserted into the channel pore and other hydrogen bonds and by hydrophobic interactions involving Phe9 and Phe23 of the toxin. The validity of the predicted structure for the complex was assessed by calculating the effects of mutating charged and polar residues of both the toxin and the channel protein, with the calculated effects correlating reasonably well with experimental data. The present study suggests a general binding mechanism, whereby proteins are pre-aligned in their diffusional encounter by long-range electrostatic attraction, and nanosecond-scale rearrangements within the initial complex then lead to a specifically bound complex.

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.

How optimal are the binding energetics of barnase and barstar?

The extracellular ribonuclease barnase and its intracellular inhibitor barstar bind fast and with high affinity. Although extensive experimental and theoretical studies have been carried out on this system, it is unclear what the relative importance of different contributions to the high affinity is and whether binding can be improved through point mutations. In this work, we first applied Poisson-Boltzmann electrostatic calculations to 65 barnase-barstar complexes with mutations in both barnase and barstar. The continuum electrostatic calculations with a van der Waals surface dielectric boundary definition result in the electrostatic interaction free energy providing the dominant contribution favoring barnase-barstar binding. The results show that the computed electrostatic binding free energy can be improved through mutations at W44/barstar and E73/barnase. Furthermore, the determinants of binding affinity were quantified by applying COMparative BINding Energy (COMBINE) analysis to derive quantitative structure-activity relationships (QSARs) for the 65 complexes. The COMBINE QSAR model highlights approximately 20 interfacial residue pairs as responsible for most of the differences in binding affinity between the mutant complexes, mainly due to electrostatic interactions. Based on the COMBINE model, together with Brownian dynamics simulations to compute diffusional association rate constants, several mutants were designed to have higher binding affinities than the wild-type proteins.

Theoretical and experimental study of the D2194G mutation in the C2 domain of coagulation factor V

Coagulation factor V (FV) is a large plasma glycoprotein with functions in both the pro- and anticoagulant pathways. In carriers of the so-called R2-FV haplotype, the FV D2194G mutation, in the C2 membrane-binding domain, is associated with low expression levels, suggesting a potential folding/stability problem. To analyze the molecular mechanisms potentially responsible for this in vitro phenotype, we used molecular dynamics (MD) and continuum electrostatic calculations. Implicit solvent simulations were performed on the x-ray structure of the wild-type C2 domain and on a model of the D2194G mutant. Because D2194 is located next to a disulfide bond (S-S bond), MD calculations were also performed on S-S bond depleted structures. D2194 is part of a salt-bridge network and investigations of the stabilizing/destabilizing role of these ionic interactions were carried out. Five mutant FV molecules were created and the expression levels measured with the aim of assessing the tolerance to amino acid changes in this region of molecule. Analysis of the MD trajectories indicated increased flexibility in some areas and energetic comparisons suggested overall destabilization of the structure due to the D2194G mutation. This substitution causes electrostatic destabilization of the domain by approximately 3 kcal/mol. Together these effects likely explain the lowered expression levels in R2-FV carriers.

The variant E233G of the RAD51D gene could be a low-penetrance allele in high-risk breast cancer families without BRCA1/2 mutations

Six SNPs have been detected in the DNA repair genes RAD51C and RAD51D, not previously characterized. The novel variant E233G in RAD51D is more highly represented in high-risk, site-specific, familial breast cancer cases that are not associated with the BRCA1/2 genes, with a frequency of 5.74% (n = 174) compared to a control population (n = 567) and another subset of breast cancer patients (n = 765) with a prevalence of around 2% only (comparison to controls, OR = 2.6, 95% CI 1.12-6.03; p < 0.021). We found that the immunohistochemical profile detected in available tumors from these patients differs slightly from those described in non-BRCA1/2 tumors. Finally, the structural prediction of the putative functional consequence of this change indicates that it can diminish protein stability and structure. This suggests a role for E233G as a low-penetrance susceptibility gene in the specific subgroup of high-risk familial breast cancer cases that are not related to BRCA1/2.

Calculating proton uptake/release and binding free energy taking into account ionization and conformation changes induced by protein-inhibitor association: Application to plasmepsin, cathepsin D and endothiapepsin-pepstatin complexes

The protein-inhibitor binding energies of enzymes are often pH dependent, and binding induces either proton uptake or proton release. The proton uptake/release and the binding energy for three complexes with available experimental data were numerically studied: pepstatin-cathepsin D, pepstatin-plasmepsin II and pepstatin-endothiapepsin. Very good agreement with the experimental data was achieved when conformational changes were taken into account. The role of the desolvation energy and the conformational changes was revealed by modeling the complex, the separated molecules in the complex conformation and the free molecules. It was shown that the conformational changes induced by the complex formation are as important for the proton transfer as the loss of solvation energy caused by the burial of interface residues. The residues responsible for the proton transfer were identified and their contribution to the proton uptake/release calculated. These residues were found to be scattered along the whole protein rather than being localized only at the active site. In the case of cathepsin D, these residues were found to be highly conserved among the cathepsin D sequences of other species. It was shown that conformation and ionization changes induced by the complex formation are critical for the correct calculation of the binding energy. Taking into account the electrostatics and the van der Waals (vdW) energies within the Boltzmann distribution of energies and allowing ionization and conformation changes to occur makes the calculated binding energy more realistic and closer to the experimental value. The interplay between electrostatic and vdW forces makes the pH dependence of the binding energy smoother, because the vdW force acts in reaction to the changes of the electrostatic energy. It was found that a small fraction of the ionizable groups remain uncharged in both the free and complexed molecules. The sequence and structural position of these groups aligns well within the three proteases, suggesting that these may have specific role.

Reversal of Negative Charges on the Surface of Escherichia coli Thioredoxin: Pockets versus Protrusions

Recent studies of proteins with reversed charged residues have demonstrated that electrostatic interactions on the surface can contribute significantly to protein stability. We have used the approach of reversing negatively charged residues using Arg to evaluate the effect of the electrostatics context on the transition temperature (T(m)), the unfolding Gibbs free energy change (DeltaG), and the unfolding enthalpy change (DeltaH). We have reversed negatively charged residues at a pocket (Asp9) and protrusions (Asp10, Asp20, Glu85), all located in interconnecting segments between elements of secondary structure on the surface of Arg73Ala Escherichia coli thioredoxin. DSC measurements indicate that reversal of Asp in a pocket (Asp9Arg/Arg73Ala, DeltaT(m) = -7.3 degrees C) produces a larger effect in thermal stability than reversal at protrusions: Asp10Arg/Arg73Ala, DeltaT(m) = -3.1 degrees C, Asp20Arg/Arg73Ala, DeltaT(m) = 2.0 degrees C, Glu85Arg/Arg73Ala, DeltaT(m) = 3.9 degrees ). The 3D structure of thioredoxin indicates that Asp20 and Glu85 have no nearby charges within 8 A, while Asp9 does not only have Asp10 as sequential neighbor, but it also forms a 5-A long-range ion pair with the solvent-exposed Lys69. Further DSC measurements indicate that neutralization of the individual charges of the ion pair Asp9-Lys69 with nonpolar residues produces a significant decrease in stability in both cases: Asp9Ala/Arg73Ala, DeltaT(m) = -3.7 degrees C, Asp9Met/Arg73Ala, DeltaT(m) = -5.5 degrees C, Lys69Leu/Arg73Ala, DeltaT(m) = -5.1 degrees C. However, thermodynamic analysis shows that reversal or neutralization of Asp9 produces a 9-15% decrease in DeltaH, while both reversal of Asp at protrusions and neutralization of Lys69 produce negligible changes. These results correlate well with the NMR analysis, which demonstrates that only the substitution of Asp9 produces extensive conformational changes and these changes occur in the surroundings of Lys69. Our results led us to suggest that reversal of a negative charge at a pocket has a larger effect on stability than a similar reversal at a protrusion and that this difference arises largely from short-range interactions with polar groups within the pocket, rather than long-range interactions with solvent-exposed charged groups.

Numerical calculations of the pH of maximal protein stability. The effect of the sequence composition and three-dimensional structure

A large number of proteins, found experimentally to have different optimum pH of maximal stability, were studied to reveal the basic principles of their preference for a particular pH. The pH-dependent free energy of folding was modeled numerically as a function of pH as well as the net charge of the protein. The optimum pH was determined in the numerical calculations as the pH of the minimum free energy of folding. The experimental data for the pH of maximal stability (experimental optimum pH) was reproducible (rmsd = 0.73). It was shown that the optimum pH results from two factors - amino acid composition and the organization of the titratable groups with the 3D structure. It was demonstrated that the optimum pH and isoelectric point could be quite different. In many cases, the optimum pH was found at a pH corresponding to a large net charge of the protein. At the same time, there was a tendency for proteins having acidic optimum pHs to have a base/acid ratio smaller than one and vice versa. The correlation between the optimum pH and base/acid ratio is significant if only buried groups are taken into account. It was shown that a protein that provides a favorable electrostatic environment for acids and disfavors the bases tends to have high optimum pH and vice versa.

GroEL stability and function. Contribution of the ionic interactions at the inter-ring contact sites

The chaperonin GroEL consists of a double ring structure made of identical subunits that display different modes of allosteric communication. The protein folding cycle requires the simultaneous positive intra-ring and negative inter-ring cooperativities of ATP binding. This ensures GroES binding to one ring and release of the ligands from the opposite one. To better characterize inter-ring allosterism, the thermal stability as well as the temperature dependence of the functional and conformational properties of wild type GroEL, a single ring mutant (SR1) and two single point mutants suppressing one interring salt bridge (E434K and E461K) were studied. The results indicate that ionic interactions at the two interring contact sites are essential to maintain the negative cooperativity for protein substrate binding and to set the protein thermostat at 39 degrees C. These electrostatic interactions contribute distinctly to the stability of the inter-ring interface and the overall protein stability, e.g. the E434K thermal inactivation curve is shifted to lower temperatures, and its unfolding temperature and activation energy are also lowered. An analysis of the ionic interactions at the inter-ring contact sites reveals that at the so called "left site" a network of electrostatic interactions involving three charged residues might be established, in contrast to what is found at the "right site" where only two oppositely charged residues interact. Our data suggest that electrostatic interactions stabilize protein-protein interfaces depending on both the number of ionic interactions and the number of residues engaged in each of these interactions. In the case of GroEL, this combination sets the thermostat of the protein so that the chaperonin distinguishes physiological from stress temperatures.

Comparison of calculation and experiment implicates significant electrostatic contributions to the binding stability of barnase and barstar

The contributions of electrostatic interactions to the binding stability of barnase and barstar were studied by the Poisson-Boltzmann model with three different protocols: a), the dielectric boundary specified as the van der Waals (vdW) surface of the protein along with a protein dielectric constant (epsilon (p)) of 4; b), the dielectric boundary specified as the molecular (i.e., solvent-exclusion (SE)) surface along with epsilon (p) = 4; and c), "SE + epsilon (p) = 20." The "vdW + epsilon (p) = 4" and "SE + epsilon (p) = 20" protocols predicted an overall electrostatic stabilization whereas the "SE + epsilon (p) = 4" protocol predicted an overall electrostatic destabilization. The "vdW + epsilon (p) = 4" protocol was most consistent with experiment. It quantitatively reproduced the observed effects of 17 mutations neutralizing charged residues lining the binding interface and the measured coupling energies of six charge pairs across the interface and reasonably rationalized the experimental ionic strength and pH dependences of the binding constant. In contrast, the "SE + epsilon (p) = 4" protocol predicted significantly larger coupling energies of charge pairs whereas the "SE + epsilon (p) = 20" protocol did not predict any pH dependence. This study calls for further scrutiny of the different Poisson-Boltzmann protocols and demonstrates potential danger in drawing conclusions on electrostatic contributions based on a particular calculation protocol.

High-resolution structures of retinol-binding protein in complex with retinol: pH-induced protein structural changes in the crystal state

The targeted delivery of non-polar ligands by binding proteins to membranes or membrane receptors involves the release of these ligands on or near the plasma membrane of target cells. Because these hydrophobic ligands are often bound inside a deep cavity of binding proteins, as shown previously for plasma retinol-binding protein (RBP), their release from these proteins might require the destabilization of the protein structure by partially denaturing conditions, such as those possibly present near plasma membranes. RBP is a plasma transport protein which delivers specifically retinol from its store sites to target cells. Here, we report the high-resolution (1.1-1.4A) crystal structures of bovine holo-RBP at five different pH values, ranging from 9 to 2. While unraveling details of the native protein structure and of the interactions with retinol at nearly atomic resolution at neutral pH, this study provides evidence for definite pH-induced modifications of several structural features of RBP. The structure most representative of the changes that holo-RBP undergoes at different pH values is that of its flexible state at pH 2. At this pH, most significant are the alteration of the arrangement of salt bridges and of the network of water molecules/H-bonds that participates in the retinol-RBP interaction, an appreciable increase of the volume of the beta-barrel cavity, a considerably higher degree of mobility of the RBP-bound ligand and of several protein regions and the disorder of a large number of solvent molecules that are ordered at neutral pH. These changes are likely to be accompanied by a modification of the pattern of charge distribution on the protein surface. All these changes, which reveal a substantially lowered conformational stability of RBP, presumably occur at the initial stages of the acidic denaturation of RBP and are possibly associated with a facilitated release of the retinol molecule from its carrier protein.

Site-specific contributions to the pH dependence of protein stability

Understanding protein stability is a significant challenge requiring characterization of interactions within both folded and unfolded states. Of these, electrostatic interactions influence ionization equilibria of acidic and basic groups and diversify their pK(a) values. The pH dependence of the thermodynamic stability (Delta G(FU)) of a protein arises as a consequence of differential pK(a) values between folded and unfolded states. Previous attempts to calculate pH-dependent contributions to stability have been limited by the lack of experimental unfolded state pK(a) values. Using recently developed NMR spectroscopic methods, we have determined residue-specific pK(a) values for a thermodynamically unstable Src homology 3 domain in both states, enabling the calculation of the pH dependence of stability based on simple analytical expressions. The calculated pH stability profile obtained agrees very well with experiment, unlike profiles derived from two current models of electrostatic interactions within unfolded states. Most importantly, per-residue contributions to the pH dependence of Delta G(FU) derived from the data provide insights into specific electrostatic interactions in both the folded and unfolded states and their roles in protein stability. These interactions include a hydrogen bond between the Asp-8 side-chain and the Lys-21 backbone amide group in the folded state, which represents a highly conserved interaction in Src homology 3 domains.

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.

Electrostatic contributions to the stability of a thermophilic cold shock protein

The thermophilic Bacillus caldolyticus cold shock protein (Bc-Csp) differs from the mesophilic Bacillus subtilis cold shock protein B (Bs-CspB) in 11 of the 66 residues. Stability measurements of Schmid and co-workers have implicated contributions of electrostatic interactions to the thermostability. To further elucidate the physical basis of the difference in stability, previously developed theoretical methods that treat electrostatic effects in both the folded and the unfolded states were used in this paper to study the effects of mutations, ionic strength, and temperature. For 27 mutations that narrow the difference in sequence between Bc-Csp and Bs-CspB, calculated changes in unfolding free energy (Delta G) and experimental results have a correlation coefficient of 0.98. Bc-Csp appears to use destabilization of the unfolded state by unfavorable charge-charge interactions as a mechanism for increasing stability. Accounting for the effects of ionic strength and temperature on the electrostatic free energies in both the folded and the unfolded states, explanations for two important experimental observations are presented. The disparate ionic strength dependences of Delta G for Bc-Csp and Bs-CspB were attributed to the difference in the total charges (-2e and -6e, respectively). A main contribution to the much higher unfolding entropy of Bs-CspB was found to come from the less favorable electrostatic interactions in the folded state. These results should provide insight for understanding the thermostability of other thermophilic proteins.

Changes in stability upon charge reversal and neutralization substitution in staphylococcal nuclease are dominated by favorable electrostatic effects

Single site mutations that reverse or neutralize a surface charge were made at 22 ionizable residues in staphylococcal nuclease. Unfolding free energies were obtained by guanidine hydrochloride denaturation. These data, in conjunction with previously obtained stabilities of the corresponding alanine mutants, unequivocally show that the dominant contribution to stability for virtually all of the wild-type side chains examined is the electrostatic effect associated with each residue's charged group. With only a few exceptions, these charges stabilize the native state, with an average loss of 0.5 kcal/mol of stability upon neutralization of a charge. When the charge is reversed, the average destabilization is doubled. Structure-based calculations of electrostatic free energy with the continuum method based on the finite difference solution to the linearized Poisson-Boltzmann equation reproduce the observed energetics when the polarizability in the protein interior is represented with a dielectric constant of 20. However, in some cases, large differences are found, giving insight into possible areas for improvement of the calculations. In particular, it appears that the assumptions made in the calculations about the absence of electrostatic interactions in the denatured state and the energetic consequences of dynamic fluctuations in the native state will have to be further explored.