The ionization of a buried glutamic acid is thermodynamically linked to the stability of Leishmania mexicana triose phosphate isomerase

Experimental structure based drug design (SBDD) applications for anti-leishmanial drugs: A paradigm shift?

Leishmaniasis is a group of neglected tropical diseases caused by at least 20 species of Leishmania protozoa, which are spread by the bite of infected sandflies. There are three main forms of the disease: cutaneous leishmaniasis (CL, the most common), visceral leishmaniasis (VL, also known as kala-azar, the most serious), and mucocutaneous leishmaniasis. One billion people live in areas endemic to leishmaniasis, with an annual estimation of 30,000 new cases of VL and more than 1 million of CL. New treatments for leishmaniasis are an urgent need, as the existing ones are inefficient, toxic, and/or expensive. We have revised the experimental structure-based drug design (SBDD) efforts applied to the discovery of new drugs against leishmaniasis. We have grouped the explored targets according to the metabolic pathways they belong to, and the key achieved advances are highlighted and evaluated. In most cases, SBDD studies follow high-throughput screening campaigns and are secondary to pharmacokinetic optimization, due to the majoritarian belief that there are few validated targets for SBDD in leishmaniasis. However, some SBDD strategies have significantly contributed to new drug candidates against leishmaniasis and a bigger number holds promise for future development.

The Role of Asn11 in Catalysis by Triosephosphate Isomerase

Four catalytic amino acids at triosephosphate isomerase (TIM) are highly conserved: N11, K13, H95, and E167. Asparagine 11 is the last of these to be characterized in mutagenesis studies. The ND2 side chain atom of N11 is hydrogen bonded to the O-1 hydroxyl of enzyme-bound dihydroxyacetone phosphate (DHAP), and it sits in an extended chain of hydrogen-bonded side chains that includes T75' from the second subunit. The N11A variants of wild-type TIM from Trypanosoma brucei brucei (TbbTIM) and Leishmania mexicana (LmTIM) undergo dissociation from the dimer to monomer under our assay conditions. Values of Kas = 8 × 103 and 1 × 106 M-1, respectively, were determined for the conversion of monomeric N11A TbbTIM and LmTIM into their homodimers. The N11A substitution at the variant of LmTIM previously stabilized by the E65Q substitution gives the N11A/E65Q variant that is stable to dissociation under our assay conditions. The X-ray crystal structure of N11A/E65Q LmTIM shows an active site that is essentially superimposable on that for wild-type TbbTIM, which also has a glutamine at position 65. A comparison of the kinetic parameters for E65Q LmTIM and N11A/E65Q LmTIM-catalyzed reactions of (R)-glyceraldehyde 3-phosphate (GAP) and (DHAP) shows that the N11A substitution results in a (13-14)-fold decrease in kcat/Km for substrate isomerization and a similar decrease in kcat for DHAP but only a 2-fold decrease in kcat for GAP.

Analysis of Biologics Molecular Descriptors towards Predictive Modelling for Protein Drug Development Using Time-Gated Raman Spectroscopy

Pharmaceutical proteins, compared to small molecular weight drugs, are relatively fragile molecules, thus necessitating monitoring protein unfolding and aggregation during production and post-marketing. Currently, many analytical techniques take offline measurements, which cannot directly assess protein folding during production and unfolding during processing and storage. In addition, several orthogonal techniques are needed during production and market surveillance. In this study, we introduce the use of time-gated Raman spectroscopy to identify molecular descriptors of protein unfolding. Raman spectroscopy can measure the unfolding of proteins in-line and in real-time without labels. Using K-means clustering and PCA analysis, we could correlate local unfolding events with traditional analytical methods. This is the first step toward predictive modeling of unfolding events of proteins during production and storage.

Parallel molecular mechanisms for enzyme temperature adaptation

The mechanisms that underly the adaptation of enzyme activities and stabilities to temperature are fundamental to our understanding of molecular evolution and how enzymes work. Here, we investigate the molecular and evolutionary mechanisms of enzyme temperature adaption, combining deep mechanistic studies with comprehensive sequence analyses of thousands of enzymes. We show that temperature adaptation in ketosteroid isomerase (KSI) arises primarily from one residue change with limited, local epistasis, and we establish the underlying physical mechanisms. This residue change occurs in diverse KSI backgrounds, suggesting parallel adaptation to temperature. We identify residues associated with organismal growth temperature across 1005 diverse bacterial enzyme families, suggesting widespread parallel adaptation to temperature. We assess the residue properties, molecular interactions, and interaction networks that appear to underly temperature adaptation.

Perdeuteration, large crystal growth and neutron data collection of Leishmania mexicana triose-phosphate isomerase E65Q variant

