Reversible Equilibrium Unfolding of Triosephosphate Isomerase from Trypanosoma cruzi in Guanidinium Hydrochloride Involves Stable Dimeric and Monomeric Intermediates

Stable monomers in the ancestral sequence reconstruction of the last opisthokont common ancestor of dimeric triosephosphate isomerase

Function and structure are strongly coupled in obligated oligomers such as Triosephosphate isomerase (TIM). In animals and fungi, TIM monomers are inactive and unstable. Previously, we used ancestral sequence reconstruction to study TIM evolution and found that before these lineages diverged, the last opisthokonta common ancestor of TIM (LOCATIM) was an obligated oligomer that resembles those of extant TIMs. Notably, calorimetric evidence indicated that ancestral TIM monomers are more structured than extant ones. To further increase confidence about the function, structure, and stability of the LOCATIM, in this work, we applied two different inference methodologies and the worst plausible case scenario for both of them, to infer four sequences of this ancestor and test the robustness of their physicochemical properties. The extensive biophysical characterization of the four reconstructed sequences of LOCATIM showed very similar hydrodynamic and spectroscopic properties, as well as ligand-binding energetics and catalytic parameters. Their 3D structures were also conserved. Although differences were observed in melting temperature, all LOCATIMs showed reversible urea-induced unfolding transitions, and for those that reached equilibrium, high conformational stability was estimated (ΔGTot = 40.6-46.2 kcal/mol). The stability of the inactive monomeric intermediates was also high (ΔGunf = 12.6-18.4 kcal/mol), resembling some protozoan TIMs rather than the unstable monomer observed in extant opisthokonts. A comparative analysis of the 3D structure of ancestral and extant TIMs shows a correlation between the higher stability of the ancestral monomers with the presence of several hydrogen bonds located in the "bottom" part of the barrel.

Structure and conformational stability of the triosephosphate isomerase from Zea mays. Comparison with the chemical unfolding pathways of other eukaryotic TIMs

We studied the structure, function and thermodynamic properties for the unfolding of the Triosephosphate isomerase (TIM) from Zea mays (ZmTIM). ZmTIM shows a catalytic efficiency close to the diffusion limit. Native ZmTIM is a dimer that dissociates upon dilution into inactive and unfolded monomers. Its thermal unfolding is irreversible with a T_m of 61.6 ± 1.4 °C and an activation energy of 383.4 ± 11.5 kJ mol^-1. The urea-induced unfolding of ZmTIM is reversible. Transitions followed by catalytic activity and spectroscopic properties are monophasic and superimposable, indicating that ZmTIM unfolds/refolds in a two-state behavior with an unfolding ΔG°_(H20) = 99.8 ± 5.3 kJ mol^-1. This contrasts with most other studied TIMs, where folding intermediates are common. The three-dimensional structure of ZmTIM was solved at 1.8 Å. A structural comparison with other eukaryotic TIMs shows a similar number of intramolecular and intermolecular interactions. Interestingly the number of interfacial water molecules found in ZmTIM is lower than those observed in most TIMs that show folding intermediates. Although with the available data, there is no clear correlation between structural properties and the number of equilibrium intermediates in the unfolding of TIM, the identification of such structural properties should increase our understanding of folding mechanisms.

Identification of the critical residues responsible for differential reactivation of the triosephosphate isomerases of two trypanosomes

The reactivation of triosephosphate isomerase (TIM) from unfolded monomers induced by guanidine hydrochloride involves different amino acids of its sequence in different stages of protein refolding. We describe a systematic mutagenesis method to find critical residues for certain physico-chemical properties of a protein. The two similar TIMs of Trypanosoma brucei and Trypanosoma cruzi have different reactivation velocities and efficiencies. We used a small number of chimeric enzymes, additive mutants and planned site-directed mutants to produce an enzyme from T. brucei with 13 mutations in its sequence, which reactivates fast and efficiently like wild-type (WT) TIM from T. cruzi, and another enzyme from T. cruzi, with 13 slightly altered mutations, which reactivated slowly and inefficiently like the WT TIM of T. brucei. Our method is a shorter alternative to random mutagenesis, saturation mutagenesis or directed evolution to find multiple amino acids critical for certain properties of proteins.

Interplay between Protein Thermal Flexibility and Kinetic Stability

