Search

HPPK (6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase) is a monomeric protein with 158 residues, which undergoes large-scale conformational changes between apo, open, and holo states responding to ligand binding for its function. It has been explored widely as an excellent target for potential antibacterial drug development. However, little is known about how conformational dynamics between the native states influences the substrate recognition and the functionality of enzymatic catalysis. Here, we report a coarse-grained triple-basin structure-based model upon ligand binding to describe such multiple-state system by the molecular dynamics simulation. With our model, we have made theoretical predictions that are in good agreement with the experimental measurements. Our results revealed the intrinsic conformational fluctuations between apo and open states without ligand binding. We found that HPPK can switch to the activated holo state upon the ordered binding of the two ligands (ATP and HP). We uncovered the underlying mechanism by which major induced fit and minor population shift pathways coexist upon ligand binding by quantitative flux analysis. Additionally, we pointed out the structural origin for the conformational changes and identified the key residues as well as contact interactions. We further explored the temperature effect on the conformational distributions and pathway weights. It gave strong support that higher temperatures promote population shift, while the induced fit pathway is always the predominant activation route of the HPPK system. These findings will provide significant insights of the mechanisms of the multistate conformational dynamics of HPPK upon ligand binding.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that, in its binary form, even bound by analogues of MgATP, loops 2 and 3 remain rather flexible; this raises questions about the putative large-scale induced-fit conformational change of the HPPK-MgATP binary complex. In this work, long-time all-atomic molecular dynamics simulations were conducted to investigate the loop dynamics in this complex. Our simulations show that, with loop 3 closed, multiple conformations of loop 2, including the open, semiopen, and closed forms, are all accessible to the binary complex. These results provide valuable structural insights into the details of conformational changes upon 6-hydroxymethyl-7,8-dihydropterin (HP) binding and biological activities of HPPK. Conformational network analysis and principal component analysis related to the loops are also discussed.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a member of the folate biosynthesis pathway found in prokaryotes and lower eukaryotes that catalyzes the pyrophosphoryl transfer from the ATP cofactor to a 6-hydroxymethyl-7,8-dihydropterin substrate. We report the chemical synthesis of a series of S-functionalized 8-mercaptoguanine (8MG) analogues as substrate site inhibitors of HPPK and quantify binding against the E. coli and S. aureus enzymes (EcHPPK and SaHPPK). The results demonstrate that analogues incorporating acetophenone-based substituents have comparable affinities for both enzymes. Preferential binding of benzyl-substituted 8MG derivatives to SaHPPK was reconciled when a cryptic pocket unique to SaHPPK was revealed by X-ray crystallography. Differential chemical shift perturbation analysis confirmed this to be a common mode of binding for this series to SaHPPK. One compound (41) displayed binding affinities of 120 nM and 1.76 μM for SaHPPK and EcHPPK, respectively, and represents a lead for the development of more potent and selective inhibitors of SaHPPK.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the first reaction in the folate biosynthetic pathway. Comparison of its X-ray and nuclear magnetic resonance structures suggests that the enzyme undergoes significant conformational change upon binding to its substrates, especially in three catalytic loops. Experimental research has shown that even when confined by crystal contacts, loops 2 and 3 remain rather flexible when the enzyme is in its apo form, raising questions about the putative large-scale induced-fit conformational change of HPPK. To investigate the loop dynamics in a crystal-free environment, we performed conventional molecular dynamics simulations of the apo-enzyme at two different temperatures (300 and 350 K). Our simulations show that the crystallographic B-factors considerably underestimate the loop dynamics; multiple conformations of loops 2 and 3, including the open, semi-open, and closed conformations that an enzyme must adopt throughout its catalytic cycle, are all accessible to the apo-enzyme. These results revise our previous view of the functional mechanism of conformational change upon MgATP binding and offer valuable structural insights into the workings of HPPK. In this paper, conformational network analysis and principal component analysis related to the loops are discussed to support the presented conclusions.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), which follows an ordered bi-bi kinetic mechanism with ATP binding to the enzyme first. HPPK undergoes dramatic conformational changes during its catalytic cycle as revealed by X-ray crystallography, and the conformational changes are essential for the enzymatic catalysis as shown by site-directed mutagenesis and biochemical and crystallographic analysis of the mutants. However, the dynamic properties of the enzyme have not been measured experimentally. Here, we report a (15)N NMR relaxation study of the dynamic properties of Escherichia coli HPPK from the apo form to the binary substrate complex with MgATP (represented by MgAMPCPP, an ATP analogue) to the Michaelis complex (ternary substrate complex) with MgATP (represented by MgAMPCPP) and HP (represented by 7,7-dimethyl-6-hydroxypterin, an HP analogue). The results show that the binding of the nucleotide to HPPK does not cause major changes in the dynamic properties of the enzyme. Whereas enzymes are often more rigid when bound to the ligand or the substrate, the internal mobility of HPPK is not reduced and is even moderately increased in the binary complex, particularly in the catalytic loops. The internal mobility of the catalytic loops is significantly quenched upon the formation of the ternary complex, but some mobility remains. The enhanced motions in the catalytic loops of the binary substrate complex may be required for the assembling of the ternary complex. On the other hand, some degrees of mobility in the catalytic loops of the ternary complex may be required for the optimal stabilization of the transition state, which may need the instantaneous adjustment and alignment of the side-chain positions of catalytic residues. Such dynamic behaviors may be characteristic of bisubstrate enzymes.

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) catalyzes the transfer of pyrophosphate from ATP to 6-hydroxymethyl-7,8-dihydropterin (HP), leading to the biosynthesis of folate cofactors. HPPK undergoes dramatic conformational changes during its catalytic cycle, and the conformational changes are essential for enzymatic catalysis. Thus, the enzyme is not only an attractive target for developing antimicrobial agents but also an excellent model system for studying the catalytic mechanism of enzymatic pyrophosphoryl transfer as well as the role of protein dynamics in enzymatic catalysis. In the present study, we report the NMR solution structures of the binary complex HPPK*MgAMPCPP and the ternary complex HPPK*MgAMPCPP*DMHP, where alpha,beta-methyleneadenosine triphosphate (AMPCPP) and 7, 7-dimethyl-6-hydroxypterin (DMHP) are the analogues of the substrates ATP and HP, respectively. The results suggest that the three catalytic loops of the binary complex of HPPK can assume multiple conformations in slow exchanges as evidenced by multiple sets of NMR signals for several residues in loops 2 and 3 and the very weak or missing NH cross-peaks for several residues in loops 1 and 3. However, the ternary complex shows only one set of NMR signals, and the cross-peak intensities are rather uniform, suggesting that the binding of the second substrate shifts the multiple conformations of the binary complex to an apparently single conformation of the ternary complex. The NMR behaviors and conformations of the binary complex HPPK*MgAMPCPP are significantly different from those of HPPK in complex with Mgbeta,gamma-methyleneadenosine triphosphate (MgAMPPCP). It is suggested that the conformational properties of the binary substrate complex HPPK*MgATP be represented by those of HPPK*MgAMPCPP, because MgAMPCPP is a better MgATP analogue for HPPK with respect to both binding affinity and bound conformation.

