Frataxin from Psychromonas ingrahamii as a model to study stability modulation within the CyaY protein family

Adaptation of life to low temperatures influences both protein stability and flexibility. Thus, proteins from psychrophilic organisms are excellent models to study relations between these properties. Here we focused on frataxin from Psychromonas ingrahamii (pFXN), an extreme psychrophilic sea ice bacterium that can grow at temperatures as low as − 12 °C. This α/β protein is highly conserved and plays a key role in iron homeostasis as an iron chaperone. In contrast to other frataxin homologs, chemical and temperature unfolding experiments showed that the thermodynamic stability of pFXN is strongly modulated by pHs: ranging from 5.5 ± 0.9 (pH 6.0) to 0.9 ± 0.3 kcal mol^− 1 (pH 8.0). This protein was crystallized and its X-ray structure solved at 1.45 Å. Comparison of B-factor profiles between Escherichia coli and P. ingrahamii frataxin variants (51% of identity) suggests that, although both proteins share the same structural features, their flexibility distribution is different. Molecular dynamics simulations showed that protonation of His44 or His67 in pFXN lowers the mobility of regions encompassing residues 20–30 and the C-terminal end, probably through favorable electrostatic interactions with residues Asp27, Glu42 and Glu99. Since the C-terminal end of the protein is critical for the stabilization of the frataxin fold, the predictions presented may be reporting on the microscopic origin of the decrease in global stability produced near neutral pH in the psychrophilic variant. We propose that suboptimal electrostatic interactions may have been an evolutionary strategy for the adaptation of frataxin flexibility and function to cold environments.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Cartilage Acidic Protein 2 a hyperthermostable,high affinity calcium-binding protein

Cartilage Acidic Protein 2 (CRTAC2) is a novel protein present from prokaryotes to vertebrates with abundant expression in the teleost fish pituitary gland and an isoform of CRTAC1, a chondrocyte marker in humans. The two proteins are non-integrins containing N-terminal integrin-like Ca^2 +-binding motifs and their structure and function remain to be assigned. Structural studies of recombinant sea bream (sb)CRTAC2 revealed it is composed of 8.8% α-helix, 33.4% β-sheet and 57.8% unordered protein. sbCRTAC2 bound Ca^2 + with high affinity (K_d = 1.46 nM) and favourable Gibbs free energy (?G = − 12.4 kcal/mol). The stoichiometry for Ca^2 + bound to sbCRTAC2 at saturation indicated six Ca^2 + ligand-binding sites exist per protein molecule. No conformational change in sbCRTAC2 occurred in the presence of Ca^2 +. Fluorescence emission revealed that the tertiary structure of the protein is hyperthermostable between 25 °C and 95 °C and the fully unfolded state is only induced by chemical denaturing (4 M GndCl). sbCRTAC has a widespread tissue distribution and is present as high molecular weight aggregates, although strong reducing conditions promote formation of the monomer. sbCRTAC2 promotes epithelial cell outgrowth in vitro suggesting it may share functional homology with mammalian CRTAC1, recently implicated in cell–cell and cell–matrix interactions.     

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

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

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.     

Charge-charge interactions influence the denatured state ensemble and contribute to protein stability

Several recent studies have shown that it is possible to increase protein stability by improving electrostatic interactions among charged groups on the surface of the folded protein. However, the stability increases are considerably smaller than predicted by a simple Coulomb's law calculation, and in some cases, a charge reversal on the surface leads to a decrease in stability when an increase was predicted. These results suggest that favorable charge-charge interactions are important in determining the denatured state ensemble, and that the free energy of the denatured state may be decreased more than that of the native state by reversing the charge of a side chain. We suggest that when the hydrophobic and hydrogen bonding interactions that stabilize the folded state are disrupted, the unfolded polypeptide chain rearranges to compact conformations with favorable long-range electrostatic interactions. These charge-charge interactions in the denatured state will reduce the net contribution of electrostatic interactions to protein stability and will help determine the denatured state ensemble. To support this idea, we show that the denatured state ensemble of ribonuclease Sa is considerably more compact at pH 7 where favorable charge-charge interactions are possible than at pH 3, where unfavorable electrostatic repulsion among the positive charges causes an expansion of the denatured state ensemble. Further support is provided by studies of the ionic strength dependence of the stability of charge-reversal mutants of ribonuclease Sa. These results may have important implications for the mechanism of protein folding.     

The energetic contribution of induced electrostatic asymmetry to DNA bending by a site-specific protein

