Tyr74 is essential for the formation, stability and function of Plasmodium falciparum triosephosphate isomerase dimer

Plasmodium falciparum triosephosphate isomerase (PfTIM) is known to be functional only as a homodimer. Although many studies have shown that the interface Cys13 plays a major role in the stability of the dimer, a few reports have demonstrated that structurally conserved Tyr74 may be essential for the stability of PfTIM dimer. To understand the role of Tyr74, we have performed molecular dynamics (MD) simulations of monomeric and dimeric PfTIM mutated to glycine and cysteine at position 74. Simulations of the monomer revealed that mutant Tyr74Gly does not produce changes in folding and stability of the monomer. Interestingly, comparison of the flexibility of Tyr74 in the monomer and dimer revealed that this residue possesses an intrinsic restricted mobility, indicating that Tyr74 is an anchor residue required for homodimerization. Tyr74 also appears to play an important role in binding by facilitating the disorder-to-order transitions of loops 1 and 3, which allows Cys13 to form favorable interactions with loop 3 and Lys12 to be locked in a favorable position for catalysis. High-temperature MD simulations of the wild-type and Tyr74Gly PfTIM dimers showed that the aromatic moiety of Tyr74 is necessary to preserve the geometry and native contacts between loops 1 and 3 at the interface of the dimer. Disulfide cross-linking between mutant Tyr74Cys and Cys13 further revealed that Tyr74 stabilizes the geometry of loop 1 (which contains the catalytic residue Lys12) and the interactions between loops 1 and 3 via aromatic-aromatic interactions with residues Phe69, Tyr101, and Phe102. Principal component analysis showed that Tyr74 is also necessary to preserve the collective motions in the dimer that contribute to the catalytic efficiency of PfTIM dimer. We conclude that Tyr74 not only plays a role in the stability of the dimer, but also participates in the dimerization process and collective motions via coupled disorder-to-order transitions of intrinsically disordered regions, necessary for efficiency in the catalytic function of PfTIM.     

Highly specific inactivation of triosephosphate isomerase from Trypanosoma cruzi

We searched for molecules that selectively inactivate homodimeric triosephosphate isomerase from Trypanosoma cruzi (TcTIM), the parasite that causes Chagas' disease. We found that some benzothiazoles inactivate the enzyme. The most potent were 3-(2-benzothiazolylthio)-propanesulfonic acid, 2-(p-aminophenyl)-6-methylbenzothiazole-7-sulfonic acid, and 2-(2-4(4-aminophenyl)benzothiazole-6-methylbenzothiazole-7-sulfonic acid. Half-maximal inactivation by these compounds was attained with 33, 56, and 8 microM, respectively; in human TIM, half-maximal inactivation required 422 microM, 3.3 mM, and 1.6 mM. In TcTIM, the effect of the benzothiazoles decreased as the concentration of the enzyme was increased. TcTIM has a cysteine (Cys 15) at the dimer interface, whereas human TIM has methionine in that position. In M15C human TIM, the benzothiazole concentrations that caused half-maximal inactivation were much lower than in the wild type. The overall findings suggest that the benzothiazoles perturb the interactions between the two subunits of TcTIM through a process in which the interface cysteine is central in their deleterious action.     

Crystal structure of glucose-6-phosphate isomerase from Thermus thermophilus HB8 showing a snapshot of active dimeric state

Glucose-6-phosphate isomerase (GPI) is a glycolytic enzyme with ill-defined oligomeric state. In order to obtain insight into the correlation between oligomerization and the catalytic function of this enzyme, the crystal structure of GPI from the extreme thermophile Thermus thermophilus HB8 (TtGPI) has been determined at 1.95 Å resolution. The crystallographic asymmetric unit contains an apparent dimer. The core fold of protomer and the interprotomer spatial arrangement of the dimer are similar to those of already reported crystal structures of other GPIs. The active site is located on the dimer interface, and putative catalytic residues are well conserved among the GPIs. These results suggest that the observed dimeric state of TtGPI in the crystal is biologically relevant and that this enzyme uses a common catalytic mechanism for the isomerase reaction. Gel-filtration chromatography, chemical cross-linking, sedimentation equilibrium by analytical ultracentrifugation, and dynamic light-scattering experiments indicate that TtGPI exists in a dynamic equilibrium between monomeric and dimeric states in solution. Several factors potentially contributing to the thermal stability of TtGPI protomer were identified: (i) a decrease in denaturation entropy by the shorter polypeptide length and by amino acid composition, including the increased number of proline residues and a higher arginine-to-lysine ratio; (ii) a larger number of ion pairs; and (iii) a reduction in cavity volume. From these results, it is suggested that transient dimer formation is sufficient for the catalytic function and that the TtGPI protomer itself has intrinsically higher thermal stability.     