Triose-phosphate isomerase (TIM) catalyses the interconversion of dihydroxyacetone phosphate and glyceraldehyde 3-phosphate. Two catalytic mechanisms have been proposed based on two reaction-intermediate analogues, 2-phosphoglycolate (2PG) and phosphoglycolohydroxamate (PGH), that have been used as mimics of the cis-enediol(ate) intermediate in several studies of TIM. The protonation states that are critical for the mechanistic interpretation of these structures are generally not visible in the X-ray structures. To resolve these questions, it is necessary to determine the hydrogen positions using neutron crystallography. Neutron crystallography requires large crystals and benefits from replacing all hydrogens with deuterium. Leishmania mexicana triose-phosphate isomerase was therefore perdeuterated and large crystals with 2PG and PGH were produced. Neutron diffraction data collected from two crystals with different volumes highlighted the importance of crystal volume, as smaller crystals required longer exposures and resulted in overall worse statistics.

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.

Substitutions of a buried glutamate residue hinder the conformational change in horse liver alcohol dehydrogenase and yield a surprising complex with endogenous 3'-Dephosphocoenzyme A

Glu-267 is highly conserved in alcohol dehydrogenases and buried as a negatively-charged residue in a loop of the NAD coenzyme binding domain. Glu-267 might have a structural role and contribute to a rate-promoting vibration that facilitates catalysis. Substitutions of Glu-267 with histidine or asparagine residues increase the dissociation constants for the coenzymes (NAD⁺ by ∼40-fold, NADH by ∼200-fold) and significantly decrease catalytic efficiencies by 16-1200-fold various substrates and substituted enzymes. The turnover numbers modestly change with the substitutions, but hydride transfer is at least partially rate-limiting for turnover for alcohol oxidation. X-ray structures of the E267H and E267 N enzymes are similar to the apoenzyme (open) conformation of the wild-type enzyme, and the substitutions are accommodated by local changes in the structure. Surprisingly, the E267H and E267 N enzymes have endogenous (from the expression in E. coli) 3'-dephosphocoenzyme A bound in the active site with the ADP moiety in the NAD binding site and the pantethiene sulfhydryl bound to the catalytic zinc. The kinetics and crystallography show that the substitutions of Glu-267 hinder the conformational change, which occurs when wild-type enzyme binds coenzymes, and affect productive binding of substrates.

A guide to the effects of a large portion of the residues of triosephosphate isomerase on catalysis, stability, druggability, and human disease

Triosephosphate isomerase (TIM) is a ubiquitous enzyme, which appeared early in evolution. TIM is responsible for obtaining net ATP from glycolysis and producing an extra pyruvate molecule for each glucose molecule, under aerobic and anaerobic conditions. It is placed in a metabolic crossroad that allows a quick balance of the triose phosphate aldolase produced by glycolysis, and is also linked to lipid metabolism through the alternation of glycerol-3-phosphate and the pentose cycle. TIM is one of the most studied enzymes with more than 199 structures deposited in the PDB. The interest for this enzyme stems from the fact that it is involved in glycolysis, but also in aging, human diseases and metabolism. TIM has been a target in the search for chemical compounds against infectious diseases and is a model to study catalytic features. Until February 2017, 62% of all residues of the protein have been studied by mutagenesis and/or using other approaches. Here, we present a detailed and comprehensive recompilation of the reported effects on TIM catalysis, stability, druggability and human disease produced by each of the amino acids studied, contributing to a better understanding of the properties of this fundamental protein. The information reviewed here shows that the role of the noncatalytic residues depend on their molecular context, the delicate balance between the short and long-range interactions in concerted action determining the properties of the protein. Each protein should be regarded as a unique entity that has evolved to be functional in the organism to which it belongs. Proteins 2017; 85:1190-1211. © 2017 Wiley Periodicals, Inc.

Structural and Genetic Studies Demonstrate Neurologic Dysfunction in Triosephosphate Isomerase Deficiency Is Associated with Impaired Synaptic Vesicle Dynamics

Triosephosphate isomerase (TPI) deficiency is a poorly understood disease characterized by hemolytic anemia, cardiomyopathy, neurologic dysfunction, and early death. TPI deficiency is one of a group of diseases known as glycolytic enzymopathies, but is unique for its severe patient neuropathology and early mortality. The disease is caused by missense mutations and dysfunction in the glycolytic enzyme, TPI. Previous studies have detailed structural and catalytic changes elicited by disease-associated TPI substitutions, and samples of patient erythrocytes have yielded insight into patient hemolytic anemia; however, the neuropathophysiology of this disease remains a mystery. This study combines structural, biochemical, and genetic approaches to demonstrate that perturbations of the TPI dimer interface are sufficient to elicit TPI deficiency neuropathogenesis. The present study demonstrates that neurologic dysfunction resulting from TPI deficiency is characterized by synaptic vesicle dysfunction, and can be attenuated with catalytically inactive TPI. Collectively, our findings are the first to identify, to our knowledge, a functional synaptic defect in TPI deficiency derived from molecular changes in the TPI dimer interface.

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.

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.

Highly perturbed pKa values in the unfolded state of hen egg white lysozyme

