Closed Conformation of the Active Site Loop of Rabbit Muscle Triosephosphate Isomerase in the Absence of Substrate: Evidence of Conformational Heterogeneity

Crystal Structure Analysis and Conformational Epitope Mutation of Triosephosphate Isomerase, a Mud Crab Allergen

The triosephosphate isomerase (TIM), Scy p 8, is a crab allergen and shows cross-reactivity in the shellfish. Here, recombinant Scy p 8 was expressed, and its crystal structure was determined at a resolution of 1.8 Å. The three-dimensional structure of Scy p 8 is primarily composed of a (β/α)₈-barrel motif prototype. Additionally, Scy p 8 showed cross-reactivity with high sequential and secondary structural identity among TIMs from shellfish species. The site-directed mutagenesis of critical amino acids of conformational epitopes was carried out, and the mutants of Trp 168 and Lys 237 to Ala reduced immunoglobulin E (IgE)-binding activity by approximately 30%, compared with wild-type TIM in an inhibition ELISA; however, it still induced basophil activation despite the interpatient variability between patients. These results can help to provide an accurate template for the analysis of the IgE binding and establish meaningful relationships between structure and allergenicity.

Characterization of human triosephosphate isomerase S-nitrosylation

Triosephosphate isomerase (TPI), the glycolytic enzyme that catalyzes the isomerization of dihydroxyacetone phosphate (DHAP) to glyceraldehyde-3-phosphate (G3P), has been frequently identified as a target of S-nitrosylation by proteomic studies. However, the effect of S-nitrosylation on its activity has only been explored in plants and algae. Here, we describe the in vitro S-nitrosylation of human TPI (hTPI), and the effect of the modification on its enzymatic parameters. NO-incorporation into the enzyme cysteine residues occurred by a time-dependent S-transnitrosylation from both, S-nitrosocysteine (CySNO) and S-nitrosoglutathione (GSNO), with CySNO being the more efficient NO-donor. Both X-ray crystal structure and mass spectrometry analyses showed that only Cys217 was S-nitrosylated. hTPI S-nitrosylation produced a 30% inhibition of the Vmax of the DHAP conversion to G3P, without affecting the Km for DHAP. This is the first study describing features of human TPI S-nitrosylation.

Biopharmaceutical Informatics: supporting biologic drug development via molecular modelling and informatics

Objectives

The purpose of this article is to introduce an emerging field called ‘Biopharmaceutical Informatics’. It describes how tools from Information technology and Molecular Biophysics can be adapted, developed and gainfully employed in discovery and development of biologic drugs.

Key Findings

The findings described here are based on literature surveys and the authors’ collective experiences in the field of biologic drug product development. A strategic framework to forecast early the hurdles faced during drug product development is weaved together and elucidated using chemical degradation as an example. Efficiency of translating biologic drug discoveries into drug products can be significantly improved by combining learnings from experimental biophysical and analytical data on the drug candidates with molecular properties computed from their sequences and structures via molecular modeling and simulations.

Summary

Biopharmaceutical Informatics seeks to promote applications of computational tools towards discovery and development of biologic drugs. When fully implemented, industry-wide, it will enable rapid materials-free developability assessments of biologic drug candidates at early stages as well as streamline drug product development activities such as commercial scale production, purification, formulation, analytical characterization, safety and in vivo performance.

Ion Mobility-Mass Spectrometry Analysis of Cross-Linked Intact Multiprotein Complexes: Enhanced Gas-Phase Stabilities and Altered Dissociation Pathways

Analysis of protein complexes by ion mobility-mass spectrometry is a valuable method for the rapid assessment of complex composition, binding stoichiometries, and structures. However, capturing labile, unknown protein assemblies directly from cells remains a challenge for the technology. Furthermore, ion mobility-mass spectrometry measurements of complexes, subcomplexes, and subunits are necessary to build complete models of intact assemblies, and such data can be difficult to acquire in a comprehensive fashion. Here, we present the use of novel mass spectrometry cleavable cross-linkers and tags to stabilize intact protein complexes for ion mobility-mass spectrometry. Our data reveal that tags and linkers bearing permanent charges are superior stabilizers relative to neutral cross-linkers, especially in the context of retaining compact forms of the assembly under a wide array of activating conditions. In addition, when cross-linked protein complexes are collisionally activated in the gas phase, a larger proportion of the product ions produced are often more compact and reflect native protein subcomplexes when compared with unmodified complexes activated in the same fashion, greatly enabling applications in structural biology.

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

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