Structure and inactivation of triosephosphate isomerase from Entamoeba histolytica

Triosephosphate isomerase (TIM) has been proposed as a target for drug design. TIMs from several parasites have a cysteine residue at the dimer interface, whose derivatization with thiol-specific reagents induces enzyme inactivation and aggregation. TIMs lacking this residue, such as human TIM, are less affected. TIM from Entamoeba histolytica (EhTIM) has the interface cysteine residue and presents more than ten insertions when compared with the enzyme from other pathogens. To gain further insight into the role that interface residues play in the stability and reactivity of these enzymes, we determined the high-resolution structure and characterized the effect of methylmethane thiosulfonate (MMTS) on the activity and conformational properties of EhTIM. The structure of this enzyme was determined at 1.5A resolution using molecular replacement, observing that the dimer is not symmetric. EhTIM is completely inactivated by MMTS, and dissociated into stable monomers that possess considerable secondary structure. Structural and spectroscopic analysis of EhTIM and comparison with TIMs from other pathogens reveal that conformational rearrangements of the interface after dissociation, as well as intramonomeric contacts formed by the inserted residues, may contribute to the unusual stability of the derivatized EhTIM monomer.     

Structural and biochemical characterization of a recombinant triosephosphate isomerase from Rhipicephalus (Boophilus) microplus

Triosephosphate isomerase (TIM) is an enzyme with a role in glycolysis and gluconeogenesis by catalyzing the interconversion between glyceraldehyde 3-phosphate and dihydroxyacetone phosphate. This enzyme has been used as a target in endoparasite drug development. In this work we cloned, expressed, purified and studied kinetic and structural characteristics of TIM from tick embryos, Rhipicephalus (Boophilus) microplus (BmTIM). The Km and Vmax of the recombinant BmTIM with glyceraldehyde 3-phosphate as substrate, were 0.47 mM and 6031 ??mol min⁻¹ mg protein⁻¹, respectively. The resolution of the diffracted crystal was estimated to be 2.4 Å and the overall data showed that BmTIM is similar to other reported dimeric TIMs. However, we found that, in comparison to other TIMs, BmTIM has the highest content of cysteine residues (nine cysteine residues per monomer). Only two cysteines could make disulfide bonds in monomers of BmTIM. Furthermore, BmTIM was highly sensitive to the action of the thiol reagents dithionitrobenzoic acid and methyl methane thiosulfonate, suggesting that there are five cysteines exposed in each dimer and that these residues could be employed in the development of species-specific inhibitors.     

Functional Roles of the Dimer-Interface Residues in Human Ornithine Decarboxylase

Ornithine decarboxylase (ODC) catalyzes the decarboxylation of ornithine to putrescine and is the rate-limiting enzyme in the polyamine biosynthesis pathway. ODC is a dimeric enzyme, and the active sites of this enzyme reside at the dimer interface. Once the enzyme dissociates, the enzyme activity is lost. In this paper, we investigated the roles of amino acid residues at the dimer interface regarding the dimerization, protein stability and/or enzyme activity of ODC. A multiple sequence alignment of ODC and its homologous protein antizyme inhibitor revealed that 5 of 9 residues (residues 165, 277, 331, 332 and 389) are divergent, whereas 4 (134, 169, 294 and 322) are conserved. Analytical ultracentrifugation analysis suggested that some dimer-interface amino acid residues contribute to formation of the dimer of ODC and that this dimerization results from the cooperativity of these interface residues. The quaternary structure of the sextuple mutant Y331S/Y389D/R277S/D332E/V322D/D134A was changed to a monomer rather than a dimer, and the K _d value of the mutant was 52.8 µM, which is over 500-fold greater than that of the wild-type ODC (ODC_WT). In addition, most interface mutants showed low but detectable or negligible enzyme activity. Therefore, the protein stability of these interface mutants was measured by differential scanning calorimetry. These results indicate that these dimer-interface residues are important for dimer formation and, as a consequence, are critical for enzyme catalysis.     

