Molecular basis of cooperativity in protein folding. IV. CORE: a general cooperative folding model

In silico Mapping of Protein Unfolding Mutations for Inherited Disease

The effect of disease-causing missense mutations on protein folding is difficult to evaluate. To understand this relationship, we developed the unfolding mutation screen (UMS) for in silico evaluation of the severity of genetic perturbations at the atomic level of protein structure. The program takes into account the protein-unfolding curve and generates propensities using calculated free energy changes for every possible missense mutation at once. These results are presented in a series of unfolding heat maps and a colored protein 3D structure to show the residues critical to the protein folding and are available for quick reference. UMS was tested with 16 crystal structures to evaluate the unfolding for 1391 mutations from the ProTherm database. Our results showed that the computational accuracy of the unfolding calculations was similar to the accuracy of previously published free energy changes but provided a better scale. Our residue identity control helps to improve protein homology models. The unfolding predictions for proteins involved in age-related macular degeneration, retinitis pigmentosa, and Leber's congenital amaurosis matched well with data from previous studies. These results suggest that UMS could be a useful tool in the analysis of genotype-to-phenotype associations and next-generation sequencing data for inherited diseases.

Evaluation of folding co-operativity of a chimeric protein based on the molecular recognition between polyproline ligands and SH3 domains

In previous work, we designed a chimeric protein, named SPCp41, to evaluate the thermodynamics of the interaction between SH3 domains and proline-rich ligands by combining thermal unfolding measurements and mutagenesis. Here, we have investigated the energetic integrity of the chain extension corresponding to the ligand sequence into the native structure, since the opposite will produce changes in the folding mechanism of the SH3 domain that may give rise to undesirable contributions to the thermodynamic parameters. We have analysed the folding-unfolding kinetics under standard conditions (50 mM phosphate pH 7). Kinetic evolutions are well described by a bi-exponential where, on top of the main kinetic phase, a low-populated slower phase appears as a consequence of cis-trans isomerisation of Pro39, as demonstrated by the influence of prolyl isomerases and by mutational analysis. There is also a burst phase possibly due to a productive formation of some helical ensembles. The main evolution, accounting for the true folding kinetics of SPCp41, can be considered as a two-state process, where the folding transition state produces essentially the same picture shown by the circular permutant S19-P20s (the 'nucleus' of the design) and the ligand will dock at the latter stages of the two-state process. Thus, all conclusions argue in favour of the effectiveness of SPCp41 to study energetic, dynamic and structural aspects of SH3-ligand interactions.

Thermodynamic and dynamic factors involved in the stability of native protein structure in amorphous solids in relation to levels of hydration

The internal, dynamical fluctuations of protein molecules exhibit many of the features typical of polymeric and bulk small molecule glass forming systems. The response of a protein's internal molecular mobility to temperature changes is similar to that of other amorphous systems, in that different types of motions freeze out at different temperatures, suggesting they exhibit the alpha-beta-modes of motion typical of polymeric glass formers. These modes of motion are attributed to the dynamic regimes that afford proteins the flexibility for function but that also develop into the large-scale collective motions that lead to unfolding. The protein dynamical transition, T(d), which has the same meaning as the T(g) value of other amorphous systems, is attributed to the temperature where protein activity is lost and the unfolding process is inhibited. This review describes how modulation of T(d) by hydration and lyoprotectants can determine the stability of protein molecules that have been processed as bulk, amorphous materials. It also examines the thermodynamic, dynamic, and molecular factors involved in stabilizing folded proteins, and the effects typical pharmaceutical processes can have on native protein structure in going from the solution state to the solid state.

Circular dichroism analysis of allergens

Recombinant allergens have gained a lot of importance lately for the diagnosis of allergic diseases and for specific immunotherapy. To characterize recombinant allergens and potential hypo-allergenic derivatives thereof circular dichroism (CD) spectroscopy is used widely. It is a convenient, fast method to assess the structural integrity of the recombinant proteins, compare them with the allergens isolated from natural sources, and to determine the effects of mutations on the structural properties. In this paper, we will describe the techniques and the most useful applications of CD spectroscopy to the field of allergy research.

Modeling protein density of states: additive hydrophobic effects are insufficient for calorimetric two-state cooperativity

