The effect of increasing the stability of non-native interactions on the folding landscape of the bacterial immunity protein Im9

How stabilising non-native interactions influence protein folding energy landscapes is currently not well understood: such interactions could speed folding by reducing the conformational search to the native state, or could slow folding by increasing ruggedness. Here, we examine the influence of non-native interactions in the folding process of the bacterial immunity protein Im9, by exploiting our ability to manipulate the stability of the intermediate and rate-limiting transition state (TS) in the folding of this protein by minor alteration of its sequence or changes in solvent conditions. By analysing the properties of these species using Phi-value analysis, and exploration of the structural properties of the TS ensemble using molecular dynamics simulations, we demonstrate the importance of non-native interactions in immunity protein folding and demonstrate that the rate-limiting step involves partial reorganisation of these interactions as the TS ensemble is traversed. Moreover, we show that increasing the contribution to stability made by non-native interactions results in an increase in Phi-values of the TS ensemble without altering its structural properties or solvent-accessible surface area. The data suggest that the immunity proteins fold on multiple, but closely related, micropathways, resulting in a heterogeneous TS ensemble that responds subtly to mutation or changes in the solvent conditions. Thus, altering the relative strength of native and non-native interactions influences the search to the native state by restricting the pathways through the folding energy landscape.     

The effect of protein stability on protein-monoglyceride interactions

We have previously shown that proteins such as beta-lactoglobulin and lysozyme insert into monoglyceride monolayers and are able to induce an L(beta) to coagel phase transition in monoglyceride bilayers. These studies gave a first indication that protein stability could be an important factor for these interactions. This study therefore aims at further investigating the potential role of protein stability on protein-monoglyceride interactions. To this end we studied the interaction of stable and destabilized alpha-lactalbumin with monostearoylglycerol. Our results show that protein stability is important for the insertion of proteins into a monostearoylglycerol monolayer, such that the lower the stability of the protein the better the protein inserts. In marked contrast to beta-lactoglobulin and lysozyme we found that destabilized alpha-lactalbumin does not induce the L(beta) to coagel phase transition in monoglyceride bilayers. We propose that this is due to an increased surface coverage by the protein which could result from the unfolding of the protein upon binding to the interface.     

Acidic conditions stabilise intermediates populated during the folding of Im7 and Im9

The helical bacterial immunity proteins Im7 and Im9 have been shown to fold via kinetic mechanisms of differing complexity, despite having 60 % sequence identity. At pH 7.0 and 10 degrees C, Im7 folds in a three-state mechanism involving an on-pathway intermediate, while Im9 folds in an apparent two-state transition. In order to examine the folding mechanisms of these proteins in more detail, the folding kinetics of both Im7 and Im9 (at 10 degrees C in 0.4 M sodium sulphate) have been examined as a function of pH. Kinetic modelling of the folding and unfolding data for Im7 between pH 5.0 and 8.0 shows that the on-pathway intermediate is stabilised by more acidic conditions, whilst the native state is destabilised. The opposing effect of pH on the stability of these states results in a significant population of the intermediate at equilibrium at pH 6.0 and below. At pH 7.0, the folding and unfolding kinetics for Im9 can be fitted adequately by a two-state model, in accord with previous results. However, under acidic conditions there is a clear change of slope in the plot of the logarithm of the folding rate constant versus denaturant concentration, consistent with the population of one or more intermediate(s) early during folding. The kinetic data for Im9 at these pH values can be fitted to a three-state model, where the intermediate ensemble is stabilised and the native state destabilised as the pH is reduced, rationalising previous results that showed that an intermediate is not observed experimentally at pH 7.0. The data suggest that intermediate formation is a general step in immunity protein folding and demonstrate that it is necessary to explore a wide range of refolding conditions in order to show that intermediates do not form in the folding of other small, single-domain proteins.     

The effect of a side chain-backbone swap on protein stability.

