Three flavonoids targeting the beta-hydroxyacyl-acyl carrier protein dehydratase from Helicobacter pylori: crystal structure characterization with enzymatic inhibition assay

Flavonoids are the major functional components of many herbal and insect preparations and demonstrate varied pharmacological functions including antibacterial activity. Here by enzymatic assay and crystal structure analysis, we studied the inhibition of three flavonoids (quercetin, apigenin, and (S)-sakuranetin) against the beta-hydroxyacyl-acyl carrier protein dehydratase from Helicobacter pylori (HpFabZ). These three flavonoids are all competitive inhibitors against HpFabZ by either binding to the entrance of substrate tunnel B (binding model A) or plugging into the tunnel C near the catalytic residues (binding model B) mainly by hydrophobic interaction and hydrogen-bond pattern. Surrounded by hydrophobic residues of HpFabZ at both positions of models A and B, the methoxy group at C-7 of (S)-sakuranetin seems to play an important role for the inhibitor's binding to HpFabZ, partly responsible for the higher inhibitory activity of (S)-sakuranetin than those of quercetin and apigenin against HpFabZ (IC(50) in microM: (S)-sakuranetin, 2.0 +/- 0.1; quercetin: 39.3 +/- 2.7; apigenin, 11.0 +/- 2.5). Our work is expected to supply useful information for understanding the potential antibacterial mechanism of flavonoids.     

Structural basis for the functional and inhibitory mechanisms of β-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The β-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of “U”. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms.     

A new beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) from Helicobacter pylori: Molecular cloning, enzymatic characterization, and structural modeling

Helicobacter pylori is a gram-negative pathogenic bacterium that causes peptic ulcer disease and gastric cancer, and studies of the related potent enzymes associated with this bacterium are urgent for the discovery of novel drug targets. In bacteria, beta-hydroxyacyl-acyl carrier protein (ACP) dehydratase (FabZ) is a potent enzyme in fatty acid biosynthesis and catalyzes the dehydration of beta-hydroxyacyl-ACP to trans-2-acyl-ACP. In this study, the cloning and enzymatic characterization of FabZ from H. pylori strain SS1 (HpFabZ) were reported, and the gene sequence of HpfabZ was deposited in the GenBank database. Enzyme dynamic analysis showed that HpFabZ had a K(m) of 82.6+/-4.3 microM toward its substrate analog crotonoyl-CoA. Dynamic light scattering and native-PAGE investigations suggested that HpFabZ exists as hexamer in native state. Enzymatic characterization and thermal-induced unfolding analysis based on circular dichroism spectral measurements indicated that HpFabZ is very stable against high temperature (90 degrees C). Such a high stability of HpFabZ was well elucidated by the strong H-bonds and hydrophobic interactions among the HpFabZ hexamer as investigated in the modeled HpFabZ hexamer structure. Our current study is hoped to provide useful information in better understanding the FabZ of H. pylori strain and further supply possible hints in the discovery of anti-bacterial compounds using HpFabZ as target.     

The structure of (3R)-hydroxyacyl-acyl carrier protein dehydratase (FabZ) from Pseudomonas aeruginosa

Type II fatty acid biosynthesis systems are essential for membrane formation in bacteria, making the constituent proteins of this pathway attractive targets for antibacterial drug discovery. The third step in the elongation cycle of the type II fatty acid biosynthesis is catalyzed by beta-hydroxyacyl-(acyl carrier protein) (ACP) dehydratase. There are two isoforms. FabZ, which catalyzes the dehydration of (3R)-hydroxyacyl-ACP to trans-2-acyl-ACP, is a universally expressed component of the bacterial type II system. FabA, the second isoform, as has more limited distribution in nature and, in addition to dehydration, also carries out the isomerization of trans-2- to cis-3-decenoyl-ACP as an essential step in unsaturated fatty acid biosynthesis. We report the structure of FabZ from the important human pathogen Pseudomonas aeruginosa at 2.5 A of resolution. PaFabZ is a hexamer (trimer of dimers) with the His/Glu catalytic dyad located within a deep, narrow tunnel formed at the dimer interface. Site-directed mutagenesis experiments showed that the obvious differences in the active site residues that distinguish the FabA and FabZ subfamilies of dehydratases do not account for the unique ability of FabA to catalyze isomerization. Because the catalytic machinery of the two enzymes is practically indistinguishable, the structural differences observed in the shape of the substrate binding channels of FabA and FabZ lead us to hypothesize that the different shapes of the tunnels control the conformation and positioning of the bound substrate, allowing FabA, but not FabZ, to catalyze the isomerization reaction.     

