How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

The tunnel region at triosephosphate isomerase (TIM)’s dimer interface, distant from its catalytic site, is a target site for certain benzothiazole derivatives that inhibit TIM’s catalytic activity in Trypanosoma cruzi, the parasite that causes Chagas disease. We performed multiple 100-ns molecular-dynamics (MD) simulations and elastic network modeling (ENM) on both apo and complex structures to shed light on the still unclear inhibitory mechanism of one such inhibitor, named bt10. Within the time frame of our MD simulations, we observed stabilization of aromatic clusters at the dimer interface and enhancement of intersubunit hydrogen bonds in the presence of bt10, which point to an allosteric effect rather than destabilization of the dimeric structure. The collective dynamics dictated by the topology of TIM is known to facilitate the closure of its catalytic loop over the active site that is critical for substrate entrance and product release. We incorporated the ligand’s effect on vibrational dynamics by applying mixed coarse-grained ENM to each one of 54,000 MD snapshots. Using this computationally efficient technique, we observed altered collective modes and positive shifts in eigenvalues due to the constraining effect of bt10 binding. Accordingly, we observed allosteric changes in the catalytic loop’s dynamics, flexibility, and correlations, as well as the solvent exposure of catalytic residues. A newly (to our knowledge) introduced technique that performs residue-based ENM scanning of TIM revealed the tunnel region as a key binding site that can alter global dynamics of the enzyme.     

Dynamic properties of the N-terminal swapped dimer of ribonuclease A

Bovine pancreatic ribonuclease (RNase A) forms two 3-dimensional domain-swapped dimers with different quaternary structures. One dimer is characterized by the swapping of the C-terminal region (C-Dimer) and presents a rather loose structure. The other dimer (N-Dimer) exhibits a very compact structure with exchange of the N-terminal helix. Here we report the results of a molecular dynamics/essential dynamics (MD/ED) study carried out on the N-Dimer. This investigation, which represents the first MD/ED analysis on a three-dimensional domain-swapped enzyme, provides information on the dynamic properties of the active site residues as well as on the global motions of the dimer subunits. In particular, the analysis of the flexibility of the active site residues agrees well with recent crystallographic and site-directed mutagenesis studies on monomeric RNase A, thus indicating that domain swapping does not affect the dynamics of the active sites. A slight but significant rearrangement of N-Dimer quaternary structure, favored by the formation of additional hydrogen bonds at subunit interface, has been observed during the MD simulation. The analysis of collective movements reveals that each subunit of the dimer retains the functional breathing motion observed for RNase A. Interestingly, the breathing motion of the two subunits is dynamically coupled, as they open and close in phase. These correlated motions indicate the presence of active site intercommunications in this dimer. On these bases, we propose a speculative mechanism that may explain negative cooperativity in systems preserving structural symmetry during the allosteric transitions.     

Active site properties of monomeric triosephosphate isomerase (monoTIM) as deduced from mutational and structural studies.

MonoTIM is a stable monomeric variant of the dimeric trypanosomal enzyme triose phosphate isomerase (TIM) with less, but significant, catalytic activity. It is known that in TIM, three residues, Lys 13 (loop 1), His 95 (loop 4), and Glu 167 (loop 6) are the crucial catalytic residues. In the wild-type TIM dimer, loop 1 and loop 4 are very rigid because of tight interactions with residues of the other subunit. Previous structural studies indicate that Lys 13 and His 95 have much increased conformational flexibility in monoTIM. Using site-directed mutagenesis, it is shown here that Lys 13 and His 95 are nevertheless essential for optimal catalysis by monoTIM: monoTIM-K13A is completely inactive, although it can still bind substrate analogues, and monoTIM-H95A is 50 times less active. The best inhibitors of wild-type TIM are phosphoglycolohydroxamate (PGH) and 2-phosphoglycolate (2PG), with KI values of 8 microM and 26 microM, respectively. The affinity of the monoTIM active site for PGH has been reduced approximately 60-fold, whereas for 2PG, only a twofold weakening of affinity is observed. The mode of binding, as determined by protein crystallographic analysis of these substrate analogues, shows that, in particular, 2PG interacts with Lys 13 and His 95 in a way similar but not identical to that observed for the wild-type enzyme. This crystallographic analysis also shows that Glu 167 has the same interactions with the substrate analogues as in the wild type. The data presented suggest that, despite the absence of the second subunit, monoTIM catalyzes the interconversion of D-glyceraldehyde-3-phosphate and dihydroxyacetone phosphate via the same mechanism as in the wild type.     