To assess the relative importance of backbone hydrogen bonding (H-bonding) vs. side chain hydrophobicity in protein structural formation, a method called side chain-backbone swap is proposed. Such a method swaps the side chain and backbone portions of certain amino acid residues, such as Asp, Glu, Asn, Gln, Lys, and Arg. Such a swap retains the sequence of a polypeptide and preserves the identity of the backbone linkage. On the other hand, the swap disrupts backbone H-bonding geometry because of the introduction of extra methylene groups into the peptide backbone. In this project, we chose the two-stranded alpha-helical coiled-coil to implement side chain-backbone swap. A pair of 36-residue peptides was designed. The two peptides have identical sequence with four residues in each heptad repeat occupied by glutamyl residues. Each glutamic acid was incorporated either as alpha-glutamyl residue (the peptide is denoted as alpha-Glu-36) or as gamma-glutamyl residue (the peptide is denoted as gamma-Glu-36). The inter-conversion between the two peptides constitutes a side chain-backbone swap. Residues constituting the hydrophobic core of the coiled-coil, however, are left unchanged. The peptide pair was characterized by circular dichroism spectroscopy, reversed-phase liquid chromatography (RPLC), and two-dimensional nuclear magnetic resonance (NMR). The results indicate that alpha-Glu-36 is a two-stranded alpha-helical coiled-coil while gamma-Glu-36 lacks stable structural elements. It is concluded that, at least for coiled-coils where hydrophobic interactions are predominantly long-range, local backbone H-bonding is a required for structural formation, consistent with a hierarchic folding mechanism. The methodological implication of side chain-backbone swap is also discussed.     

Configurational entropy modulates the mechanical stability of protein GB1

Configurational entropy plays important roles in defining the thermodynamic stability as well as the folding/unfolding kinetics of proteins. Here we combine single-molecule atomic force microscopy and protein engineering techniques to directly examine the role of configurational entropy in the mechanical unfolding kinetics and mechanical stability of proteins. We used a small protein, GB1, as a model system and constructed four mutants that elongate loop 2 of GB1 by 2, 5, 24 and 46 flexible residues, respectively. These loop elongation mutants fold properly as determined by far-UV circular dichroism spectroscopy, suggesting that loop 2 is well tolerant of loop insertions without affecting GB1′s native structure. Our single-molecule atomic force microscopy results reveal that loop elongation decreases the mechanical stability of GB1 and accelerates the mechanical unfolding kinetics. These results can be explained by the loss of configurational entropy upon closing an unstructured flexible loop using classical polymer theory, highlighting the important role of loop regions in the mechanical unfolding of proteins. This study not only demonstrates a general approach to investigating the structural deformation of the loop regions in mechanical unfolding transition state, but also provides the foundation to use configurational entropy as an effective means to modulate the mechanical stability of proteins, which is of critical importance towards engineering artificial elastomeric proteins with tailored nanomechanical properties.     

Ligand-binding alters the mechanical stability of muscle protein titin domains

Abstract

Communicated by Ramaswamy H. Sarma     

The effects of drug complexation on the stability and conformation of human serum albumin

We report different analytical methods used to study the effects of 3\'-azido-3\'-deoxythymidine, aspirin, taxol, cisplatin, atrazine, 2,4-dichlorophenoxyacetic, biogenic polyamines, chlorophyll, chlorophyllin, poly(ethylene glycol), vanadyl cation, vanadate anion, cobalt-hexamine cation, and As2O3, on the stability and secondary structure of human serum albumin (HSA) in aqueous solution, using capillary electrophoresis, Fourier transform infrared, ultraviolet visible, and circular dichroism (CD) spectroscopic methods. The concentrations of HSA used were 4% to 2% or 0.6 to 0.3 mM, while different ligand concentrations were 1 microM to 1 mM. Structural data showed drugs are mostly located along the polypeptide chains with both specific and nonspecific interactions. The stability of drug-protein complexes were in the order K(VO(2+)) 1.2 x 10(8) M(-1) > K(AZT) 1.9 x 10(6) M(-)1 > K(PEG) 4.1 x 10(5) M(-1) > K(atrazine) 3.5 x 10(4) M(-1) > K(chlorophyll) 2.9 x 10(4) M(-1) > K2,4-D 2.5 x 10(4) M-1 > K(spermine) 1.7 x 10(4) M(-1) > K(taxol) 1.43 x 10(4) M(-1) > K(Co(3+)) > 1.1 x 10(4) M(-1) > K(aspirin) 1.04 x 10(4)i(-1) > K(chlorophyllin) 7.0 x 10(3) M(-1) > K(VO(3)(-)) 6.0 x 103 M(-1) > K(spermidine) 5.4 x 10(3) M(-1) > K(putrescine) 3.9 x 10(3) M(-1) > K(As(2)O(3)) 2.2 x 10(3) M(-1)> K(cisplatin) 1.2 x 10(2) M(-1). The protein conformation was altered (infrared and CD results) with major reduction of alpha-helix from 60 to 55% (free HSA) to 49 to 40% and increase of beta-structure from 22 to 15% (free HSA) to 33 to 23% in the drug-protein complexes. The alterations of protein secondary structure are attributed to a partial unfolding of HSA on drug complexation.     