Crystal structure of dimeric FabZ of Plasmodium falciparum reveals conformational switching to active hexamers by peptide flips

The crystal structure of beta-hydroxyacyl acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) has been determined at a resolution of 2.4 A. PfFabZ has been found to exist as a homodimer (d-PfFabZ) in the crystals of the present study in contrast to the reported hexameric form (h-PfFabZ) which is a trimer of dimers crystallized in a different condition. The catalytic sites of this enzyme are located in deep narrow tunnel-shaped pockets formed at the dimer interface. A histidine residue from one subunit of the dimer and a glutamate residue from the other subunit lining the tunnel form the catalytic dyad in the reported crystal structures. While the position of glutamate remains unaltered in the crystal structure of d-PfFabZ compared to that in h-PfFabZ, the histidine residue takes up an entirely different conformation and moves away from the tunnel leading to a His-Phe cis-trans peptide flip at the histidine residue. In addition, a loop in the vicinity has been observed to undergo a similar flip at a Tyr-Pro peptide bond. These alterations not only prevent the formation of a hexamer but also distort the active site geometry resulting in a dimeric form of FabZ that is incapable of substrate binding. The dimeric state and an altered catalytic site architecture make d-PfFabZ distinctly different from the FabZ structures described so far. Dynamic light scattering and size exclusion chromatographic studies clearly indicate a pH-related switching of the dimers to active hexamers.     

Conformational Changes in the tryptophan synthase from a hyperthermophile upon alpha2beta2 complex formation: crystal structure of the complex

The three-dimensional structure of the bifunctional tryptophan synthase alpha(2)beta(2) complex from Pyrococcus furiosus was determined by crystallographic analysis. This crystal structure, with the structures of an alpha subunit monomer and a beta(2) subunit dimer that have already been reported, is the first structural set in which changes in structure that occur upon the association of the individual tryptophan synthase subunits were observed. To elucidate the structural basis of the stimulation of the enzymatic activity of each of the alpha and beta(2) subunits upon alpha(2)beta(2) complex formation, the conformational changes due to complex formation were analyzed in detail compared with the structures of the alpha monomer and beta(2) subunit dimer. The major conformational changes due to complex formation occurred in the region correlated with the catalytic function of the enzyme as follows. (1) Structural changes in the beta subunit were greater than those in the alpha subunit. (2) Large movements of A46 and L165 in the alpha subunit due to complex formation caused a more open conformation favoring the entry of the substrate at the alpha active site. (3) The major changes in the beta subunit were the broadening of a long tunnel through which the alpha subunit product (indole) is transferred to the beta active site and the opening of an entrance at the beta active site. (4) The changes in the conformations of both the alpha and beta subunits due to complex formation contributed to the stabilization of the subunit association, which is critical for the stimulation of the enzymatic activities.     

Characterization of isoflavonoids as inhibitors of β-hydroxyacyl-acyl carrier protein dehydratase (FabZ) from Moraxella catarrhalis: Kinetics,spectroscopic, thermodynamics and in silico studies

Backgroud

β-hydroxyacyl-acyl carrier protein dehydratase (FabZ) is an essential component of type II fatty acid biosynthesis (FAS II) pathway in bacteria. It performs dehydration of β-hydroxyacyl-ACP to trans-2-acyl-ACP in the elongation cycle of the FAS II pathway. FabZ is ubiquitously expressed and has uniform distribution, which makes FabZ an excellent target for developing novel drugs against pathogenic bacteria.

Crystal structure of red sea bream transglutaminase

