Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Computational and Experimental Characterization of dVHL Establish a Drosophila Model of VHL Syndrome

The von Hippel-Lindau (VHL) cancer syndrome is associated with mutations in the VHL gene. The pVHL protein is involved in response to changes in oxygen availability as part of an E3-ligase that targets the Hypoxia-Inducible Factor for degradation. pVHL has a molten globule configuration with marginal thermodynamic stability. The cancer-associated mutations further destabilize it. The Drosophila homolog, dVHL, has relatively low sequence similarity to pVHL, and is also involved in regulating HIF1-α. Using in silico, in vitro and in vivo approaches we demonstrate high similarity between the structure and function of dVHL and pVHL. These proteins have a similar fold, secondary and tertiary structures, as well as thermodynamic stability. Key functional residues in dVHL are evolutionary conserved. This structural homology underlies functional similarity of both proteins, evident by their ability to bind their reciprocal partner proteins, and by the observation that transgenic pVHL can fully maintain normal dVHL-HIF1-α downstream pathways in flies. This novel transgenic Drosophila model is thus useful for studying the VHL syndrome, and for testing drug candidates to treat it.     

Fold change in evolution of protein structures

Typically, protein spatial structures are more conserved in evolution than amino acid sequences. However, the recent explosion of sequence and structure information accompanied by the development of powerful computational methods led to the accumulation of examples of homologous proteins with globally distinct structures. Significant sequence conservation, local structural resemblance, and functional similarity strongly indicate evolutionary relationships between these proteins despite pronounced structural differences at the fold level. Several mechanisms such as insertions/deletions/substitutions, circular permutations, and rearrangements in beta-sheet topologies account for the majority of detected structural irregularities. The existence of evolutionarily related proteins that possess different folds brings new challenges to the homology modeling techniques and the structure classification strategies and offers new opportunities for protein design in experimental studies.     

The role of zinc in the stability of the marginally stable IscU scaffold protein

Understanding the factors that determine protein stability is interesting because it directly reflects the evolutionary pressure coming from function and environment. Here, we have combined experimental and computational methods to study the stability of IscU, a bacterial scaffold protein highly conserved in most organisms and an essential component of the iron–sulfur cluster biogenesis pathway. We demonstrate that the effect of zinc and its consequence strongly depend on the sample history. IscU is a marginally stable protein at low ionic strength to the point that undergoes cold denaturation at around −8°C with a corresponding dramatic decrease of enthalpy, which is consistent with the fluxional nature of the protein. Presence of constitutively bound zinc appreciably stabilizes the IscU fold, whereas it may cause protein aggregation when zinc is added back posthumously. We discuss how zinc coordination can be achieved by different side chains spatially available and all competent for tetrahedral coordination. The individual absence of some of these residues can be largely compensated by small local rearrangements of the others. We discuss the potential importance of our findings in vitro for the function in vivo of the protein.     

Thermodynamic properties distinguish human mitochondrial aspartyl-tRNA synthetase from bacterial homolog with same 3D architecture

In the mammalian mitochondrial translation apparatus, the proteins and their partner RNAs are coded by two genomes. The proteins are nuclear-encoded and resemble their homologs, whereas the RNAs coming from the rapidly evolving mitochondrial genome have lost critical structural information. This raises the question of molecular adaptation of these proteins to their peculiar partner RNAs. The crystal structure of the homodimeric bacterial-type human mitochondrial aspartyl-tRNA synthetase (DRS) confirmed a 3D architecture close to that of Escherichia coli DRS. However, the mitochondrial enzyme distinguishes by an enlarged catalytic groove, a more electropositive surface potential and an alternate interaction network at the subunits interface. It also presented a thermal stability reduced by as much as 12°C. Isothermal titration calorimetry analyses revealed that the affinity of the mitochondrial enzyme for cognate and non-cognate tRNAs is one order of magnitude higher, but with different enthalpy and entropy contributions. They further indicated that both enzymes bind an adenylate analog by a cooperative allosteric mechanism with different thermodynamic contributions. The larger flexibility of the mitochondrial synthetase with respect to the bacterial enzyme, in combination with a preserved architecture, may represent an evolutionary process, allowing nuclear-encoded proteins to cooperate with degenerated organelle RNAs.     