Crystal structure of the C-terminal domain of a flagellar hook-capping protein from Xanthomonas campestris

The crystal structure of the C-terminal domain of a hook-capping protein FlgD from the plant pathogen Xanthomonas campestris (Xc) has been determined to a resolution of ca 2.5 Å using X-ray crystallography. The monomer of whole FlgD comprises 221 amino acids with a molecular mass of 22.7 kDa, but the flexible N-terminus is cleaved for up to 75 residues during crystallization. The final structure of the C-terminal domain reveals a novel hybrid comprising a tudor-like domain interdigitated with a fibronectin type III domain. The C-terminal domain of XcFlgD forms three types of dimers in the crystal. In agreement with this, analytical ultracentrifugation and gel filtration experiments reveal that they form a stable dimer in solution. From these results, we propose that the Xc flagellar hook cap protein FlgD comprises two individual domains, a flexible N-terminal domain that cannot be detected in the current study and a stable C-terminal domain that forms a stable dimer.     

Quantum chemical investigation of the intra- and intermolecular proton transfer reactions and hydrogen bonding interactions in 4-amino-5-(2-hydroxyphenyl)-2H-1,2,4-triazole-3(4H)-thione

The intramolecular thione-thiol tautomerism and intermolecular double proton transfer reaction of the hydrogen-bonded thione and thiol dimers in the title triazole compound were studied at the B3LYP level of theory using 6?311++G(d,p) basis function. The influence of the solvent on the single and double proton transfer reactions was examined in three solvents (chloroform, methanol and water) using the polarizable continuum model (PCM) approximation. The computational results show that the thione tautomer is the most stable isomer with a very high tautomeric energy barrier both in the gas phase and in solution phase, indicating a quite disfavored process. The solvent effect is found to be sizable with increasing polarity. In the double proton transfer reaction, the thione dimer is found to be more stable than thiol dimer both in the gas phase and in solution phase. The energetic and thermodynamic parameters of the double proton transfer process show that the double proton exchange from thione dimer to thiol dimer is thermodynamically unfavored. However, the exchange from thiol dimer to thione dimer for the gas phase and water phase seems to be feasible with a low barrier height and with a negative value in enthalpy and free energy changes. In addition, the hydrogen bonding interactions were analyzed in the gas phase regarding their geometries and energies. It is found that all complex formations are enthalpically favored, and the stability of the H-bonds comes in the order of S1—H2···N2 > N2—H2···S1 > N3—H3B···O1. Finally, non-linear optical properties were carried out at the same calculation level in the gas phase.

A rationally designed monomeric variant of anthranilate phosphoribosyltransferase from Sulfolobus solfataricus is as active as the dimeric wild-type enzyme but less thermostable

The anthranilate phosphoribosyltransferase from Sulfolobus solfataricus (ssAnPRT) forms a homodimer with a hydrophobic subunit interface. To elucidate the role of oligomerisation for catalytic activity and thermal stability of the enzyme, we loosened the dimer by replacing two apolar interface residues with negatively charged residues (mutations I36E and M47D). The purified double mutant I36E+M47D formed a monomer with wild-type catalytic activity but reduced thermal stability. The single mutants I36E and M47D were present in a monomer-dimer equilibrium with dissociation constants of about 1 μM and 20 μM, respectively, which were calculated from the concentration-dependence of their heat inactivation kinetics. The monomeric form of M47D, which is populated at low subunit concentrations, was as thermolabile as monomeric I36E+M47D. Likewise, the dimeric form of I36E, which was populated at high subunit concentrations, was as thermostable as dimeric wild-type ssAnPRT. These findings show that the increased stability of wild-type ssAnPRT compared to the I36E+M47D double mutant is not caused by the amino acid exchanges per se but by the higher intrinsic stability of the dimer compared to the monomer. In accordance with the negligible effect of the mutations on catalytic activity and stability, the X-ray structure of M47D contains only minor local perturbations at the dimer interface. We conclude that the monomeric double mutant resembles the individual wild-type subunits, and that ssAnPRT is a dimer for stability but not for activity reasons.     

Dimer Formation of a Stabilized Gβ1 Variant: A Structural and Energetic Analysis