Kinetic stability is a key parameter to comprehend protein behavior and it plays a central role to understand how evolution has reached the balance between function and stability in cell-relevant timescales. Using an approach that includes simulations, protein engineering, and calorimetry, we show that there is a clear correlation between kinetic stability determined by differential scanning calorimetry and protein thermal flexibility obtained from a novel method based on temperature-induced unfolding molecular dynamics simulations. Thermal flexibility quantitatively measures the increment of the conformational space available to the protein when energy in provided. The (β/α)₈ barrel fold of two closely related by evolution triosephosphate isomerases from two trypanosomes are used as model systems. The kinetic stability-thermal flexibility correlation has predictive power for the studied proteins, suggesting that the strategy and methodology discussed here might be applied to other proteins in biotechnological developments, evolutionary studies, and the design of protein based therapeutics.

Structural insights from a novel invertebrate triosephosphate isomerase from Litopenaeus vannamei

Triosephosphate isomerase (TIM; EC 5.3.1.1) is a key enzyme involved in glycolysis and gluconeogenesis. Glycolysis is one of the most regulated metabolic pathways, however little is known about the structural mechanisms for its regulation in non-model organisms, like crustaceans. To understand the structure and function of this enzyme in invertebrates, we obtained the crystal structure of triosephosphate isomerase from the marine Pacific whiteleg shrimp (Litopenaeus vannamei, LvTIM) in complex with its inhibitor 2-phosphogyceric acid (2-PG) at 1.7Å resolution. LvTIM assembles as a homodimer with residues 166-176 covering the active site and residue Glu166 interacting with the inhibitor. We found that LvTIM is the least stable TIM characterized to date, with the lowest range of melting temperatures, and with the lowest activation enthalpy associated with the thermal unfolding process reported. In TIMs dimer stabilization is maintained by an interaction of loop 3 by a set of hydrophobic contacts between subunits. Within these contacts, the side chain of a hydrophobic residue of one subunit fits into a cavity created by a set of hydrophobic residues in the neighboring subunit, via a "ball and socket" interaction. LvTIM presents a Cys47 at the "ball" inter-subunit contact indicating that the character of this residue is responsible for the decrease in dimer stability. Mutational studies show that this residue plays a role in dimer stability but is not a solely determinant for dimer formation.

Substrate-Induced Dimerization of Engineered Monomeric Variants of Triosephosphate Isomerase from Trichomonas vaginalis

The dimeric nature of triosephosphate isomerases (TIMs) is maintained by an extensive surface area interface of more than 1600 Ų. TIMs from Trichomonas vaginalis (TvTIM) are held in their dimeric state by two mechanisms: a ball and socket interaction of residue 45 of one subunit that fits into the hydrophobic pocket of the complementary subunit and by swapping of loop 3 between subunits. TvTIMs differ from other TIMs in their unfolding energetics. In TvTIMs the energy necessary to unfold a monomer is greater than the energy necessary to dissociate the dimer. Herein we found that the character of residue I45 controls the dimer-monomer equilibrium in TvTIMs. Unfolding experiments employing monomeric and dimeric mutants led us to conclude that dimeric TvTIMs unfold following a four state model denaturation process whereas monomeric TvTIMs follow a three state model. In contrast to other monomeric TIMs, monomeric variants of TvTIM1 are stable and unexpectedly one of them (I45A) is only 29-fold less active than wild-type TvTIM1. The high enzymatic activity of monomeric TvTIMs contrast with the marginal catalytic activity of diverse monomeric TIMs variants. The stability of the monomeric variants of TvTIM1 and the use of cross-linking and analytical ultracentrifugation experiments permit us to understand the differences between the catalytic activities of TvTIMs and other marginally active monomeric TIMs. As TvTIMs do not unfold upon dimer dissociation, herein we found that the high enzymatic activity of monomeric TvTIM variants is explained by the formation of catalytic dimeric competent species assisted by substrate binding.

Reversibility and two state behaviour in the thermal unfolding of oligomeric TIM barrel proteins

The reversible thermal unfolding of oligomeric TIM barrels results from a delicate balance of physicochemical properties related to the sequence, the native and unfolded states and the transition between them.

Triosephosphate isomerase I170V alters catalytic site, enhances stability and induces pathology in a Drosophila model of TPI deficiency