Rapid folding with and without populated intermediates in the homologous four-helix proteins Im7 and Im9

The kinetics and thermodynamics of the folding of the homologous four-helix proteins Im7 and Im9 have been characterised at pH 7.0 and 10 degrees C. These proteins are 60 % identical in sequence and have the same three-dimensional structure, yet appear to fold by different kinetic mechanisms. The logarithm of the folding and unfolding rates of Im9 change linearly as a function of urea concentration and fit well to an equation describing a two-state mechanism (with a folding rate of 1500 s-1, an unfolding rate of 0. 01 s-1, and a highly compact transition state that has approximately 95 % of the native surface area buried). By contrast, there is clear evidence for the population of an intermediate during the refolding of Im7, as indicated by a change in the urea dependence of the folding rate and the presence of a significant burst phase amplitude in the refolding kinetics. Under stabilising conditions (0.25 M Na2SO4, pH 7.0 and 10 degrees C) the folding of Im9 remains two-state, whilst under similar conditions (0.4 M Na2SO4, pH 7.0 and 10 degrees C) the intermediate populated during Im7 refolding is significantly stabilised (KUI=125). Equilibrium denaturation experiments, under the conditions used in the kinetic measurements, show that Im7 is significantly less stable than Im9 (DeltaDeltaG 9.3 kJ/mol) and the DeltaG and m values determined accord with those obtained from the fit to the kinetic data. The results show, therefore, that the population of an intermediate in the refolding of the immunity protein structure is defined by the precise amino acid sequence rather than the global stability of the protein. We discuss the possibility that the intermediate of Im7 is populated due to differences in helix propensity in Im7 and Im9 and the relevance of these data to the folding of helical proteins in general.     

The effect of the 540-kilodalton actin cross-linking protein, actin-binding protein, on the mechanical properties of F-actin

This study describes the effect of actin-binding protein derived from rabbit lung macrophages on the mechanical properties of F-actin. The dynamic storage modulus, G'(omega), and loss modulus, G"(omega) of F-actin, at concentrations from 1 to 4 mg/ml, in the absence or presence of actin-binding protein at molar ratios to actin of 1:1000 to 1:125, were measured at frequencies ranging from 3 X 10(-3) to 0.5 Hz. Actin-binding protein increased the dynamic moduli of F-actin, but this increase was much greater as either the actin-binding protein/actin ratio or the total protein concentration increased. Moreover, there was a convergence of the values of G' and G" at high frequencies for F-actin which became more prominent upon the addition of actin-binding protein. The value of the modulus obtained by an extrapolation of these data to actin concentrations similar to that found in the cell cortex was close to values which have been obtained by direct measurements. The addition of actin-binding protein to an F-actin solution enabled it to reach an equilibrium strain following the application of a stress, in contrast to pure F-actin. These data allow a more rigorous definition of the "sol" to "gel" transition and suggest that the cross-linking of actin filaments by actin-binding protein leads to the formation of a network structure whose underlying mechanism of mechanical behavior is short range intrafilament bending in contrast to the classical rubber network.     

Hydrophobic effect on the stability and folding of a hyperthermophilic protein