A well-established experimental criterion for two-state thermodynamic cooperativity in protein folding is that the van't Hoff enthalpy DeltaH(vH) around the transition midpoint is equal, or very nearly so, to the calorimetric enthalpy DeltaH(cal) of the entire transition. This condition is satisfied by many small proteins. We use simple lattice models to provide a statistical mechanical framework to elucidate how this calorimetric two-state picture may be reconciled with the hierarchical multistate scenario emerging from recent hydrogen exchange experiments. We investigate the feasibility of using inverse Laplace transforms to recover the underlying density of states (i.e., enthalpy distribution) from calorimetric data. We find that the constraint imposed by DeltaH(vH)/DeltaH(cal) approximately 1 on densities of states of proteins is often more stringent than other "two-state" criteria proposed in recent theoretical studies. In conjunction with reasonable assumptions, the calorimetric two-state condition implies a narrow distribution of denatured-state enthalpies relative to the overall enthalpy difference between the native and the denatured conformations. This requirement does not always correlate with simple definitions of "sharpness" of a transition and has important ramifications for theoretical modeling. We find that protein models that assume capillarity cooperativity can exhibit overall calorimetric two-state-like behaviors. However, common heteropolymer models based on additive hydrophobic-like interactions, including highly specific two-dimensional Gō models, fail to produce proteinlike DeltaH(vH)/DeltaH(cal) approximately 1. A simple model is constructed to illustrate a proposed scenario in which physically plausible local and nonlocal cooperative terms, which mimic helical cooperativity and environment-dependent hydrogen bonding strength, can lead to thermodynamic behaviors closer to experiment. Our results suggest that proteinlike thermodynamic cooperativity may require a cooperative interplay between local and nonlocal interactions. The prospect of using calorimetric data to constrain Z-scores of knowledge-based potentials is discussed.

Thermozymes

Enzymes synthesized by thermophiles (organisms with optimal growth temperatures > 60 degrees C) and hyperthermophiles (optimal growth temperatures > 80 degrees C) are typically thermostable (resistant to irreversible inactivation at high temperatures) and thermophilic (optimally active at high temperatures, i.e., > 60 degrees C). These enzymes, called thermozymes, share catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, thermozymes usually retain their thermal properties, suggesting that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, and crystal structure comparisons indicate that thermozymes are, indeed, very similar to mesophilic enzymes. No obvious sequence or structural features account for enzyme thermostability and thermophilicity. Thermostability and thermophilicity molecular mechanisms are varied, differing from enzyme to enzyme. Thermostability and thermophilicity are usually caused by the accumulation of numerous subtle sequence differences. This review concentrates on the mechanisms involved in enzyme thermostability and thermophilicity. Their relationships with protein rigidity and flexibility and with protein folding and unfolding are discussed. Intrinsic stabilizing forces (e.g., salt bridges, hydrogen bonds, hydrophobic interactions) and extrinsic stabilizing factors are examined. Finally, thermozymes' potential as catalysts for industrial processes and specialty uses are discussed, and lines of development (through new applications, and protein engineering) are also proposed.

The conformational equilibrium of human growth hormone

The structural stability of recombinant human growth hormone (rhGH) has been studied by differential scanning calorimetry, circular dichroism and by following the tyrosine and histidine chemical shifts in the 1H NMR spectrum. These studies demonstrate that the folding/unfolding equilibrium of rhGH involves a partially folded dimeric intermediate. The formation of this dimeric intermediate is a reversible process. At acid pH (pH 3) the conformational equilibrium is reversible even at high protein concentrations (10 mg/ml). At neutral pH reversibility is observed only at low protein concentrations (<0.5 mg/ml). The free energy of this intermediate conformation is only approximately 3 kcal/mol apart from the native state indicating that the conformational equilibrium can be effectively modulated by changes in solvent composition or physical conditions. According to the spectroscopic and thermodynamic results, the formation of the dimeric intermediate occurs without a major loss in helical content and is driven by the formation of substantial hydrophobic contacts between two partially folded molecules. A thermodynamic model that accounts quantitatively for the experimental data has been developed. These studies demonstrate that partially folded conformations of certain proteins are able to form stoichiometric complexes, and that the formation of these complexes provide a significant source of stabilizing Gibbs energy for conformational states that, otherwise, will be characterized by extremely unfavorable free energies.

Consistency in structural energetics of protein folding and peptide recognition

We report a new free energy decomposition that includes structure-derived atomic contact energies for the desolvation component, and show that it applies equally well to the analysis of single-domain protein folding and to the binding of flexible peptides to proteins. Specifically, we selected the 17 single-domain proteins for which the three-dimensional structures and thermodynamic unfolding free energies are available. By calculating all terms except the backbone conformational entropy change and comparing the result to the experimentally measured free energy, we estimated that the mean entropy gain by the backbone chain upon unfolding (delta Sbb) is 5.3 cal/K per mole of residue, and that the average backbone entropy for glycine is 6.7 cal/K. Both numbers are in close agreement with recent estimates made by entirely different methods, suggesting a promising degree of consistency between data obtained from disparate sources. In addition, a quantitative analysis of the folding free energy indicates that the unfavorable backbone entropy for each of the proteins is balanced predominantly by favorable backbone interactions. Finally, because the binding of flexible peptides to receptors is physically similar to folding, the free energy function should, in principle, be equally applicable to flexible docking. By combining atomic contact energies, electrostatics, and sequence-dependent backbone entropy, we calculated a priori the free energy changes associated with the binding of four different peptides to HLA-A2, 1 MHC molecule and found agreement with experiment to within 10% without parameter adjustment.

