Similarity in Shape Dictates Signature Intrinsic Dynamics Despite No Functional Conservation in TIM Barrel Enzymes.

The conservation of the intrinsic dynamics of proteins emerges as we attempt to understand the relationship between sequence, structure and functional conservation. We characterise the conservation of such dynamics in a case where the structure is conserved but function differs greatly. The triosephosphate isomerase barrel fold (TBF), renowned for its 8 ß-strand-α-helix repeats that close to form a barrel, is one of the most diverse and abundant folds found in known protein structures. Proteins with this fold have diverse enzymatic functions spanning five of six Enzyme Commission classes, and we have picked five different superfamily candidates for our analysis using elastic network models. We find that the overall shape is a large determinant in the similarity of the intrinsic dynamics, regardless of function. In particular, the ß-barrel core is highly rigid, while the α-helices that flank the ß-strands have greater relative mobility, allowing for the many possibilities for placement of catalytic residues. We find that these elements correlate with each other via the loops that link them, as opposed to being directly correlated. We are also able to analyse the types of motions encoded by the normal mode vectors of the α-helices. We suggest that the global conservation of the intrinsic dynamics in the TBF contributes greatly to its success as an enzymatic scaffold both through evolution and enzyme design.

Insights into the fold organization of TIM barrel from interaction energy based structure networks.

There are many well-known examples of proteins with low sequence similarity, adopting the same structural fold. This aspect of sequence-structure relationship has been extensively studied both experimentally and theoretically, however with limited success. Most of the studies consider remote homology or "sequence conservation" as the basis for their understanding. Recently "interaction energy" based network formalism (Protein Energy Networks (PENs)) was developed to understand the determinants of protein structures. In this paper we have used these PENs to investigate the common non-covalent interactions and their collective features which stabilize the TIM barrel fold. We have also developed a method of aligning PENs in order to understand the spatial conservation of interactions in the fold. We have identified key common interactions responsible for the conservation of the TIM fold, despite high sequence dissimilarity. For instance, the central beta barrel of the TIM fold is stabilized by long-range high energy electrostatic interactions and low-energy contiguous vdW interactions in certain families. The other interfaces like the helix-sheet or the helix-helix seem to be devoid of any high energy conserved interactions. Conserved interactions in the loop regions around the catalytic site of the TIM fold have also been identified, pointing out their significance in both structural and functional evolution. Based on these investigations, we have developed a novel network based phylogenetic analysis for remote homologues, which can perform better than sequence based phylogeny. Such an analysis is more meaningful from both structural and functional evolutionary perspective. We believe that the information obtained through the "interaction conservation" viewpoint and the subsequently developed method of structure network alignment, can shed new light in the fields of fold organization and de novo computational protein design.

Role of Protein in Fungal Biomineralization of Copper Carbonate Nanoparticles.

Biomineralization processes are of key importance in the biogeochemical cycling of metals and other elements by microorganisms, and several studies have highlighted the potential applications of nanoparticle synthesis via biomineralization. The roles played by proteins in the transformation and biologically induced biomineralization of metals by microorganisms is not well understood, despite the interactions of protein and nanoparticles at mineral interfaces attracting much interest in various emerging fields for novel biomaterial synthesis. Here, we have elucidated the association and involvement of fungal proteins in the formation of biogenic copper carbonate nanoparticles (CuNPs) using a carbonate-enriched biomass-free ureolytic fungal culture supernatant. Proteomic analysis was conducted that identified the major proteins present in the culture supernatant. Of the proteins identified, triosephosphate isomerase (TPI) exhibited a strong affinity to the CuNPs, and the impact of purified TPI on CuNP formation was studied in detail. The combined use of scanning electron microscopy (SEM) and confocal laser scanning microscopy (CLSM) confirmed that TPI played an important role in controlling the morphology and structure of the nanomaterials. Fourier transform infrared spectroscopy (FTIR) was applied to examine conformational changes of the proteins to further clarity the interaction mechanisms with CuNPs during biomineralization. Such analyses revealed unfolding of proteins on the mineral surface and an increase in ß sheets within the protein structure. These results extend understanding of how microbial systems can influence biomineral formation through protein secretion, the mechanisms involved in formation of complex protein/inorganic systems, and provide useful guidelines for the synthesis of inorganic-protein based nanomaterials.

Probing the flexibility of large conformational changes in protein structures through local perturbations.

