Unique GMP-binding site in Mycobacterium tuberculosis guanosine monophosphate kinase

Development of an HPLC-based guanosine monophosphate kinase assay and application to Plasmodium vivax guanylate kinase

The development of a high-performance liquid chromatography (HPLC)-based method, for guanosine monophosphate kinase activity assays, is presented. The method uses the intrinsic UV absorption (at 260 nm) of substrates and products of the enzymatic reaction (GMP, ATP, ADP and GDP) to unambiguously determine percent conversion of substrate into product. It uses a commercially available C18 column which can separate reaction samples by elution under isocratic conditions in 12 min per run. The kinetics of the forward reaction catalyzed by Plasmodium vivax guanylate kinase (PvGK), a potential drug target against malaria, was determined. The relative concentrations of the two substrates (GMP and ATP) have a distinct effect on reaction velocity. Kinetic analyses showed the PvGK-catalyzed reaction to be associated with atypical kinetics, where substrate inhibition kinetics and non-Michaelis-Menten (sigmoidal) kinetics were found with respect to GMP and ATP, respectively. Additionally, the method was used in inhibition assays to screen twenty fragment-like compounds. The assays were robust and reproducible, with a signal window of 3.8 and a Z’ factor of 0.6. For the best inhibitor, an IC₅₀ curve was generated.

•

Simple HPLC separation of nucleotides involved in the guanylate kinase reaction.

•

Direct and unambiguous determination of percent conversion of substrate into product.

•

Successful application to Plasmodium vivax guanylate kinase (PvGK) activity studies.

•

Reaction catalyzed by PvGK found to be associated with atypical kinetics.

•

Robust and reproducible inhibition assay for compound screening.

Solution structure and functional investigation of human guanylate kinase reveals allosteric networking and a crucial role for the enzyme in cancer

Human guanylate kinase (hGMPK) is the only known enzyme responsible for cellular GDP production, making it essential for cellular viability and proliferation. Moreover, hGMPK has been assigned a critical role in metabolic activation of antiviral and antineoplastic nucleoside-analog prodrugs. Given that hGMPK is indispensable for producing the nucleotide building blocks of DNA, RNA, and cGMP and that cancer cells possess elevated GTP levels, it is surprising that a detailed structural and functional characterization of hGMPK is lacking. Here, we present the first high-resolution structure of hGMPK in the apo form, determined with NMR spectroscopy. The structure revealed that hGMPK consists of three distinct regions designated as the LID, GMP-binding (GMP-BD), and CORE domains and is in an open configuration that is nucleotide binding-competent. We also demonstrate that nonsynonymous single-nucleotide variants (nsSNVs) of the hGMPK CORE domain distant from the nucleotide-binding site of this domain modulate enzymatic activity without significantly affecting hGMPK's structure. Finally, we show that knocking down the hGMPK gene in lung adenocarcinoma cell lines decreases cellular viability, proliferation, and clonogenic potential while not altering the proliferation of immortalized, noncancerous human peripheral airway cells. Taken together, our results provide an important step toward establishing hGMPK as a potential biomolecular target, from both an orthosteric (ligand-binding sites) and allosteric (location of CORE domain-located nsSNVs) standpoint.

Dissipation at the angstrom scale: Probing the surface and interior of an enzyme

Pursuing a materials science approach to understanding the deformability of enzymes, we introduce measurements of the phase of the mechanical response function within the nanorheology paradigm. Driven conformational motion of the enzyme is dissipative as characterized by the phase measurements. The dissipation originates both from the surface hydration layer and the interior of the molecule, probed by examining the effect of point mutations on the mechanics. We also document changes in the mechanics of the enzyme examined, guanylate kinase, upon binding its four substrates. GMP binding stiffens the molecule, ATP and ADP binding softens it, while there is no clear mechanical signature of GDP binding. A hyperactive two-Gly mutant is found to possibly trade specificity for speed. Global deformations of enzymes are shown to be dependent on both hydration layer and polypeptide chain dynamics.

1H, 13C and 15N resonance assignment of human guanylate kinase

