Dynamic origins of substrate promiscuity in bacterial galactokinases.

Galactokinase catalyses the ATP-dependent phosphorylation of galactose and structurally related sugars. The enzyme has attracted interest as a potential biocatalyst for the production of sugar 1-phosphates and several attempts have been made to broaden its specificity. In general, bacterial galactokinases have wider substrate ranges than mammalian ones. The enzymes from Escherichia coli and Lactococcus lactis have received particular attention and a number of variants with increased promiscuity have been identified. Here, we present a molecular dynamics study designed to investigate the molecular causes of the wider substrate ranges of these enzymes and their variants with particular reference to protein mobility. Some regions close to the active site of the enzyme have different structures in the bacterial enzymes compared to the human one. Alterations known to increase the substrate range (e.g. Y371H in the E. coli enzyme), tend to alter the conformation of a key α-helical region (residues 216-232 in the E. coli enzyme). The equivalent helix in the human enzyme has previously been predicted to be altered in variants which affect catalytic activity or protein stability. This helix appears to be a key region in galactokinases from a range of species and may represent an interesting target for future attempts to broaden the specificity of galactokinases.

Modulation of the mobility of a key region in human galactokinase: Impacts on catalysis and stability.

Galactokinase catalyses the phosphorylation of α-d-galactose and some structurally related monosaccharides. The enzyme is of interest due to its potential as a biocatalyst for the production of sugar 1-phosphates and due to its involvement in the inherited metabolic disease type II galactosemia. It has been previously shown that a region (residues 231-245) in human galactokinase often has altered mobility when active site residues are varied. We hypothesised that the reverse may be true and that designing changes to this region might affect the functioning of the active site of the enzyme. Focussing on four residues (Leu-231, Gln-242, Glu-244 and Glu-245) we conducted molecular dynamics simulations to explore the effects of changing these residues to glycine or serine. In most cases the variations resulted in local changes to the 231-245 region and global changes to the root mean squared fluctuation (RMSF) of the protein. The four serine variants were expressed as recombinant proteins. All had altered steady state enzyme kinetic parameters with α-d-galactose as a substrate. However, these changes were generally less than ten-fold in magnitude. Changes were also observed with 2-deoxy-α-d-galactose, α-d-galactosamine and α-d-talose as substrates, including (in some cases) loss of detectable activity, suggesting that these variations can tune the specificity of the enzyme. This study demonstrates that activity and specificity of human galactokinase can be modulated by variations designed to affect active site flexibility. It is likely that this principle can be generalised to other enzymes.

A galactokinase-like protein from the liver fluke Fasciola hepatica.

Galactokinase catalyses the ATP-dependent phosphorylation of galactose. A galactokinase-like sequence was identified in a Fasciola hepatica EST library. Recombinant expression of the corresponding protein in Escherichia coli resulted in a protein of approximately 50 kDa. The protein is monomeric, like galactokinases from higher animals, yeasts and some bacteria. The protein has no detectable enzymatic activity with galactose or N-acetylgalactosamine as a substrate. However, it does bind to ATP. Molecular modelling predicted that the protein adopts a similar fold to galactokinase and other GHMP kinases. However, a key loop in the active site was identified which may influence the lack of activity. Sequence analysis strongly suggested that this protein (and other proteins annotated as "galactokinase" in the trematodes Schistosoma mansoni and Clonorchis sinensis) are closer to N-acetylgalactosamine kinases. No other galactokinase-like sequences appear to be present in the genomes of these three species. This raises the intriguing possibility that these (and possibly other) trematodes are unable to catabolise galactose through the Leloir pathway due to the lack of a functional galactokinase.

Improving the Activity and Stability of Human Galactokinase for Therapeutic and Biotechnological Applications.

Galactokinase catalyses the site- and stereospecific phosphorylation of α-d-galactose. As such it has attracted interest as a biocatalyst for the introduction of phosphate groups into monosaccharides. However, attempts to broaden the substrate range of human galactokinase have generally resulted in substantially reduced activity. The enzyme also has biotechnological potential in enzyme replacement therapy (ERT) for type II galactosaemia. The return-to-consensus approach can be used to identify residues that can be altered to increase protein stability and enzyme activity. This approach identified six residues of potential interest in human galactokinase. Some of the single consensus variants (M60V, D268E, A334S and G373S) increased the catalytic turnover of the enzyme, but none resulted in improved stability. When all six changes were introduced into the protein (M60V/M180V/D268E/A334S/R366Q/G373S), thermal stability was increased. Molecular dynamics simulations suggested that these changes altered the protein's conformation at key sites. The number of salt bridges and hydrogen bonds was also increased. Combining the six consensus variations with Y379W (a variant with greater substrate promiscuity) increased the stability of this variant and its turnover towards some substrates. Thus, the six consensus variants can be used to stabilise catalytically interesting variants of human galactokinase and might also be useful if the protein were to be used in ERT.

Modulating Mobility: a Paradigm for Protein Engineering?

Proteins are highly mobile structures. In addition to gross conformational changes occurring on, for example, ligand binding, they are also subject to constant thermal motion. The mobility of a protein varies through its structure and can be modulated by ligand binding and other events. It is becoming increasingly clear that this mobility plays an important role in key functions of proteins including catalysis, allostery, cooperativity, and regulation. Thus, in addition to an optimum structure, proteins most likely also require an optimal dynamic state. Alteration of this dynamic state through protein engineering will affect protein function. A dramatic example of this is seen in some inherited metabolic diseases where alternation of residues distant from the active site affects the mobility of the protein and impairs function. We postulate that using molecular dynamics simulations, experimental data or a combination of the two, it should be possible to engineer the mobility of active sites. This may be useful in, for example, increasing the promiscuity of enzymes. Thus, a paradigm for protein engineering is suggested in which the mobility of the active site is rationally modified. This might be combined with more "traditional" approaches such as altering functional groups in the active site.

Insight into the mechanism of galactokinase: Role of a critical glutamate residue and helix/coil transitions.

Galactokinase, the enzyme which catalyses the first committed step in the Leloir pathway, has attracted interest due to its potential as a biocatalyst and as a possible drug target in the treatment of type I galactosemia. The mechanism of the enzyme is not fully elucidated. Molecular dynamics (MD) simulations of galactokinase with the active site residues Arg-37 and Asp-186 altered predicted that two regions (residues 174-179 and 231-240) had different dynamics as a consequence. Interestingly, the same two regions were also affected by alterations in Arg-105, Glu-174 and Arg-228. These three residues were identified as important in catalysis in previous computational studies on human galactokinase. Alteration of Arg-105 to methionine resulted in a modest reduction in activity with little change in stability. When Arg-228 was changed to methionine, the enzyme's interaction with both ATP and galactose was affected. This variant was significantly less stable than the wild-type protein. Changing Glu-174 to glutamine (but not to aspartate) resulted in no detectable activity and a less stable enzyme. Overall, these combined in silico and in vitro studies demonstrate the importance of a negative charge at position 174 and highlight the critical role of the dynamics in to key regions of the protein. We postulate that these regions may be critical for mediating the enzyme's structure and function.