Energetics of peptide recognition by the second PDZ domain of human protein tyrosine phosphatase 1E

Formation of protein-protein assemblies is essential in maintaining cell structure and function. Conservation of structural motifs and binding sites is the result of evolutionary pressure for solutions compatible with both molecular economy and regulation. PDZ domains are a typical example: A conserved fold governs specificity toward recognition of C-terminal protein sequences by small sequential and/or structural deviations within a canonical binding mode. The energetic principles underlying the strength and specificity of PDZ-protein interactions are practically unknown. We use the second PDZ domain (PDZ2) of the human protein tyrosine phosphatase (hPTP1E) as a model to study the energetics of peptide binding to a class I PDZ domain. Calorimetric experiments reveal the enthalpy, entropy, and heat capacity changes accompanying PDZ2 binding to the C-terminal pentadecapeptide derived from the guanine nucleotide exchange factor RA-GEF2. Association is driven by favorable enthalpy and entropy changes below 18 degrees C. Above that temperature the entropy change opposes complex formation. Structure-based predictions poorly reproduce the observed thermodynamic profile of the PDZ-peptide complex. On the basis of MD simulations and experimental findings by others we suggest that changes in the dynamics of the PDZ domain upon peptide binding make a large contribution to the observed thermodynamic parameters. Possible impacts of subtle, ligand-induced structural "stiffening" of PDZ domains are discussed. In our hands, the C-terminal segment of the tumor suppressor APC binds much less tightly to PDZ2 than what has been proposed earlier from surface plasmon resonance experiments.     

How Do Thermophilic Proteins and Proteomes Withstand High Temperature?

We attempt to understand the origin of enhanced stability in thermophilic proteins by analyzing thermodynamic data for 116 proteins, the largest data set achieved to date. We compute changes in entropy and enthalpy at the convergence temperature where different driving forces are maximally decoupled, in contrast to the majority of previous studies that were performed at the melting temperature. We find, on average, that the gain in enthalpy upon folding is lower in thermophiles than in mesophiles, whereas the loss in entropy upon folding is higher in mesophiles than in thermophiles. This implies that entropic stabilization may be responsible for the high melting temperature, and hints at residual structure or compactness of the denatured state in thermophiles. We find a similar trend by analyzing a homologous set of proteins classified based only on the optimum growth temperature of the organisms from which they were extracted. We find that the folding free energy at the temperature of maximal stability is significantly more favorable in thermophiles than in mesophiles, whereas the maximal stability temperature itself is similar between these two classes. Furthermore, we extend the thermodynamic analysis to model the entire proteome. The results explain the high optimal growth temperature in thermophilic organisms and are in excellent quantitative agreement with full thermal growth rate data obtained in a dozen thermophilic and mesophilic organisms.     

Two-dimensional differential scanning calorimetry: simultaneous resolution of intrinsic protein structural energetics and ligand binding interactions by global linkage analysis.

A general theoretical development for the design and analysis of two-dimensional thermal stability surfaces of proteins is presented. The surfaces are generated from multiple excess heat capacity profiles ( vs T) obtained at varying concentrations of an interacting ligand. The energetics of both the intrinsic protein stability and the protein-ligand interaction are simultaneously resolved by employing statistical thermodynamic models in global linkage analysis. This formalism allows resolution of the intrinsic protein folding-unfolding parameters (enthalpy, entropy, and heat capacity changes) as well as the ligand interaction parameters (binding stoichiometry, enthalpy, entropy, and heat capacity changes). The theory has been applied to the case of ribonuclease A and its interaction with cytidine-2'-monophosphate. The accuracy of the thermodynamic parameters obtained by this approach compares within error with those parameters that can be obtained by direct measurements.     

Structural/functional assignment of unknown bacteriophage T4 proteins by iterative database searches