The majority of pK(a) values in protein unfolded states are close to the amino acid model pK(a) values, thus reflecting the weak intramolecular interactions present in the unfolded ensemble of most proteins. We have carried out thermal denaturation measurements on the WT and eight mutants of HEWL from pH 1.5 to pH 11.0 to examine the unfolded state pK(a) values and the pH dependence of protein stability for this enzyme. The availability of accurate pK(a) values for the folded state of HEWL and separate measurements of mutant-induced effects on the folded state pK(a) values, allows us to estimate the pK(a) values of seven acidic residues in the unfolded state of HEWL. Asp-48 and Asp-66 display pK(a) values of 2.9 and 3.1 in our analysis, thus representing the most depressed unfolded state pK(a) values observed to date. We observe a strong correlation between the folded state pK(a) values and the unfolded state pK(a) values of HEWL, thus suggesting that the unfolded state of HEWL possesses a large degree of native state characteristics.

Triosephosphate isomerase: a highly evolved biocatalyst

Triosephosphate isomerase (TIM) is a perfectly evolved enzyme which very fast interconverts dihydroxyacetone phosphate and D: -glyceraldehyde-3-phosphate. Its catalytic site is at the dimer interface, but the four catalytic residues, Asn11, Lys13, His95 and Glu167, are from the same subunit. Glu167 is the catalytic base. An important feature of the TIM active site is the concerted closure of loop-6 and loop-7 on ligand binding, shielding the catalytic site from bulk solvent. The buried active site stabilises the enediolate intermediate. The catalytic residue Glu167 is at the beginning of loop-6. On closure of loop-6, the Glu167 carboxylate moiety moves approximately 2 Å to the substrate. The dynamic properties of the Glu167 side chain in the enzyme substrate complex are a key feature of the proton shuttling mechanism. Two proton shuttling mechanisms, the classical and the criss-cross mechanism, are responsible for the interconversion of the substrates of this enolising enzyme.

The pK(a) values of acidic and basic residues buried at the same internal location in a protein are governed by different factors

The pK(a) values of internal ionizable groups are usually very different from the normal pK(a) values of ionizable groups in water. To examine the molecular determinants of pK(a) values of internal groups, we compared the properties of Lys, Asp, and Glu at internal position 38 in staphylococcal nuclease. Lys38 titrates with a normal or elevated pK(a), whereas Asp38 and Glu38 titrate with elevated pK(a) values of 7.0 and 7.2, respectively. In the structure of the L38K variant, the buried amino group of the Lys38 side chain makes an ion pair with Glu122, whereas in the structure of the L38E variant, the buried carboxyl group of Glu38 interacts with two backbone amides and has several nearby carboxyl oxygen atoms. Previously, we showed that the pK(a) of Lys38 is normal owing to structural reorganization and water penetration concomitant with ionization of the Lys side chain. In contrast, the pK(a) values of Asp38 and Glu38 are perturbed significantly owing to an imbalance between favorable polar interactions and unfavorable contributions from dehydration and from Coulomb interactions with surface carboxylic groups. Their ionization is also coupled to subtle structural reorganization. These results illustrate the complex interplay between local polarity, Coulomb interactions, and structural reorganization as determinants of pK(a) values of internal groups in proteins. This study suggests that improvements to computational methods for pK(a) calculations will require explicit treatment of the conformational reorganization that can occur when internal groups ionize.

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.

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.

Pressure and Denaturants in the Unfolding of Triosephosphate Isomerase: The Monomeric Intermediates of the Enzymes from Saccharomyces cerevisiae and Entamoeba histolytica

Triosephosphate isomerase (TIM) is a dimeric enzyme formed by two identical (beta/alpha)8 barrels. In this work, we compare the unfolding and refolding of the TIMs from Entamoeba histolytica (EhTIM) and baker's yeast (yTIM). A monomeric intermediate was detected in the GdnHCl-induced unfolding of EhTIM. The thermodynamic, spectroscopic, catalytic, and hydrodynamic properties of this intermediate were found to be very similar to those previously described for a monomeric intermediate of yTIM observed in GdnHCl. Monomer unfolding was reversible for both TIMs; however, the dissociation step was reversible in yTIM and irreversible in EhTIM. Monomer unfolding induced by high pressure in the presence of GdnHCl was a reversible process. DeltaGUnf, DeltaVUnf, and P1/2 were obtained for the 0.7-1.2 M GdnHCl range. The linear extrapolation of these thermodynamic parameters to the absence of denaturant showed the same values for both intermediates. The DeltaVUnfH2O values calculated for EhTIM and yTIM monomeric intermediates are the same within experimental error (-57 +/- 10 and -76 +/- 14 mL/mol, respectively). These DeltaVUnf H2O values are smaller than those reported for the unfolding of monomeric proteins of similar size, suggesting that TIM intermediates are only partially hydrated. |DeltaVUnf| increased with denaturant concentration; this behavior is probably related to structural changes in the unfolded state induced by GdnHCl and pressure. From the thermodynamic parameters that were obtained, it is predicted that in the absence of denaturants, pressure would induce monomer unfolding (P1/2 approximately 140 MPa) prior to dimer dissociation (P1/2 approximately 580 MPa). Therefore, dimerization prevents the pressure unfolding of the monomer.

Redesigning protein pKa values