Human guanylate kinase (hGMPK) is a critical enzyme that, in addition to phosphorylating its physiological substrate (d)GMP, catalyzes the second phosphorylation step in the conversion of anti-viral and anti-cancer nucleoside analogs to their corresponding active nucleoside analog triphosphates. Until now, a high-resolution structure of hGMPK is unavailable and thus, we studied free hGMPK by NMR and assigned the chemical shift resonances of backbone and side chain ¹H, ¹³C, and ¹⁵N nuclei as a first step towards the enzyme's structural and mechanistic analysis with atomic resolution.

Eukaryotic-type serine/threonine kinase mediated phosphorylation at Thr169 perturbs mycobacterial guanylate kinase activity

Guanylate kinase is an essential and conserved enzyme in nucleotide biosynthetic pathway that transfers phosphoryl group of ATP to GMP for yielding GDP. Here, we report the phosphorylation of guanylate kinase from Mycobacterium tuberculosis (mGmk) by eukaryotic-type Ser/Thr kinase, PknA. Mass spectrometric studies identified Thr¹⁰¹ and Thr¹⁶⁹ as phosphorylatable residues in mGmk. To evaluate the significance of phosphorylation in these threonines, two point (T101A and T169A) and one double (T101A-T169A) mutants were generated. The kinase assay with these mutant proteins revealed the major contribution of Thr¹⁶⁹ compared with Thr¹⁰¹ in the phosphorylation of mGmk. Kinetic analysis indicated that p-mGmk was deficient in its enzymatic activity compared with that of its un-phosphorylated counterpart. Surprisingly, its phosphoablated (T169A) as well as phosphomimic (T169E) variants exhibited decreased activity as was observed with p-mGmk. Structural analysis suggested that phosphorylation of Thr¹⁶⁹ might affect its interaction with Arg¹⁶⁶, which is crucial for the functioning of mGmk. In fact, the R166A and R166K mutant proteins displayed a drastic decrease in enzymatic activity compared with that of the wild-type mGmk. Molecular dynamics (MD) studies of mGmk revealed that upon phosphorylation of Thr¹⁶⁹, the interactions of Arg¹⁶⁵/Arg¹⁶⁶ with Glu¹⁵⁸, Asp¹²¹ and residues of the loop in GMP-binding domain are perturbed. Taken together, our results illuminate the mechanistic insights into phosphorylation-mediated modulation of the catalytic activity of mGmk.

Structural studies of a hyperthermophilic thymidylate kinase enzyme reveal conformational substates along the reaction coordinate

Thymidylate kinase (TMK) is a key enzyme which plays an important role in DNA synthesis. It belongs to the family of nucleoside monophosphate kinases, several of which undergo structure-encoded conformational changes to perform their function. However, the absence of three-dimensional structures for all the different reaction intermediates of a single TMK homolog hinders a clear understanding of its functional mechanism. We herein report the different conformational states along the reaction coordinate of a hyperthermophilic TMK from Aquifex aeolicus, determined via X-ray diffraction and further validated through normal-mode studies. The analyses implicate an arginine residue in the Lid region in catalysis, which was confirmed through site-directed mutagenesis and subsequent enzyme assays on the wild-type protein and mutants. Furthermore, the enzyme was found to exhibit broad specificity toward phosphate group acceptor nucleotides. Our comprehensive analyses of the conformational landscape of TMK, together with associated biochemical experiments, provide insights into the mechanistic details of TMK-driven catalysis, for example, the order of substrate binding and the reaction mechanism for phosphate transfer. Such a study has utility in the design of potent inhibitors for these enzymes.

DATABASE

Structural data are available in the PDB under the accession numbers 2PBR, 4S2E, 5H5B, 5XAI, 4S35, 5XB2, 5H56, 5XB3, 5H5K, 5XB5, and 5XBH.

Strain analysis of protein structures and low dimensionality of mechanical allosteric couplings