A novel inhibitor of Mammalian triosephosphate isomerase found by an in silico approach

Triosephosphate isomerase (TIM) is an essential, highly conserved component of glycolysis. Tumors are often dependent on glycolysis for energy and metabolite production (the Warburg effect). Glycolysis inhibitors thus show promise as cancer treatments. TIM inhibition, unlike inhibition of other glycolysis enzymes, also produces toxic methylglyoxal targeted to regions of high glycolysis, an effect that might also be therapeutically useful. Thus TIM is an attractive drug target. A total of 338,562 lead-like molecules were analyzed computationally to find TIM inhibitors by an efficient “double screen” approach. The first fragment-sized compounds were studied using structure-based virtual screening to identify binding motifs for mammalian TIM. Subsequently, larger compounds, filtered to meet the binding criteria developed in the first analysis, were ranked using a second round of structure-based virtual screening. A compound was found that inhibited mammalian TIM in vitro in the micromolar range. Docking and molecular dynamics (MD) suggested that the inhibitor made hydrogen bond contacts with TIM catalytic residues. In addition, hydrophobic contacts were made throughout the binding site. All predicted inhibitor-TIM interactions involved TIM residues that were highly conserved. The discovered compound may provide a scaffold for elaboration of other inhibitors.

Inhibition of triosephosphate isomerase by phosphoenolpyruvate in the feedback-regulation of glycolysis

The inhibition of triosephosphate isomerase (TPI) in glycolysis by the pyruvate kinase (PK) substrate phosphoenolpyruvate (PEP) results in a newly discovered feedback loop that counters oxidative stress in cancer and actively respiring cells. The mechanism underlying this inhibition is illuminated by the co-crystal structure of TPI with bound PEP at 1.6 Å resolution, and by mutational studies guided by the crystallographic results. PEP is bound to the catalytic pocket of TPI and occludes substrate, which accounts for the observation that PEP competitively inhibits the interconversion of glyceraldehyde-3-phosphate and dihydroxyacetone phosphate. Replacing an isoleucine residue located in the catalytic pocket of TPI with valine or threonine altered binding of substrates and PEP, reducing TPI activity in vitro and in vivo. Confirming a TPI-mediated activation of the pentose phosphate pathway (PPP), transgenic yeast cells expressing these TPI mutations accumulate greater levels of PPP intermediates and have altered stress resistance, mimicking the activation of the PK-TPI feedback loop. These results support a model in which glycolytic regulation requires direct catalytic inhibition of TPI by the pyruvate kinase substrate PEP, mediating a protective metabolic self-reconfiguration of central metabolism under conditions of oxidative stress.

High-resolution crystal structure and redox properties of chloroplastic triosephosphate isomerase from Chlamydomonas reinhardtii

Triosephosphate isomerase (TPI) catalyzes the interconversion of glyceraldehyde-3-phosphate to dihydroxyacetone phosphate. Photosynthetic organisms generally contain two isoforms of TPI located in both cytoplasm and chloroplasts. While the cytoplasmic TPI is involved in the glycolysis, the chloroplastic isoform participates in the Calvin-Benson cycle, a key photosynthetic process responsible for carbon fixation. Compared with its cytoplasmic counterpart, the functional features of chloroplastic TPI have been poorly investigated and its three-dimensional structure has not been solved. Recently, several studies proposed TPI as a potential target of different redox modifications including dithiol/disulfide interchanges, glutathionylation, and nitrosylation. However, neither the effects on protein activity nor the molecular mechanisms underlying these redox modifications have been investigated. Here, we have produced recombinantly and purified TPI from the unicellular green alga Chlamydomonas reinhardtii (Cr). The biochemical properties of the enzyme were delineated and its crystallographic structure was determined at a resolution of 1.1 Å. CrTPI is a homodimer with subunits containing the typical (β/α)8-barrel fold. Although no evidence for TRX regulation was obtained, CrTPI was found to undergo glutathionylation by oxidized glutathione and trans-nitrosylation by nitrosoglutathione, confirming its sensitivity to multiple redox modifications.

Structural and functional perturbation of Giardia lamblia triosephosphate isomerase by modification of a non-catalytic, non-conserved region

BACKGROUND