Triosephosphate isomerase (TPI) is a glycolytic enzyme which homodimerizes for full catalytic activity. Mutations of the TPI gene elicit a disease known as TPI Deficiency, a glycolytic enzymopathy noted for its unique severity of neurological symptoms. Evidence suggests that TPI Deficiency pathogenesis may be due to conformational changes of the protein, likely affecting dimerization and protein stability. In this report, we genetically and physically characterize a human disease-associated TPI mutation caused by an I170V substitution. Human TPI(I170V) elicits behavioral abnormalities in Drosophila. An examination of hTPI(I170V) enzyme kinetics revealed this substitution reduced catalytic turnover, while assessments of thermal stability demonstrated an increase in enzyme stability. The crystal structure of the homodimeric I170V mutant reveals changes in the geometry of critical residues within the catalytic pocket. Collectively these data reveal new observations of the structural and kinetic determinants of TPI Deficiency pathology, providing new insights into disease pathogenesis.

The E104D mutation increases the susceptibility of human triosephosphate isomerase to proteolysis. Asymmetric cleavage of the two monomers of the homodimeric enzyme

The deficiency of human triosephosphate isomerase (HsTIM) generates neurological alterations, cardiomyopathy and premature death. The mutation E104D is the most frequent cause of the disease. Although the wild type and mutant exhibit similar kinetic parameters, it has been shown that the E104D substitution induces perturbation of an interfacial water network that, in turn, reduces the association constant between subunits promoting enzyme inactivation. To gain further insight into the effects of the mutation on the structure, stability and function of the enzyme, we measured the sensitivity of recombinant E104D mutant and wild type HsTIM to limited proteolysis. The mutation increases the susceptibility to proteolysis as consequence of the loss of rigidity of its overall 3-D structure. Unexpectedly, it was observed that proteolysis of wild type HsTIM generated two different stable nicked dimers. One was formed in relatively short times of incubation with proteinase K; as shown by spectrometric and crystallographic data, it corresponded to a dimer containing a nicked monomer and an intact monomer. The formation of the other nicked species requires relatively long incubation times with proteinase K and corresponds to a dimer with two clipped subunits. The first species retains 50% of the original activity, whereas the second species is inactive. Collectively, we found that the E104D mutant is highly susceptible to proteolysis, which in all likelihood contributes to the pathogenesis of enzymopathy. In addition, the proteolysis data on wild type HsTIM illustrate an asymmetric conduct of the two monomers.

Structural and thermodynamic folding characterization of triosephosphate isomerases from Trichomonas vaginalis reveals the role of destabilizing mutations following gene duplication

We report the structures and thermodynamic analysis of the unfolding of two triosephosphate isomerases (TvTIM1 and TvTIM2) from Trichomonas vaginalis. Both isoforms differ by the character of four amino acids: E/Q 18, I/V 24, I/V 45, and P/A 239. Despite the high sequence and structural similarities between both isoforms, they display substantial differences in their stabilities. TvTIM1 (E18, I24, I45, and P239) is more stable and less dissociable than TvTIM2 (Q18, V24, V45, and A239). We postulate that the identities of residues 24 and 45 are responsible for the differences in monomer stability and dimer dissociability, respectively. The structural difference between both amino acids is one methyl group. In TvTIMs, residue 24 is involved in packing α-helix 1 against α-helix 2 of each monomer and residue 45 is located at the center of the dimer interface forming a "ball and socket" interplay with a hydrophobic cavity. The mutation of valine at position 45 for an alanine in TvTIM2 produces a protein that migrates as a monomer by gel filtration. A comparison with known TIM structures indicates that this kind of interplay is a conserved feature that stabilizes dimeric TIM structures. In addition, TvTIMs are located in the cytoplasm and in the membrane. As TvTIM2 is an easily dissociable dimer, the dual localization of TvTIMs may be related to the acquisition of a moonlighting activity of monomeric TvTIM2. To our knowledge, this is the simplest example of how a single amino acid substitution can provide alternative function to a TIM barrel protein.

Different catalytic properties of two highly homologous triosephosphate isomerase monomers