In many proteins, especially allosteric proteins that communicate regulatory states from allosteric to active sites, structural deformations are functionally important. To understand these deformations, dynamical experiments are ideal but challenging. Using static structural information, although more limited than dynamical analysis, is much more accessible. Underused for protein analysis, strain is the natural quantity for studying local deformations. We calculate strain tensor fields for proteins deformed by ligands or thermal fluctuations using crystal and NMR structure ensembles. Strains-primarily shears-show deformations around binding sites. These deformations can be induced solely by ligand binding at distant allosteric sites. Shears reveal quasi-2D paths of mechanical coupling between allosteric and active sites that may constitute a widespread mechanism of allostery. We argue that strain-particularly shear-is the most appropriate quantity for analysis of local protein deformations. This analysis can reveal mechanical and biological properties of many proteins.

Nucleotide Metabolism and DNA Replication

The development and application of a highly versatile suite of tools for mycobacterial genetics, coupled with widespread use of "omics" approaches to elucidate the structure, function, and regulation of mycobacterial proteins, has led to spectacular advances in our understanding of the metabolism and physiology of mycobacteria. In this article, we provide an update on nucleotide metabolism and DNA replication in mycobacteria, highlighting key findings from the past 10 to 15 years. In the first section, we focus on nucleotide metabolism, ranging from the biosynthesis, salvage, and interconversion of purine and pyrimidine ribonucleotides to the formation of deoxyribonucleotides. The second part of the article is devoted to DNA replication, with a focus on replication initiation and elongation, as well as DNA unwinding. We provide an overview of replication fidelity and mutation rates in mycobacteria and summarize evidence suggesting that DNA replication occurs during states of low metabolic activity, and conclude by suggesting directions for future research to address key outstanding questions. Although this article focuses primarily on observations from Mycobacterium tuberculosis, it is interspersed, where appropriate, with insights from, and comparisons with, other mycobacterial species as well as better characterized bacterial models such as Escherichia coli. Finally, a common theme underlying almost all studies of mycobacterial metabolism is the potential to identify and validate functions or pathways that can be exploited for tuberculosis drug discovery. In this context, we have specifically highlighted those processes in mycobacterial DNA replication that might satisfy this critical requirement.

Purification and characterization of guanylate kinase, a nucleoside monophosphate kinase of Brugia malayi

Guanylate kinase, a nucleoside monophosphate kinase of Brugia malayi which is involved in reversible transfer of phosphate groups from ATP to GMP, was cloned, expressed and characterized. The native molecular mass of BmGK was found to be 45 kDa as determined by size exclusion chromatography and glutaraldehyde cross-linking which revealed that the protein is homodimer in nature. This is a unique characteristic among known eukaryotic GKs. GMP and ATP served as the most effective phosphate acceptor and donor, respectively. Recombinant BmGK utilized both GMP and dGMP, as substrates showing Km values of 30 and 38 μm, respectively. Free Mg+2 (un-complexed to ATP) and GTP play a regulatory role in catalysis of BmGK. The enzyme showed higher catalytic efficiency as compared with human GK and showed ternary complex (BmGK-GMP-ATP) formation with sequential substrate binding. The secondary structure of BmGK consisted of 45% α-helices, 18% β-sheets as revealed by CD analysis. Homology modelling and docking with GMP revealed conserved substrate binding residues with slight differences. Differences in kinetic properties and oligomerization of BmGK compared with human GK can provide the way for design of parasite-specific inhibitors.

Enzyme closure and nucleotide binding structurally lock guanylate kinase

We investigate the conformational dynamics and mechanical properties of guanylate kinase (GK) using a multiscale approach combining high-resolution atomistic molecular dynamics and low-resolution Brownian dynamics simulations. The GK enzyme is subject to large conformational changes, leading from an open to a closed form, which are further influenced by the presence of nucleotides. As suggested by recent work on simple coarse-grained models of apo-GK, we primarily focus on GK's closure mechanism with the aim to establish a detailed picture of the hierarchy and chronology of structural events essential for the enzymatic reaction. We have investigated open-versus-closed, apo-versus-holo, and substrate-versus-product-loaded forms of the GK enzyme. Bound ligands significantly modulate the mechanical and dynamical properties of GK and rigidity profiles of open and closed states hint at functionally important differences. Our data emphasizes the role of magnesium, highlights a water channel permitting active site hydration, and reveals a structural lock that stabilizes the closed form of the enzyme.