We have previously proposed triosephosphate isomerase of Giardia lamblia (GlTIM) as a target for rational drug design against giardiasis, one of the most common parasitic infections in humans. Since the enzyme exists in the parasite and the host, selective inhibition is a major challenge because essential regions that could be considered molecular targets are highly conserved. Previous biochemical evidence showed that chemical modification of the non-conserved non-catalytic cysteine 222 (C222) inactivates specifically GlTIM. The inactivation correlates with the physicochemical properties of the modifying agent: addition of a non-polar, small chemical group at C222 reduces the enzyme activity by one half, whereas negatively charged, large chemical groups cause full inactivation.

RESULTS

In this work we used mutagenesis to extend our understanding of the functional and structural effects triggered by modification of C222. To this end, six GlTIM C222 mutants with side chains having diverse physicochemical characteristics were characterized. We found that the polarity, charge and volume of the side chain in the mutant amino acid differentially alter the activity, the affinity, the stability and the structure of the enzyme. The data show that mutagenesis of C222 mimics the effects of chemical modification. The crystallographic structure of C222D GlTIM shows the disruptive effects of introducing a negative charge at position 222: the mutation perturbs loop 7, a region of the enzyme whose interactions with the catalytic loop 6 are essential for TIM stability, ligand binding and catalysis. The amino acid sequence of TIM in phylogenetic diverse groups indicates that C222 and its surrounding residues are poorly conserved, supporting the proposal that this region is a good target for specific drug design.

CONCLUSIONS

The results demonstrate that it is possible to inhibit species-specifically a ubiquitous, structurally highly conserved enzyme by modification of a non-conserved, non-catalytic residue through long-range perturbation of essential regions.

Novel c-di-GMP recognition modes of the mouse innate immune adaptor protein STING

The mammalian ER protein STING (stimulator of interferon genes; also known as MITA, ERIS, MPYS or TMEM173) is an adaptor protein that links the detection of cytosolic dsDNA to the activation of TANK-binding kinase 1 (TBK1) and its downstream transcription factor interferon regulatory factor 3 (IFN3). Recently, STING itself has been found to be the direct receptor of bacterial c-di-GMP, and crystal structures of several human STING C-terminal domain (STING-CTD) dimers in the apo form or in complex with c-di-GMP have been published. Here, a novel set of structures of mouse STING-CTD (mSTING(137-344)) in apo and complex forms determined from crystals obtained under different crystallization conditions are reported. These novel closed-form structures exhibited considerable differences from previously reported open-form human STING-CTD structures. The novel mSTING structures feature extensive interactions between the two monomers, a unique asymmetric c-di-GMP molecule with one guanine base in an unusual syn conformation that is well accommodated in the dimeric interface with many direct specific interactions and two unexpected equivalent secondary peripheral c-di-GMP binding sites. Replacement of the amino acids crucial for specific c-di-GMP binding in mSTING significantly changes the ITC titration profiles and reduces the IFN-β reporter luciferase activity. Taken together, these results reveal a more stable c-di-GMP binding mode of STING proteins that could serve as a template for rational drug design to stimulate interferon production by mammalian cells.

Crystal structures of triosephosphate isomerase from methicillin resistant Staphylococcus aureus MRSA252 provide structural insights into novel modes of ligand binding and unique conformations of catalytic loop

Staphylococcus aureus is one of the most dreaded pathogens worldwide and emergence of notorious antibiotic resistant strains have further exacerbated the present scenario. The glycolytic enzyme, triosephosphate isomerase (TIM) is one of the cell envelope proteins of the coccus and is involved in biofilm formation. It also plays an instrumental role in adherence and invasion of the bacteria into the host cell. To structurally characterize this important enzyme and analyze it's interaction with different inhibitors, substrate and transition state analogues, the present article describes several crystal structures of SaTIM alone and in complex with different ligands: glycerol-3-phosphate (G3P), glycerol-2-phosphate (G2P), 3-phosphoglyceric acid (3PG) and 2-phosphoglyceric acid (2PG). Unique conformations of the catalytic loop 6 (L6) has been observed in the different complexes. It is found to be in "almost closed" conformation in both subunits of the structure complexed to G3P. However L6 adopts the open conformation in presence of G2P and 2PG. The preference of the conformation of the catalytic loop can be correlated with the position of the phosphate group in the ligand. Novel modes of binding have been observed for G2P and 3PG for the very first time. The triose moiety is oriented away from the catalytic residues and occupies an entirely different position in some subunits. A completely new binding site for phosphate has also been identified in the complex with 2PG which differs substantially from the conventional phosphate binding site of the ligand in the crystal structures of TIM determined so far.