Solution NMR and computer simulation studies of active site loop motion in triosephosphate isomerase

Solution NMR spin relaxation experiments and classical MD simulations are used to study the dynamics of triosephosphate isomerase (TIM) in complex with glycerol 3-phosphate (G3P). Three regions in TIM exhibit conformational transitions on the micros-ms time scale as detected by chemical exchange broadening effects in NMR spectroscopy: residue Lys 84 on helix C, located at the dimeric interface; active site loop 6; and helix G. The results indicate that the conformational exchange process affecting the residues of loop 6 is the correlated opening and closing of the loop. Distinct processes are responsible for the chemical exchange linebroadening observed in the other regions of TIM. MD simulations confirm that motions of individual residues within the active site loop are correlated and suggest that the chemical exchange processes observed for residues in helix G arise from transitions between 3(10)- and alpha-helical structures. The results of the joint NMR and MD study provide global insight into the role of conformational dynamic processes in the function of TIM.     

Molecular dynamics simulations of "loop closing" in the enzyme triose phosphate isomerase

We present molecular dynamics simulations on the active site region of dimeric triose phosphate isomerase (TIM) using the co-ordinates of native chicken muscle TIM as a starting point and performing simulations with no substrate, with dihydroxyacetone phosphate (DHAP), the natural substrate, and with dihydroxyacetone sulfate (DHAS), a substrate analog. Whereas most of the protein moves less than 1 A during the simulation, some residues in the active site loop move more than 8 A during the 10.5 picoseconds of dynamics for each of the simulations. Most interestingly, the nature of the loop motion depends on the substrate, with the largest motion found in the presence of DHAP, and only in the presence of DHAP does the loop move to "close off" the active site pocket. The final structure found for the DHAP-chicken TIM complex is qualitatively similar to that described by Alber et al. for DHAP-yeast TIM. Simulations on the monomeric protein gives insight into why the molecule is active only as a dimer.     

Coupling between Catalytic Loop Motions and Enzyme Global Dynamics

Catalytic loop motions facilitate substrate recognition and binding in many enzymes. While these motions appear to be highly flexible, their functional significance suggests that structure-encoded preferences may play a role in selecting particular mechanisms of motions. We performed an extensive study on a set of enzymes to assess whether the collective/global dynamics, as predicted by elastic network models (ENMs), facilitates or even defines the local motions undergone by functional loops. Our dataset includes a total of 117 crystal structures for ten enzymes of different sizes and oligomerization states. Each enzyme contains a specific functional/catalytic loop (10–21 residues long) that closes over the active site during catalysis. Principal component analysis (PCA) of the available crystal structures (including apo and ligand-bound forms) for each enzyme revealed the dominant conformational changes taking place in these loops upon substrate binding. These experimentally observed loop reconfigurations are shown to be predominantly driven by energetically favored modes of motion intrinsically accessible to the enzyme in the absence of its substrate. The analysis suggests that robust global modes cooperatively defined by the overall enzyme architecture also entail local components that assist in suitable opening/closure of the catalytic loop over the active site.     

Structures of the “open” and “closed” state of trypanosomal triosephosphate isomerase,as observed in a new crystal form: Implications for the reaction mechanism