Functional modes and residue flexibility control the anisotropic response of guanylate kinase to mechanical stress

The coupling between the mechanical properties of enzymes and their biological activity is a well-established feature that has been the object of numerous experimental and theoretical works. In particular, recent experiments show that enzymatic function can be modulated anisotropically by mechanical stress. We study such phenomena using a method for investigating local flexibility on the residue scale that combines a reduced protein representation with Brownian dynamics simulations. We performed calculations on the enzyme guanylate kinase to study its mechanical response when submitted to anisotropic deformations. The resulting modifications of the protein's rigidity profile can be related to the changes in substrate binding affinity observed experimentally. Further analysis of the principal components of motion of the trajectories shows how the application of a mechanical constraint on the protein can disrupt its dynamics, thus leading to a decrease of the enzyme's catalytic rate. Eventually, a systematic probe of the protein surface led to the prediction of potential hotspots where the application of an external constraint would produce a large functional response both from the mechanical and dynamical points of view. Such enzyme-engineering approaches open the possibility to tune catalytic function by varying selected external forces.

Structural and functional characterization of the Mycobacterium tuberculosis uridine monophosphate kinase: insights into the allosteric regulation

Nucleoside Monophosphate Kinases (NMPKs) family are key enzymes in nucleotide metabolism. Bacterial UMPKs depart from the main superfamily of NMPKs. Having no eukaryotic counterparts they represent attractive therapeutic targets. They are regulated by GTP and UTP, while showing different mechanisms in Gram(+), Gram(–) and archaeal bacteria. In this work, we have characterized the mycobacterial UMPK (UMPKmt) combining enzymatic and structural investigations with site-directed mutagenesis. UMPKmt exhibits cooperativity toward ATP and an allosteric regulation by GTP and UTP. The crystal structure of the complex of UMPKmt with GTP solved at 2.5 Å, was merely identical to the modelled apo-form, in agreement with SAXS experiments. Only a small stretch of residues was affected upon nucleotide binding, pointing out the role of macromolecular dynamics rather than major structural changes in the allosteric regulation of bacterial UMPKs. We further probe allosteric regulation by site-directed mutagenesis. In particular, a key residue involved in the allosteric regulation of this enzyme was identified.

Binding dynamics and energetic insight into the molecular forces driving nucleotide binding by guanylate kinase

Plasmodium deoxyguanylate pathways are an attractive area of investigation for future metabolic and drug discovery studies due to their unique substrate specificities. We investigated the energetic contribution to guanylate kinase substrate binding and the forces underlying ligand recognition. In the range from 20 to 35°C, the thermodynamic profiles displayed marked decrease in binding enthalpy, while the free energy of binding showed little changes. GMP produced a large binding heat capacity change of -356 cal mol(-1) K(-1), indicating considerable conformational changes upon ligand binding. Interestingly, the calculated ΔCp was -32 cal mol(-1) K(-1), indicating that the accessible surface area is not the central change in substrate binding, and that other entropic forces, including conformational changes, are more predominant. The thermodynamic signature for GMP is inconsistent with rigid-body association, while dGMP showed more or less rigid-body association. These binding profiles explain the poor catalytic efficiency and low affinity for dGMP compared with GMP. At low temperature, the ligands bind to the receptor site under the effect of hydrophobic forces. Interestingly, by increasing the temperature, the entropic forces gradually vanish and proceed to a nonfavorable contribution, and the interaction occurs mainly through bonding, electrostatic forces, and van der Waals interactions.

Complex molecular assemblies at hand via interactive simulations