GalaxyWEB server for protein structure prediction and refinement

Three-dimensional protein structures provide invaluable information for understanding and regulating biological functions of proteins. The GalaxyWEB server predicts protein structure from sequence by template-based modeling and refines loop or terminus regions by ab initio modeling. This web server is based on the method tested in CASP9 (9th Critical Assessment of techniques for protein Structure Prediction) as ‘Seok-server’, which was assessed to be among top performing template-based modeling servers. The method generates reliable core structures from multiple templates and re-builds unreliable loops or termini by using an optimization-based refinement method. In addition to structure prediction, a user can also submit a refinement only job by providing a starting model structure and locations of loops or termini to refine. The web server can be freely accessed at http://galaxy.seoklab.org/.

Initiation of the reaction of deamidation in triosephosphate isomerase: investigations by means of molecular dynamics simulations

Deamidation of asparagine is the spontaneous degradation of this residue into aspartic acid. The kinetics of this slow reaction is mainly dependent on the nature of the adjacent amino acid that follows asparagine in a peptide or protein primary sequence. In the homodimer triosephosphate isomerase (TPI), there are two main deamidation sites per subunit: Asn15-Gly16 and Asn71-Gly72 for which deamidation dynamics are known to be interrelated. In this study, we investigate the initiation of the deamidation reaction in TPI by means of molecular dynamics. Simulations based on classical AMBER force field are performed in a 60 to 90 ns time scale for six distinct samples. Conformational changes, desolvation effects, and hydrogen bond networks are analyzed to interpret the experimental findings and previous quantum mechanical (QM) results. Results that are based on desolvation analysis clarify the assignments in the literature about the different behaviors of two deamidating sites in TPI. Conformational analysis supports findings suggested by QM studies: the most favorable reaction mechanism is the one that yields to succinimide intermediate via one or two step routes. The mechanism leading to the succinimide intermediate most likely involves the formation of a tetrahedral intermediate that is formed either directly from asparagine or via a side chain tautomer intermediate. In all cases, surrounding water molecules are present to assist the reaction.

Cooperativity in monomeric enzymes with single ligand-binding sites

Cooperativity is widespread in biology. It empowers a variety of regulatory mechanisms and impacts both the kinetic and thermodynamic properties of macromolecular systems. Traditionally, cooperativity is viewed as requiring the participation of multiple, spatially distinct binding sites that communicate via ligand-induced structural rearrangements; however, cooperativity requires neither multiple ligand binding events nor multimeric assemblies. An underappreciated manifestation of cooperativity has been observed in the non-Michaelis-Menten kinetic response of certain monomeric enzymes that possess only a single ligand-binding site. In this review, we present an overview of kinetic cooperativity in monomeric enzymes. We discuss the primary mechanisms postulated to give rise to monomeric cooperativity and highlight modern experimental methods that could offer new insights into the nature of this phenomenon. We conclude with an updated list of single subunit enzymes that are suspected of displaying cooperativity, and a discussion of the biological significance of this unique kinetic response.

Ionic strength dependence of F-actin and glycolytic enzyme associations: a Brownian dynamics simulations approach

The association of glycolytic enzymes with F-actin is proposed to be one mechanism by which these enzymes are compartmentalized, and, as a result, may possibly play important roles for: regulation of the glycolytic pathway, potential substrate channeling, and increasing glycolytic flux. Historically, in vitro experiments have shown that many enzyme/actin interactions are dependent on ionic strength. Herein, Brownian dynamics (BD) examines how ionic strength impacts the energetics of the association of F-actin with the glycolytic enzymes: lactate dehydrogenase (LDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), fructose-1,6-bisphosphate aldolase (aldolase), and triose phosphate isomerase (TPI). The BD simulations are steered by electrostatics calculated by Poisson-Boltzmann theory. The BD results confirm experimental observations that the degree of association diminishes as ionic strength increases but also suggest that these interactions are significant, at physiological ionic strengths. Furthermore, BD agrees with experiments that muscle LDH, aldolase, and GAPDH interact significantly with F-actin whereas TPI does not. BD indicates similarities in binding regions for aldolase and LDH among the different species investigated. Furthermore, the residues responsible for salt bridge formation in stable complexes persist as ionic strength increases. This suggests the importance of the residues determined for these binary complexes and specificity of the interactions. That these interactions are conserved across species, and there appears to be a general trend among the enzymes, support the importance of these enzyme-F-actin interactions in creating initial complexes critical for compartmentation.

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.