The structure of trypanosomal triosephosphate isomerase (TIM)has been solved at a resolution of 2.1Å in a new crystal form grown at pH 8.8 from PEG6000. In this new crystal form (space group C2, cell dimensions 94.8 Å, 48.3 Å, 131.0 Å, 90.0°, 100.3°, 90.0°), TIM is present in a ligand-free state. The asymmetric unit consists of two TIM subunits. Each of these subunits is part of a dimer which is sitting on a crystallographic twofold axis, such that the crystal packing is formed from two TIM dimers in two distinct environments. The two constituent monomers of a given dimer are, therefore, crystallographically equivalent. In the ligand-free state of TIM in this crystal form, the two types of dimer are very similar in structure, with the flexible loops in the “Open” conformation. For one dimer (termed molecule-1), the flexible loop (loop-6) is involved in crystal contacts. Crystals of this type have been used in soaking experiments with 0.4 M ammonium sulphate (studied at 2.4 Å resolution), and with 40 μM phosphoglycolohydroxamate (studied at 2.5 Å resolution). It is found that transfer to 0.4 M ammonuum sulphate (equal to 80 times the Ki of sulphate for TIM), gives rise to significant sulphate binding at the active site of one dimer (termed molecule-2), and less significant binding at the active site of the other. In neither dimer does sulphate induce a “closed” conformation. In a mother liquor containing 40 μM phosphoglycolohydroxamate (equal to 10 times the Ki of phosphoglycolohydroxamate for TIM), an inhibitor molecule binds at the active site of only that dimer of which the flexible loop is free from crystal contacts (molecule-2). In this dimer, it induces a closed conformation. These three structures are compared and discussed with respect to the mode of binding of ligand in the active site as well as with respect to the conformational changes resulting from ligand binding. © 1993 Wiley-Liss, Inc.     

Comparison of the structures and the crystal contacts of trypanosomal triosephosphate isomerase in four different crystal forms.

Triosephosphate isomerase (TIM) is a dimeric enzyme consisting of 2 identical subunits. Trypanosomal TIM can be crystallized in 4 different spacegroups: P2(1)2(1)2(1), C2(big cell), C2(small cell), and P1. The P1 crystal form only grows in the presence of 1.4 M DMSO; there are 2 DMSO binding sites per subunit. The structures have been refined at a resolution of 1.83 A, 2.10 A, 2.13 A, and 1.80 A, respectively. In the 4 different spacegroups the TIM subunit can be observed in the context of 7 different crystallographic environments. In the C2 cells, the dimer 2-fold axis coincides with a crystallographic 2-fold axis. The similarities and differences of the 7 subunits are discussed. In 6 subunits the flexible loop (loop 6) is open, whereas in the P2(1)2(1)2(1) cell, the flexible loop of subunit 2 is in an almost closed conformation. The crystal contacts in the 4 different crystal forms are predominantly generated by polar residues in loops. A statistical analysis of the residues involved in crystal contacts shows that, in particular, serines are frequently involved in these interactions; 19% of the exposed serines are involved in crystal contacts.     

The crystal structure of an engineered monomeric triosephosphate isomerase, monoTIM: the correct modelling of an eight-residue loop

BACKGROUND: The triosephosphate isomerase (TIM) fold is found in several different classes of enzymes, most of which are oligomers; TIM itself always functions as a very tight dimer. It has recently been shown that a monomeric form of TIM ('monoTIM') can be constructed by replacing a 15-residue interface loop, loop-3, with an eight-residue fragment; modelling suggests that this should result in a short strain-free turn, resulting in the subsequent helix, helix-A3, having an additional turn at its amino terminus. RESULTS: The crystal structure of monoTIM shows that it retains the characteristic TIM-barrel (betaalpha)8-fold and that the new loop has a structure very close to that predicted. Two other interface loops, loop-1 and loop-4, which contain the active site residues Lys13 and His95, respectively, show significant changes in structure in monoTIM compared with dimeric wild-type TIM. CONCLUSION: The observed structural differences between monoTIM and wild-type TIM indicate that the dimeric appearance of TIM determines the location and conformation of two of the four catalytic residues.     

Comparison of the refined crystal structures of liganded and unliganded chicken, yeast and trypanosomal triosephosphate isomerase.