DNA bending can be promoted by reducing the net negative electrostatic potential around phosphates on one face of the DNA, such that electrostatic repulsion among phosphates on the opposite face drives bending toward the less negative surface. To provide the first assessment of energetic contribution to DNA bending when electrostatic asymmetry is induced by a site-specific DNA binding protein, we manipulated the electrostatics in the EcoRV endonuclease-DNA complex by mutation of cationic side chains that contact DNA phosphates and/or by replacement of a selected phosphate in each strand with uncharged methylphosphonate. Reducing the net negative charge at two symmetrically located phosphates on the concave DNA face contributes − 2.3 kcal mol^− 1 to − 0.9 kcal mol^− 1 (depending on position) to complex formation. In contrast, reducing negative charge on the opposing convex face produces a penalty of + 1.3 kcal mol^− 1. Förster resonance energy transfer experiments show that the extent of axial DNA bending (about 50°) is little affected in modified complexes, implying that modification affects the energetic cost but not the extent of DNA bending. Kinetic studies show that the favorable effects of induced electrostatic asymmetry on equilibrium binding derive primarily from a reduced rate of complex dissociation, suggesting stabilization of the specific complex between protein and markedly bent DNA. A smaller increase in the association rate may suggest that the DNA in the initial encounter complex is mildly bent. The data imply that protein-induced electrostatic asymmetry makes a significant contribution to DNA bending but is not itself sufficient to drive full bending in the specific EcoRV-DNA complex.     

The courtship of proteins: Understanding the encounter complex

The formation of protein complexes involves an encounter complex, in which proteins show few specific interactions and assume many orientations. Recent kinetic and structural studies have shed light on this elusive state. It is generally dominated by electrostatic interactions, although hydrophobic interactions can play a role. During the encounter phase the proteins remain largely solvated. In extreme cases, the proteins only form an encounter complex, and in many other complexes, the encounter state constitutes a significant amount (5% or more), indicating that the energy difference between encounter and productive complexes is small. Thus, the encounter complex represents an essential part of protein complexes.     

Residual structure in unfolded proteins

The denatured state ensemble (DSE) of unfolded proteins, once considered to be well-modeled by an energetically featureless random coil, is now well-known to contain flickering elements of residual structure. The position and nature of DSE residual structure may provide clues toward deciphering the protein folding code. This review focuses on recent advances in our understanding of the nature of DSE collapse under folding conditions, the quantification of the stability of residual structure in the DSE, the determination of the location and types of residues involved in thermodynamically significant residual structure and advances in detection of long-range interactions in the DSE.     

pK(a) values for the unfolded state under native conditions explain the pH-dependent stability of PGB1

Understanding the role of electrostatics in protein stability requires knowledge of these interactions in both the folded and unfolded states. Electrostatic interactions can be probed experimentally by characterizing ionization equilibria of titrating groups, parameterized as pK_a values. However, pK_a values of the unfolded state are rarely accessible under native conditions, where the unfolded state has a very low population. Here, we report pK_a values under nondenaturing conditions for two unfolded fragments of the protein G B1 domain that mimic the unfolded state of the intact protein. pK_a values were determined for carboxyl groups by monitoring their pH-dependent ¹³C chemical shifts. Monte Carlo simulations using a Gaussian chain model provide corrections for changes in electrostatic interactions that arise from fragmentation of the protein. Most pK_a values for the unfolded state agree well with model values, but some residues show significant perturbations that can be rationalized by local electrostatic interactions. The pH-dependent stability was calculated from the experimental pK_a values of the folded and unfolded states and compared to experimental stability data. The use of experimental pK_a values for the unfolded state results in significantly improved agreement with experimental data, as compared to calculations based on model data alone.     

Characterization of Ca and phosphocholine interactions with C-reactive protein using a surface plasmon resonance biosensor