It is assumed that amino acid sequence differences in highly homologous enzymes would be found at the peripheral level, subtle changes that would not necessarily affect catalysis. Here, we demonstrate that, using the same set of mutations at the level of the interface loop 3, the activity of a triosephosphate isomerase monomeric enzyme is ten times higher than that of a homologous enzyme with 74% identity and 86% similarity, whereas the activity of the native, dimeric enzymes is essentially the same. This is an example of how the dimeric biological unit evolved to compensate for the intrinsic differences found at the monomeric species level. Biophysical techniques of size exclusion chromatography, dynamic light scattering, X-ray crystallography, fluorescence and circular dichroism, as well as denaturation/renaturation assays with guanidinium hydrochloride and ANS binding, allowed us to fully characterize the properties of the new monomer.

Between-species variation in the kinetic stability of TIM proteins linked to solvation-barrier free energies

Theoretical, computational, and experimental studies have suggested the existence of solvation barriers in protein unfolding and denaturation processes. These barriers are related to the finite size of water molecules and can be envisioned as arising from the asynchrony between water penetration and breakup of internal interactions. Solvation barriers have been proposed to play roles in protein cooperativity and kinetic stability; therefore, they may be expected to be subject to natural selection. We study the thermal denaturation, in the presence and in the absence of chemical denaturants, of triosephosphate isomerases (TIMs) from three different species: Trypanosoma cruzi, Trypanosoma brucei, and Leishmania mexicana. In all cases, denaturation was irreversible and kinetically controlled. Surprisingly, however, we found large differences between the kinetic denaturation parameters, with T. cruzi TIM showing a much larger activation energy value (and, consequently, much lower room-temperature, extrapolated denaturation rates). This disparity cannot be accounted for by variations in the degree of exposure to solvent in transition states (as measured by kinetic urea m values) and is, therefore, to be attributed mainly to differences in solvation-barrier contributions. This was supported by structure-energetics analyses of the transition states and by application of a novel procedure to estimate from experimental data the solvation-barrier impact at the entropy and free-energy levels. These analyses were actually performed with an extended protein set (including six small proteins plus seven variants of lipase from Thermomyces lanuginosus and spanning a wide range of activation parameters), allowing us to delineate the general trends of the solvation-barrier contributions. Overall, this work supports that proteins sharing the same structure and function but belonging to different organisms may show widely different solvation barriers, possibly as a result of different levels of the selection pressure associated with cooperativity, kinetic stability, and related factors.

Identification of a dimeric intermediate in the unfolding pathway for the calcium-binding protein S100B

The S100 proteins comprise 25 calcium-signalling members of the EF-hand protein family. Unlike typical EF-hand signalling proteins such as calmodulin and troponin-C, the S100 proteins are dimeric, forming both homo- and heterodimers in vivo. One member of this family, S100B, is a homodimeric protein shown to control the assembly of several cytoskeletal proteins and regulate phosphorylation events in a calcium-sensitive manner. Calcium binding to S100B causes a conformational change involving movement of helix III in the second calcium-binding site (EF2) that exposes a hydrophobic surface enabling interactions with other proteins such as tubulin and Ndr kinase. In several S100 proteins, calcium binding also stabilizes dimerization compared to the calcium-free states. In this work, we have examined the guanidine hydrochloride (GuHCl)-induced unfolding of dimeric calcium-free S100B. A series of tryptophan substitutions near the dimer interface and the EF2 calcium-binding site were studied by fluorescence spectroscopy and showed biphasic unfolding curves. The presence of a plateau near 1.5 M GuHCl showed the presence of an intermediate that had a greater exposed hydrophobic surface area compared to the native dimer based on increased 4,4-dianilino-1,1'-binaphthyl-5,5'-disulfonic acid fluorescence. Furthermore, (1)H-(15)N heteronuclear single quantum coherence analyses as a function of GuHCl showed significant chemical shift changes in regions near the EF1 calcium-binding loop and between the linker and C-terminus of helix IV. Together these observations show that calcium-free S100B unfolds via a dimeric intermediate.

Unraveling the mechanisms of tryptophan fluorescence quenching in the triosephosphate isomerase from Giardia lamblia