The refined crystal structures of chicken, yeast and trypanosomal triosephosphate isomerase (TIM) have been compared. TIM is known to exist in an "open" (unliganded) and "closed" (liganded) conformation. For chicken TIM only the refined open structure is available, whereas for yeast TIM and trypanosomal TIM refined structures of both the open and the closed structure have been used for this study. Comparison of these structures shows that the open structures of chicken TIM, yeast TIM and trypanosomal TIM are essentially identical. Also it is shown that the closed structures of yeast TIM and trypanosomal TIM are essentially identical. The conformational difference between the open and closed structures concerns a major shift (7 A) in loop-6. Minor shifts are observed in the two adjacent loops, loop-5 (1 A) and loop-7 (1 A). The pairwise comparison of the three different TIM barrels shows that the 105C alpha atoms of the core superimpose within 0.9 A. The sequences of these three TIMs have a pairwise sequence identity of approximately 50%. The residues that line the active site are 100% conserved. The residues interacting with each other across the dimer interface show extensive variability, but the direct hydrogen bonds between the two subunits are well conserved. The orientation of the two monomers with respect to each other is almost identical in the three different TIM structures. There are 56 (22%) conserved residues out of approximately 250 residues in 13 sequences. The functions of most of these conserved residues can be understood from the available open and closed structures of the three different TIMs. Some of these residues are quite far from the active site. For example, at a distance of 19 A from the active site there is a conserved saltbridge interaction between residues at the C-terminal ends of alpha-helix-6 and alpha-helix-7. This anchoring contrasts with the large conformational flexibility of loop-6 and loop-7 near the N termini of these helices. The flexibility of loop-6 is facilitated by a conserved large empty cavity near the N terminus of alpha-helix-6, which exists only in the open conformation.     

The essential dynamics of Cu, Zn superoxide dismutase: suggestion of intersubunit communication.

A 300-ps molecular dynamics simulation of the whole Cu, Zn superoxide dismutase dimer has been carried out in water, and the trajectory has been analyzed by the essential dynamics method. The results indicate that the motion is defined by few preferred directions identified by the first four to six eigenvectors and that the motion of the two monomers at each instant is not symmetrical. The vectors symmetrical to the eigenvectors are significantly sampled, suggesting that, on average, the motions of the two subunits will exchange. Large intra- and intersubunit motions involving different subdomains of the protein are observed. A mechanical coupling between the two subunits is also suggested, because displacements of the loops surrounding the active site in one monomer are correlated with the motion of parts of the second toward the intersubunit interface.     

Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Dynamic requirements for a functional protein hinge

The enzyme triosephosphate isomerase (TIM) is a model of catalytic efficiency. The 11 residue loop 6 at the TIM active site plays a major role in this enzymatic prowess. The loop moves between open and closed states, which facilitate substrate access and catalysis, respectively. The N and C-terminal hinges of loop 6 control this motion. Here, we detail flexibility requirements for hinges in a comparative solution NMR study of wild-type (WT) TIM and a quintuple mutant (PGG/GGG). The latter contained glycine substitutions in the N-terminal hinge at Val167 and Trp168, which follow the essential Pro166, and in the C-terminal hinge at Lys174, Thr175, and Ala176. Previous work demonstrated that PGG/GGG has a tenfold higher Km value and 10(3)-fold reduced k(cat) relative to WT with either d-glyceraldehyde 3-phosphate or dihyrdroxyacetone phosphate as substrate. Our NMR results explain this in terms of altered loop-6 dynamics in PGG/GGG. In the mutant, loop 6 exhibits conformational heterogeneity with corresponding motional rates <750 s(-1) that are an order of magnitude slower than the natural WT loop 6 motion. At the same time, nanosecond timescale motions of loop 6 are greatly enhanced in the mutant relative to WT. These differences from WT behavior occur in both apo PGG/GGG and in the form bound to the reaction-intermediate analog, 2-phosphoglycolate (2-PGA). In addition, as indicated by 1H, 15N and 13CO chemical-shifts, the glycine substitutions diminished the enzyme's response to ligand, and induced structural perturbations in apo and 2-PGA-bound forms of TIM that are atypical of WT. These data show that PGG/GGG exists in multiple conformations that are not fully competent for ligand binding or catalysis. These experiments elucidate an important principle of catalytic hinge design in proteins: structural rigidity is essential for focused motional freedom of active-site loops.     

Closed conformation of the active site loop of rabbit muscle triosephosphate isomerase in the absence of substrate: evidence of conformational heterogeneity