Among the total of 274 orfs within bacteriophage T4, only half have been reasonably well characterized, and the functions of the rest have remained obscure. In order to predict the molecular functions of the orfs, a position-specific iterated (PSI)-BLAST search of bacteriophage T4 against the sequence database of known 3D structures was carried out. PSI-BLAST is one of the most powerful iterative sequence search methods using multiple sequence alignment, with the ability to detect many more proteins with distant homology than standard pairwise methods. The 3D structures of proteins are considered to be better preserved than the sequences, and the detected distantly homologous proteins are likely to possess highly similar 3D structures. Thirteen orfs of phage T4, whose homologues were not detected by standard pairwise methods, were found to have significantly homologous counterparts by this method. The plausibility of the results was confirmed by checking whether important residues at substrate/ligand-binding sites were conserved. Among them, two orfs, vs.1 and e.1, which are similar to Escherichia coli lytic enzyme and MutT protein, respectively, had not been studied previously. Also, gp rIIA, a rapid lysis protein, whose gene structure had been intensively studied during the development of molecular biology in the 1950s and yet whose molecular function remains unknown, has an N-terminal domain that is significantly similar to the N-terminal region of the heat shock protein Hsp90.     

Why water-soluble, compact, globular proteins have similar specific enthalpies of unfolding at 110 degrees C.

The changes in free energy, enthalpy, and entropy of unfolding have been measured for many water-soluble, compact, globular proteins by a number of workers. In principle, a wide range in stability could be achieved by proteins, as measured by the free energy of unfolding; in practice, evolution only allows a narrow range in this quantity. Proteins are only marginally stable at room temperature for many possible reasons, including ensuring that folding is reversible and polypeptide chains are not trapped in incorrectly folded structures. Many of these proteins have approximately the same values of enthalpy of unfolding around 110 degrees C. We show here that this arises because the change in entropy of unfolding at room temperature and the change in heat capacity on unfolding, which governs the temperature variation of the enthalpy and entropy, both vary with the magnitude of the hydrophobic effect in the protein. As all these proteins have evolved to achieve similar stabilities at room temperature, the enthalpy of unfolding will also vary with the size of the hydrophobic effect in the protein. A consequence of this is that curves of the specific unfolding enthalpy against temperature for different proteins intersect around 110 degrees C. A similar conclusion, on the basis of similar melting points rather than similar free energies of unfolding, has been reached independently by Baldwin and Muller (R. L. Baldwin, personal communication).     

Structure determination of archaea-specific ribosomal protein L46a reveals a novel protein fold

Three archaea-specific ribosomal proteins recently identified show no sequence homology with other known proteins. Here we determined the structure of L46a, the most conserved one among the three proteins, from Sulfolobus solfataricus P2 using NMR spectroscopy. The structure presents a twisted β-sheet formed by the N-terminal part and two helices at the C-terminus. The L46a structure has a positively charged surface which is conserved in the L46a protein family and is the potential rRNA-binding site. Searching homologous structures in Protein Data Bank revealed that the structure of L46a represents a novel protein fold. The backbone dynamics identified by NMR relaxation experiments reveal significant flexibility at the rRNA binding surface. The potential position of L46a on the ribosome was proposed by fitting the structure into a previous electron microscopy map of the ribosomal 50S subunit, which indicated that L46a contacts to domain I of 23S rRNA near a multifunctional ribosomal protein L7ae.     

Predicting the effect of a point mutation on a protein fold: the villin and advillin headpieces and their Pro62Ala mutants

Homology modeling of unknown proteins is based on the assumption that highly similar sequences are likely to share the same fold. However, this does not provide any information on the stability of a given fold, which is ultimately determined by the subtle interplay of enthalpic and entropic contributions. Herein it is shown that ab initio atomistic simulations can be used to predict the effect of point mutations on the stability of a protein fold. The calculations indicate that the fold stabilities of two proteins of similar sequence and identical fold, the villin and advillin C-terminal headpiece fragments, are different and that the same P62A point mutation has a dramatic effect on the fold of villin but a minor one on that of advillin. These predictions were subsequently validated by NMR and CD experiments.     