Studying complex molecular assemblies interactively is becoming an increasingly appealing approach to molecular modeling. Here we focus on interactive molecular dynamics (IMD) as a textbook example for interactive simulation methods. Such simulations can be useful in exploring and generating hypotheses about the structural and mechanical aspects of biomolecular interactions. For the first time, we carry out low-resolution coarse-grain IMD simulations. Such simplified modeling methods currently appear to be more suitable for interactive experiments and represent a well-balanced compromise between an important gain in computational speed versus a moderate loss in modeling accuracy compared to higher resolution all-atom simulations. This is particularly useful for initial exploration and hypothesis development for rare molecular interaction events. We evaluate which applications are currently feasible using molecular assemblies from 1900 to over 300,000 particles. Three biochemical systems are discussed: the guanylate kinase (GK) enzyme, the outer membrane protease T and the soluble N-ethylmaleimide-sensitive factor attachment protein receptors complex involved in membrane fusion. We induce large conformational changes, carry out interactive docking experiments, probe lipid-protein interactions and are able to sense the mechanical properties of a molecular model. Furthermore, such interactive simulations facilitate exploration of modeling parameters for method improvement. For the purpose of these simulations, we have developed a freely available software library called MDDriver. It uses the IMD protocol from NAMD and facilitates the implementation and application of interactive simulations. With MDDriver it becomes very easy to render any particle-based molecular simulation engine interactive. Here we use its implementation in the Gromacs software as an example.

Molecular cloning, expression, characterization and mutation of Plasmodium falciparum guanylate kinase

The present work describes cloning, expression, purification, characterization, and mutation of Plasmodium falciparum guanylate kinase (PlasmoDB ID PFI1420w). Amino-acid sequence alignment revealed important differences especially in K42-V51, Y73-A77, and F100-L110, which include residues important for kinase activity, and at helix 3, which is important for domain movements. The catalytic efficiency for dGMP was 22-fold lower than that for GMP, whose value is the lowest among known guanylate kinases. dGMP was found to a competitive inhibitor for GMP with K(i)=0.148 mM and a mixed-type inhibitor with regard to ATP with measured K(i)=0.4 mM. The specificity constant (K(cat)/K(m)) of the four examined mutants varied for natural substrate GMP/dGMP, indicating the involvement of different mechanisms in substrate recognition and subsequent loop-domain movement. These results show that P. falciparum guanylate kinase is structurally and biochemically distinct from other guanylate kinases and could be a possible target in drug development.

Structural and functional consequences of single amino acid substitutions in the pyrimidine base binding pocket of Escherichia coli CMP kinase

Bacterial CMP kinases are specific for CMP and dCMP, whereas the related eukaryotic NMP kinase phosphorylates CMP and UMP with similar efficiency. To explain these differences in structural terms, we investigated the contribution of four key amino acids interacting with the pyrimidine ring of CMP (Ser36, Asp132, Arg110 and Arg188) to the stability, catalysis and substrate specificity of Escherichia coli CMP kinase. In contrast to eukaryotic UMP/CMP kinases, which interact with the nucleobase via one or two water molecules, bacterial CMP kinase has a narrower NMP-binding pocket and a hydrogen-bonding network involving the pyrimidine moiety specific for the cytosine nucleobase. The side chains of Arg110 and Ser36 cannot establish hydrogen bonds with UMP, and their substitution by hydrophobic amino acids simultaneously affects the K(m) of CMP/dCMP and the k(cat) value. Substitution of Ser for Asp132 results in a moderate decrease in stability without significant changes in K(m) value for CMP and dCMP. Replacement of Arg188 with Met does not affect enzyme stability but dramatically decreases the k(cat)/K(m) ratio compared with wild-type enzyme. This effect might be explained by opening of the enzyme/nucleotide complex, so that the sugar no longer interacts with Asp185. The reaction rate for different modified CMP kinases with ATP as a variable substrate indicated that none of changes induced by these amino acid substitutions was 'propagated' to the ATP subsite. This 'modular' behavior of E. coli CMP kinase is unique in comparison with other NMP kinases.

Crystal structure and nucleotide binding of the Thermus thermophilus RNA helicase Hera N-terminal domain