In the native state several proteins exhibit a quenching of fluorescence of their tryptophans. We studied triosephosphate isomerase from Giardia lamblia (GlTIM) to dissect the mechanisms that account for the quenching of fluorescence of its Trp. GlTIM contains four Trp per monomer (Trp75, Trp162, Trp173, and Trp196) distributed throughout the 3D structure. The fluorescence of the denatured enzyme is 3-fold higher than that of native GlTIM. To ascertain the origin of this phenomenon, single and triple mutants of Trp per Phe were made. The intrinsic fluorescence was determined, and the data were interpreted on the basis of the crystal structure of the enzyme. Our data show that the fluorescence of all Trp residues is quenched through two different mechanisms. In one, fluorescence is quenched by aromatic-aromatic interactions due to the proximity and orientation of the indole groups of Trp196 and Trp162. The magnitude of the quenching of fluorescence in Trp162 is higher than in the other three Trp. Fluorescence quenching is also due to energy transfer to the charged residues that surround Trp 75, 173 and 196. Further analysis of the fluorescence of GlTIM showed that, among TIMs from other parasites, Trp at position 12 exhibits rather unique properties.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

Thermal-unfolding Reaction of Triosephosphate Isomerase from Trypanosoma cruzi

Thermal denaturation of triosephosphate isomerase from Trypanosoma cruzi was studied by circular dicrhoism and fluorescence spectroscopies. The unfolding transition was found to be highly irreversible even at the very early stages of the reaction. Kinetic studies, allowed us to identify consecutive reactions. Firstly, only the tryptophan environment is altered. Next, changes on the secondary structure and hydrophobic surface exposure measured by 1-anilino-8-naphthalenesulfonate (ANS) binding were observed. Further conformational changes imply additional modifications on the secondary and tertiary structures and release of the hydrophobic dye leading to the formation of the unfolded state that is prone to aggregate.

Thermodynamic and kinetic characterization of the association of triosephosphate isomerase: the role of diffusion

It is known that diffusion plays a central role in the folding of small monomeric proteins and in the rigid-body association of proteins, however, the role of diffusion in the association of the folding intermediates of oligomeric proteins has been scarcely explored. In this work, catalytic activity and fluorescence measurements were used to study the effect of viscosity in the unfolding and refolding of the homodimeric enzyme triosephosphate isomerase from Saccharamyces cerevisiae. Two transitions were found by equilibrium and kinetic experiments, suggesting a three-state model with a monomeric intermediate. Glycerol barely affects DeltaG(0)(fold) whereas DeltaG(0)(assoc) becomes more favourable in the presence of the cosolvent. From 0 to 60% (v/v) glycerol, the association rate constant showed a near unitary dependence on solvent viscosity. However, at higher glycerol concentrations deviations from Kramers theory were observed. The dissociation rate constant showed a viscosity effect much higher than one. This may be related to secondary effects such as short-range glycerol-induced repulsion between monomers. Nevertheless, after comparison under isostability conditions, a slope near one was also observed for the dissociation rate. These results strongly suggest that the bimolecular association producing the native dimer is limited by diffusional events of the polypeptide chains through the solvent.

Biophysical characterization of the enzyme I of the Streptomyces coelicolor phosphoenolpyruvate:sugar phosphotransferase system

The first protein in the bacterial phosphoenolpyruvate (PEP):sugar phosphotransferase system is the homodimeric 60-kDa enzyme I (EI), which autophosphorylates in the presence of PEP and Mg2+. The conformational stability and structure of the EI from Streptomyces coelicolor, EI(sc), were explored in the absence and in the presence of its effectors by using several biophysical probes (namely, fluorescence, far-ultraviolet circular dichroism, Fourier transform infrared spectroscopy (FTIR), and differential scanning calorimetry) and computational approaches. The structure of EI(sc) was obtained by homology modeling of the isolated N- and C-terminal domains of other EI proteins. The experimental results indicate that at physiological pH, the dimeric EI(sc) had a well-folded structure; however, at low pH, EI(sc) showed a partially unfolded state with the features of a molten globule, as suggested by fluorescence, far-ultraviolet circular dichroism, FTIR, and 8-anilino-1-naphthalene-sulfonic acid binding. The thermal stability of EI(sc), in the absence of PEP and Mg2+, was maximal at pH 7. The presence of PEP and Mg2+ did not change substantially the secondary structure of the protein, as indicated by FTIR measurements. However, quenching experiments and proteolysis patterns suggest conformational changes in the presence of PEP; furthermore, the thermal stability of EI(sc) was modified depending on the effector added. Our approach suggests that thermodynamical analysis might reveal subtle conformational changes.