The ability to re-engineer enzymatic pH-activity profiles is of importance for industrial applications of enzymes. We theoretically explore the feasibility of re-engineering enzymatic pH-activity profiles by changing active site pK(a) values using point mutations. We calculate the maximum achievable DeltapK(a) values for 141 target titratable groups in seven enzymes by introducing conservative net-charge altering point mutations. We examine the importance of the number of mutations introduced, their distance from the target titratable group, and the characteristics of the target group itself. The results show that multiple mutations at 10A can change pK(a) values up to two units, but that the introduction of a requirement to keep other pK(a) values constant reduces the magnitude of the achievable DeltapK(a). The algorithm presented shows a good correlation with existing experimental data and is available for download and via a web server at http://enzyme.ucd.ie/pKD.

Fast empirical pKa prediction by Ewald summation

pK(a) calculations for macromolecules are normally performed by solving the Poisson-Boltzmann equation, accounting for the different dielectric constants of solvent and solute, as well as the ionic strength. Despite the large number of successful applications, there are some situations where the current algorithms are not suitable: (1) large scale, high-throughput analysis which requires calculations to be completed within a fraction of a second, e.g. when permanently monitoring pK(a) shifts during a molecular dynamics simulation; (2) prediction of pK(a)s in periodic boundaries, e.g. when reconstructing entire protein crystal unit cells from PDB files, including the correct protonation patterns at experimental pH. Such in silico crystals are needed by 'self-parameterizing' molecular dynamics force fields like YASARA YAMBER, that optimize their parameters while energy-minimizing high-resolution protein crystals. To address both problems, we define an empirical equation that expresses the pK(a) as a function of electrostatic potential, hydrogen bonds and accessible surface area. The electrostatic potential is evaluated by Ewald summation, which captures periodic crystal environments and the uncertainty in atom positions using Gaussian charge densities. The empirical proportionality constants are derived from 217 experimentally determined pK(a)s, and despite its simplicity, this pK(a) calculation method reaches a high overall jack-knifed accuracy, and is fast enough to be used during a molecular dynamics simulation. A reliable null-model to judge pK(a) prediction accuracies is also presented.

pKD: re-designing protein pKa values

The pK_a values in proteins govern the pH-dependence of protein stability and enzymatic activity. A large number of mutagenesis experiments have been carried out in the last three decades to re-engineer the pH-activity and pH-stability profile of enzymes and proteins. We have developed the pKD webserver (), which predicts sets of point mutations that will change the pK_a values of a set of target residues in a given direction, thus allowing for targeted re-design of the pH-dependent characteristics of proteins. The server provides the user with an interactive experience for re-designing pK_a values by pre-calculating ΔpK_a values from all feasible point mutations. Design solutions are found in less than 10 min for a typical design job for a medium-sized protein. Mutant ΔpK_a values calculated by the pKD web server are in close agreement with those produced by comparing results from full-fledged pK_a calculation methods.

Structural Tolerance versus Functional Intolerance to Mutation of Hydrophobic Core Residues Surrounding Cavities in a Parvovirus Capsid

The structural and functional relevance of amino acid residues surrounding cavities within the hydrophobic core of the protein subunits that form the capsid of parvoviruses has been investigated. Several of the evolutionarily conserved, hydrophobic residues that delimit these cavities in the capsid of the minute virus of mice were replaced by other hydrophobic residues that would affect the size and/or shape of the cavity. When four or more methylene-sized groups were introduced, or six or more groups removed, capsid assembly was drastically impaired. In contrast, the introduction or removal of up to three groups had no significant effect on capsid assembly or thermostability. However, many of these mutations affected a capsid conformational transition needed for viral infectivity. Replacement of some polar residues around the largest cavity showed that capsid assembly requires a carboxylate buried within this cavity, but both aspartate and glutamate are structurally accepted. Again, only the aspartate allowed the production of infectious viruses, because of a specific role in encapsidation of the viral genome. These observations provide evidence of a remarkable structural tolerance to mutation of the hydrophobic core of the protein subunits in a viral capsid, and of an involvement of core residues and internal cavities in capsid functions needed for infectivity.

PPD v1.0--an integrated, web-accessible database of experimentally determined protein pKa values

The Protein pK(a) Database (PPD) v1.0 provides a compendium of protein residue-specific ionization equilibria (pK(a) values), as collated from the primary literature, in the form of a web-accessible postgreSQL relational database. Ionizable residues play key roles in the molecular mechanisms that underlie many biological phenomena, including protein folding and enzyme catalysis. The PPD serves as a general protein pK(a) archive and as a source of data that allows for the development and improvement of pK(a) prediction systems. The database is accessed through an HTML interface, which offers two fast, efficient search methods: an amino acid-based query and a Basic Local Alignment Search Tool search. Entries also give details of experimental techniques and links to other key databases, such as National Center for Biotechnology Information and the Protein Data Bank, providing the user with considerable background information. The database can be found at the following URL: http://www.jenner.ac.uk/PPD.

Hydrogen bonding is the prime determinant of carboxyl pKa values at the N-termini of α-helices