A single disulfide bond restores thermodynamic and proteolytic stability to an extensively mutated protein

The potential for engineering stable proteins with multiple amino acid substitutions was explored. Eleven lysine, five methionine, two tryptophan, one glycine, and three threonine substitutions were simultaneously made in barley chymotrypsin inhibitor-2 (CI-2) to substantially improve the essential amino acid content of the protein. These substitutions were chosen based on the three-dimensional structure of CI-2 and an alignment of homologous sequences. The initial engineered protein folded into a wild-type-like structure, but had a free energy of unfolding of only 2.2 kcal/mol, considerably less than the wild-type value of 7.5 kcal/mol. Restoration of the lysine mutation at position 67 to the wild-type arginine increased the free energy of unfolding to 3.1 kcal/mol. Subsequent cysteine substitutions at positions 22 and 82 resulted in disulfide bond formation and a protein with nearly wild-type thermodynamic stability (7.0 kcal/mol). None of the engineered proteins retained inhibitory activity against chymotrypsin or elastase, and all had substantially reduced inhibitory activity against subtilisin. The proteolytic stabilities of the proteins correlated with their thermodynamic stabilities. Reduction of the disulfide bond resulted in substantial loss of both thermodynamic and proteolytic stabilities, confirming that the disulfide bond, and not merely the cysteine substitutions, was responsible for the increased stability. We conclude that it is possible to replace over a third of the residues in CI-2 with minimal disruption of stability and structural integrity.     

Positive selection dictates the choice between kinetic and thermodynamic protein folding and stability in subtilases

Subtilisin E (SbtE) is a member of the ubiquitous superfamily of serine proteases called subtilases and serves as a model for understanding propeptide-mediated protein folding mechanisms. Unlike most proteins that adopt thermodynamically stable conformations, the native state of SbtE is trapped into a kinetically stable conformation. While kinetic stability offers distinct functional advantages to the native state, the constraints that dictate the selection between kinetic and thermodynamic folding and stability remain unknown. Using highly conserved subtilases, we demonstrate that adaptive evolution of sequence dictates selection of folding pathways. Intracellular and extracellular serine proteases (ISPs and ESPs, respectively) constitute two subfamilies within the family of subtilases that have highly conserved sequences, structures, and catalytic activities. Our studies on the folding pathways of subtilisin E (SbtE), an ESP, and its homologue intracellular serine protease 1 (ISP1), an ISP, show that although topology, contact order, and hydrophobicity that drive protein folding reactions are conserved, ISP1 and SbtE fold through significantly different pathways and kinetics. While SbtE absolutely requires the propeptide to fold into a kinetically trapped conformer, ISP1 folds to a thermodynamically stable state more than 1 million times faster and independent of a propeptide. Furthermore, kinetics establish that ISP1 and SbtE fold through different intermediate states. An evolutionary analysis of folding constraints in subtilases suggests that observed differences in folding pathways may be mediated through positive selection of specific residues that map mostly onto the protein surface. Together, our results demonstrate that closely related subtilases can fold through distinct pathways and mechanisms, and suggest that fine sequence details can dictate the choice between kinetic and thermodynamic folding and stability.     

A Thermodynamic Atlas of Proteomes Reveals Energetic Innovation across the Tree of Life

Protein stability is a fundamental molecular property enabling organisms to adapt to their biological niches. How this is facilitated and whether there are kingdom specific or more general universal strategies are unknown. A principal obstacle to addressing this issue is that the vast majority of proteins lack annotation, specifically thermodynamic annotation, beyond the amino acid and chromosome information derived from genome sequencing. To address this gap and facilitate future investigation into large-scale patterns of protein stability and dynamics within and between organisms, we applied a unique ensemble-based thermodynamic characterization of protein folds to a substantial portion of extant sequenced genomes. Using this approach, we compiled a database resource focused on the position-specific variation in protein stability. Interrogation of the database reveals: 1) domains of life exhibit distinguishing thermodynamic features, with eukaryotes particularly different from both archaea and bacteria; 2) the optimal growth temperature of an organism is proportional to the average apolar enthalpy of its proteome; 3) intrinsic disorder content is also proportional to the apolar enthalpy (but unexpectedly not the predicted stability at 25 °C); and 4) secondary structure and global stability information of individual proteins is extractable. We hypothesize that wider access to residue-specific thermodynamic information of proteomes will result in deeper understanding of mechanisms driving functional adaptation and protein evolution. Our database is free for download at https://afc-science.github.io/thermo-env-atlas/ (last accessed January 18, 2022).     