Protein conformational changes and dynamic behavior are fundamental for such processes as catalysis, regulation, and substrate recognition. Although protein dynamics have been successfully explored in computer simulation, there is an intermediate-scale of motions that has proven difficult to simulate - the motion of individual segments or domains that move independently of the body the protein. Here, we introduce a molecular-dynamics perturbation method, the Rotamerically Induced Perturbation (RIP), which can generate large, coherent motions of structural elements in picoseconds by applying large torsional perturbations to individual sidechains. Despite the large-scale motions, secondary structure elements remain intact without the need for applying backbone positional restraints. Owing to its computational efficiency, RIP can be applied to every residue in a protein, producing a global map of deformability. This map is remarkably sparse, with the dominant sites of deformation generally found on the protein surface. The global map can be used to identify loops and helices that are less tightly bound to the protein and thus are likely sites of dynamic modulation that may have important functional consequences. Additionally, they identify individual residues that have the potential to drive large-scale coherent conformational change. Applying RIP to two well-studied proteins, Dihdydrofolate Reductase and Triosephosphate Isomerase, which possess functionally-relevant mobile loops that fluctuate on the microsecond/millisecond timescale, the RIP deformation map identifies and recapitulates the flexibility of these elements. In contrast, the RIP deformation map of alpha-lytic protease, a kinetically stable protein, results in a map with no significant deformations. In the N-terminal domain of HSP90, the RIP deformation map clearly identifies the ligand-binding lid as a highly flexible region capable of large conformational changes. In the Estrogen Receptor ligand-binding domain, the RIP deformation map is quite sparse except for one large conformational change involving Helix-12, which is the structural element that allosterically links ligand binding to receptor activation. RIP analysis has the potential to discover sites of functional conformational changes and the linchpin residues critical in determining these conformational states.

EMatch: discovery of high resolution structural homologues of protein domains in intermediate resolution cryo-EM maps.

Cryo-EM has become an increasingly powerful technique for elucidating the structure, dynamics, and function of large flexible macromolecule assemblies that cannot be determined at atomic resolution. However, due to the relatively low resolution of cryo-EM data, a major challenge is to identify components of complexes appearing in cryo-EM maps. Here, we describe EMatch, a novel integrated approach for recognizing structural homologues of protein domains present in a 6-10 A resolution cryo-EM map and constructing a quasi-atomic structural model of their assembly. The method is highly efficient and has been successfully validated on various simulated data. The strength of the method is demonstrated by a domain assembly of an experimental cryo-EM map of native GroEL at 6 A resolution.

Refined 1.83 A structure of trypanosomal triosephosphate isomerase crystallized in the presence of 2.4 M-ammonium sulphate. A comparison with the structure of the trypanosomal triosephosphate isomerase-glycerol-3-phosphate complex.

Triosephosphate isomerase (TIM) is a dimeric glycolytic enzyme. TIM from Trypanosoma brucei brucei has been crystallized at pH 7.0 in 2.4 M-ammonium sulphate. The well-diffracting crystals have one dimer per asymmetric unit. The structure has been refined at 1.83 A resolution with an R-factor of 18.3% for all data between 6 A and 1.83 A (37,568 reflections). The model consists of 3778 protein atoms and 297 solvent atoms. Subunit 1 is involved in considerably more crystal contacts than subunit 2. Correlated with these differences in crystal packing is the observation that only in the active site of subunit 2 is a sulphate ion bound. Furthermore, significant differences with respect to structure and flexibility are observed in three loops near the active site. In particular, there is a 7 A positional difference of the tip of the flexible loop (loop 6) when comparing subunit 1 and subunit 2. Also, the neighbouring loops (loop 5 and loop 7) have significantly different conformations and flexibility. In subunit 1, loop 6 is in an "open" conformation, in subunit 2, loop 6 is in an "almost closed" conformation. Only in the presence of a phosphate-containing ligand, such as glycerol-3-phosphate, does loop 6 take up the "closed" conformation. Loop 6 and loop 7 (and also to some extent loop 5) are rather flexible in the almost closed conformation, but well defined in the open and closed conformations. The closing of loop 6 (167 to 180), as observed in the almost closed conformation, slightly changes the main-chain conformation of the catalytic glutamate, Glu167, leading to a change of the chi 1 angle of this residue from approximately -60 degrees to approximately 60 degrees and the weakening of the hydrogen bonds between its polar side-chain atoms and Ser96. In the closed conformation, in the presence of glycerol-3-phosphate, the main-chain atoms of Glu167 remain in the same position as in the almost closed conformation, but the side-chain has rotated around the CA-CB bond changing chi 1 from approximately 60 degrees to approximately -60 degrees. In this new position the hydrogen bonding to Ser96 is completely lost and also a water-mediated salt bridge between OE2(Glu167) and NE(Arg99) is lost. Comparison of the two independently refined subunits, showed that the root-mean-square deviation for all 249 CA atoms is 0.9 A; for the CA atoms of the beta-strands this is only 0.2 A. The average B-factor for all subunit 1 and subunit 2 atoms is 20 A2 and 25 A2, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)

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.

Structure and Stability of the Dimeric Triosephosphate Isomerase from the Thermophilic Archaeon Thermoplasma acidophilum.