The interactions between Ca²⁺ and C-reactive protein (CRP) have been characterized using a surface plasmon resonance (SPR) biosensor. The protein was immobilized on a sensor chip, and increasing concentrations of Ca²⁺ or phosphocholine were injected. Binding of Ca²⁺ induced a 10-fold higher signal than expected from the molecular weight of Ca²⁺. It was interpreted to result from the conformational change that occurs on binding of Ca²⁺. Two sites with different characteristics were distinguished: a high-affinity site with K_D = 0.03 mM and a low-affinity site with K_D = 5.45 mM. The pH dependencies of the two Ca²⁺ interactions were different and enabled the assignment of the different sites in the three-dimensional structure of CRP. There was no evidence for cooperativity in the phosphocholine interaction, which had K_D = 5 μM at 10 mM Ca²⁺. SPR biosensors can clearly detect and quantify the binding of very small molecules or ions to immobilized proteins despite the theoretically very low signals expected on binding, provided that significant conformational changes are involved. Both the interactions and the conformational changes can be characterized. The data have important implications for the understanding of the function of CRP and suggest that Ca²⁺ is an efficient regulator under physiological conditions.     

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the β-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.     

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₂- group), than to the stability of a large protein, VlsE (1.6 ± 0.3 kcal/mol per -CH₂- 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₂- 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.     

Ionic-strength-dependent effects in protein folding: analysis of rate equilibrium free-energy relationships and their interpretation

The traditional approach to studying protein folding involves applying a perturbation, usually denaturant or mutation, and determining the effect upon the free energy of folding, DeltaG0, and the activation free energy, DeltaG(not equal). Data collected as a function of the perturbation can be used to construct rate equilibrium free-energy relationships, which report on the development of interactions in the transition state for folding. We examine the use of the ionic-strength-dependent rate equilibrium free-energy relationship in protein folding using the N-terminal domain of L9, a small alpha-beta protein, as a model system. Folding is two-state for the range of ionic strength examined, 0.045-1.52 M. The plot of DeltaG(not equal) versus DeltaG0 is linear (r2= 0.918), with a slope equal to 0.45. The relatively low value of the slope indicates that the ionic-strength-dependent interactions are modestly developed in the transition state. The slope is, however, greater than that of a plot of DeltaG(not equal) versus DeltaG0 constructed by varying pH, thus demonstrating directly that ionic-strength-dependent studies probe more than simple electrostatic interactions. Potential transition movement was probed by analysis of the denaturant, ionic strength cross-interaction parameters. The values are small but nonzero and positive, suggesting a small shift of the transition state toward the native state as the protein is destabilized, i.e., Hammond behavior. The complications that arise in the interpretation of ionic-strength-dependent rate equilibrium free-energy relationships are discussed, and it is concluded that the ionic-strength-dependent studies do not provide a reliable indicator of the role of electrostatic interactions. Complications include incomplete screening of electrostatic interactions, specific ion binding, Hofmeister effects, and the potential presence of electrostatic interactions in the denatured state ensemble.     

Propensities of Aromatic Amino Acids versus Leucine and Proline to Induce Residual Structure in the Denatured-State Ensemble of Iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X = Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.     

Structural and mechanistic exploration of acid resistance: kinetic stability facilitates evolution of extremophilic behavior

Kinetically stable proteins are unique in that their stability is determined solely by kinetic barriers rather than by thermodynamic equilibria. To better understand how kinetic stability promotes protein survival under extreme environmental conditions, we analyzed the unfolding behavior and determined the structure of Nocardiopsis alba Protease A (NAPase), an acid-resistant, kinetically stable protease, and compared these results with a neutrophilic homolog, α-lytic protease (αLP). Although NAPase and αLP have the same number of acid-titratable residues, kinetic studies revealed that the height of the unfolding free energy barrier for NAPase is less sensitive to acid than that of αLP, thereby accounting for NAPase's improved tolerance of low pH. A comparison of the αLP and NAPase structures identified multiple salt-bridges in the domain interface of αLP that were relocated to outer regions of NAPase, suggesting a novel mechanism of acid stability in which acid-sensitive electrostatic interactions are rearranged to similarly affect the energetics of both the native state and the unfolding transition state. An acid-stable variant of αLP in which a single interdomain salt-bridge is replaced with a corresponding intradomain NAPase salt-bridge shows a dramatic > 15-fold increase in acid resistance, providing further evidence for this hypothesis. These observations also led to a general model of the unfolding transition state structure for αLP protease family members in which the two domains separate from each other while remaining relatively intact themselves. These results illustrate the remarkable utility of kinetic stability as an evolutionary tool for developing longevity over a broad range of harsh conditions.     

A logical OR redundancy within the Asx-Pro-Asx-Gly type I beta-turn motif

Turn secondary structure is essential to the formation of globular protein architecture. Turn structures are, however, much more complex than either α-helix or β-sheet, and the thermodynamics and folding kinetics are poorly understood. Type I β-turns are the most common type of reverse turn, and they exhibit a statistical consensus sequence of Asx-Pro-Asx-Gly (where Asx is Asp or Asn). A comprehensive series of individual and combined Asx mutations has been constructed within three separate type I 3:5 G1 bulge β-turns in human fibroblast growth factor-1, and their effects on structure, stability, and folding have been determined. The results show a fundamental logical OR relationship between the Asx residues in the motif, involving H-bond interactions with main-chain amides within the turn. These interactions can be modulated by additional interactions with residues adjacent to the turn at positions i + 4 and i + 6. The results show that the Asx residues in the turn motif make a substantial contribution to the overall stability of the protein, and the Asx logical OR relationship defines a redundant system that can compensate for deleterious point mutations. The results also show that the stability of the turn is unlikely to be the prime determinant of formation of turn structure in the folding transition state.