Contribution of hydrogen bonds to protein stability.

Our goal was to gain a better understanding of the contribution of the burial of polar groups and their hydrogen bonds to the conformational stability of proteins. We measured the change in stability, Δ(ΔG), for a series of hydrogen bonding mutants in four proteins: villin headpiece subdomain (VHP) containing 36 residues, a surface protein from Borrelia burgdorferi (VlsE) containing 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa (RNase Sa) and T1 (RNase T1). Crystal structures were determined for three of the hydrogen bonding mutants of RNase Sa: S24A, Y51F, and T95A. The structures are very similar to wild type RNase Sa and the hydrogen bonding partners form intermolecular hydrogen bonds to water in all three mutants. We compare our results with previous studies of similar mutants in other proteins and reach the following conclusions. (1) Hydrogen bonds contribute favorably to protein stability. (2) The contribution of hydrogen bonds to protein stability is strongly context dependent. (3) Hydrogen bonds by side chains and peptide groups make similar contributions to protein stability. (4) Polar group burial can make a favorable contribution to protein stability even if the polar groups are not hydrogen bonded. (5) The contribution of hydrogen bonds to protein stability is similar for VHP, a small protein, and VlsE, a large protein.

Cytotoxicity of RNase Sa to the acute myeloid leukemia Kasumi-1 cells depends on the net charge.

The majority of known cytotoxic RNases are basic proteins which destroy intracellular RNA. Cationization of RNases is considered to be an effective strategy for strengthening their antitumor properties. We constructed a set of RNase Sa variants consisting of charge reversal mutants, charge neutralization mutants, and variants with positively charged cluster at the N-terminus. All constructs retain a high level of catalytic activity and differ in net charge. Using acute myeloid leukemia cells Kasumi-1 we have shown that (i) cytotoxicity of RNase Sa mutants is linearly enhanced by cationization, (ii) the ability of cytotoxic mutants to induce cell death is caused by induction of apoptosis and (iii) localization of positive charge on N-terminus does not contribute to RNase Sa cytotoxicity. Capacity to induce apoptosis in malignant cells and the absence of necrotic effects make the RNase Sa mutants with high positive charge a suitable anti-cancer agent.

Toward a molecular understanding of protein solubility: increased negative surface charge correlates with increased solubility.

Protein solubility is a problem for many protein chemists, including structural biologists and developers of protein pharmaceuticals. Knowledge about how intrinsic factors influence solubility is limited due to the difficulty of obtaining quantitative solubility measurements. Solubility measurements in buffer alone are difficult to reproduce, because gels or supersaturated solutions often form, making it impossible to determine solubility values for many proteins. Protein precipitants can be used to obtain comparative solubility measurements and, in some cases, estimations of solubility in buffer alone. Protein precipitants fall into three broad classes: salts, long-chain polymers, and organic solvents. Here, we compare the use of representatives from two classes of precipitants, ammonium sulfate and polyethylene glycol 8000, by measuring the solubility of seven proteins. We find that increased negative surface charge correlates strongly with increased protein solubility and may be due to strong binding of water by the acidic amino acids. We also find that the solubility results obtained for the two different precipitants agree closely with each other, suggesting that the two precipitants probe similar properties that are relevant to solubility in buffer alone.

Contribution of hydrophobic interactions to protein stability.

Our goal was to gain a better understanding of the contribution of hydrophobic interactions to protein stability. We measured the change in conformational stability, Δ(ΔG), for hydrophobic mutants of four proteins: villin headpiece subdomain (VHP) with 36 residues, a surface protein from Borrelia burgdorferi (VlsE) with 341 residues, and two proteins previously studied in our laboratory, ribonucleases Sa and T1. We compared our results with those of previous studies and reached the following conclusions: (1) Hydrophobic interactions contribute less to the stability of a small protein, VHP (0.6±0.3 kcal/mol per -CH(2)- group), than to the stability of a large protein, VlsE (1.6±0.3 kcal/mol per -CH(2)- group). (2) Hydrophobic interactions make the major contribution to the stability of VHP (40 kcal/mol) and the major contributors are (in kilocalories per mole) Phe18 (3.9), Met13 (3.1), Phe7 (2.9), Phe11 (2.7), and Leu21 (2.7). (3) Based on the Δ(ΔG) values for 148 hydrophobic mutants in 13 proteins, burying a -CH(2)- group on folding contributes, on average, 1.1±0.5 kcal/mol to protein stability. (4) The experimental Δ(ΔG) values for aliphatic side chains (Ala, Val, Ile, and Leu) are in good agreement with their ΔG(tr) values from water to cyclohexane. (5) For 22 proteins with 36 to 534 residues, hydrophobic interactions contribute 60±4% and hydrogen bonds contribute 40±4% to protein stability. (6) Conformational entropy contributes about 2.4 kcal/mol per residue to protein instability. The globular conformation of proteins is stabilized predominantly by hydrophobic interactions.