The active site loop of triosephosphate isomerase (TIM) exhibits a hinged-lid motion, alternating between the two well defined "open" and "closed" conformations. Until now the closed conformation had only been observed in protein complexes with substrate analogues. Here, we present the first rabbit muscle apo TIM structure, refined to 1.5A resolution, in which the active site loop is either in the open or in the closed conformation in different subunits of the enzyme. In the closed conformation described here, the lid loop residues participate in stabilizing hydrogen bonds characteristic of holo TIM structures, whereas chemical interactions observed in the open loop conformation are similar to those found in the apo structures of TIM. In the closed conformation, a number of water molecules are observed at the projected ligand atom positions that are hydrogen bonded to the active site residues. Additives used during crystallization (DMSO and Tris molecules and magnesium atoms) were modeled in the electron density maps. However, no specific binding of these molecules is observed at, or close to, the active site and the lid loop. To further investigate this unusual closed conformation of the apo enzyme, two more rabbit muscle TIM structures, one in the same and another in a different crystal form, were determined. These structures present the open lid conformation only, indicating that the closed conformation cannot be explained by crystal contact effects. To rationalize why the active site loop is closed in the absence of ligand in one of the subunits, extensive comparison with previously solved TIM structures was carried out, supported by the bulk of available experimental information about enzyme kinetics and reaction mechanism of TIM. The observation of both open and closed lid conformations in TIM crystals might be related to a persistent conformational heterogeneity of this protein in solution.     

Toward an understanding of the role of dynamics on enzymatic catalysis in lactate dehydrogenase

The motions of key residues at the substrate binding site of lactate dehydrogenase (LDH) were probed on the 10 ns to 10 ms time scale using laser-induced temperature-jump relaxation spectroscopy employing both UV fluorescence and isotope-edited IR absorption spectroscopy as structural probes. The dynamics of the mobile loop, which closes over the active site and is important for catalysis and binding, were characterized by studies of the inhibitor oxamate binding to the LDH/NADH binary complex monitoring the changes in emission of bound NADH. The bound NAD-pyruvate adduct, whose pyruvate moiety likely interacts with the same residues that interact with pyruvate in its ternary complex with LDH, served as a probe for any relative motions of active site residues against the substrate. The frequencies of its C=O stretch and -COO(-) antisymmetric stretch shift substantially should any relative motion of the polar moieties at the active site (His-195, Asp-168, Arg-109, and Arg-171) occur. The dynamics associated with loop closure are observed to involve several steps with motions from 1 to 300 microms. Apart from the "melting" of a few residues on the protein's surface, no kinetics were observed on any time scale in experiments of the bound NAD-pyr adduct although the measurements were made with a high degree of accuracy, even for final temperatures close to the unfolding transition of the protein. This is contrary to simple physical considerations and models. These results show that, once a productive protein/substrate complex is formed, the binding pocket is very rigid with very little, if any, motion apart from the mobile loop. The results also show that loop opening involves concomitant movement of the substrate out of the binding pocket.     

Comparison of backbone dynamics of monomeric and domain-swapped stefin A

Three-dimensional domain swapping has been observed in increasing number of proteins and has been implicated in the initial stages of protein aggregation, including that of the cystatins. Stefin A folds as a monomer under native conditions, while under some denaturing conditions domain-swapped dimer is formed. We have determined the backbone dynamics of the monomeric and domain-swapped dimeric forms of stefin A by (15)N relaxation using a model-free approach. The overall correlation times of the molecules were determined to be 4.6 +/- 0.1 ns and 9.2 +/- 0.2 ns for the monomer and the dimer, respectively. In the monomer, decreased order parameters indicate an increased mobility for the N-terminal trunk, the first and the second binding loops. At the opposite side of the molecule, the loop connecting the alpha-helix with strand B, the beginning of strand B and the loop connecting strands C and D show increased localized mobility. In the domain-swapped dimer, a distinctive feature of the structure is the concatenation of strands B and C into a single long beta-strand. The newly formed linker region between strands B and C, which substitutes for the first binding loop in the monomer, has order parameters typical for the remainder of the beta-strands. Thus, the interaction between subunits that occurs on domain-swapping has consequences for the dynamics of the protein at long-range from the site of conformational change, where an increased rigidity in the newly formed linker region is accompanied by an increased mobility of loops remote from that site.     