Experimentally determined mean pK(a) values of carboxyl residues located at the N-termini of alpha-helices are lower than their overall mean values. Here, we perform three types of analyses to account for this phenomenon. We estimate the magnitude of the helix macrodipole to determine its potential role in lowering carboxyl pK(a) values at the N-termini. No correlation between the magnitude of the macrodipole and the pK(a) values is observed. Using the pK(a) program propKa we compare the molecular surroundings of 18 N-termini carboxyl residues versus 233 protein carboxyl groups from a previously studied database. Although pK(a) lowering interactions at the N-termini are similar in nature to those encountered in other protein regions, pK(a) lowering backbone and side-chain hydrogen bonds appear in greater number at the N-termini. For both Asp and Glu, there are about 0.5 more hydrogen bonds per residue at the N-termini than in other protein regions, which can be used to explain their lower than average pK(a) values. Using a QM-based pK(a) prediction model, we investigate the chemical environment of the two lowest Asp and the two lowest Glu pK(a) values at the N-termini so as to quantify the effect of various pK(a) determinants. We show that local interactions suffice to account for the acidity of carboxyl residues at the N-termini. The effect of the helix dipole on carboxyl pK(a) values, if any, is marginal. Backbone amide hydrogen bonds constitute the single biggest contributor to the lowest carboxyl pK(a) values at the N-termini. Their estimated pK(a) lowering effects range from about 1.0 to 1.9 pK(a) units.

Asp79 makes a large, unfavorable contribution to the stability of RNase Sa

The two most buried carboxyl groups in ribonuclease Sa (RNase Sa) are Asp33 (99% buried; pK 2.4) and Asp79 (85% buried; pK 7.4). Above these pK values, the stability of the D33A variant is 6kcal/mol less than wild-type RNase Sa, and the stability of the D79A variant is 3.3kcal/mol greater than wild-type RNase Sa. The key structural difference between the carboxyl groups is that Asp33 forms three intramolecular hydrogen bonds, and Asp79 forms no intramolecular hydrogen bond. Here, we focus on Asp79 and describe studies of 11 Asp79 variants. Most of the variants were at least 2kcal/mol more stable than wild-type RNase Sa, and the most interesting was D79F. At pH 3, below the pK of Asp79, RNase Sa is 0.3kcal/mol more stable than the D79F variant. At pH 8.5, above the pK of Asp79, RNase Sa is 3.7kcal/mol less stable than the D79F variant. The unfavorable contribution of Asp79 to the stability appears to result from the Born self-energy of burying the charge and, more importantly, from unfavorable charge-charge interactions. To counteract the effect of the negative charge on Asp79, we prepared the Q94K variant and the crystal structure showed that the amino group of the Lys formed a hydrogen-bonded ion pair (distance, 2.71A; angle, 100 degrees ) with the carboxyl group of Asp79. The stability of the Q94K variant was about the same as the wild-type at pH 3, where Asp79 is uncharged, but 1kcal/mol greater than that of wild-type RNase Sa at pH 8.5, where Asp79 is charged. Differences in hydrophobicity, steric strain, Born self-energy, and electrostatic interactions all appear to contribute to the range of stabilities observed in the variants. When it is possible, replacing buried, non-hydrogen bonded, ionizable side-chains with non-polar side-chains is an excellent means of increasing protein stability.

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

The reversible guanidinium hydrochloride-induced unfolding of Trypanosoma cruzi triosephosphate isomerase (TcTIM) was characterized under equilibrium conditions. The catalytic activity was followed as a native homodimeric functional probe. Circular dichroism, intrinsic fluorescence, and size-exclusion chromatography were used as secondary, tertiary, and quaternary structural probes, respectively. The change in ANS fluorescence intensity with increasing denaturant concentrations was also determined. The results show that two stable intermediates exist in the transition from the homodimeric native enzyme to the unfolded monomers: one (N(2*)) is a slightly more expanded, non-native, and active dimer, and the other is a partially expanded monomer (M) that binds ANS. Spectroscopic and activity data were used to reach a thermodynamic characterization. The results indicate that the Gibbs free energies for the partial reactions are 4.5 (N(2) <==> N(2*)), 65.8 (N(2*) <==> 2M), and 17.8 kJ/mol (M <==> U). It appears that TcTIM monomers are more stable than those found for other TIM species (except yeast TIM), where monomer stability is only marginal. These results are compared with those for the guanidinium hydrochloride-induced denaturation of TIM from different species, where despite the functional and three-dimensional similarities, a remarkable heterogeneity exists in the unfolding pathways.

New inhibitors of rabbit muscle triose-phosphate isomerase

We describe the synthesis and evaluation of three new competitive inhibitors of triose-phosphate isomerase. One of them (phosphoglycoloamidoxime: K(i) = 4.5 microM) is among the best reversible inhibitors so far reported for this enzyme.

Potential use of sugar binding proteins in reactors for regeneration of CO2 fixation acceptor D-Ribulose-1,5-bisphosphate