Increasing protein stability: importance of DeltaC(p) and the denatured state.

Increasing the conformational stability of proteins is an important goal for both basic research and industrial applications. In vitro selection has been used successfully to increase protein stability, but more often site-directed mutagenesis is used to optimize the various forces that contribute to protein stability. In previous studies, we showed that improving electrostatic interactions on the protein surface and improving the beta-turn sequences were good general strategies for increasing protein stability, and used them to increase the stability of RNase Sa. By incorporating seven of these mutations in RNase Sa, we increased the stability by 5.3 kcal/mol. Adding one more mutation, D79F, gave a total increase in stability of 7.7 kcal/mol, and a melting temperature 28 degrees C higher than the wild-type enzyme. Surprisingly, the D79F mutation lowers the change in heat capacity for folding, DeltaC(p), by 0.6 kcal/mol/K. This suggests that this mutation stabilizes structure in the denatured state ensemble. We made other mutants that give some insight into the structure present in the denatured state. Finally, the thermodynamics of folding of these stabilized variants of RNase Sa are compared with those observed for proteins from thermophiles.

Factors that influence helical preferences for singly charged gas-phase peptide ions: the effects of multiple potential charge-carrying sites.

Ion mobility-mass spectrometry is used to investigate the structure(s) of a series of model peptide [M + H](+) ions to better understand how intrinsic properties affect structure in low dielectric environments. The influence of peptide length, amino acid sequence, and composition on gas-phase structure is examined for a series of model peptides that have been previously studied in solution. Collision cross sections for the [M + H](+) ions of Ac-(AAKAA)(n)Y-NH(2) (n = 3-6) and Ac-Y(AEAAKA)(n)F-NH(2) (n = 2-5) are reported and correlated with candidate structures generated using molecular modeling techniques. The [M + H](+) ions of the AAKAA peptide series each exhibit a single, dominant ion mobility arrival time distribution (ATD) which correlates to partial helical structures, whereas the [M + H](+) ions of the AEAAKA ion series are composed of ATDs which correlate to charge-solvated globules (i.e., the charge is coordinated or solvated by polar peptide functional groups). These data raise numerous questions concerning intrinsic properties (amino acid sequence and composition as well as charge location) that dictate gas-phase peptide ion structure, which may reflect trends for peptide ion structure in low dielectric environments, such as transmembrane segments.

Increasing protein stability by improving beta-turns.

Our goal was to gain a better understanding of how protein stability can be increased by improving beta-turns. We studied 22 beta-turns in nine proteins with 66-370 residues by replacing other residues with proline and glycine and measuring the stability. These two residues are statistically preferred in some beta-turn positions. We studied: Cold shock protein B (CspB), Histidine-containing phosphocarrier protein, Ubiquitin, Ribonucleases Sa2, Sa3, T1, and HI, Tryptophan synthetase alpha-subunit, and Maltose binding protein. Of the 15 single proline mutations, 11 increased stability (Average = 0.8 +/- 0.3; Range = 0.3-1.5 kcal/mol), and the stabilizing effect of double proline mutants was additive. On the basis of this and our previous work, we conclude that proteins can generally be stabilized by replacing nonproline residues with proline residues at the i + 1 position of Type I and II beta-turns and at the i position in Type II beta-turns. Other turn positions can sometimes be used if the phi angle is near -60 degrees for the residue replaced. It is important that the side chain of the residue replaced is less than 50% buried. Identical substitutions in beta-turns in related proteins give similar results. Proline substitutions increase stability mainly by decreasing the entropy of the denatured state. In contrast, the large, diverse group of proteins considered here had almost no residues in beta-turns that could be replaced by Gly to increase protein stability. Improving beta-turns by substituting Pro residues is a generally useful way of increasing protein stability.

Protein ionizable groups: pK values and their contribution to protein stability and solubility.

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

Determining the conformational stability of a protein using urea denaturation curves.

The stability of globular proteins is an important factor in determining their usefulness in basic research and medicine. A number of environmental factors contribute to the conformational stability of a protein, including pH, temperature, and ionic strength. In addition, variants of proteins may show remarkable differences in stability from their wild-type form. In this chapter, we describe the method and analysis of urea denaturation curves to determine the conformational stability of a protein. This involves relatively simple experiments that can be done in a typical biochemistry laboratory, especially when using ordinary spectroscopic techniques to follow unfolding.