Ribonuclease HII from hyperthermophile Thermococcus kodakaraensis (Tk-RNase HII) is a kinetically robust monomeric protein. The conformational stability and folding kinetics of Tk-RNase HII were measured for nine mutant proteins in which a buried larger hydrophobic side chain is replaced by a smaller one (Leu/Ile to Ala). The mutant proteins were destabilized by 8.9 to 22.0 kJ mol^− 1 as compared with the wild-type protein. The removal of each -CH₂- group burial decreased the stability by 5.1 kJ mol^− 1 on average in the mutant proteins of Tk-RNase HII examined. This is comparable with the value of 5.3 kJ mol^− 1 obtained from experiments for proteins from organisms growing at moderate temperature. We conclude that the hydrophobic residues buried inside protein molecules contribute to the stabilization of hyperthermophilic proteins to a similar extent as proteins at normal temperature. In the folding experiments, the mutant proteins of Tk-RNase HII examined exhibited faster unfolding compared with the wild-type protein. These results indicate that the buried hydrophobic residues strongly contribute to the kinetic robustness of Tk-RNase HII. This is the first report that provides a practical cause of slow unfolding of hyperthermostable proteins.     

Mutational and structural-based analyses of the osmolyte effect on protein stability

It is known that several naturally occurring substances known as osmolytes increase the conformational stability of proteins. Bolen and co-worker proposed the osmophobic theory, which asserts the osmolyte effect occurs because of an unfavorable interaction of osmolytes mainly with the protein backbone, based on the results on the transfer Gibbs energy of amino acids (Deltag) [Bolen and Baskakov (2001) J. Mol. Biol. 310, 955-963]. In this paper, we report the effect of sarcosine on the conformational stability (DeltaG) of RNase Sa (96 residues and one disulfide bond) and four mutant proteins. The thermal denaturation curves for RNase Sa in sarcosine fitted a two-state model on nonlinear least-squares analysis. All the RNase Sa proteins were stabilized by sarcosine. For example, the increase in stability of the wild-type protein in 4 M sarcosine due to the osmolyte effect (Delta(o)DeltaG) is 3.2 kcal/mol. Mutational analysis of the osmolyte effect indicated that the changed Delta(o)DeltaG values upon mutation (Delta(m)Delta(o)DeltaG), as estimated from the Deltag values, are similar to the experimental values. Structural-based analysis of the osmolyte effect was also performed using model denatured structures: (a) a fully extended model (single chain) with no disulfide bond, (b) two-part, unfolded models (two chains) with a disulfide bond constructed through molecular dynamic (MD) simulation, and (c) a two-part, folded model (two chains). The two-part, unfolded models were expected to be more suitable as denatured structures. The Delta(o)DeltaG values calculated using the two-part, unfolded models were more consistent with experimental values than those calculated using the fully extended and two-part, folded models. This suggests that MD simulation is useful for testing denatured structures. These results indicate that the osmophobic theory can explain the osmolyte effect on protein stability.     

Atypical effect of salts on the thermodynamic stability of human prion protein

Prion diseases are associated with the conversion of cellular prion protein, PrPC, into a misfolded oligomeric form, PrPSc. Previous studies indicate that salts promote conformational conversion of the recombinant prion protein into a PrPSc-like form. To gain insight into the mechanism of this effect, here we have studied the influence of a number of salts (sodium sulfate, sodium fluoride, sodium acetate, and sodium chloride) on the thermodynamic stability of the recombinant human prion protein. Chemical unfolding studies in urea show that at low concentrations (below approximately 50 mm), all salts tested significantly reduced the thermodynamic stability of the protein. This highly unusual response to salts was observed for both the full-length prion protein as well as the N-truncated fragments huPrP90-231 and huPrP122-231. At higher salt concentrations, the destabilizing effect was gradually reversed, and salts behaved according to their ranking in the Hofmeister series. The present data indicate that electrostatic interactions play an unusually important role in the stability of the prion protein. The abnormal effect of salts is likely because of the ion-induced destabilization of salt bridges (Asp144-Arg148 and/or Asp147-Arg151) in the extremely hydrophilic helix 1. Contrary to previous suggestions, this effect is not due to the interaction of ions with the glycine-rich flexible N-terminal region of the prion protein. The results of this study suggest that ionic species present in the cellular environment may control the PrPC to PrPSc conversion by modulating the thermodynamic stability of the native PrPC isoform.     

Switching two-state to three-state kinetics in the helical protein Im9 via the optimisation of stabilising non-native interactions by design