Sugar binding proteins and binders of intermediate sugar metabolites derived from microbes are increasingly being used as reagents in new and expanding areas of biotechnology. The fixation of carbon dioxide at emission source has recently emerged as a technology with potentially significant implications for environmental biotechnology. Carbon dioxide is fixed onto a five carbon sugar D-ribulose-1,5-bisphosphate. We present a review of enzymatic and non-enzymatic binding proteins, for 3-phosphoglycerate (3PGA), 3-phosphoglyceraldehyde (3PGAL), dihydroxyacetone phosphate (DHAP), xylulose-5-phosphate (X5P) and ribulose-1,5-bisphosphate (RuBP) which could be potentially used in reactors regenerating RuBP from 3PGA. A series of reactors combined in a linear fashion has been previously shown to convert 3-PGA, (the product of fixed CO₂ on RuBP as starting material) into RuBP (Bhattacharya et al., 2004; Bhattacharya, 2001). This was the basis for designing reactors harboring enzyme complexes/mixtures instead of linear combination of single-enzyme reactors for conversion of 3PGA into RuBP. Specific sugars in such enzyme-complex harboring reactors requires removal at key steps and fed to different reactors necessitating reversible sugar binders. In this review we present an account of existing microbial sugar binding proteins and their potential utility in these operations.

The determinants of carboxyl pKa values in turkey ovomucoid third domain

A computational methodology for protein pK(a) predictions, based on ab initio quantum mechanical treatment of part of the protein and linear Poisson-Boltzmann equation treatment of the bulk solvent, is presented. The method is used to predict and interpret the pK(a) values of the five carboxyl residues (Asp7, Glu10, Glu19, Asp27, and Glu43) in the serine protease inhibitor turkey ovomucoid third domain. All the predicted pK(a) values are within 0.5 pH units of experiment, with a root-mean-square deviation of 0.31 pH units. We show that the decreased pK(a) values observed for some of the residues are primarily due to hydrogen bonds to the carboxyl oxygens. Hydrogen bonds involving amide protons are shown to be particularly important, and the effect of hydrogen bonding is shown to be nonadditive. Hydrophobic effects are also shown to be important in raising the pK(a). Interactions with charged residues are shown to have relatively little effect on the carboxyl pK(a) values in this protein, in general agreement with experiment.

Dissecting Structural and Electrostatic Interactions of Charged Groups in α-Sarcin. An NMR Study of Some Mutants Involving the Catalytic Residues

The cytotoxic ribonuclease alpha-sarcin is the best characterized member of the ribotoxin family. Ribotoxins share a common structural core, catalytic residues, and active site topology with members of the broader family of nontoxic microbial extracellular RNases. They are, however, much more specific in their biological action. To shed light on the highly specific alpha-sarcin activity, we have evaluated the structural and electrostatic interactions of its charged groups, by combining the structural and pK(a) characterization by NMR of several variants with theoretical calculations based on the Tanford-Kirkwood and Poisson-Boltzmann models. The NMR data reveal that the global conformation of wild-type alpha-sarcin is preserved in the H50Q, E96Q, H137Q, and H50/137Q variants, and that His137 is involved in an H-bond that is crucial in maintaining the active site structure and in reinforcing the stability of the enzyme. The loss of this H-bond in the H137Q and H50/137Q variants modifies the local structure of the active site. The pK(a) values of active site groups H50, E96, and H137 in the four variants have been determined by two-dimensional NMR. The catalytic dyad of E96 and H137 is not sensitive to charge replacements, since their pK(a) values vary less than +/-0.3 pH unit with respect to those of the wild type. On the contrary, the pK(a) of His50 undergoes drastic changes when compared to its value in the intact protein. These amount to an increase of 0.5 pH unit or a decrease of 1.1 pH units depending on whether a positive or negative charge is substituted at the active site. The main determinants of the pK(a) values of most of the charged groups in alpha-sarcin have been established by considering the NMR results in conjunction with those derived from theoretical pK(a) calculations. With regard to the active site residues, the H50 pK(a) is chiefly influenced by electrostatic interactions with E96 and H137, whereas the effect of the low dielectric constant and the interaction with R121 appear to be the main determinants of the altered pK(a) value of E96 and H137. Charge-charge interactions and an increased level of burial perturb the pK(a) values of the active site residues of alpha-sarcin, which can account for its reduced ribonucleolytic activity and its high specificity.

Thermodynamic characterization of yeast triosephosphate isomerase refolding: insights into the interplay between function and stability as reasons for the oligomeric nature of the enzyme

The reasons underlying the oligomeric nature of some proteins such as triosephosphate isomerase (TIM) are unclear. It has been proposed that this enzyme is an oligomer, mainly because of its stability rather than for functional reasons. To address this issue, the reversible denaturation and renaturation of the homodimeric TIM from baker's yeast ( Saccharomyces cerevisiae ) induced by guanidinium chloride and urea have been characterized by spectroscopic, functional and hydrodynamic techniques. The unfolding and refolding of this enzyme are not coincident after 'conventional' equilibrium times. Unfolding experiments did not reach equilibrium, owing to a very slow dissociation and/or unfolding process. By contrast, equilibrium was reached in the refolding direction. The simplest equilibrium pathway compatible with the obtained data was found to be a three-state process involving an inactive and expanded monomer. The Gibbs energy changes for monomer folding (delta G (0)(fold) = -16.6+/-0.7 kJ x mol(-1)) and monomer association (delta G (0)(assoc) = -70.3+/-1.1 kJ x mol(-1)) were calculated from data obtained in the two denaturants. From an analysis of the present data and data from the literature on the stability of TIM from different species and for other beta/alpha barrels, and model simulations on the effect of stability in the catalytic activity of the enzyme, it is concluded that the low stability of the monomers is neither the only, nor the main, cause for the dimeric nature of TIM. There is interplay between function and stability.