Conformational differences between the wild type and V30M mutant transthyretin modulate its binding to genistein: implications to tetramer stability and ligand-binding

Transthyretin (TTR) is a tetrameric beta-sheet-rich transporter protein directly involved in human amyloid diseases. It was recently found that the isoflavone genistein (GEN) potently inhibits TTR amyloid fibril formation (Green et al., 2005) and is therefore a promising candidate for TTR amyloidosis treatment. Here we used structural and biophysical approaches to characterize genistein binding to the wild type (TTRwt) and to its most frequent amyloidogenic variant, the V30M mutant. In a dose-dependent manner, genistein elicited considerable increases in both mutant and TTRwt stability as demonstrated by high hydrostatic pressure (HHP) and acid-mediated dissociation/denaturation assays. TTR:GEN crystal complexes and isothermal titration calorimetry (ITC) experiments showed that the binding mechanisms of genistein to the TTRwt and to V30M are different and are dependent on apoTTR structure conformations. Furthermore, we could also identify potential allosteric movements caused by genistein binding to the wild type TTR that explains, at least in part, the frequently observed negatively cooperative process between the two sites of TTRwt when binding ligands. These findings show that TTR mutants may present different ligand recognition and therefore are of value in ligand design for inhibiting TTR amyloidosis.

A coleopteran triosephosphate isomerase: X-ray structure and phylogenetic impact of insect sequences

A coleopteran triosephosphate isomerase (TIM) from Tenebrio molitor (yellow mealworm beetle) was recombinantly expressed in Escherichia coli and characterized with respect to thermal stability, kinetic parameters and oligomeric state. The enzyme was successfully crystallized and the structure determined by X-ray analysis to 2.0 A resolution. This is the first example of an invertebrate TIM. We compare structural features with known structures of TIMs from microorganisms, plants and vertebrates, and discuss the utility of the Tenebrio TIM sequence, together with several newly sequenced insect TIMs, for molecular phylogenetic analysis.

Expression, purification, crystallization and preliminary X-ray diffraction studies of triosephosphate isomerase from methicillin-resistant Staphylococcus aureus (MRSA252)

Triosephosphate isomerase from methicillin-resistant Staphylococcus aureus (MRSA252) was cloned in pQE30 vector, overexpressed in Escherichia coli M15 (pREP4) cells and purified to homogeneity. The protein was crystallized from 1.6 M trisodium citrate dihydrate pH 6.5 using the hanging-drop vapour-diffusion method. The crystals belonged to space group P4(3)2(1)2, with unit-cell parameters a = b = 79.15, c = 174.27 A. X-ray diffraction data were collected and processed to a maximum resolution of 1.9 A. The presence of two molecules in the asymmetric unit gave a Matthews coefficient (V(M)) of 2.64 A(3) Da(-1), with a solvent content of 53.63%.

Conformational selection in silico: loop latching motions and ligand binding in enzymes

Ligand binding frequently induces significant conformational changes in a protein receptor. Understanding and predicting such conformational changes represent an important challenge for computational biology, including applications to structure-based drug design. We describe an approach to this problem based on the assumption that the holo state is at least transiently populated in the absence of a ligand; this hypothesis has been referred to as "conformational selection." Here, we apply a method that tests this hypothesis on a challenging class of ligand-induced conformational changes, which we refer to as loop latching: the closing of a loop around an active site that sequesters the ligand from solvent. The method uses a combination of replica exchange molecular dynamics and a loop prediction algorithm to generate low-energy loop structures, and docking to select the conformation appropriate for binding a particular ligand. On a test set of six proteins, it yields loop structures including hololike conformations, generally below 2 A RMSD from the liganded structure, for loops that span up to 15 residues. Docking serves as a stringent test of the predictions. In five of the six cases, the predicted loop conformations improve the ranks of cognate ligands relative to using the apo structure, although the results remain, in most cases, significantly worse than using a holo structure. The poses of the cognate ligands are correct in four of the six test cases, while they are correct for five of the six using a holo structure.

Slow proton transfer from the hydrogen-labelled carboxylic acid side chain (Glu-165) of triosephosphate isomerase to imidazole buffer in D2O