The crystal structure of the tissue-type transglutaminase from red sea bream liver (fish-derived transglutaminase, FTG) has been determined at 2.5-A resolution using the molecular replacement method, based on the crystal structure of human blood coagulation factor XIII, which is a transglutaminase zymogen. The model contains 666 residues of a total of 695 residues, 382 water molecules, and 1 sulfate ion. FTG consists of four domains, and its overall and active site structures are similar to those of human factor XIII. However, significant structural differences are observed in both the acyl donor and acyl acceptor binding sites, which account for the difference in substrate preferences. The active site of the enzyme is inaccessible to the solvent, because the catalytic Cys-272 hydrogen-bonds to Tyr-515, which is thought to be displaced upon acyl donor binding to FTG. It is postulated that the binding of an inappropriate substrate to FTG would lead to inactivation of the enzyme because of the formation of a new disulfide bridge between Cys-272 and the adjacent Cys-333 immediately after the displacement of Tyr-515. Considering the mutational studies previously reported on the tissue-type transglutaminases, we propose that Cys-333 and Tyr-515 are important in strictly controlling the enzymatic activity of FTG.     

Crystal structure of Mycobacterium tuberculosis MenB,a key enzyme in vitamin K2 biosynthesis

Bacterial enzymes of the menaquinone (Vitamin K2) pathway are potential drug targets because they lack human homologs. MenB, 1,4-dihydroxy-2-naphthoyl-CoA synthase, the fourth enzyme in the biosynthetic pathway leading from chorismate to menaquinone, catalyzes the conversion of O-succinylbenzoyl-CoA (OSB-CoA) to 1,4-dihydroxy-2-naphthoyl-CoA (DHNA-CoA). Based on our interest in developing novel tuberculosis chemotherapeutics, we have solved the structures of MenB from Mycobacterium tuberculosis and its complex with acetoacetyl-coenzyme A at 1.8 and 2.3 A resolution, respectively. Like other members of the crotonase superfamily, MenB folds as an (alpha3)2 hexamer, but its fold is distinct in that the C terminus crosses the trimer-trimer interface, forming a flexible part of the active site within the opposing trimer. The highly conserved active site of MenB contains a deep pocket lined by Asp-192, Tyr-287, and hydrophobic residues. Mutagenesis shows that Asp-192 and Tyr-287 are essential for enzymatic catalysis. We postulate a catalytic mechanism in which MenB enables proton transfer within the substrate to yield an oxyanion as the initial step in catalysis. Knowledge of the active site geometry and characterization of the catalytic mechanism of MenB will aid in identifying new inhibitors for this potential drug target.     

Analysis of conserved active site residues in monoamine oxidase A and B and their three-dimensional molecular modeling

Monoamine oxidase (MAO) is a key enzyme responsible for the degradation of serotonin, norepinephrine, dopamine, and phenylethylamine. It is an outer membrane mitochondrial enzyme existing in two isoforms, A and B. We have recently generated 14 site-directed mutants of human MAO A and B, and we found that four key amino acids, Lys-305, Trp-397, Tyr-407, and Tyr-444, in MAO A and their corresponding amino acids in MAO B, Lys-296, Trp-388, Tyr-398, and Tyr-435, play important roles in MAO catalytic activity. Based on the polyamine oxidase three-dimensional crystal structure, it is suggested that Lys-305, Trp-397, and Tyr-407 in MAO A and Lys-296, Trp-388, and Tyr-398 in MAO B may be involved in the non-covalent binding to FAD. Tyr-407 and Tyr-444 in MAO A (Tyr-398 and Tyr-435 in MAO B) may form an aromatic sandwich that stabilizes the substrate binding. Asp-132 in MAO A (Asp-123 in MAO B) located at the entrance of the U-shaped substrate-binding site has no effect on MAO A nor MAO B catalytic activity. The similar impact of analogous mutants in MAO A and MAO B suggests that these amino acids have the same function in both isoenzymes. Three-dimensional modeling of MAO A and B using polyamine oxidase as template suggests that the overall tertiary structure and the active sites of MAO A and B may be similar.     