Crystal structure of triosephosphate isomerase complexed with 2-phosphoglycolate at 0.83-A resolution

The atomic resolution structure of Leishmania mexicana triosephosphate isomerase complexed with 2-phosphoglycolate shows that this transition state analogue is bound in two conformations. Also for the side chain of the catalytic glutamate, Glu(167), two conformations are observed. In both conformations, a very short hydrogen bond exists between the carboxylate group of the ligand and the catalytic glutamate. The distance between O11 of PGA and Oepsilon2 of Glu(167) is 2.61 and 2.55 A for the major and minor conformations, respectively. In either conformation, Oepsilon1 of Glu(167) is hydrogen-bonded to a water network connecting the side chain with bulk solvent. This network also occurs in two mutually exclusive arrangements. Despite the structural disorder in the active site, the C termini of the beta strands that construct the active site display the least anisotropy compared with the rest of the protein. The loops following these beta strands display various degrees of anisotropy, with the tip of the dimer interface loop 3 having very low anisotropy and the C-terminal region of the active site loop 6 having the highest anisotropy. The pyrrolidine ring of Pro(168) at the N-terminal region of loop 6 is in a strained planar conformation to facilitate loop opening and product release.

On the evaluation and optimization of protein X-ray structures for pKa calculations

The calculation of the physical properties of a protein from its X-ray structure is of importance in virtually every aspect of modern biology. Although computational algorithms have been developed for calculating everything from the dynamics of a protein to its binding specificity, only limited information is available on the ability of these methods to give accurate results when used with a particular X-ray structure. We examine the ability of a pKa calculation algorithm to predict the proton-donating residue in the catalytic mechanism of hen egg white lysozyme. We examine the correlation between the ability of the pKa calculation method to obtain the correct result and the overall characteristics of 41 X-ray structures such as crystallization conditions, resolution, and the output of structure validation software. We furthermore examine the ability of energy minimizations (EM), molecular dynamics (MD) simulations, and structure-perturbation methods to optimize the X-ray structures such that these give correct results with the pKa calculation algorithm. We propose a set of criteria for identifying the proton donor in a catalytic mechanism, and demonstrate that the application of these criteria give highly accurate prediction results when using unmodified X-ray structures. More specifically, we are able to successfully identify the proton donor in 85% of the X-ray structures when excluding structures with crystal contacts near the active site. Neither the use of the overall characteristics of the X-ray structures nor the optimization of the structure by EM, MD, or other methods improves the results of the pKa calculation algorithm. We discuss these results and their implications for the design of structure-based energy calculation algorithms in general.

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.

Inhibition of plasmodium falciparum triose-phosphate isomerase by chemical modification of an interface cysteine. Electrospray ionization mass spectrometric analysis of differential cysteine reactivities

Plasmodium falciparum triose-phosphate isomerase, a homodimeric enzyme, contains four cysteine residues at positions 13, 126, 196, and 217 per subunit. Among these, Cys-13 is present at the dimer interface and is replaced by methionine in the corresponding human enzyme. We have investigated the effect of sulfhydryl labeling on the parasite enzyme, with a view toward developing selective covalent inhibitors by targeting the interface cysteine residue. Differential labeling of the cysteine residues by iodoacetic acid and iodoacetamide has been followed by electrospray ionization mass spectrometry and positions of the labels determined by analysis of tryptic fragments. The rates of labeling follows the order Cys-196 > Cys-13 Cys-217/Cys-126, which correlates well with surface accessibility calculations based on the enzyme crystal structure. Iodoacetic acid labeling leads to a soluble, largely inactive enzyme, whereas IAM labeling leads to precipitation. Carboxyl methylation of Cys-13 results in formation of monomeric species detectable by gel filtration. Studies with an engineered C13D mutant permitted elucidation of the effects of introducing a negative charge at the interface. The C13D mutant exhibits a reduced stability to denaturants and 7-fold reduction in the enzymatic activity even under the concentrations in which dimeric species are observed.

Subunit interface mutation disrupting an aromatic cluster in Plasmodium falciparum triosephosphate isomerase: effect on dimer stability