The catalytic base at the active site of triosephosphate isomerase (TIM) was labelled with -H by abstraction of a proton from substrate d-glyceraldehyde 3-phosphate to form an enzyme-bound enediol(ate) in D2O solvent. The partitioning of this labelled enzyme between intramolecular transfer of -H to form dihydroxyacetone phosphate (DHAP), and irreversible exchange with -D from solvent was examined by determining the yields of H- and D-labelled products by 1H NMR spectroscopy. The yield of hydrogen-labelled product DHAP remains constant as the concentration of the basic form of imidazole buffer is increased from 0.014 to 0.56 M. This shows that the active site of free TIM, which has an open conformation needed to allow substrate binding, adopts a closed conformation at the enediolate-complex intermediate where the catalytic side chain is sequestered from interaction with imidazole dissolved in D2O.

Comparative analysis of structural and dynamic properties of the loaded and unloaded hemophore HasA: functional implications

A heme-acquisition system present in several Gram-negative bacteria requires the secretion of hemophores. These extracellular carrier proteins capture heme and deliver it to specific outer membrane receptors. The Serratia marcescens HasA hemophore is a monodomain protein that binds heme with a very high affinity. Its alpha/beta structure, as that of its binding pocket, has no common features with other iron- or heme-binding proteins. Heme is held by two loops L1 and L2 and coordinated to iron by an unusual ligand pair, H32/Y75. Two independent regions of the hemophore beta-sheet are involved in HasA-HasR receptor interaction. Here, we report the 3-D NMR structure of apoHasA and the backbone dynamics of both loaded and unloaded hemophore. While the overall structure of HasA is very similar in the apo and holo forms, the hemophore presents a transition from an open to a closed form upon ligand binding, through a large movement, of up to 30 A, of loop L1 bearing H32. Comparison of loaded and unloaded HasA dynamics on different time scales reveals striking flexibility changes in the binding pocket. We propose a mechanism by which these structural and dynamic features provide the dual function of heme binding and release to the HasR receptor.

Dynamics, Hydration, and Motional Averaging of a Loop-Gated Artificial Protein Cavity: The W191G Mutant of Cytochrome c Peroxidase in Water as Revealed by Molecular Dynamics Simulations

Five molecular dynamics simulations of the W191G cavity mutant of cytochrome c peroxidase in explicit water reveal distinct dynamic and hydration behavior depending on the closed or open state of the flexible loop gating the cavity, the binding of (K+ or small molecule) cations, and the system temperature. The conformational spaces sampled by the loop region and by the cavity significantly reduce upon binding. The largest ordering factor on water dynamics is the presence of the K+ ion occupying the gated cavity. Considerable water exchange occurs for the open-gate cavity when no ligand or cation is bound. In all cases, good correspondence is found between the calculated (ensemble-averaged) location of water molecules and the water sites determined by X-ray crystallography experiments. However, our simulations suggest that these sites do not necessarily correspond to the presence of bound water molecules. In fact, individual water molecules may repeatedly exchange within the cavity volume yet occupy on average these water sites. Four major conclusions emerge. First, it seems misleading to interpret the conformation of protein loop regions in terms of single dominant structures. Second, our simulations support the general picture of Pro 190 cis-trans isomerization as a determinant of the loop-opening mechanism. Third, receptor flexibility is fundamental for ligand binding and molecular recognition, and our results suggest its importance for the docking of small compounds to the artificial cavity. Fourth, after validation against the available experimental data, molecular dynamics simulations can be used to characterize the dynamics and exchange of water molecules and ions, providing atomic level and time-dependent information otherwise inaccessible to experiments.

Functional role of the conserved active site proline of triosephosphate isomerase

The importance of the fully conserved active site proline, Pro168, for the reaction mechanism of triosephosphate isomerase (TIM) has been investigated by studying the enzymatic and crystallographic properties of the P168A variant of trypanosomal TIM. In TIM, Pro168 follows the key catalytic residue Glu167, situated at the beginning of the flexible active site loop (loop 6). Turnover numbers of the P168A variant for its substrates are reduced approximately 50-fold, whereas the Km values are approximately 2 times lower. The affinity of the P168A variant for the transition state analogue 2-phosphoglycolate (2PG) is reduced 5-fold. The crystal structures of unliganded and liganded (2PG) P168A show that the phosphate moiety of 2PG is bound similarly as in wild-type TIM, whereas the interactions of the carboxylic acid moiety with the side chain of the catalytic Glu167 differ. The unique properties of the proline side chain at position 168 are required to transmit ligand binding to the conformational change of Glu167: the side chain of Glu167 flips from the inactive swung-out to the active swung-in conformation on ligand binding in wild-type TIM, whereas in the mutant this conformational change does not occur. Further structural comparisons show that in the wild-type enzyme the concerted movement of loop 6 and loop 7 from unliganded-open to liganded-closed appears to be facilitated by the interactions of the phosphate moiety with loop 7. Apparently, the rotation of 90 degrees of the Gly211-Gly212 peptide plane of loop 7 plays a key role in this concerted movement.