In previous work, a strongly stabilized variant of the β1 domain of streptococcal protein G (Gβ1) was obtained by an in vitro selection method. This variant, termed Gβ1-M2, contains the four substitutions E15V, T16L, T18I, and N37L. Here we elucidated the molecular basis of the observed strong stabilizations. The contributions of these four residues were analyzed individually and in various combinations, additional selections with focused Gβ1 gene libraries were performed, and the crystal structure of Gβ1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gβ1 by contributions of between 1.6 and 6.0 kJ mol^− 1 (at 70 °C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gβ1. They also increase the tendency to form dimers, as shown by dependence on the concentration of apparent molecular mass in analytical ultracentrifugation, by concentration-dependent stability, and by a strongly increased van't Hoff enthalpy of unfolding. The 0.88-Å crystal structure of Gβ1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gβ1-M2 molecules via six intermolecular hydrogen bonds between the two β strands 2 and 2′ and an adjacent self-complementary hydrophobic surface area, which is created by the T16L and N37L substitutions and a large 120° rotation of the Tyr33 side chain. This removal of hydrophilic groups and the malleability of the created hydrophobic surface provide the basis for the dimer formation of stabilized Gβ1 variants.     

Thermodynamic, enzymatic and structural effects of removing a salt bridge at the base of loop 4 in (pro)caspase-3

Interactions between loops 2, 2′ and 4, known as the loop bundle, stabilize the active site of caspase-3. Loop 4 (L4) is of particular interest due to its location between the active site and the dimer interface. We have disrupted a salt bridge between K242 and E246 at the base of L4 to determine its role in overall conformational stability and in maintaining the active site environment. Stability measurements show that only the K242A single mutant decreases stability of the dimer, whereas both single mutants and the double mutant demonstrate much lower activity compared to wild-type caspase-3. Structural studies of the caspase-3 variants show the involvement of K242 in hydrophobic interactions that stabilize helix 5, near the dimer interface, and the role of E246 appears to be to neutralize the positive charge of K242 within the hydrophobic cluster. Overall, the results suggest E246 and K242 are important in procaspase-3 for their interaction with neighboring residues, not with one another. Conversely, formation of the K242–E246 salt bridge in caspase-3 is needed for an accurate, stable conformation of loop L4 and proper active site formation in the mature enzyme.     

Engineering of an intersubunit disulphide bridge in glutathione reductase from Escherichia coli

By site-directed mutagenesis, Thr-75 was converted to Cys-75 in the glutathione reductase (EC 1.6.4.2) of Escherichia coli. This led to the spontaneous formation of an intersubunit disulphide bridge across the 2-fold axis of the dimeric enzyme. The disulphide bridge had no deleterious effect on the catalytic activity, but nor did it increase the thermal stability of the enzyme, possibly because of local conformational flexibility on the dimer interface. The T75C mutant, like the wild-type enzyme, was inactivated by NADPH, proving that this inactivation cannot be due to simple dissociation of the dimer.     

Gas phase dimerization of neuropeptide head activator analogs useful for the noncovalent constraint of peptides

We have used electrospray mass spectrometry to examine the dimerization of EFLIVKS, a reversed sequence analog of part of the neuropeptide head activator, and other similar analogs. Observation of the EFLIVKS gas phase dimer is concentration-dependent, with a half-saturation concentration for relative dimer formation of 7.8 microM, similar to that of SKVILFE of 12 microM. The lowest energy conformers from quenched molecular dynamics simulations suggest EFLIVKS may dimerize in the gas phase by formation of multiple ion pairs across the dimer interface. Alanine-scan mutants also dimerize in the gas phase, with replacement of the interior residues FLIVK diminishing dimerization. The concentration-dependence of the EFLIVKS circular dichroism spectrum at pH 7.5 suggests the existence of different conformation states at different concentrations, but does not provide evidence supporting the saturable dimer formation in solution. Different analogs of EFLIVKS, when fused to each end of a 18mer unfolded peptide, induce solution structures with T(m)s of 42-50 degrees C. These peptides and analogs may thus be useful for the noncovalent constraint of peptides and peptide library members used in cellular screens.     

PCR-based methods for <Emphasis Type="Italic">CDA</Emphasis> K27Q and A70T genotyping: genotypes and alleles distribution in a central Italy population

Cytidine deaminase (CDA) is a pyrimidine salvage pathway enzyme that catalyzes the hydrolytic deamination of cytidine and deoxycytidine to their corresponding uracil nucleosides. CDA also catalyzes the inactivation of some chemotherapeutic nucleoside analogues such as cytosine arabinoside and gemcitabine. CDA 79A > C (K27Q, rs2072671) and 208G > A (A70T, rs60369023) were found to be associated either with clinical outcomes as well as with pharmacokinetics and toxicity of drugs administered to different subsets of patients. In this paper we reported two PCR-based methods for CDA 79A > C (K27Q) and 208G > A (A70T) genotyping and tested their feasibility using DNA extracted from whole blood as well as from buccal swabs. The aim of this study was also to assess the distribution of genotypic variants in a central Italy population. The allele frequencies were 56.3% (K*) and 43.7% (Q*) for K27Q and 100% (A*) and 0% (T*) for A70T. The genotype frequencies were 32.8% (K*/K*), 46.9% (K*/Q*) and 20.3% (Q*/Q*) for K27Q. The genotype frequencies did not deviate from Hardy–Weinberg equilibrium. The results were compared with those of other reported populations. They showed marked ethnic group differences.     

Long-range interactions in the dimer interface of ornithine decarboxylase are important for enzyme function.

Ornithine decarboxylase (ODC) is a pyridoxal 5'-phosphate dependent enzyme that catalyzes the first committed step in the biosynthesis of polyamines. ODC is a proven drug target for the treatment of African sleeping sickness. The enzyme is an obligate homodimer, and the two identical active sites are formed at the dimer interface. Alanine scanning mutagenesis of dimer interface residues in Trypanosoma brucei ODC was undertaken to determine the energetic contribution of these residues to subunit association. Twenty-three mutant enzymes were analyzed by analytical ultracentrifugation, and none of the mutations were found to cause a greater than 1 kcal/mol decrease in dimer stability. These data suggest that the energetics of the interaction may be distributed across the interface. Most significantly, many of the mutations had large effects (DeltaDeltaG kcat/Km > 2.5 kcal/mol) on the catalytic efficiency of the enzyme. Residues that affected activity included those in or near the substrate binding site but also a number of residues that are distant (15-20 A) from this site. These data provide evidence that long-range energetic coupling of interface residues to the active site is essential for enzyme function, even though structural changes upon ligand binding to wild-type ODC are limited to local conformational changes in the active site. The ODC dimer interface appears to be optimized for catalytic function and not for dimer stability. Thus, small molecules directed to the ODC interfaces could impact biological function without having to overcome the difficult energetic barrier of dissociating the interacting partners.     

Missense variant in TPI1 (Arg189Gln) causes neurologic deficits through structural changes in the triosephosphate isomerase catalytic site and reduced enzyme levels in vivo

Mutations in the gene triosephosphate isomerase (TPI) lead to a severe multisystem condition that is characterized by hemolytic anemia, a weakened immune system, and significant neurologic symptoms such as seizures, distal neuropathy, and intellectual disability. No effective therapy is available. Here we report a compound heterozygous patient with a novel TPI pathogenic variant (NM_000365.5:c.569G>A:p.(Arg189Gln)) in combination with the common (NM_000365.5:c.315G>C:p.(Glu104Asp)) allele. We characterized the novel variant by mutating the homologous Arg in Drosophila using a genomic engineering system, demonstrating that missense mutations at this position cause a strong loss of function. Compound heterozygote animals were generated and exhibit motor behavioural deficits and markedly reduced protein levels. Furthermore, examinations of the TPI^Arg189Gln/TPI^Glu104Asp patient fibroblasts confirmed the reduction of TPI levels, suggesting that Arg189Gln may also affect the stability of the protein. The Arg189 residue participates in two salt bridges on the backside of the TPI enzyme dimer, and we reveal that a mutation at this position alters the coordination of the substrate-binding site and important catalytic residues. Collectively, these data reveal a new human pathogenic variant associated with TPI deficiency, identify the Arg189 salt bridge as critical for organizing the catalytic site of the TPI enzyme, and demonstrates that reduced TPI levels are associated with human TPI deficiency. These findings advance our understanding of the molecular pathogenesis of the disease, and suggest new therapeutic avenues for pre-clinical trials.     

Chemical reactivity of brome mosaic virus capsid protein

Viral particles are biological machines that have evolved to package, protect, and deliver the viral genome into the host via regulated conformational changes of virions. We have developed a procedure to modify lysine residues with S-methylthioacetimidate across the pH range from 5.5 to 8.5. Lysine residues that are not completely modified are involved in tertiary or quaternary structural interactions, and their extent of modification can be quantified as a function of pH. This procedure was applied to the pH-dependent structural transitions of brome mosaic virus (BMV). As the reaction pH increases from 5.5 to 8.5, the average number of modified lysine residues in the BMV capsid protein increases from 6 to 12, correlating well with the known pH-dependent swelling behavior of BMV virions. The extent of reaction of each of the capsid protein's lysine residues has been quantified at eight pH values using coupled liquid chromatography-tandem mass spectrometry. Each lysine can be assigned to one of three structural classes identified by inspection of the BMV virion crystal structure. Several lysine residues display reactivity that indicates their involvement in dynamic interactions that are not obvious in the crystal structure. The influence of several capsid protein mutants on the pH-dependent structural transition of BMV has also been investigated. Mutant H75Q exhibits an altered swelling transition accompanying solution pH increases. The H75Q capsids show increased reactivity at lysine residues 64 and 130, residues distal from the dimer interface occupied by H75, across the entire pH range.     

Stabilizing interactions in the dimer interface of alpha-subunit in Escherichia coli RNA polymerase: a graph spectral and point mutation study

The formation of alpha(2) dimer in Escherichia coli core RNA polymerase (RNAP) is thought to be the first step toward the assembly of the functional enzyme. A large number of evidences indicate that the alpha-subunit dimerizes through its N-terminal domain (NTD). The crystal structures of the alpha-subunit NTD and that of a homologous Thermus aquaticus core RNAP are known. To identify the stabilizing interactions in the dimer interface of the alpha-NTD of E. coli RNAP, we identified side-chain clusters by using the crystal structure coordinates of E. coli alpha-NTD. A graph spectral algorithm was used to identify side-chain clusters. This algorithm considers the global nonbonded side-chain interactions of the residues for the clustering procedure and is unique in identifying residues that make the largest number of interactions among the residues that form clusters in a very quantitative way. By using this algorithm, a nine-residue cluster consisting of polar and hydrophobic residues was identified in the subunit interface adjacent to the hydrophobic core. The residues forming the cluster are relatively rigid regions of the interface, as measured by the thermal factors of the residues. Most of the cluster residues in the E. coli enzyme were topologically and sequentially conserved in the T. aquaticus RNAP crystal structure. Residues 35F and 46I were predicted to be important in the stability of the alpha-dimer interface, with 35F forming the center of the cluster. The predictions were tested by isolating single-point mutants alpha-F35A and alpha-I46S on the dimer interface, which were found to disrupt dimerization. Thus, the identified cluster at the edge of the dimer interface seems to be a vital component in stabilizing the alpha-NTD.     

NDM-1 Zn1-binding residue His116 plays critical roles in antibiotic hydrolysis

Bacteria expressing NDM-1 have been labeled as superbugs because it confers upon them resistance to a broad range of β-lactam antibiotics. The enzyme has a di?zinc active centre, with the Zn2 site extensively studied. The roles of active-site Zn1 ligand residues are, however, still not fully understood. We carried out structure-function studies using the mutants, H116A, H116N, and H116Q. Zinc content analysis showed that Zn1 binding was weakened by 40 to 60% in the H116 mutants. The enzymatic-activity studies showed that the lower hydrolysis rates were mainly caused by their weaker substrate binding. The catalytic efficiency (k_cat/K_m) of the mutants followed the order: WT > > H116Q (decreased by 4–20 fold) > H116A (decreased by 20–700 fold) ≥ H116N (decreased by 6–800 fold). The maximum effect was observed on H116N against penicillin G, whereas ampicillin was not hydrolyzed at all. The fold-increase of K_m values, which informs the weakening of substrate binding, were: H116A by 5–45 fold; H116N by 6–100 fold; H116Q by 2–10 fold. Molecular dynamics simulations suggested that the Zn1 site mutations affected the positions of Zn2 and the bridging hydroxide, by 0.8 to 1.2 Å, with the largest changes of ~1.5 Å observed on Zn2 ligand C221. A native hydrogen bond between H118 and D236 was disrupted in the H116N and H116Q mutants, which led to increased flexibility of loop 10. Consequently, residue N233 was no longer maintained at an optimal position for substrate binding. H116 connected loop 7 across Zn1 to loop 10, thereby contributed to the overall integrity. This work revealed that the H116-Zn1 interaction plays a critical role in defining the substrate-binding site. From these results, it can be inferred that inhibition strategies targeting the zinc ions may be a new direction for drug development.