A summary of the measured pK values of the ionizable groups in folded proteins.

We tabulated 541 measured pK values reported in the literature for the Asp, Glu, His, Cys, Tyr, and Lys side chains, and the C and N termini of 78 folded proteins. The majority of these values are for the Asp, Glu, and His side chains. The average pK values are Asp 3.5 +/- 1.2 (139); Glu 4.2 +/- 0.9 (153); His 6.6 +/- 1.0 (131); Cys 6.8 +/- 2.7 (25); Tyr 10.3 +/- 1.2 (20); Lys 10.5 +/- 1.1 (35); C-terminus 3.3 +/- 0.8 (22) and N-terminus 7.7 +/- 0.5 (16). We compare these results with the measured pK values of these groups in alanine pentapeptides, and comment on our overall findings.

Solvent denaturation of proteins and interpretations of the m value.

The stability of globular proteins is important in medicine, proteomics, and basic research. The conformational stability of the folded state can be determined experimentally by analyzing urea, guanidinium chloride, and thermal denaturation curves. Solvent denaturation curves in particular may give useful information about a protein such as the existence of domains or the presence of stable folding intermediates. The linear extrapolation method (LEM) for analyzing solvent denaturation curves gives the parameter m, which is a measure of the dependence of ΔG on denaturant concentration. There is much recent interest in the m value as it relates to the change in accessible surface area of a protein when it unfolds and what it may reveal about the denatured states of proteins.

Measuring and increasing protein solubility.

High concentration protein delivery is difficult to achieve for several protein pharmaceuticals due to low solubility. In this review, we discuss different types of low protein solubility, including low in vitro solubility, which is relevant to the formulation of protein pharmaceuticals. We also discuss different methods of measuring protein solubility with an emphasis on the method of inducing amorphous precipitation using ammonium sulfate. Finally, we discuss strategies for increasing protein solubility, including site-directed mutagenesis. Evidence from solubility-changing mutations in the literature indicate that some hydrophilic residues (aspartic acid, glutamic acid, and serine) contribute significantly more favorably to protein solubility than other hydrophilic residues (asparagine, glutamine, threonine, lysine, and arginine). These findings should prove useful especially in cases where protein structure is not known. In these cases, instead of targeting hydrophobic residues that are often buried, one could target hydrophilic residues that do not contribute favorably to protein solubility and replace them with hydrophilic residues that contribute more favorably.

RNase-induced apoptosis: fate of calcium-activated potassium channels.

The connection between the action of microbial RNases and Ca2+-activated K+ (KCa) channels was investigated in human embryo kidney cells HEKhSK4 artificially expressing the channels. These channels protected HEKhSK4 cells from apoptosis induced by binase and 5K charge reversal mutant of RNase Sa. After the first 24h, potassium current increased without increase in intracellular Ca2+, and mitochondrial potential remained high. After 72 h, the concentration of calcium increased and mitochondria lost their potential. Whole-cell recordings of membrane currents through KCa channels in RNase-treated cells demonstrated a biphasic pattern: initially their activity in cell population increased, peaked at 24h, and then gradually decreased. In each individual cell we observed either an increase of the amplitude of KCa current, or a complete shutdown of the channels. The activity of KCa channels could be restored by removing RNases from the media. Based on this pattern and especially its timing, we hypothesize that toxic RNases downregulate KCa channels at the level of transcription or translation. Our results indicate that new anticancer agents could be created on the basis of microbial RNases targeting KCa channels.

Tryptophan fluorescence reveals the presence of long-range interactions in the denatured state of ribonuclease Sa.

Characterizing the denatured state ensemble is crucial to understanding protein stability and the mechanism of protein folding. The aim of this research was to see if fluorescence could be used to gain new information on the denatured state ensemble. Ribonuclease Sa (RNase Sa) contains no Trp residues. We made five variants of RNase Sa by adding Trp residues at locations where they are found in other members of the microbial ribonuclease family. To better understand the protein denatured state, we also studied the fluorescence properties of the following peptides: N-acetyl-Trp-amide (NATA), N-acetyl-Ala-Trp-Ala-amide (AWA), N-acetyl-Ala-Ala-Trp-Ala-Ala-amide (AAWAA), and the five pentapeptides with the same sequence as the Trp substitution sites in RNase Sa. The major conclusions are: 1), the wavelength of maximum fluorescence intensity, lambda(max), does not differ significantly for the peptides and the denatured proteins; 2), the fluorescence intensity at lambda(max), I(F), differs significantly for the five Trp containing variants of RNase Sa; 3), the I(F) differences for the denatured proteins are mirrored in the peptides, showing that the short-range effects giving rise to the I(F) differences in the peptides are also present in the proteins; 4) the I(F) values for the denatured proteins are more than 30% greater than for the peptides, showing the presence of long-range effects in the proteins; 5), fluorescence quenching of Trp by acrylamide and iodide is more than 50% greater in the peptides than in the denatured proteins, showing that long-range effects limit the accessibility of the quenchers to the Trp side chains in the proteins; and 6), these results show that nonlocal effects in the denatured states of proteins influence Trp fluorescence and accessibility significantly.