Thermoplasma acidophilum is a thermophilic archaeon that uses both non-phosphorylative Entner-Doudoroff (ED) pathway and Embden-Meyerhof-Parnas (EMP) pathway for glucose degradation. While triosephosphate isomerase (TPI), a well-known glycolytic enzyme, is not involved in the ED pathway in T. acidophilum, it has been considered to play an important role in the EMP pathway. Here, we report crystal structures of apo- and glycerol-3-phosphate-bound TPI from T. acidophilum (TaTPI). TaTPI adopts the canonical TIM-barrel fold with eight α-helices and parallel eight ß-strands. Although TaTPI shares ~30% sequence identity to other TPIs from thermophilic species that adopt tetrameric conformation for enzymatic activity in their harsh physiological environments, TaTPI exists as a dimer in solution. We confirmed the dimeric conformation of TaTPI by analytical ultracentrifugation and size-exclusion chromatography. Helix 5 as well as helix 4 of thermostable tetrameric TPIs have been known to play crucial roles in oligomerization, forming a hydrophobic interface. However, TaTPI contains unique charged-amino acid residues in the helix 5 and adopts dimer conformation. TaTPI exhibits the apparent Td value of 74.6°C and maintains its overall structure with some changes in the secondary structure contents at extremely acidic conditions (pH 1-2). Based on our structural and biophysical analyses of TaTPI, more compact structure of the protomer with reduced length of loops and certain patches on the surface could account for the robust nature of Thermoplasma acidophilum TPI.

An interface point-mutation variant of triosephosphate isomerase is compactly folded and monomeric at low protein concentrations.

Wild-type trypanosomal triosephosphate isomerase (wtTIM) is a very tight dimer. The interface residue His-47 of wtTIM has been mutated into an asparagine. Ultracentrifugation data show that this variant (H47N) only dimerises at protein concentrations above 3 mg/ml. H47N has been characterised at a protein concentration where it is predominantly a monomer. Circular dichroism measurements in the near-UV and far-UV show that this monomer is a compactly folded protein with secondary structure similar as in wtTIM. The thermal stability of the monomeric H47N is decreased compared to wtTIM: temperature gradient gel electrophoresis (TGGE) measurements give Tm-values of 41 degrees C for wtTIM, whereas the Tm-value for the monomeric form of H47N is approximately 7 degrees C lower.

Simulations of the folding pathway of triose phosphate isomerase-type alpha/beta barrel proteins.

Simulations of the folding pathways of two large alpha/beta proteins, the alpha subunit of tryptophan synthase and triose phosphate isomerase, are reported using the knight's walk lattice model of globular proteins and Monte Carlo dynamics. Starting from randomly generated unfolded states and with no assumptions regarding the nature of the folding intermediates, for the tryptophan synthase subunit these simulations predict, in agreement with experiment, the existence and location of a stable equilibrium intermediate comprised of six beta strands on the amino terminus of the molecule. For the case of triose phosphate isomerase, the simulations predict that both amino- and carboxyl-terminal intermediates should be observed. In a significant modification of previous lattice models, this model includes a full heavy atom side chain description and is capable of representing native conformations at the level of 2.5- to 3-A rms deviation for the C alpha positions, as compared to the crystal structure. With a well-balanced compromise between accuracy of the protein description and the computer requirements necessary to perform simulations spanning biologically significant amounts of time, the lattice model described here brings the possibility of studying important biological processes to present-day computers.

Inhibition of triosephosphate isomerase from Trypanosoma brucei with cyclic hexapeptides.

Two series of oligopeptides have been synthesized. Their effects on the activity of purified triosephosphate isomerase from Trypanosoma brucei and various other organisms have been studied. Using detailed three-dimensional structure information, the first series consisted of both cyclic and linear hydrophilic peptides that were designed to mimic the beta turns of the subunit interface loops of the trypanosome triosephosphate isomerase dimer. None of these exerted any inhibitory effect. The second series consisted of more hydrophobic cyclic peptides, originally designed to inhibit a hepatic transport system. Several of these were very effective in inhibiting the trypanosome triosephosphate isomerase, but not the homologous enzymes from rabbit, dog, yeast or Escherichia coli. The most active peptide, cyclo[-Trp-Phe-D-Pro-Phe-Phe-Lys(Z)-], exerted 50% inhibitory activity at a concentration of 3 microM. The nature of the inhibitory action of one of these compounds cyclo[-Trp-Tyr(OSO3Na)-D-Pro-Phe-Thr(OSO3Na)-Lys(Z)-] was studied in more detail. Its inhibition was noncompetitive and reversible and more than one peptide was able to bind/active site.

A helix-turn-strand structural motif common in alpha-beta proteins.

By exhaustive structural comparisons, we have found that about one-third of the alpha-helix-turn-beta-strand polypeptides in alpha-beta barrel domains share a common structural motif. The chief characteristics of this motif are that first, the geometry of the turn between the alpha-helix and the beta-strand is somewhat constrained, and second, the beta-strand contains a hydrophobic patch that fits into a hydrophobic pocket on the alpha-helix. The geometry of the turn does not seem to be a major determinant of the alpha-beta unit, because the turns vary in length from four to six residues. However, the motif does not occur when there are few constraints on the geometry of the turn-for instance, when the turns between the alpha-helix and the beta-strands are very long. It also occurs much less frequently in flat-sheet alpha-beta proteins, where the topology is much less regular and the amount of twist on the sheet varies considerably more than in the barrel proteins. The motif may be one of the basic building blocks from which alpha-beta barrels are constructed.