The four-helix protein Im7 folds through an on-pathway intermediate at pH 7.0 and 10 degrees C. By contrast, under these conditions there is no evidence for a populated intermediate in the folding of its more stable homologue, Im9, even in the presence of 0.4 M sodium sulphate. Previous studies using phi-value analysis have shown that the Im7 intermediate is misfolded, in that three of its four native helices are formed, but are docked in a non-native manner. Using knowledge of the structure of the intermediate of Im7, we have used rational design to stabilise an intermediate formed during the folding of Im9 by the introduction of specific stabilising interactions at positions known to stabilise the Im7 folding intermediate through non-native interactions. We show that the redesigned Im9 sequence folds with three-state kinetics at pH 7.0 and have used phi-value analysis to demonstrate that this species resembles the misfolded intermediate populated during Im7 folding. The redesigned Im9 sequence folds 20-fold faster than the wild-type protein under conditions in which folding is two-state. The data show that intermediate formation is an important feature of folding, even for small proteins such as Im9 for which these partially folded states do not become significantly populated. In addition, they show that the introduction of stabilising interactions can lead to rapid refolding, even when the contacts introduced are non-native.     

The effect of elevated temperature on the complexation of am3+ with chloride

The temperature dependence of the formation of the 1:1 Am³⁺-Cl^- solution complex is reported. Stability constants were determined as a function of temperature using a solvent extraction technique. The stability constant for the formation of the AmCl²⁺ complex is measurably higher at 50‡C compared to values determined under ambient conditions. From these studies, we have estimated the enthalpy of complexation in this temperature range. A small, positive enthalpy was observed. The results are discussed in the context of previous studies.     

The effect of Glu75 of staphylococcal nuclease on enzyme activity, protein stability and protein unfolding.

Staphylococcal nuclease mutants, E57G and E75G, were generated. A comparison of the kinetic parameters both for mutants and wild-type protein shows that the Michaelis constants (Km) were almost identical for the wild-type protein and E57G mutant. An approximately 30-fold decrease in Km compared with the wild-type protein was observed for the E75G mutant. The turnover numbers for the enzyme (kcat) were higher with both the wild-type protein and the E57G mutant (3.88 +/- 0.21 x 103 s-1 and 3.71 +/- 0.28 x 103 s-1) than with the E75G mutant (3.04 +/- 0.02 x 102 s-1). The results of thermal denaturation with differential scanning microcalorimetry indicate that the excess calorimetric enthalpy of denaturations, DeltaHcal, was almost identical for the wild-type protein and E57G mutant (84.1 +/- 6.2 kcal.mol-1 and 79.3 +/- 7.1 kcal.mol-1, respectively). An approximately twofold decrease in DeltaHcal compared with the wild-type protein was observed for the E75G mutant (42.7 +/- 5.5 kcal.mol-1). These outcomes imply that Glu at position 75 plays a significant role in maintaining enzyme activity and protein stability. Further study of the unfolding of the wild-type protein and E75G mutant was conducted by using time-resolved fluorescence with a picosecond laser pulse. Two fluorescent lifetimes were found in the subnanosecond time range. The faster lifetime (tau2) did not generally vary with either pH or the concentration of guanidinium hydrochloride (GdmHCl) in the wild-type protein and the E75G mutant. The slow lifetime (tau1), however, did vary with these parameters and was faster as the protein is unfolded by either pH or GdmHCl denaturation. The midpoints of the transition for tau1 are pH 3.5 and 5.8 for the wild-type protein and E75G mutant, respectively, and the GdmHCl concentrations are 1.1 m and 0.6 m for the wild-type protein and E75G mutant, respectively. Parallel steady-state fluorescence measurements have also been carried out and the results are in general agreement with the time-resolved fluorescence experiments, indicating that Glu at position 75 plays an important role in protein unfolding.     

Using protein design to dissect the effect of charged residues on metal binding and protein stability