Peptide sequence and conformation strongly influence tryptophan fluorescence.

This article probes the denatured state ensemble of ribonuclease Sa (RNase Sa) using fluorescence. To interpret the results obtained with RNase Sa, it is essential that we gain a better understanding of the fluorescence properties of tryptophan (Trp) in peptides. We describe studies of N-acetyl-L-tryptophanamide (NATA), a tripeptide: AWA, and six pentapeptides: AAWAA, WVSGT, GYWHE, HEWTV, EAWQE, and DYWTG. The latter five peptides have the same sequence as those surrounding the Trp residues studied in RNase Sa. The fluorescence emission spectra, the fluorescence lifetimes, and the fluorescence quenching by acrylamide and iodide were measured in concentrated solutions of urea and guanidine hydrochloride. Excited-state electron transfer from the indole ring of Trp to the carbonyl groups of peptide bonds is thought to be the most important mechanism for intramolecular quenching of Trp fluorescence. We find the maximum fluorescence intensities vary from 49,000 for NATA with two carbonyls, to 24,400 for AWA with four carbonyls, to 28,500 for AAWAA with six carbonyls. This suggests that the four carbonyls of AWA are better able to quench Trp fluorescence than the six carbonyls of AAWAA, and this must reflect a difference in the conformations of the peptides. For the pentapeptides, EAWQE has a fluorescence intensity that is more than 50% greater than DYWTG, showing that the amino acid sequence influences the fluorescence intensity either directly through side-chain quenching and/or indirectly through an influence on the conformational ensemble of the peptides. Our results show that peptides are generally better models for the Trp residues in proteins than NATA. Finally, our results emphasize that we have much to learn about Trp fluorescence even in simple compounds.

Increasing protein conformational stability by optimizing beta-turn sequence.

Protein conformational stability is an important concern in many fields. Here, we describe a strategy for significantly increasing conformational stability by optimizing beta-turn sequence. Proline and glycine residues are statistically preferred at several beta-turn positions, presumably because their unique side-chains contribute favorably to conformational stability in certain beta-turn positions. However, beta-turn sequences often deviate from preferred proline or preferred glycine. Therefore, our strategy involves replacing non-proline and non-glycine beta-turn residues with preferred proline or preferred glycine residues. Here, we develop guidelines for selecting appropriate mutations, and present results for five mutations (S31P, S42G, S48P, T76P, and Q77G) that significantly increase the conformational stability of RNase Sa. The increases in stability ranged from 0.7 kcal/mol to 1.3 kcal/mol. The strategy was successful in overlapping or isolated beta-turns, at buried (up to 50%) or completely exposed sites, and at relatively flexible or inflexible sites. Considering the significant number of beta-turn residues in every globular protein and the frequent deviation of beta-turn sequences from preferred proline and preferred glycine residues, this simple, efficient strategy will be useful for increasing the conformational stability of proteins.

Amino acid contribution to protein solubility: Asp, Glu, and Ser contribute more favorably than the other hydrophilic amino acids in RNase Sa.

Poor protein solubility is a common problem in high-resolution structural studies, formulation of protein pharmaceuticals, and biochemical characterization of proteins. One popular strategy to improve protein solubility is to use site-directed mutagenesis to make hydrophobic to hydrophilic mutations on the protein surface. However, a systematic investigation of the relative contributions of all 20 amino acids to protein solubility has not been done. Here, 20 variants at the completely solvent-exposed position 76 of ribonuclease (RNase) Sa are made to compare the contributions of each amino acid. Stability measurements were also made for these variants, which occur at the i+1 position of a type II beta-turn. Solubility measurements in ammonium sulfate solutions were made at high positive net charge, low net charge, and high negative net charge. Surprisingly, there was a wide range of contributions to protein solubility even among the hydrophilic amino acids. The results suggest that aspartic acid, glutamic acid, and serine contribute significantly more favorably than the other hydrophilic amino acids especially at high net charge. Therefore, to increase protein solubility, asparagine, glutamine, or threonine should be replaced with aspartic acid, glutamic acid or serine.

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.