Crystal Structures of a Populus tomentosa 4-Coumarate:CoA Ligase Shed Light on Its Enzymatic Mechanisms

4-Coumaric acid:CoA ligase (4CL) is the central enzyme of the plant-specific phenylpropanoid pathway. It catalyzes the synthesis of hydroxycinnamate-CoA thioesters, the precursors of lignin and other important phenylpropanoids, in two-step reactions involving the formation of hydroxycinnamate-AMP anhydride and then the nucleophilic substitution of AMP by CoA. In this study, we determined the crystal structures of Populus tomentosa 4CL1 in the unmodified (apo) form and in forms complexed with AMP and adenosine 5′-(3-(4-hydroxyphenyl)propyl)phosphate (APP), an intermediate analog, at 2.4, 2.5, and 1.9 Å resolution, respectively. 4CL1 consists of two globular domains connected by a flexible linker region. The larger N-domain contains a substrate binding pocket, while the C-domain contains catalytic residues. Upon binding of APP, the C-domain rotates 81° relative to the N-domain. The crystal structure of 4CL1-APP reveals its substrate binding pocket. We identified residues essential for catalytic activities (Lys-438, Gln-443, and Lys-523) and substrate binding (Tyr-236, Gly-306, Gly-331, Pro-337, and Val-338) based on their crystal structures and by means of mutagenesis and enzymatic activity studies. We also demonstrated that the size of the binding pocket is the most important factor in determining the substrate specificities of 4CL1. These findings shed light on the enzymatic mechanisms of 4CLs and provide a solid foundation for the bioengineering of these enzymes.     

Crystal structures of c-Src reveal features of its autoinhibitory mechanism.

Src family kinases are maintained in an assembled, inactive conformation by intramolecular interactions of their SH2 and SH3 domains. Full catalytic activity requires release of these restraints as well as phosphorylation of Tyr-416 in the activation loop. In previous structures of inactive Src kinases, Tyr-416 and flanking residues are disordered. We report here four additional c-Src structures in which this segment adopts an ordered but inhibitory conformation. The ordered activation loop forms an alpha helix that stabilizes the inactive conformation of the kinase domain, blocks the peptide substrate-binding site, and prevents Tyr-416 phosphorylation. Disassembly of the regulatory domains, induced by SH2 or SH3 ligands, or by dephosphorylation of Tyr-527, could lead to exposure and phosphorylation of Tyr-416.     

Biological and structural characterization of Trypanosoma cruzi phosphodiesterase C and Implications for design of parasite selective inhibitors

Trypanosoma cruzi phosphodiesterase C (TcrPDEC) is a potential new drug target for the treatment of Chagas disease but has not been well studied. This study reports the enzymatic properties of various kinetoplastid PDECs and the crystal structures of the unliganded TcrPDEC1 catalytic domain and its complex with an inhibitor. Mutations of PDEC during the course of evolution led to inactivation of PDEC in Trypanosoma brucei/Trypanosoma evansi/Trypanosoma congolense, whereas the enzyme is active in all other kinetoplastids. The TcrPDEC1 catalytic domain hydrolyzes both cAMP and cGMP with a K(m) of 23.8 μm and a k(cat) of 31 s(-1) for cAMP and a K(m) of 99.1 μm and a k(cat) of 17 s(-1) for cGMP, thus confirming its dual specificity. The crystal structures show that the N-terminal fragment wraps around the TcrPDEC catalytic domain and may thus regulate its enzymatic activity via direct interactions with the active site residues. A PDE5 selective inhibitor that has an IC(50) of 230 nm for TcrPDEC1 binds to TcrPDEC1 in an orientation opposite to that of sildenafil. This observation, together with the screen of the inhibitory potency of human PDE inhibitors against TcrPDEC, implies that the scaffold of some human PDE inhibitors might be used as the starting model for design of parasite PDE inhibitors. The structural study also identified a unique parasite pocket that neighbors the active site and may thus be valuable for the design of parasite-specific inhibitors.     