Disulfide bridges in the mesophilic triosephosphate isomerase from Giardia lamblia are related to oligomerization and activity

Triosephosphate isomerase from the mesophile Giardia lamblia (GlTIM) is the only known TIM with natural disulfide bridges. We previously found that oxidized and reduced thiol states of GlTIM are involved in the interconversion between native dimers and higher oligomeric species, and in the regulation of enzymatic activity. Here, we found that trophozoites and cysts have different oligomeric species of GlTIM and complexes of GlTIM with other proteins. Our data indicate that the internal milieu of G. lamblia is favorable for the formation of disulfide bonds. Enzyme mutants of the three most solvent exposed Cys of GlTIM (C202A, C222A, and C228A) were prepared to ascertain their contribution to oligomerization and activity. The data show that the establishment of a disulfide bridge between two C202 of two dimeric GlTIMs accounts for multimerization. In addition, we found that the establishment of an intramonomeric disulfide bond between C222 and C228 abolishes catalysis. Multimerization and inactivation are both reversed by reducing conditions. The 3D structure of the C202A GlTIM was solved at 2.1 A resolution, showing that the environment of the C202 is prone to hydrophobic interactions. Molecular dynamics of an in silico model of GlTIM when the intramonomeric disulfide bond is formed, showed that S216 is displaced 4.6 A from its original position, causing loss of hydrogen bonds with residues of the active-site loop. This suggests that this change perturb the conformational state that aligns the catalytic center with the substrate, inducing enzyme inactivation.

Diabetes mellitus due to misfolding of a beta-cell transcription factor: stereospecific frustration of a Schellman motif in HNF-1alpha

Maturity-onset diabetes of the young (MODY3), a monogenic form of type II diabetes mellitus, results most commonly from mutations in hepatocyte nuclear factor 1alpha (HNF-1alpha). Diabetes-associated mutation G20R perturbs the dimerization domain of HNF-1alpha, an intertwined four-helix bundle. In the wild-type structure G20 participates in a Schellman motif to cap an alpha-helix; its dihedral angles lie in the right side of the Ramachandran plot (alpha(L) region; phi 97 degrees). Substitutions G20R and G20A lead to dimeric molten globules of low stability, suggesting that the impaired function of the diabetes-associated transcription factor is due in large part to a main-chain perturbation rather than to specific features of the Arg side-chain. This hypothesis is supported by the enhanced stability of non-standard analogues containing D-Ala or D-Ser at position 20. The crystal structure of the D-Ala20 analogue, determined to a resolution of 1.4 A, is essentially identical to the wild-type structure in the same crystal form. The mean root-mean-square deviation between equivalent C(alpha) atoms (residues 5-28) is 0.3 A; (phi, psi) angles of D-Ala20 are the same as those of G20 in the wild-type structure. Whereas the side-chain of A20 or R20 would be expected to clash with the preceding carbonyl oxygen (thus accounting for its frustrated energy landscape), the side-chain of D-Ala20 projects into solvent without perturbation of the Schellman motif. Calorimetric studies indicate that the increased stability of the D-Ala20 analogue (DeltaDeltaG(u) 1.5 kcal/mol) is entropic in origin, consistent with a conformational bias toward native-like conformations in the unfolded state. Studies of multiple substitutions at G20 and neighboring positions highlight the essential contributions of a glycine-specific tight turn and adjoining inter-subunit side-chain hydrogen bonds to the stability and architectural specificity of the intertwined dimer. Comparison of L- and D amino acid substitutions thus provides an example of the stereospecific control of an energy landscape by a helix-capping residue.

Long-lived conformational isomerism of protein dimers: the role of the free energy of subunit association