DEAD box RNA helicases use the energy of ATP hydrolysis to unwind double-stranded RNA regions or to disrupt RNA/protein complexes. A minimal RNA helicase comprises nine conserved motifs distributed over two RecA-like domains. The N-terminal domain contains all motifs involved in nucleotide binding, namely the Q-motif, the DEAD box, and the P-loop, as well as the SAT motif, which has been implicated in the coordination of ATP hydrolysis and RNA unwinding. We present here the crystal structure of the N-terminal domain of the Thermus thermophilus RNA helicase Hera in complex with adenosine monophosphate (AMP). Upon binding of AMP the P-loop adopts a partially collapsed or half-open conformation that is still connected to the DEAD box motif, and the DEAD box in turn is linked to the SAT motif via hydrogen bonds. This network of interactions communicates changes in the P-loop conformation to distant parts of the helicase. The affinity of AMP is comparable to that of ADP and ATP, substantiating that the binding energy from additional phosphate moieties is directly converted into conformational changes of the entire helicase. Importantly, the N-terminal Hera domain forms a dimer in the crystal similar to that seen in another thermophilic prokaryote. It is possible that this mode of dimerization represents the prototypic architecture in RNA helicases of thermophilic origin.

Crystal structures of GMP kinase in complex with ganciclovir monophosphate and Ap5G

Guanosine monophosphate kinases (GMPK), by catalyzing the phosphorylation of GMP or dGMP, are of dual potential in assisting the activation of anti-viral prodrugs or as candidates for antibiotic strategies. Human GMPK is an obligate step for the activation of acyclic guanosine analogs, such as ganciclovir, which necessitate efficient phosphorylation, while GMPK from bacterial pathogens, in which this enzyme is essential, are potential targets for therapeutic inhibition. Here we analyze these two aspects of GMPK activity with the crystal structures of Escherichia coli GMPK in complex with ganciclovir-monophosphate (GCV-MP) and with a bi-substrate inhibitor, Ap5G. GCV-MP binds as GMP to the GMP-binding domain, which is identical in E. coli and human GMPKs, but unlike the natural substrate fails to stabilize the closed, catalytically-competent conformation of this domain. Comparison with GMP- and GDP-bound GMPK structures identifies the 2'hydroxyl of the ribose moiety as responsible for hooking the GMP-binding domain onto the CORE domain. Absence of this hydroxyl in GCV-MP impairs the stabilization of the active conformation, and explains why GCV-MP is phosphorylated less efficiently than GMP, but as efficiently as dGMP. In contrast, Ap5G is an efficient inhibitor of GMPK. The crystal structure shows that Ap5G locks an incompletely closed conformation of the enzyme, in which the adenine moiety is located outside its expected binding site. Instead, it binds at a subunit interface that is unique to the bacterial enzyme, which is in equilibrium between a dimeric and an hexameric form in solution. This suggests that inhibitors could be designed to bind at this interface such as to prevent nucleotide-induced domain closure. Altogether, these complexes point to domain motions as critical components to be evaluated in therapeutic strategies targeting NMP kinases, with opposite effects depending on whether efficient phosphorylation or inhibition is being sought after.

Calorimetric and crystallographic analysis of the oligomeric structure of Escherichia coli GMP kinase

Guanosine monophosphate kinases (GMPKs), which catalyze the phosphorylation of GMP and dGMP to their diphosphate form, have been characterized as monomeric enzymes in eukaryotes and prokaryotes. Here, we report that GMPK from Escherichia coli (ecGMPK) assembles in solution and in the crystal as several different oligomers. Thermodynamic analysis of ecGMPK using differential scanning calorimetry shows that the enzyme is in equilibrium between a dimer and higher order oligomers, whose relative amounts depend on protein concentration, ionic strength, and the presence of ATP. Crystallographic structures of ecGMPK in the apo, GMP and GDP-bound forms were solved at 3.2A, 2.9A and 2.4A resolution, respectively. ecGMPK forms a hexamer with D3 symmetry in all crystal forms, in which the two nucleotide-binding domains are able to undergo closure comparable to that of monomeric GMPKs. The 2-fold and 3-fold interfaces involve a 20-residue C-terminal extension and a sequence signature, respectively, that are missing from monomeric eukaryotic GMPKs, explaining why ecGMPK forms oligomers. These signatures are found in GMPKs from proteobacteria, some of which are human pathogens. GMPKs from these bacteria are thus likely to form the same quaternary structures. The shift of the thermodynamic equilibrium towards the dimer at low ecGMPK concentration together with the observation that inter-subunit interactions partially occlude the ATP-binding site in the hexameric structure suggest that the dimer may be the active species at physiological enzyme concentration.