Deletion mutagenesis, biochemical, and X-ray crystallographic studies have shown that loop 3 of Escherichia coli 6-hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is required for the assembly of the active center, plays an important role in the stabilization of the ternary complex of HPPK with MgATP and 6-hydroxymethyl-7,8-dihydropterin (HP), and is essential for catalysis. Whether the critical functional importance of loop 3 is due to the interactions between residues R84 and W89 and the two substrates has been addressed by site-directed mutagenesis, biochemical, and X-ray crystallographic studies. Substitution of R84 with alanine causes little changes in the dissociation constants and kinetic constants of the HPPK-catalyzed reaction, indicating that R84 is not important for either substrate binding or catalysis. Substitution of W89 with alanine increases the K(d) for the binding of MgATP by a factor of 3, whereas the K(d) for HP increases by a factor of 6, which is due to the increase in the dissociation rate constant. The W89A mutation decreases the rate constant for the chemical step of the forward reaction by a factor of 15 and the rate constant for the chemical step of the reverse reaction by a factor of 25. The biochemical results of the W89A mutation indicate that W89 contributes somewhat to the binding of HP and more significantly to the chemical step. The crystal structures of W89A show that W89A has different conformations in loops 2 and 3, but the critical catalytic residues are positioned for catalysis. When these results are taken together, they suggest that the critical functional importance of loop 3 is not due to the interactions of the R84 guanidinium group or the W89 indole ring with the substrates.