- Open Access
Structural analysis on mutation residues and interfacial water molecules for human TIM disease understanding
© Li et al.; licensee BioMed Central Ltd. 2013
Published: 22 October 2013
Human triosephosphate isomerase (HsTIM) deficiency is a genetic disease caused often by the pathogenic mutation E104D. This mutation, located at the side of an abnormally large cluster of water in the inter-subunit interface, reduces the thermostability of the enzyme. Why and how these water molecules are directly related to the excessive thermolability of the mutant have not been investigated in structural biology.
This work compares the structure of the E104D mutant with its wild type counterparts. It is found that the water topology in the dimer interface of HsTIM is atypical, having a "wet-core-dry-rim" distribution with 16 water molecules tightly packed in a small deep region surrounded by 22 residues including GLU104. These water molecules are co-conserved with their surrounding residues in non-archaeal TIMs (dimers) but not conserved across archaeal TIMs (tetramers), indicating their importance in preserving the overall quaternary structure. As the structural permutation induced by the mutation is not significant, we hypothesize that the excessive thermolability of the E104D mutant is attributed to the easy propagation of atoms' flexibility from the surface into the core via the large cluster of water. It is indeed found that the B factor increment in the wet region is higher than other regions, and, more importantly, the B factor increment in the wet region is maintained in the deeply buried core. Molecular dynamics simulations revealed that for the mutant structure at normal temperature, a clear increase of the root-mean-square deviation is observed for the wet region contacting with the large cluster of interfacial water. Such increase is not observed for other interfacial regions or the whole protein. This clearly suggests that, in the E104D mutant, the large water cluster is responsible for the subunit interface flexibility and overall thermolability, and it ultimately leads to the deficiency of this enzyme.
Our study reveals that a large cluster of water buried in protein interfaces is fragile and high-maintenance, closely related to the structure, function and evolution of the whole protein.
Human triosephosphate isomerase (TIM) deficiency, first reported by Schneider et al. in 1965 , is a genetic disease caused by the dysfunction of TIM. Clinical phenotypic characteristics of this disease include chronic hemolytic anemia and progressive neuromuscular disorder, which can eventually lead to early childhood death. In aged people, induced dysfunction of TIM is related to the neurodegenerative Alzheimer's disease [2, 3]. Diseased cells of patients with TIM deficiency usually exhibit reduced TIM activity and a high level of the substrate DHAP. There are also many misfolded TIM proteins in diseased cells which accumulate to form large protein aggregates directly responsible for the neurodegenerative disorder . A more recent study shows that reduced activity of TIM can lead to an oxidized redox state, making the subject sensitive to oxidative stress, and this explains why the disease phenotype progresses gradually .
More than 10 mutations  have been observed in Homo sapiens TIM (HsTIM) that causes the TIM deficiency, such as E104D , V231M , F240L , C41Y , and I170V . The mutation E104D, a glutamate at position 104 substituted by an aspartate, occurs most often of these mutations. E104D was first studied by Daar et al. in 1986  by comparing two unrelated TIM patients. An important finding is that the E104D mutant is thermolabile--that is, when the temperature goes up slightly, the enzyme is easily subject to destruction or great change, and then the enzyme loses its activity quickly. Following these pioneering studies [1, 7], many other patients with the E104D mutation from several families across the world have been investigated [8, 11–13]. All these patients are believed to originate from a common ancestor dated back to more than 1000 years ago . In 1994, a recombinant HsTIM structure was solved at 2.8 Å resolution , and the local residue organization of GLU104 was unveiled. However, the conclusion--the E104D mutation changes the structure of the active sites through a chain of perturbations--is not exactly correct. Rather, in 2006, Ralser et al. concluded that TIM deficiency is caused by the altered dimerization but not the inactivity of the enzyme . In 2008, the structure of E104D HsTIM mutant was solved at a better resolution 1.85 Å . Comparing with the wild type structure, it was found that there is no significant change in the active site region, and found that the catalytic activity of the mutant enzyme is at the same magnitude as the native enzyme. However, the structural comparison of the mutated site between the wild type and the mutant reveals that there is a perturbation in the organization of the interfacial water molecules. Despite of these conflicting statements in the past research, the common understanding on the pathogenesis of the mutation E104D shares the following points: (i) There is an alteration of the binding of the two subunits in the E104D mutant ; (ii) The altered binding of the two subunits harms thermostability of the protein ; and (iii) the excessive thermolability causes the dysfunction of the enzyme . These are wet-lab results and are hence valid. The only gap yet to be filled in the chain of knowledge explaining the pathogenesis of the mutation E104D is to find out how exactly the E104D alters the binding of the two subunits. It is already known that the interfacial water molecules play an important role in the binding of the two subunits . However, the literature work does not answer why and how the interfacial water molecules contribute to the excessive thermolability of the mutant. In fact, it is unknown whether or not the presence of the large water cluster is responsible for the thermolability.
In this work, we conduct a comparative analysis on the wild type HsTIM and the E104D mutant's structures to understand the pathogenesis of this mutation and address three specific questions:
why the mutation E104D can alter the binding affinity of the two subunits,
how the mutation E104D introduces excessive thermolability to the enzyme, and
what role the interfacial water molecules play.
This work goes beyond the scope of structural perturbation studies. We hypothesize that the large cluster of water in region A transmits the atoms' mobility from the surface deeply into the core of the interface, and that this mobility transmission can harm the thermostability of the subunit binding and the function of the protein. As the B factor of atoms is closely related to the thermolability of a protein [18, 19], the elevated normalized B factors are used to understand the increased flexibility of the interfacial water and other atoms clustered near the mutated site (i.e., water in region A).
To validate our hypothesis of the role of water cluster in region A, molecular dynamics simulations were run for the wild type and the mutant HsTIM. The difference in the root-mean-square deviation (RMSD) between the two structures over the simulation is only observed in region A. For other interfacial regions and the whole protein, the RMSD does not change much. Thus we conclude that it is the water molecules in region A that are directly responsible for the excessive thermolability of the enzyme.
As highlights, our pathogenesis study for the mutation E104D is to affirm: (i) this mutation can lead to an excessive thermolability of the enzyme, and (ii) the binding between the two subunits is harmed by the increased flexibility of the interfacial atoms, which are amplified by the large cluster of buried water near the mutation site. In other words, the abnormal organizational topology of water molecules in the HsTIM dimer interface is the key to understanding the fundamental mechanism of the HsTIM deficiency. Our results can be generalized to understand that large clusters of water molecules in protein binding interface are very high-maintenance, and they can greatly affect the overall structure, function and evolution of the protein.
Protein structures used in this study
TIM structures and subunit interfaces used in this study.
The water distribution topology of HsTIM dimer interface is compared with a data set of homodimeric interfaces. They are selected from 206 obligate interfaces collected from a few previous studies [21–24]. Obligate protein interfaces are considered as a homodimeric interface if the two interacting protein chains are identical. Some interfaces with less than 10 water molecules are removed as their rWBL and gini coefficient (defined later) are of large variance. Finally, 91 homodimers are left and they are used in a comparative analysis to illustrate the abnormality of HsTIM dimer interface. The information of these 91 interfaces and their properties are available in Additional file 1.
Definitions of protein-water-protein tripartite interface, burial level and gini coefficient
We consider only heavy atoms in this study. Water molecules immobilized in a protein/protein complex are considered as a part of the protein/protein complex rather than a part of bulk solvent. Water molecules with a larger than 10 Å2 solvent accessible surface area (SASA) are defined as exposed as a part of bulk solvent and they are removed iteratively from a structure until no water molecules are exposed in the structure. The remaining structure including those buried water molecules are considered as the protein complex and are used in the subsequent processes in our work.
For each protein complex structure, an atomic contact graph is defined with atoms as its nodes and atomic contacts as its edges. Atomic contacts are defined by combining the Voronoi diagram and an atom-atom distance. The distance threshold for an atomic contact between two atoms is the sum of their radii plus the diameter of a water molecule 2.75 Å.
Protein-water-protein tripartite interface
We define the interface between two proteins as a tripartite graph, which is a subgraph of the atomic contact graph of the whole protein complex. The nodes of an interface tripartite graph consist of three sets of atoms: the oxygen atoms of the interfacial water molecules and the atoms from the two interacting partners that have atomic contact with interfacial water or with the interacting partner. A water molecule is defined as interfacial water if it has atomic contacts with both sides. The edges in the tripartite graph are the atomic contacts among the three sets. In this way the atom level interface is defined as a set of atoms. We also define any residue that has at least one atom in the interface as an interfacial residue.
Burial level and relative water burial level
The burial level of an atom a (denoted as BL a ) in a protein complex is calculated in the atomic contact graph. It is the length of the shortest path to the nearest exposed atom-an atom with an SASA larger than 10 Å2.
where I is the set of atoms in the interface, I W is the oxygen atoms of the water molecules in interface, and |I| is the cardinality of set I, which is the number of atoms in this interface. According to this definition, a large rWBL indicates that the water molecules are deeply buried.
Gini coefficient of water distribution
Gini coefficient is widely used to describe the unevenness of a distribution. Usually it is used to measure the inequality of income or wealth. Its value is between 0 and 1. A lower gini coefficient indicates a more uneven income distribution of a country. Gini coefficient is defined based on the Lorenz curve . Lorenz curve is a plot of the cumulative share of population ordered from the lowest income to the highest income versus the cumulative share of the income. If the incomes are absolutely equal, the Lorenz curve will be a straight line, which is called line of equality. The gini coefficient is defined as the area between the Lorenz curve and the line of equality divided by the area under the line of equality.
Structural alignment, water correspondence and B factor normalization
As the sequences of the wild type HsTIM and E104D mutant differ only from each other in one position, we perform structural alignment of the two structures by superimposing them using Pymol . The algorithm is based on the amino acid sequences of the two structures. Once the aligned wild type and mutant structures are obtained, the 25 water molecules in wild type are searched in the mutant structure to determine whether it reappears or not. If two interfacial water molecules, one in wild type and the other in mutant, are within 1.0 Å distance from each other, and they are mutually the closest water molecule to each other, we say the water in the wild type structure reappears in the mutant.
where B is the original B factor reported in the PDB record, and and σ a are the mean and standard deviation of the B factors of non water atoms which are at least 15 Å away from the C γ of the mutated residue--GLU104. We do not consider water here because the water information quality is correlated with the resolution  and, maybe, other environmental issues. We also exclude the atoms near the mutated site from calculating mean and standard deviation as the mutation is expected to change the B factors of atoms near it. In this way we compare the B factor of an atom in wild type and mutant with non water atoms far from the mutation site as references.
Molecular dynamics simulation
Two molecular dynamics simulations were run for the wild type [PDB: 2JK2] and the E104D mutant [PDB: 2VOM] HsTIMs. Buried water molecules were included in the initial structures. For both simulations, the protein structure was solvated in a water box where a minimum of 10 Å distance was kept between the protein and the boundary. Charges were neutralized by either Cl- or Na+ ions. Solvation and ionization were performed using the the VMD software .
For both simulations, CHARMM22 parameter set  with CMAP correction  was used, and the step size was set to 2 fs. The solvated and ionized system was minimized for 1000 steps and simulated for 2500000 steps (5 ns). Initial velocities were set with human body temperature (310 K). Langevin piston pressure control was used to control the system pressure at 1 atm. Periodic boundary conditions were used, and a threshold cutoff of 12 Å was set for non-bonded interactions. Particle mesh Ewald method  was used to calculate long-range electrostatic interactions. The simulations was carried out with the NAMD software .
We first present results to illustrate how exactly atypical the HsTIM dimer interface is, followed by results on the evolution of the water and residues in the interface of TIMs. We then present a hypothesis to understand the pathogenesis of the mutation E104D, and validate the hypothesis by molecular dynamics simulations.
The abnormal hydration in HsTIM dimer interfaces
Figure 1(A) shows one side of the interface and all of the 25 interfacial water molecules (in blue). We divide the interface residues and water molecules into four regions for detailed analysis: (i) region A which is a wet region consisting of 16 water molecules and their surrounding 11 residues (in magenta and black color) from each side including the pathogenic mutation site 104 (in black color), (ii) region B which is another wet region with 5 water molecules and 6 residues (in orange color) from each side, (iii) region C (in red color) which is the functional region of the enzyme, containing 15 residues from each side including three of the four active site residues (ASN11, LYS13 and HIS95), and (iv) region D (in yellow color) which is a region in the interface containing residues and water that are not in region A, B or C. As the interface is homodimeric and symmetric, the residues in these 4 regions can be doubled actually.
Specifically in region A, the residue GLU104, whose mutation into aspartate causes TIM deficiency, is highlighted in black as shown in Figures 1(A) and 1(B). It can be seen that this residue is located at the waterside and reaches its side chain deeply into the water cluster in region A to "lock" this large amount of water molecules inside.
The total 25 water molecules in this interface account for 4.6% of the entire interfacial atoms. This wetness is similar to the average wetness of obligate interfaces . However, it can be seen that these water molecules are distributed in an unusual way with a very special topology. First, most of these water molecules are deeply buried and clustered near the core of the interface. Second, these water molecules are unevenly distributed, mainly in two wet regions. Normally, water molecules in protein binding interfaces are organized in a "dry-core-wet-rim" manner with most water molecules placed near the rim and the wetness goes down progressively from rim to core . This interface does not follow such a trend at all. In particular, the average burial level of the water molecules is 2.32, yielding an extraordinarily high rWBL (the relative burial level of water in an interface, see Materials and Methods) of 1.282 with regard to the average burial level of all the interfacial atoms 1.81. This kind of large water cluster is also unusual in terms of binding free energy, as immobilizing water molecules at the core of a protein interface is energetically expensive .
Evolutionary studies on HsTIM dimer interface hydration
There are 21 species' wild type triosephosphate isomerase structures that have been solved by X-ray crystallography and the data is deposited in the protein data bank  as of 2011. One of them (PDB: 1BTM) does not have solvent information, so it is excluded from our analysis. The 20 TIM structures are listed in Table 1.
As shown in Figure 3(A), region A is wet in all of the eukaryotic and bacterial TIMs, indicating a large cluster of water is present in all non-archaeal TIMs. In Figure 4, it seems that three bacterial TIM dimer interfaces (in E. coli, M. marina and T. maritima TIMs) are dry in region A. However, these three structures are solved at very poor resolutions (Table 1), which usually implies a under-reporting of water information . Meanwhile, as in the other non-archaeal TIMs, large inter-protein cavities are observed in these structures, which may accommodate water molecules. Thus, we believe these three interfaces are as wet as the other non-bacterial TIM dimer interfaces. Region A of the archaeal TIMs does not have as much water as that of the non-archaeal organisms, as can be seen from Figure 3(A) and Figure 4. Furthermore, three positions (positions 102, 104 and 112 according to HsTIM numbering) in the eukaryotic and bacterial TIM dimer interfaces are not in the corresponding interface in the archaeal TIMs. Positions 104 (where the E104D mutation occurs) and 112 are fully conserved in eukaryotic and bacterial organisms, always with position 104 as glutamate and position 112 as lysine. In the archaeal TIMs, position 104 changes to aspartate. Recall that, the E104D mutation in HsTIM, where the glutamate at position 104 is substituted by aspartate, is pathogenic. Also, as can be seen from Figure 4, in the archaeal TIMs (the last three panels), positions 102 (colored in magenta, left bottom) and 104 (in black) are separated from the other residues in region A (large area colored in magenta). Position 112 does not even exist in the archaeal TIMs.
These observations affirm that the large cluster of water in Region A is essential for the eukaryotic and bacterial TIMs. It contributes positively to the binding affinity, because if the interface loses the large cluster of water (as in the archaeal TIMs), the enzyme must adopt another quaternary structure to stabilize the enzyme. In evolution, the tertiary structures of non-archaeal TIMs are not ready to oligomerize into tetramers or other higher order oligomers. So, excluding such a large cluster of water from the interface without harming the function of the enzyme seems too hard--once a mutation disturbs the water topology in this region (such as E104D), it would significantly destabilize the enzyme and further turn off its functionality. Thus safely excluding these water molecules requires simultaneous mutation of many residues, which is too difficult in natural evolution. That is the reason why the large water cluster is necessary and is well maintained when TIMs are in the form of dimers. Also, to maintain such a large pocket of water, its nearby residues must be conserved.
Region B seems to be wet in almost all of the 20 TIMs. Region C's water contacting profile is not very conserved. In human it is dry, but in some other species, such as C. elegans and H. Pylori, it is wet. The hydration profile of region D is not very conserved, either. Combining these hydration profiles and conservation scores of the four regions, we can conclude that
region A: both the water and residues are conserved in the bacterial and eukaryotic organisms;
region B: only water is conserved;
region C: only residues are conserved; and
region D: neither water nor residue is conserved.
Comparison between wild type and E104D mutant: why this mutation is pathogenic
Comparison between HsTIM wild type and E104D mutant dimer interfaces.
The overall structural perturbation by the mutation is not significant
Propagating atoms' mobility by water in region A: a hypothesis to explain the excessive thermolability of the mutant
As shown in the previous section, the overall structural change in the interfacial atoms and water molecules are not significant, and the overall structure of the interface is maintained in the E104D mutant. Thus the structural perturbation theories [7, 15] cannot explain why the mutation is pathogenic. Recall that the mutation only makes the enzyme thermolabile, and the mutant is still functional at normal temperature. We thus hypothesize that the extra-large cluster of water molecules in region A propagates the flexibility of interfacial atoms deeply down to the core of the interface, adding excessive thermolability to the binding of the two subunits. As monomers of TIM is not functional, the thermolability of the TIM function is thus observed.
To get a better understanding of this hypothesis, we investigated the B factor of interfacial atoms in the two structures. The B factor of an atom indicates the level of thermal motion of the atom around its average position. It is closely related to the thermolability of a protein. It was found by Reetz's group that the thermostability/thermolability of a protein can be manipulated by using B factor as a guide [18, 19]. Thus, B factor can also be used to study the thermostability/thermolability of the wild type and mutated HsTIM.
In order to better understand the role of the large water cluster in region A in enhancing the interface mobility after the mutation, the average change in normalized B factor of interfacial non water atoms is plotted at different burial levels (the burial levels in the wild type enzyme); see Figure 7(B). First of all, the mobility of an interface atom tends to increase no matter where the atom is located or how deep it is buried, as all the average changes are larger than 0. Second, the mobility of atoms in region A increases the most comparing with the mobility increment of atoms in the other three regions. More importantly, the increase in mobility propagates itself deeply into the core of region A, while in other three regions the mobility increment drops quickly when burial level goes up. When the region has more water, the mobility increment decrease slower when the burial level goes higher, and that is why region B also have a higher mobility increment than the two dry regions (regions C and D), though it is much drier than region A.
The mobility change of the interfacial water molecules and atoms indicates that after the mutation, water and atoms in region A are more flexible, and the flexibility increment flows to those deeply buried atoms in this region through the large water cluster. This change is far weaker in other dry regions.
Molecular dynamics simulations: validating the role of interfacial water
Regional RMSD of the subunit binding interfaces is compared in Figures 8(B) and 8(C). For interfacial regions B, C and D, the RMSD in the mutant is also similar to that in the wild type structure. Interestingly, a remarkable difference of RMSD is observed in region A, where the RMSD in the mutant is in general higher than that in the wild type, indicating that excessive flexibility is gained after the mutation. Meanwhile, at normal temperature, the excessive flexibility in region A does not reach other regions that is not in contact with the large cluster of water. This observation clearly shows that it is the large cluster of water that is responsible for the thermolability of the mutant as region A in the subunit interface is the only region that is more flexible after the mutation.
From the abovementioned analyses, we can conclude that the large cluster of water in region A is responsible for the excessive thermolability of the mutant. Although this work only analyzed the mutation E104D, the concept can be generalized to understand the high-maintenance of large water clusters in proteins. In HsTIM, we believe that this large cluster of water may also play important role in the pathogenesis of other point mutations in this enzyme and this water cluster could also be one of the reasons why there are so many pathogenic point mutations in this enzyme .
An investigation of a mutation that causes HsTIM deficiency, E104D, is conducted. The HsTIM dimer interface is abnormally hydrated with a very strange water distribution pattern--wet-core-dry-rim. The water molecules are mainly clustered compactly in a region with the residue GLU104 aside. This residue, along with several other residues in this region, is highly conserved when this region is wet. Comparing the wild type and the E104D mutant structures, no significant structural change was observed. The overall structure of the protein, including most of the large cluster of water, is maintained after the mutation. We hypothesize that, in the mutant, the water molecules in the subunit interface introduce the excessive thermolability to the protein by propagating atoms' flexibility into the core of the interface. This hypothesize was supported by the fact that atoms near the large cluster of water have a larger B factor increment than those in other interfacial regions. The hypothesize was further validated by using molecular dynamics simulations. We showed that the interfacial region near the large cluster of water was the only region that had a higher RMSD in the mutant than in the wild type.
This work was supported in part by two Singapore MOE Tier-2 grants (T208B2203 and MOE2009-T2-2-004), and supported in part by an ARC Discovery Project (DP130102124).
Publication of this article was funded by a UTS 2013 Early Career Research Grant.
This article has been published as part of BMC Bioinformatics Volume 14 Supplement 16, 2013: Twelfth International Conference on Bioinformatics (InCoB2013): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/14/S16.
- Schneider AS, Valentine WN, Hattori M, Heins HL: Hereditary Hemolytic Anemia with Triosephosphate Isomerase Deficiency. N Engl J Med. 1965, 272 (5): 229-235. 10.1056/NEJM196502042720503.View ArticlePubMedGoogle Scholar
- Guix FX, Ill-Raga G, Bravo R, Nakaya T, de Fabritiis G, Coma M, Miscione GP, Villa-Freixa J, Suzuki T, Fernandez-Busquets X, Valverde MA, de Strooper B, Munoz FJ: Amyloid-dependent triosephosphate isomerase nitrotyrosination induces glycation and tau fibrillation. Brain. 2009, 132: 1335-1345. 10.1093/brain/awp023.View ArticlePubMedGoogle Scholar
- Tajes M, Guivernau B, Ramos-Fernández E, Bosch-Morató M, Palomer E, Guix FX, Muñoz FJ: The patho-physiology of triose phosphate isomerase dysfunction in Alzheimer's disease. Histol Histopathol. 2013, 28: 43-51.PubMedGoogle Scholar
- Orosz F, Oláh J, Ovádi J: Triosephosphate isomerase deficiency: Facts and doubts. IUBMB Life. 2006, 58 (12): 703-715. 10.1080/15216540601115960.View ArticlePubMedGoogle Scholar
- Hrizo SL, Fisher IJ, Long DR, Hutton JA, Liu Z, Palladino MJ: Early mitochondrial dysfunction leads to altered redox chemistry underlying pathogenesis of TPI deficiency. Neurobiology of Disease. 2013,Google Scholar
- Orosz F, Oláh J, Ovádi J: Triosephosphate isomerase deficiency: new insights into an enigmatic disease. Biochim. Biophys. Acta. 2009, 1792: 1168-1174. 10.1016/j.bbadis.2009.09.012.View ArticlePubMedGoogle Scholar
- Daar IO, Artymiuk PJ, Phillips DC, Maquat LE: Human triose-phosphate isomerase deficiency: a single amino acid substitution results in a thermolabile enzyme. Proc Natl Acad Sci USA. 1986, 83 (20): 7903-7907. 10.1073/pnas.83.20.7903.PubMed CentralView ArticlePubMedGoogle Scholar
- Neubauer BA, Pekrun A, Beauregard SW, Lakomek M, Schroeter W: Relation between genetic defect, altered protein structure, and enzyme function in triose-phosphate isomerase (TPI) deficiency. Eur J Pediatr. 1992, 232-232. Suppl 151Google Scholar
- Hollan S, Fujii H, Hirono A, Hirono K, Karro H, Miwa S, Harsanyi V, Gyodi E, Inselt-Kovacs M: Hereditary triosephosphate isomerase (TPI) deficiency: two severely affected brothers one with and one without neurological symptoms. Hum Genet. 1993, 92 (5): 486-490. 10.1007/BF00216456.View ArticlePubMedGoogle Scholar
- Arya R, Lalloz MRA, Bellingham AJ, Layton DM: Evidence for founder effect of the glu104asp substitution and identification of new mutations in triosephosphate isomerase deficiency. Human Mutation. 1997, 10 (4): 290-294. 10.1002/(SICI)1098-1004(1997)10:4<290::AID-HUMU4>3.0.CO;2-L.View ArticlePubMedGoogle Scholar
- Pekrun A, Neubauer BA, Eber SW, Lakomek M, Seidel H, Schroter W: Triosephosphate isomerase deficiency: biochemical and molecular genetic analysis for prenatal diagnosis. Clin Genet. 1995, 47 (4): 175-179.View ArticlePubMedGoogle Scholar
- Schneider A, Westwood B, Yim C, Prchal J, Berkow R, Labotka R, Warrier R, Beutler E: Triosephosphate isomerase deficiency: repetitive occurrence of point mutation in amino acid 104 in multiple apparently unrelated families. Am J Hematol. 1995, 50 (4): 263-268. 10.1002/ajh.2830500407.View ArticlePubMedGoogle Scholar
- Linarello RE, Shetty AK, Thomas T, Warrier RP: Triosephosphate isomerase deficiency in a child with congenital hemolytic anemia and severe hypotonia. Pediatr Hematol Oncol. 1998, 15 (6): 553-556. 10.3109/08880019809018318.View ArticlePubMedGoogle Scholar
- Schneider A, Westwood B, Yim C, Cohen-Solal M, Rosa R, Labotka R, Eber S, Wolf R, Lammi A, Beutler E: The 1591C mutation in triosephosphate isomerase (TPI) deficiency. Tightly linked polymorphisms and a common haplotype in all known families. Blood Cells Mol Dis. 1996, 22 (12): 115-125.View ArticlePubMedGoogle Scholar
- Mande SC, Mainfroid V, Kalk KH, Goraj K, Martial JA, Hol WG: Crystal structure of recombinant human triosephosphate isomerase at 2.8 Å resolution. Triosephosphate isomerase-related human genetic disorders and comparison with the trypanosomal enzyme. Protein Sci. 1994, 3 (5): 810-821.PubMed CentralView ArticlePubMedGoogle Scholar
- Ralser M, Heeren G, Breitenbach M, Lehrach H, Krobitsch S: Triose Phosphate Isomerase Deficiency Is Caused by Altered Dimerization--not Catalytic Inactivity--of the Mutant Enzymes. PLoS ONE. 2006, 1: e30-10.1371/journal.pone.0000030.PubMed CentralView ArticlePubMedGoogle Scholar
- Rodriguez-Almazan C, Arreola R, Rodriguez-Larrea D, Aguirre-Lopez B, de Gomez-Puyou MT, Perez-Montfort R, Costas M, Gomez-Puyou A, Torres-Larios A: Structural Basis of Human Triosephosphate Isomerase Deficiency: MUTATION E104D IS RELATED TO ALTERATIONS OF A CONSERVED WATER NETWORK AT THE DIMER INTERFACE. J Biol Chem. 2008, 283 (34): 23254-23263. 10.1074/jbc.M802145200.View ArticlePubMedGoogle Scholar
- Reetz MT, Carballeira JD, Vogel A: Iterative Saturation Mutagenesis on the Basis of B Factors as a Strategy for Increasing Protein Thermostability. Angew Chem Int Ed. 2006, 45 (46): 7745-7751. 10.1002/anie.200602795.View ArticleGoogle Scholar
- Reetz MT, Soni P, Fernández L: Knowledge-guided laboratory evolution of protein thermolability. Biotechnol Bioeng. 2009, 102 (6): 1712-1717. 10.1002/bit.22202.View ArticlePubMedGoogle Scholar
- Wierenga R, Kapetaniou E, Venkatesan R: Triosephosphate isomerase: a highly evolved biocatalyst. Cell Mol Life Sci. 2010, 67: 3961-3982. 10.1007/s00018-010-0473-9.View ArticlePubMedGoogle Scholar
- Ponstingl H, Henrick K, Thornton JM: Discriminating between homodimeric and monomeric proteins in the crystalline state. Proteins Struct Funct Bioinf. 2000, 41: 47-57. 10.1002/1097-0134(20001001)41:1<47::AID-PROT80>3.0.CO;2-8.View ArticleGoogle Scholar
- Bahadur RP, Chakrabarti P, Rodier F, Janin J: A dissection of specific and non-specific protein-protein interfaces. J Mol Biol. 2004, 336 (4): 943-955. 10.1016/j.jmb.2003.12.073.View ArticlePubMedGoogle Scholar
- Mintseris J, Weng Z: Structure, function, and evolution of transient and obligate protein-protein interactions. Proc Natl Acad Sci USA. 2005, 102 (31): 10930-10935. 10.1073/pnas.0502667102.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhu H, Domingues F, Sommer I, Lengauer T: NOXclass: prediction of protein-protein interaction types. BMC Bioinf. 2006, 7: 27-10.1186/1471-2105-7-27.View ArticleGoogle Scholar
- Lorenz MO: Methods of Measuring the Concentration of Wealth. Publications of the American Statistical Association. 1905, 9 (70): 209-219. 10.2307/2276207.View ArticleGoogle Scholar
- Schrödinger LLC: The PyMOL Molecular Graphics System, Version 1.3r1. 2010Google Scholar
- Carugo O, Bordo D: How many water molecules can be detected by protein crystallography?. Acta Crystallogr D Biol Crystallogr. 1999, 55 (Pt 2): 479-483.View ArticlePubMedGoogle Scholar
- Humphrey W, Dalke A, Schulten K: VMD: Visual molecular dynamics. Journal of Molecular Graphics. 1996, 14: 33-38. 10.1016/0263-7855(96)00018-5.View ArticlePubMedGoogle Scholar
- MacKerell AD, Bashford D, Dunbrack RL, Evanseck JD, Field MJ, Fischer S, Gao J, Guo H, Ha S, Joseph-McCarthy D, Kuchnir L, Kuczera K, Lau FTK, Mattos C, Michnick S, Ngo T, Nguyen DT, Prodhom B, Reiher WE, Roux B, Schlenkrich M, Smith JC, Stote R, Straub J, Watanabe M, Wiórkiewicz-Kuczera J, Yin D, Karplus M: All-Atom Empirical Potential for Molecular Modeling and Dynamics Studies of Proteins. The Journal of Physical Chemistry B. 1998, 102 (18): 3586-3616. 10.1021/jp973084f.View ArticlePubMedGoogle Scholar
- Mackerell AD: Empirical force fields for biological macromolecules: Overview and issues. Journal of Computational Chemistry. 2004, 25 (13): 1584-1604. 10.1002/jcc.20082.View ArticlePubMedGoogle Scholar
- Darden T, York D, Pedersen L: Particle mesh Ewald: An N.log(N) method for Ewald sums in large systems. The Journal of Chemical Physics. 1993, 98 (12): 10089-10092. 10.1063/1.464397.View ArticleGoogle Scholar
- Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, Villa E, Chipot C, Skeel RD, Kalé L, Schulten K: Scalable molecular dynamics with NAMD. J Comput Chem. 2005, 26 (16): 1781-1802. 10.1002/jcc.20289.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Z, He Y, Wong L, Li J: Progressive dry-core-wet-rim hydration trend in a nested-ring topology of protein binding interfaces. BMC Bioinf. 2012, 13: 51+-10.1186/1471-2105-13-51.View ArticleGoogle Scholar
- Fischer S, Verma CS: Binding of buried structural water increases the flexibility of proteins. Proc Natl Acad Sci USA. 1999, 96 (17): 9613-9615. 10.1073/pnas.96.17.9613.PubMed CentralView ArticlePubMedGoogle Scholar
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235-242. 10.1093/nar/28.1.235.PubMed CentralView ArticlePubMedGoogle Scholar
- Keeling PJ, Doolittle WF: Evidence that eukaryotic triosephosphate isomerase is of alpha-proteobacterial origin. Proc Natl Acad Sci USA. 1997, 94 (4): 1270-1275. 10.1073/pnas.94.4.1270.PubMed CentralView ArticlePubMedGoogle Scholar
- Ashkenazy H, Erez E, Martz E, Pupko T, Ben-Tal N: ConSurf 2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res. 2010, 38 (suppl 2): W529-W533.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones S, Thornton JM: Principles of protein-protein interactions. Proc Natl Acad Sci USA. 1996, 93: 13-20. 10.1073/pnas.93.1.13.PubMed CentralView ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.