Characterization of ligand binding of a soluble human insulin-like growth factor I receptor variant suggests a ligand-induced conformational change

Details of the signal transduction mechanisms of the tyrosine kinase family of growth factor receptors remain elusive. In this work, we describe an extensive study of kinetic and thermodynamic aspects of growth factor binding to a soluble extracellular human insulin-like growth factor-I receptor (sIGF-IR) variant. The extracellular receptor domains were produced fused to an IgG-binding protein domain (Z) in transfected human 293 cells as a correctly processed secreted alpha-beta'-Z dimer. The receptor was purified using IgG affinity chromatography, rendering a pure and homogenous protein in yields from 1 to 5 mg/liter of conditioned cell media. Biosensor technology (BIAcore) was applied to measure the insulin-like growth factor-I (IGF-I), des(1-3)IGF-I, insulin-like growth factor-II, and insulin ligand binding rate constants to the immobilized IGF-IR-Z. The association equilibrium constant, Ka, for the IGF-I interaction is determined to 2.8 x 10(8) M-1 (25 degrees C). Microcalorimetric titrations on IGF-I/IGF-IR-Z were performed at three different temperatures (15, 25, and 37 degrees C) and in two different buffer systems at 25 degrees C. From these measurements, equilibrium constants for the 1:1 (IGF-I:(alpha-beta'-Z)2) receptor complex in solution are deduced to 0.96 x 10(8) M-1 (25 degrees C). The determined heat capacity change for the process is large and negative, -0.51 kcal (K mol)-1. Further, the entropy change (DeltaS) at 25 degrees C is large and negative. Far- and near-UV circular dichroism measurements display significant changes over the entire wavelength range upon binding of IGF-I to IGF-IR-Z. These data are all consistent with a significant change in structure of the system upon IGF-I binding.

Structure-based statistical thermodynamic analysis of T4 lysozyme mutants: structural mapping of cooperative interactions

The recent development of a structural parameterization of the energetics of protein folding has permitted the incorporation of the functions that describe the enthalpy, entropy and heat capacity changes, i.e. the individual components of the Gibbs energy, into a statistical thermodynamic formalism that describes the distribution of conformational states under equilibrium conditions. The goal of this approach is to construct with the computer a large ensemble of conformational states, and then to derive the most probable population distribution, i.e. the distribution of states that best accounts for a wide array of experimental observables. This analysis has been applied to four different mutants of T4 lysozyme (S44A, S44G, V131A, V131G). It is shown that the structural parameterization predicts well the stability of the protein and the effects of the mutations. The entire set of folding constants per residue has been calculated for the four mutants. In all cases, the effect of the mutations propagates beyond the mutation site itself through sequence and three-dimensional space. This phenomenon occurs despite the fact that the mutations are at solvent-exposed locations and do not directly affect other interactions in the protein. These results suggest that single amino acid mutations at solvent-exposed locations, or other locations that cause a minimal perturbation, can be used to identify the extent of cooperative interactions. The magnitude and extent of these effects and the accuracy of the algorithm can be tested by means of NMR-detected hydrogen exchange.

Side-chain conformational entropy in protein folding

An important, but often neglected, contribution to the thermodynamics of protein folding is the loss of entropy that results from restricting the number of accessible side-chain conformers in the native structure. Conformational entropy changes can be found by comparing the number of accessible rotamers in the unfolded and folded states, or by estimating fusion entropies. Comparison of several sets of results using different techniques shows that the mean conformational free energy change (T delta S) is 1 kcal.mol-1 per side chain or 0.5 kcal.mol-1 per bond. Changes in vibrational entropy appear to be negligible compared to the entropy change resulting from the loss of accessible rotamers. Side-chain entropies can help rationalize alpha-helix propensities, predict protein/inhibitor complex structures, and account for the distribution of side chains on the protein surface or interior.

Forces and factors that contribute to the structural stability of membrane proteins

While a considerable amount of literature deals with the structural energetics of water-soluble proteins, relatively little is known about the forces that determine the stability of membrane proteins. Similarly, only a few membrane protein structures are known at atomic resolution, although new structures have recently been described. In this article, we review the current knowledge about the structural features of membrane proteins. We then proceed to summarize the existing literature regarding the thermal stability of bacteriorhodopsin, cytochrome-c oxidase, the band 3 protein, Photosystem II and porins. We conclude that a fundamental difference between soluble and membrane proteins is the high thermal stability of intrabilayer secondary structure elements in membrane proteins. This property manifests itself as incomplete unfolding, and is reflected in the observed low enthalpies of denaturation of most membrane proteins. By contrast, the extramembranous parts of membrane proteins may behave much like soluble proteins. A brief general account of thermodynamics factors that contribute to the stability of water soluble and membrane proteins is presented.