A mutation at the dimer interface of Plasmodium falciparum triosephosphate isomerase (PfTIM) was created by mutating a tyrosine residue at position 74, at the subunit interface, to glycine. Tyr74 is a critical residue, forming a part of an aromatic cluster at the interface. The resultant mutant, Y74G, was found to have considerably reduced stability compared with the wild-type protein (TIMWT). The mutant was found to be much less stable to denaturing agents such as urea and guanidinium chloride. Fluorescence and circular dichroism studies revealed that the Y74G mutant and TIMWT have similar spectroscopic properties, suggestive of similar folded structures. Further, the Y74G mutant also exhibited a concentration-dependent loss of enzymatic activity over the range 0.1-10 microM. In contrast, the wild-type enzyme did not show a concentration dependence of activity in this range. Fluorescence quenching of intrinsic tryptophan emission was much more efficient in case of Y74G than TIMWT, suggestive of greater exposure of Trp11, which lies adjacent to the dimer interface. Analytical gel filtration studies revealed that in Y74G, monomeric and dimeric species are in dynamic equilibrium, with the former predominating at low protein concentration. Spectroscopic studies established that the monomeric form of the mutant is largely folded. Low concentrations of urea also drive the equilibrium towards the monomeric form. These studies suggest that the replacement of tyrosine with a small residue at the interface of triosephosphate isomerase weakens the subunit-subunit interactions, giving rise to structured, but enzymatically inactive, monomers at low protein concentration.

Unfolding of triosephosphate isomerase from Trypanosoma brucei: identification of intermediates and insight into the denaturation pathway using tryptophan mutants

The unfolding of triosephosphate isomerase (TIM) from Trypanosoma brucei (TbTIM) induced by guanidine hydrochloride (GdnHCl) was characterized. In contrast to other TIMs, where unfolding is a two or three state process, TbTIM showed two intermediates. The solvent exposure of different regions of the protein in the unfolding process was characterized spectroscopically with mutant proteins in which tryptophans (W) were changed to phenlylalanines (F). The midpoints of the transitions measured by circular dichroism, intrinsic fluorescence, and catalytic activity, as well as the increase in 1-aniline 8-naphthalene sulfonate fluorescence, show that the native state was destabilized in the W12F and W12F/W193F mutants, relative to the wild-type enzyme. Using the hydrodynamic profile for the unfolding of a monomeric TbTIM mutant (RMM0-1TIM) measured by size-exclusion chromatography as a standard, we determined the association state of these intermediates: D*, a partially expanded dimer, and M*, a partially expanded monomeric intermediate. High-molecular-weight aggregates were also detected. At concentrations over 2.0 M GdnHCl, the hydrodynamic properties of TbTIM and RMM0-1TIM are the same, suggesting that the dimeric intermediate dissociates and the unfolding proceeds through the denaturation of an expanded monomeric intermediate. The analysis of the denaturation process of the TbTIM mutants suggests a sequence for the gradual exposure of W residues: initially the expansion of the native dimer to form D* affects the environments of W12 and W159. The dissociation of D* to M* and further unfolding of M* to U induces the exposure of W170. The role of protein concentration in the formation of intermediates and aggregates is discussed considering the irreversibility of this unfolding process.

Structural determinants for ligand binding and catalysis of triosephosphate isomerase

The crystal structure of leishmania triosephosphate isomerase (TIM) complexed with 2-(N-formyl-N-hydroxy)-aminoethyl phosphonate (IPP) highlights the importance of Asn11 for binding and catalysis. IPP is an analogue of the substrate D-glyceraldehyde-3-phosphate, and it is observed to bind with its aldehyde oxygen in an oxyanion hole formed by ND2 of Asn11 and NE2 of His95. Comparison of the mode of binding of IPP and the transition state analogue phosphoglycolohydroxamate (PGH) suggests that the Glu167 side chain, as well as the triose part of the substrate, adopt different conformations as the catalysed reaction proceeds. Comparison of the TIM-IPP and the TIM-PGH structures with other liganded and unliganded structures also highlights the conformational flexibility of the ligand and the active site, as well as the conserved mode of ligand binding.

Modeling, mutagenesis, and structural studies on the fully conserved phosphate-binding loop (loop 8) of triosephosphate isomerase: toward a new substrate specificity

Loop 8 (residues 232-242) in triosephosphate isomerase (TIM) is a highly conserved loop that forms a tight binding pocket for the phosphate moiety of the substrate. Its sequence includes the fully conserved, solvent-exposed Leu238. The tight phosphate-binding pocket explains the high substrate specificity of TIM being limited to the in vivo substrates dihydroxyacetone-phosphate and D-glyceraldehyde-3-phosphate. Here we use the monomeric variant of trypanosomal TIM for exploring the structural consequences of shortening this loop. The mutagenesis, guided by extensive modeling calculations and followed up by crystallographic characterization, is aimed at widening the phosphate-binding pocket and, consequently, changing the substrate specificity. Two new variants were characterized. The crystal structures of these variants indicate that in monomeric forms of TIM, the Leu238 side-chain is nicely buried in a hydrophobic cluster. Monomeric forms of wild-type dimeric TIM are known to exist transiently as folding intermediates; our structural analysis suggests that in this monomeric form, Leu238 of loop 8 also adopts this completely buried conformation, which explains its full conservation across the evolution. The much wider phosphate-binding pocket of the new variant allows for the development of a new TIM variant with a different substrate specificity.