The association of protein subunits to form N-mers (N >or= 3) does not follow the dependence on the law of mass action predicted by the classical thermodynamic description used for the equilibrium of association of small molecules. For those anomalous cases, a so-called deterministic model has been previously proposed. The latter model was based on the empirical observation that the dynamics of subunit exchange between protein oligomers can be very slow, leading to the existence of long-lived conformational isomers and to a persistently heterogeneous ensemble of oligomers in solution. Contrary to the expectation for a protein dimer, we have recently shown that the subunit association of triosephosphate isomerase (TIM) could also be described as a deterministic process and that long-lived conformational isomers of TIM could be isolated in solution. Here we show that a), observation of hysteresis in pressure dissociation curves is an additional indicator of deterministic behavior; b), the extent of deviation from the classical thermodynamic behavior correlates with the free-energy change of subunit association; and c), experimental manipulation of the free energy of subunit association through the addition of a subdenaturing concentration of a chaotropic agent restores the concentration dependence of subunit association of TIM. A model that explains these features and its biological relevance is discussed.

Gated binding of ligands to HIV-1 protease: Brownian dynamics simulations in a coarse-grained model

The internal motions of proteins may serve as a "gate" in some systems, which controls ligand-protein association. This study applies Brownian dynamics simulations in a coarse-grained model to study the gated association rate constants of HIV-1 proteases and drugs. The computed gated association rate constants of three protease mutants, G48V/V82A/I84V/L90M, G48V, and L90M with three drugs, amprenavir, indinavir, and saquinavir, yield good agreements with experiments. The work shows that the flap dynamics leads to "slow gating". The simulations suggest that the flap flexibility and the opening frequency of the wild-type, the G48V and L90M mutants are similar, but the flaps of the variant G48V/V82A/I84V/L90M open less frequently, resulting in a lower gated rate constant. The developed methodology is fast and provides an efficient way to predict the gated association rate constants for various protease mutants and ligands.

Coupling of structural fluctuations to deamidation reaction in triosephosphate isomerase by Gaussian network model

We study the structural fluctuations of triosephosphate isomerase (TIM) by an elastic model, namely, the Gaussian network model (GNM), to identify a network of coupled motions in the allosteric communication between its deamidation and catalytic sites, and the promoting motions for the deamidation activity. For this, three TIM structures have been studied: one crystal structure and two model structures designed to describe different putative models for the deamidation reaction taking place at the subunit interface. The structural fluctuations have been mapped on the functional properties; then the differences in the fluctuations between the two models in relation to the deamidation reaction have been considered. The results demonstrate that the qualitative picture of the mean-square fluctuations and the correlations between the fluctuations are similar in both, but the differences may affect the observed barrier height of the deamidation reaction. The higher packing density at regions close to deamidation sites, reflected by the high-frequency fluctuating residues in the respective regions, the stronger positive correlation between the fluctuations of the deamidation sites, and enhanced positive correlation of the primary deamidation site with the extended vicinity of the catalytic region on the juxtaposed unit promote the probability of the deamidation reaction. The results in general emphasize the importance of structural fluctuations in enzyme reactions, as well as proposing the present methodology as a plausible approach for studies on the network of coupled promoting motions in protein functions.

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.

Structures of unliganded and inhibitor complexes of W168F, a Loop6 hinge mutant of Plasmodium falciparum triosephosphate isomerase: observation of an intermediate position of loop6

The enzymatic reaction of triosephosphate isomerase (TIM) is controlled by the movement of a loop (loop6, residues 166-176). Crystal structures of TIMs from a variety of sources have revealed that the loop6, which is in an open conformation in the unliganded enzyme, adopts a closed conformation in inhibitor complexes. In contrast, structures with loop open conformation are obtained in most of the complexes of TIM from the malarial parasite Plasmodium falciparum (PfTIM). W168 is a conserved N-terminal hinge residue, involved in different sets of interactions in the "open" and "closed" forms of loop6. The role of W168 in determining the loop conformation was examined by structural studies on the mutant W168F and its complexes with ligands. The three-dimensional structures of unliganded mutant (1.8 A) and complexes with sulfate (2.8 A) and glycerol-2-phosphate (G2P) (2.8 A) have been determined. Loop6 was found disordered in these structures, reflecting the importance of W168 in stabilizing either the open or the closed states. Critical sequence differences between the Plasmodium enzyme and other TIMs may influence the equilibrium between the closed and open forms. Examination of the environment of the loop6 shows that its propensity for the open or the closed forms is influenced not only by Phe96 as suggested earlier, but also by Asn233, which occurs in the vicinity of the active site. This residue is Gly in the other TIM sequences and probably plays a crucial role in the mode of ligand binding, which in turn affects the loop opening/closing process in PfTIM.