Crystallographic, kinetic, and spectroscopic study of the first ligninolytic peroxidase presenting a catalytic tyrosine

Trametes cervina lignin peroxidase (LiP) is a unique enzyme lacking the catalytic tryptophan strictly conserved in all other LiPs and versatile peroxidases (more than 30 sequences available). Recombinant T. cervina LiP and site-directed variants were investigated by crystallographic, kinetic, and spectroscopic techniques. The crystal structure shows three substrate oxidation site candidates involving His-170, Asp-146, and Tyr-181. Steady-state kinetics for oxidation of veratryl alcohol (the typical LiP substrate) by variants at the above three residues reveals a crucial role of Tyr-181 in LiP activity. Moreover, assays with ferrocytochrome c show that its ability to oxidize large molecules (a requisite property for oxidation of the lignin polymer) originates in Tyr-181. This residue is also involved in the oxidation of 1,4-dimethoxybenzene, a reaction initiated by the one-electron abstraction with formation of substrate cation radical, as described for the well known Phanerochaete chrysosporium LiP. Detailed spectroscopic and kinetic investigations, including low temperature EPR, show that the porphyrin radical in the two-electron activated T. cervina LiP is unstable and rapidly receives one electron from Tyr-181, forming a catalytic protein radical, which is identified as an H-bonded neutral tyrosyl radical. The crystal structure reveals a partially exposed location of Tyr-181, compatible with its catalytic role, and several neighbor residues probably contributing to catalysis: (i) by enabling substrate recognition by aromatic interactions; (ii) by acting as proton acceptor/donor from Tyr-181 or H-bonding the radical form; and (iii) by providing the acidic environment that would facilitate oxidation. This is the first structure-function study of the only ligninolytic peroxidase described to date that has a catalytic tyrosine.     

Molecular basis for the binding of competitive inhibitors of maize polyamine oxidase

Maize polyamine oxidase (MPAO), the only member of the polyamine oxidase (PAO) family whose three-dimensional structure is known, is characterized by a 30 A long U-shaped catalytic tunnel located between the substrate binding domain and the FAD. To shed light on the MPAO ligand binding mode, we studied the inhibition properties of linear diamines, agmatine, prenylagmatine (G3), G3 analogues, and guazatine, and analyzed the structural determinants of their biological activity. Linear diamines competitively inhibited MPAO, with the inhibitory activity increasing as a function of the number of methylene groups. With regard to the guanidino competitive inhibitors, including agmatine, G3, and G3 analogues, the presence of a hydrophobic substituent constitutes the principal factor influencing MPAO inhibition, as the addition of a hydrophobic substituent to the guanidino group of both G3 and G3 analogues greatly increases the inhibitory activity. Moreover, results obtained by a molecular modeling procedure indicated that in their preferred orientation, G3 analogues point the ammonium group toward the narrow entrance of the tunnel, while the terminal hydrophobic group is located within the large entrance. The high binding affinity for MPAO exhibited by G3 and G3 analogues bearing a prenyl group as a substituent on the guanidino moiety is in agreement with the observation that the prenyl group binds in a well-defined hydrophobic pocket, mainly formed by aromatic residues. Finally, docking simulations performed with the charged and uncharged forms of MPAO inhibitors indicate that the stereoelectronic properties of the MPAO active site are consistent with the binding of inhibitors in the protonated form.     

Catalytic features of the botulinum neurotoxin A light chain revealed by high resolution structure of an inhibitory peptide complex

