In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

The prototypical tryptophan synthase is a stable heterotetrameric alpha-betabeta-alpha complex. The constituting TrpA and TrpB1 subunits, which are encoded by neighboring genes in the trp operon, activate each other in a bi-directional manner. Recently, a novel class of TrpB2 proteins has been identified, whose members contain additional amino acids that might sterically prevent complex formation with TrpA. To test this hypothesis, we characterized the TrpA and TrpB proteins from Sulfolobus solfataricus. This hyperthermophilic archaeon does not contain a TrpB1 protein but instead contains two TrpB2 homologues that are encoded within (TrpB2i) and outside (TrpB2o) the trp operon. We find that TrpB2i and TrpA form a weak and transient complex during catalysis, with a uni-directional activation of TrpA by TrpB2i. In contrast, TrpB2o and TrpA do not form a detectable complex. These results suggest a model for the evolution of the tryptophan synthase in which TrpB2o, TrpB2i, and TrpB1 reflect the stepwise increase of TrpB affinity for TrpA and the refinement of functional subunit interaction, concomitant with the co-localization of the encoding genes in the trp operon.