Life at low temperatures: is disorder the driving force?

The thermodynamic characterization of various biological systems from psychrophiles points to a larger entropic contribution when compared to the corresponding mesophilic or (hyper) thermophilic counterparts, either at the level of the macromolecules (thermodynamic and kinetic stabilities) or of their function (ligand binding, catalytic activity). It is suggested here that in an environment characterized by a low heat content (enthalpy) and at temperatures that strongly slowdown molecular motions, the cold-adapted biological systems rely on a larger disorder to maintain macromolecular dynamics and function. Such pre-eminent involvement of entropy is observed in the experimental results and, from a macroscopic point of view, is also reflected for instance by the steric hindrances introduced by cis-unsaturated and branched lipids to maintain membrane fluidity, by the loose conformation of psychrophilic proteins or by the local destabilization of tRNA by dihydrouridine in psychrophilic bacteria.     

A study on the enthalpy-entropy compensation in protein unfolding

A large number of thermodynamic data including the free energy, enthalpy, entropy, and heat capacity changes were collected for the denaturation of various proteins. Regression indicated that remarkable enthalpy-entropy compensation occurred in protein unfolding, which meant that the change in enthalpy was almost compensated by a corresponding change in entropy resulting in a smaller net free energy change. This behavior was proposed to result from the water molecule reorganization, which contributed significantly to the enthalpy and entropy changes but little to the free energy change in protein unfolding. It turned out that the enthalpy-entropy compensation could provide novel insights into the problem of enthalpy and entropy convergence in protein unfolding.     

Thermodynamic studies of the core histones: stability of the octamer subunits is not altered by removal of their terminal domains.

We have investigated the role of the labile terminal domains of the core histones on the stability of the subunits of the protein core of the nucleosome by studying the thermodynamic behavior of the products of limited trypsin digestion of these subunits. The thermal stabilities of the truncated H2A-H2B dimer and the truncated (H3-H4)/(H3-H4)(2) system were studied by high-sensitivity differential scanning calorimetry and circular dichroism spectroscopy. The thermal denaturation of the truncated H2A-H2B dimer at pH 6.0 and low ionic strength is centered at 47.3 degrees C. The corresponding enthalpy change is 35 kcal/mol of 11.5 kDa monomer unit, and the heat capacity change upon unfolding is 1.2 kcal/(K mol of 11.5 kDa monomer unit). At pH 4.5 and low ionic strength, the truncated (H3-H4)/(H3-H4)(2) system, like its full-length counterpart, is quantitatively dissociated into two truncated H3-H4 dimers. The thermal denaturation of the truncated H3-H4 dimer is characterized by the presence of a single calorimetric peak centered at 60 degrees C. The enthalpy change is 25 kcal/mol of 10 kDa monomer unit, and the change in heat capacity upon unfolding is 0.5 kcal/(K mol of 10 kDa monomer unit). The thermal stabilities of both types of truncated dimers exhibit salt and pH dependencies similar to those of the full-length proteins. Finally, like their full-length counterparts, both truncated core histone dimers undergo thermal denaturation as highly cooperative units, without the involvement of any significant population of melting intermediates. Therefore, removal of the histone "tails" does not generally affect the thermodynamic behavior of the subunits of the core histone complex, indicating that the more centrally located regions of the histone fold and the extra-fold structured elements are primarily responsible for their stability and responses to parameters of their chemical microenvironment.