The Clostridium botulinum neurotoxin serotype A light chain (BoNT/A-LC) is a Zn(II)-dependent metalloprotease that blocks the release of acetylcholine at the neuromuscular junction by cleaving SNAP-25, one of the SNARE proteins required for exocytosis. Because of the potential for use of the toxin in bioterrorism and the increasingly widespread application of the toxin in the medical field, there is significant interest in the development of small-molecule inhibitors of the metalloprotease. Efforts to design such inhibitors have not benefited from knowledge of how peptides bind to the active site since the enzyme-peptide structures available previously either were not occupied in the vicinity of the catalytic Zn(II) ion or did not represent the product of SNAP-25 substrate cleavage. Herein we report the 1.4 A-resolution X-ray crystal structure of a complex between the BoNT/A-LC and the inhibitory peptide N-Ac-CRATKML, the first structure of the light chain with an inhibitory peptide bound at the catalytic Zn(II) ion. The peptide is bound with the Cys S gamma atom coordinating the metal ion. Surprisingly, the cysteine sulfur is oxidized to the sulfenic acid form. Given the unstable nature of this species in solution, is it likely that oxidation occurs on the enzyme. In addition to the peptide-bound structure, we report two structures of the unliganded light chain with and without the Zn(II) cofactor bound at 1.25 and 1.20 A resolution, respectively. The two structures are nearly identical, confirming that the Zn(II) ion plays a purely catalytic role. Additionally, the structure of the Zn(II)-bound uncomplexed enzyme allows identification of the catalytic water molecule and a second water molecule that occupies the same position as the peptidic oxygen in the tetrahedral intermediate. This observation suggests that the enzyme active site is prearranged to stabilize the tetrahedral intermediate of the protease reaction.     

Dimerization Interface of 3-Hydroxyacyl-CoA Dehydrogenase Tunes the Formation of Its Catalytic Intermediate

3-hydroxyacyl-CoA dehydrogenase (HAD, EC 1.1.1.35) is a homodimeric enzyme localized in the mitochondrial matrix, which catalyzes the third step in fatty acid β-oxidation. The crystal structures of human HAD and subsequent complexes with cofactor/substrate enabled better understanding of HAD catalytic mechanism. However, numerous human diseases were found related to mutations at HAD dimerization interface that is away from the catalytic pocket. The role of HAD dimerization in its catalytic activity needs to be elucidated. Here, we solved the crystal structure of Caenorhabditis elegans HAD (cHAD) that is highly conserved to human HAD. Even though the cHAD mutants (R204A, Y209A and R204A/Y209A) with attenuated interactions on the dimerization interface still maintain a dimerization form, their enzymatic activities significantly decrease compared to that of the wild type. Such reduced activities are in consistency with the reduced ratios of the catalytic intermediate formation. Further molecular dynamics simulations results reveal that the alteration of the dimerization interface will increase the fluctuation of a distal region (a.a. 60–80) that plays an important role in the substrate binding. The increased fluctuation decreases the stability of the catalytic intermediate formation, and therefore the enzymatic activity is attenuated. Our study reveals the molecular mechanism about the essential role of the HAD dimerization interface in its catalytic activity via allosteric effects.     

Domain swapping between Enterococcus faecalis FabN and FabZ proteins localizes the structural determinants for isomerase activity

Anaerobic unsaturated fatty acid synthesis in bacteria occurs through the introduction of a double bond into the growing acyl chain. In the Escherichia coli model system, FabA catalyzes both the dehydration of beta-hydroxydecanoyl-ACP and the isomerization of trans-2-decenoyl-ACP to cis-3-decenoyl-ACP as the essential step. A second dehydratase, FabZ, functions in acyl chain elongation but cannot carry out the isomerization reaction. Enterococcus faecalis has two highly related FabZ homologs. One of these, termed EfFabN, carries out the isomerization reaction in vivo, whereas the other, EfFabZ, does not (Wang, H., and Cronan, J. E. (2004) J. Biol. Chem. 279, 34489-34495). We carried out a series of domain swapping and mutagenesis experiments coupled with in vitro biochemical analyses to define the structural feature(s) that specify the catalytic properties of these two enzymes. Substitution of the beta3 and beta4 strands of EfFabZ with the corresponding strands from EfFabN was necessary and sufficient to convert EfFabZ into an isomerase. These data are consistent with the hypothesis that the isomerase potential of beta-hydroxyacyl-ACP dehydratases is determined by the properties of the beta-sheets that dictate the orientation of the central alpha-helix and thus the shape of the substrate binding tunnel rather than the catalytic machinery at the active site.     

How an Inhibitor Bound to Subunit Interface Alters Triosephosphate Isomerase Dynamics

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