Forces and factors that contribute to the structural stability of membrane proteins

While a considerable amount of literature deals with the structural energetics of water-soluble proteins, relatively little is known about the forces that determine the stability of membrane proteins. Similarly, only a few membrane protein structures are known at atomic resolution, although new structures have recently been described. In this article, we review the current knowledge about the structural features of membrane proteins. We then proceed to summarize the existing literature regarding the thermal stability of bacteriorhodopsin, cytochrome-c oxidase, the band 3 protein, Photosystem II and porins. We conclude that a fundamental difference between soluble and membrane proteins is the high thermal stability of intrabilayer secondary structure elements in membrane proteins. This property manifests itself as incomplete unfolding, and is reflected in the observed low enthalpies of denaturation of most membrane proteins. By contrast, the extramembranous parts of membrane proteins may behave much like soluble proteins. A brief general account of thermodynamics factors that contribute to the stability of water soluble and membrane proteins is presented.

Configurational effects in antibody-antigen interactions studied by microcalorimetry

In this paper we study the binding of two monoclonal antibodies, E3 and E8, to cytochrome c using high-sensitivity isothermal titration calorimetry. We combine the calorimetric results with empirical calculations which relate changes in heat capacity to changes in entropy which arise from the hydrophobic effect. The change in heat capacity for binding E3 is -350 +/- 60 cal K-1 mol-1 while for E8 it is -165 +/- 40 cal K-1 mol-1. This result indicates that the hydrophobic effect makes a much larger contribution for E3 than for E8. Since the total entropy change at 25 degrees C is very similar for both antibodies, it follows that the configurational entropy cost for binding E3 is much larger than for binding E8 (-77 +/- 15 vs. -34 +/- 11 cal K-1 mol-1). These results illustrate a case of entropy compensation in which the cost of restricting conformational degrees of freedom is to a large extent compensated by solvent release. We also show that the thermodynamic data can be used to make estimates of the surface area changes that occur upon binding. The results of the present study are consistent with previous hydrogen-deuterium exchange data, detected using 2D NMR, on the two antibody-antigen interactions. The NMR study indicated that protection from exchange is limited to the binding epitope for E8, but extends beyond the epitope for E3.(ABSTRACT TRUNCATED AT 250 WORDS)

Thermodynamic prediction of structural determinants of the molten globule state of barnase

Recently, it has been demonstrated that the enthalpy and heat capacity changes for protein folding/unfolding can be predicted rather accurately from the crystallographic or NMR solution structure of a protein. (K.P. Murphy, V. Bhakuni, D. Xie and E. Freire, Mol. Biol. 227 (1992) 293-306.) Under some conditions proteins do not unfold completely, giving rise to states in which the molecule remains in a compact configuration after denaturation. These compact denatured or molten globule states retain a hydrophobic core, exhibit residual structure and a compactness close to that of the native state. This phenomenon is reflected in the thermodynamics of the process. By using the structural parametrization of the energetics, it is possible to develop an algorithm aimed at selecting partly folded states that conform to the experimental thermodynamic constraints of the molten globule. We have applied our molten globule search algorithm to the globular protein barnase. This approach has allowed a structure based selection of a unique family of structural states that satisfy the experimental criteria of the molten globule. The prediction of the molten globule search algorithm indicates that the first helix together with most of the beta-sheet structure (beta 2, beta 3-5) and loop 5 constitute the main determinants of the molten globule intermediate, in agreement with the NMR data. These results open the prospect for an automated search of the structural determinants of the molten globule state of proteins and suggest that solvation parameters can be effectively used to probe structural states of proteins.

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the "compact denatured" or "molten globule" state. Recently is has been shown that alpha-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P.L. Biochemistry 33:1889-1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of alpha-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the beta sheet and the D helix were found to be unfolded. At temperatures below 45 degrees C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of alpha-lactalbumin [980 cal (mol.K)-1]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data.

Entropy in biological binding processes: estimation of translational entropy loss

The loss of translational degrees of freedom makes an important, unfavorable contribution to the free energy of binding. Examination of experimental values suggest that calculation of this entropy using the Sackur-Tetrode equation produces largely overestimated values. Better agreement is obtained using the cratic entropy. Theoretical considerations suggest that the volumes available for the movement of a ligand in solution and in a complex are rather similar, suggesting also that the cratic entropy provides the best estimate of the loss of translational entropy.