Crystal structure of coagulation factor IX-binding protein from habu snake venom at 2.6 A: implication of central loop swapping based on deletion in the linker region.

Coagulation factor IX-binding protein (IX-bp) isolated from the venom of the habu snake (Trimeresurus flavoviridis) is a disulfide-linked heterodimer consisting of homologous subunits A and B. The structure of IX-bp has been solved by X-ray crystallography at 2.6 A resolution to a crystallographic R -value of 0.181. The main-chain fold of each subunit is homologous to the carbohydrate-recognition domain of C-type lectins (C-type CRDs) except for the extended central loop. The structure is almost identical with that of factors IX and X-binding protein (IX/X-bp) as expected from the high level of amino acid sequence homology. The functional difference in ligand recognition from IX/X-bp must reside in the amino acid differences. A continuity of different amino acid residues located from the C-terminal of the second alpha-helix to the following loop forms the local conformational difference in this region between the two proteins. This loop participates in the formation of the concave surface between the two subunits, the putative binding site for the Gla-domain (gamma-carboxyglutamic acid-containing domain) of the coagulation factors. Another difference between the two proteins is in the relative disposition of subunits A and B. When the B subunits are superimposed, about a 6 degrees rotation is required for the superposition of the A subunits. A calcium ion links the second alpha-helix region to the C-terminal tail in each subunit and helps to stabilize the structure for Gla-domain binding. The interface created by the central loop swapping in the dimer IX-bp is almost identical with that seen within the monomeric C-type CRDs. This dimer forms as the result of the amino acid deletion in the linker region of the central loop of the original C-type lectins. Such a dimerization disrupts the lectin active site and creates a Gla-domain binding site, imparting functional diversity.     

A Physical Picture of Protein Dynamics and Conformational Changes

A physical model is reviewed which explains different aspects of protein dynamics consistently. At low temperatures, the molecules are frozen in conformational substates. Their average energy is 3/2RT. Solid-state vibrations occur on a time scale of femtoseconds to nanoseconds. Above a characteristic temperature, often called the dynamical transition temperature, slow modes of motions can be observed occurring on a time scale between about 140 and 1 ns. These motions are overdamped, quasidiffusive, and involve collective motions of segments of the size of an α-helix. Molecules performing these types of motion are in the “flexible state”. This state is reached by thermal activation. It is shown that these motions are essential for conformational relaxation. Based on this picture, a new approach is proposed to understand conformational changes. It connects structural fluctuations and conformational transitions.     

The critical role of the loops of triosephosphate isomerase for its oligomerization,dynamics, and functionality

Triosephosphate isomerase (TIM) catalyzes the reaction to convert dihydroxyacetone phosphate into glyceraldehyde 3‐phosphate, and vice versa. In most organisms, its functional oligomeric state is a homodimer; however, tetramer formation in hyperthermophiles is required for functional activity. The tetrameric TIM structure also provides added stability to the structure, enabling it to function at more extreme temperatures. We apply Principal Component Analysis to find that the TIM structure space is clearly divided into two groups—the open and the closed TIM structures. The distribution of the structures in the open set is much sparser than that in the closed set, showing a greater conformational diversity of the open structures. We also apply the Elastic Network Model to four different TIM structures—an engineered monomeric structure, a dimeric structure from a mesophile—Trypanosoma brucei, and two tetrameric structures from hyperthermophiles Thermotoga maritima and Pyrococcus woesei. We find that dimerization not only stabilizes the structures, it also enhances their functional dynamics. Moreover, tetramerization of the hyperthermophilic structures increases their functional loop dynamics, enabling them to function in the destabilizing environment of extreme temperatures. Computations also show that the functional loop motions, especially loops 6 and 7, are highly coordinated. In summary, our computations reveal the underlying mechanism of the allosteric regulation of the functional loops of the TIM structures, and show that tetramerization of the structure as found in the hyperthermophilic organisms is required to maintain the coordination of the functional loops at a level similar to that in the dimeric mesophilic structure.