Ca2+ controls biological processes by interacting with proteins with different affinities, which are largely influenced by the electrostatic interaction from the local negatively charged ligand residues in the coordination sphere. We have developed a general strategy for rationally designing stable Ca2+- and Ln3+-binding proteins that retain the native folding of the host protein. Domain 1 of cluster differentiation 2 (CD2) is the host for the two designed proteins in this study. We investigate the effect of local charge on Ca2+-binding affinity based on the folding properties and metal-binding affinities of the two proteins that have similarly located Ca2+-binding sites with two shared ligand positions. While mutation and Ca2+ binding do not alter the native structure of the protein, Ca2+ binding specifically induced changes around the designed Ca2+-binding site. The designed protein with a -5 charge at the binding sphere displays a 14-, 20-, and 12-fold increase in the binding affinity for Ca2+, Tb3+, and La3+, respectively, compared to the designed protein with a -3 charge, which suggests that higher local charges are preferred for both Ca2+ and Ln3+ binding. The localized charged residues significantly decrease the thermal stability of the designed protein with a -5 charge, which has a T(m) of 41 degrees C. Wild-type CD2 has a T(m) of 61 degrees C, which is similar to the designed protein with a -3 charge. This decrease is partially restored by Ca2+ binding. The effect on the protein stability is modulated by the environment and the secondary structure locations of the charged mutations. Our study demonstrates the capability and power of protein design in unveiling key determinants to Ca2+-binding affinity without the complexities of the global conformational changes, cooperativity, and multibinding process found in most natural Ca2+-binding proteins.     

Contribution of the hydrophobic effect to globular protein stability.

The decrease in conformational stability, delta(delta G), has been measured for 72 aliphatic side-chain mutants from four proteins in which a larger side-chain is replaced by a smaller side-chain so that steric effects are minimal. When these delta(delta G) values are corrected to the same accessibility, namely 100% buried, then the following -delta(delta G) values per -CH2- group (in kcal/mol) are obtained: Ile----Val (1.26), Ala (1.26), Gly (1.26); Leu----Ala (1.16), Gly (1.21); Val----Ala (1.23), Gly (1.53). The average of these values is 1.27(+/- 0.07) kcal/mol. The 72 individual values range from 0 to 2.4 kcal/mol with an average value of 1.27(+/- 0.51) (standard deviation) kcal/mol. When the delta Gtr values from n-octanol to water are corrected for the difference in volume between the solutes and the solvents, the average value for the same substitutions is 1.25(+/- 0.05) kcal/mol. This suggests that proteins gain 1.3(+/- 0.5) kcal/mol in stability for each -CH2- group buried in folding, and, furthermore, that the volume corrected delta Gtr values for n-octanol for the amino acid side-chains provide good estimates of the contribution of the hydrophobic effect to globular protein stability.     

The fluorescent monomeric protein Kusabira Orange. Pressure effect on its structure and stability

The structure and stability of the fluorescent protein monomeric Kusabira Orange (mKO), a GFP-like protein, was studied under different pressure levels and in different chemical environments. At different pH values (between pH 7.4 and pH 4.0) and under a pressure up to 600 MPa (at 25 °C), mKO did not show significant fluorescence spectral changes, indicating a structural stability of the protein. In more extreme chemical conditions (at pH 4.0 in the presence of 0.8 M guanidine hydrochloride), a marked reduction of mKO fluorescence intensity emission was observed at pressures above 300 MPa. This fluorescence emission quenching may be due to the loss of the intermolecular bonds and, consequently, to the destructuration of the mKO chromophore structure. Since the electrostatic and hydrophobic interactions as well as the salt bridges present in proteins are usually perturbed under high pressure, the reduction of mKO fluorescence intensity emission is associated to the perturbation of the protein salt bridges network.     

The effect of heme on the conformational stability of micro-myoglobin

Micro-myoglobin, the isolated heme-binding subdomain of myoglobin, is a valuable model system for the investigation of heme recognition and binding by proteins, and provides an example of protein folding induced by cofactor binding. Theoretical studies by molecular dynamics simulations on apo- and holo-micro-myoglobin show that, by contrast with the case of the full-length wild-type protein and in agreement with earlier experimental evidence, the apo-protein is not stably folded in a native-like conformation. With the cofactor bound, however, the protein fragment maintains its folded conformation over 1.5 ns in molecular dynamics simulations. Further inspection of the model structures reveals that the role of heme in stabilizing the folded state is not only a result of its direct interactions with binding residues (His93, Arg45 and Lys96), but also derives from its shielding effect on a long-range electrostatic interaction between Arg45 and Asp60, which, in the molecular dynamics simulations, apparently triggers the unfolding process of apo-micro-myoglobin.