Thermodynamic linkages in rabbit muscle pyruvate kinase: kinetic, equilibrium, and structural studies

The mechanism of allosteric regulation of rabbit muscle pyruvate kinase (PK) was examined in the presence of the allosteric inhibitor phenylalanine (Phe). Steady-state kinetic, equilibrium binding, and structural studies were conducted to provide a broad data base to establish a reasonable model for the interactions. Phe was shown to induce apparent cooperativity in the steady-state kinetic measurements at pH 7.5 and 23 degrees C. The apparent Km for phosphoenolpyruvate was shown to increase with increasing Phe concentrations. These results imply that Phe reduces the affinity of PK for phosphoenolpyruvate. This conclusion was substantiated by equilibrium binding studies which yielded association constants of phosphoenolpyruvate as a function of Phe concentration. The binding constant of Phe was also determined at pH 7.0 and 23 degrees C. The effect of ligands on the hydrodynamic properties of PK was monitored by difference sedimentation velocity, sedimentation velocity, and equilibrium experiments. The results showed that PK remains tetrameric both in the presence and in the absence of Phe. However, Phe induces a small decrease in the sedimentation coefficient of the enzyme; hence, it suggests a loosening of the protein structure. The accessibility of the sulfhydryl residues of the enzyme also increases in the presence of Phe. Furthermore, the Phe-induced conformational change was approximately 90% complete when only 25% of the binding sites were saturated. This result suggested that the regulatory behavior of PK might satisfactorily be described by the two-state model of Monod-Wyman-Changeux [Monod, J., Wyman, J., & Changeux, J.-P. (1965) J. Mol. Biol. 12, 88-118].     

Modulation of allostery of pyruvate kinase by shifting of an ensemble of microstates

Since the introduction of the concepts of allostery about four decades ago, much advancement has been made in elucidating the structure-function correlation in allostery. However, there are still a number of issues that remain unresolved. In this review we used mammalian pyruvate kinase (PK) as a model system to understand the role of protein dynamics in modulating cooperativity. PK has a triosephosphate isomerase (TIM)(α/β)₈ barrel structural motif. PK is an ideal system to address basic questions regarding regulatory mechanisms about this common (α/β)₈ structural motif. The simplest model accounting for all of the solution thermodynamic and kinetic data on ligand-enzyme interactions involves two conformational states, inactive E^T and active E^R. These conformational states are represented by domain movements. Further studies provide the first evidence for a differential effect of ligand binding on the dynamics of the structural elements, not major secondary structural changes. These data are consistent with our model that allosteric regulation of PK is the consequence of perturbation of the distribution of an ensemble of states in which the inactive E^T and active E^R represent the two extreme end states. Sequence differences and ligands can modulate the distribution of states leading to alterations of functions. The future work includes: defining the network of functionally connected residues; elucidating the chemical principles governing the sequence differences which affect functions; and probing the nature of mutations on the stability of the secondary structural elements, which in turn modulate allostery.     

Phosphoenolpyruvate and Mg Binding to Pyruvate Kinase Monitored by Infrared Spectroscopy

Structural changes in rabbit muscle pyruvate kinase (PK) induced by phosphoenolpyruvate (PEP) and Mg²⁺ binding were studied by attenuated total reflection Fourier transform infrared spectroscopy in combination with a dialysis accessory. The experiments indicated a largely preserved secondary structure upon PEP and Mg²⁺ binding but also revealed small backbone conformational changes of PK involving all types of secondary structure. To assess the effect of the protein environment on the bound PEP, we assigned and evaluated the infrared absorption bands of bound PEP. These were identified using 2,3-¹³C₂-labeled PEP. We obtained the following assignments: 1589 cm⁻¹ (antisymmetric carboxylate stretching vibration); 1415 cm⁻¹ (symmetric carboxylate stretching vibration); 1214 cm⁻¹ (C-O stretching vibration); 1124 and 1110 cm⁻¹ (asymmetric PO₃^2- stretching vibrations); and 967 cm⁻¹ (symmetric PO₃^2- stretching vibration). The corresponding band positions in solution are 1567, 1407, 1229, 1107, and 974 cm⁻¹. The differences for bound and free PEP indicate specific interactions between ligand and protein. Quantification of the interactions with the phosphate group indicated that the enzyme environment has little influence on the P-O bond strengths, and that the bridging P-O bond, which is broken in the catalytic reaction, is weakened by <3%. Thus, there is only little distortion toward a dissociative transition state of the phosphate transfer reaction when PEP binds to PK. Therefore, our results are in line with an associative transition state. Carboxylate absorption bands indicated a maximal shortening of the length of the shorter C-O bond by 1.3 pm. PEP bound to PK in the presence of the monovalent ion Na⁺ exhibited the same band positions as in the presence of K⁺, indicating very similar interaction strengths between ligand and protein in both cases.     

高粱叶片PEP羧化酶的溶液构象——近紫外圆二色性光谱研究

用近紫外CD光谱技术追踪了PEP羧化酶与各种配基的相互作用。底物PEP、必需金属离子Mg~( )、PEP-Mg~( )以及效应剂G6P、Gly、G6P-Gly,均可引起高粱叶片PEP羧化酶近紫外CD光谱各不相同的变化。这表明高粱叶片PEP羧化酶分子构象有较大的灵活性,不同的配基与酶相互作用可引起酶分子不同的构象变化,因而使酶分子表现出催化功能、调节特性、必需氨基酸残基的化学反应性以及稳定性诸方面的差异。     

Pyruvate kinase from Calicophoron ijimai trematodes and its potential inhibition by anthelmintic preparations

It was shown that pyruvate kinase (PK) in the supernatant fraction from Calicophoron ijimai is able to regulate the direction of metabolic flow at glucose break down from phosphoenolpyruvate (PEP) level. The enzyme for activity required substrate, dinucleotides, cations K+ and Mn++. The activity with Mg++ as divalent cation is low. The addition of fructose-1.6-diphosphate (FDP) did not affect the enzyme activity with Mn++, however, increased the affinity for PEP. The velocity of Mg++ activated reaction increased by 8.2 times in the presence of FDP. PK in C. ijimai is sensitive to ATP inhibition, being weakly inhibited by malate. L-alanine did not influence on the enzyme activity. The effect of some anthelminthic preparations on the PK activity was shown.     

Partial purification and characterization of pyruvate kinase from the plant fraction of soybean root nodules

Pyruvate kinase (PK, EC 2.7.1.40) was partially purified from the plant cytosolic fraction of N₂-fixing soybean ( Glycine max [L.] Merr.) root nodules. The partially purified PK preparation was completely free of contamination by phospho enol pyruvate carboxylase (PEPC, EC 4.1.1.31), the other major phospho enol pyruvate (PEP)-utilizing enzyme in legume root nodules. Latency experiments with sonicated nodule extracts showed that Bradyrhizobium japonicum bacteroids do not express either PK or PEPC activity in symbiosis. In contrast, free-living B. japonicum bacteria expressed PK activity, but not PEPC activity. Antibodies specific for the cytosolic isoform of PK from castor bean endosperm cross-reacted with a 52-kDa polypeptide in the partially purified PK preparation. At the optimal assay pH (pH 8.0 for PEPC and pH 6.9 for PK) and in the absence of malate, PEPC activity in crude nodule extracts was 2.6 times the corresponding PK activity. This would tend to favour PEP metabolism by PEPC over PEP metabolism by PK. However, at pH 7.0 in the presence of 5 m M malate, PEPC activity was strongly inhibited, but PK activity was unaffected. Thus, we propose that PK and PEPC activity in legume root nodules may be coordinately regulated by fluctuations in malate concentration in the plant cytosolic fraction of the bacteroid-containing cells. Reduced uptake of malate by the bacteroids, as a result of reduced rates of N₂ fixation, may favour PEP metabolism by PK over PEP metabolism by PEPC.     

Dynamic cross-talk among remote binding sites: the molecular basis for unusual synergistic allostery

Allosteric regulation of protein function is critical for metabolic control. Binding of allosteric effectors elicits a functional change in a remote ligand binding site on a protein by altering the equilibrium between different forms in the protein ensemble. 3-Deoxy-d-arabino-heptulosonate 7-phosphate synthase (DAH7PS) catalyzes the first step in the shikimate pathway, which is responsible for the biosynthesis of aromatic amino acids Trp, Phe, and Tyr. Feedback regulation by the aromatic amino acids is important for controlling the cellular levels of the aromatic amino acids, and many organisms have two or more DAH7PS isozymes that show differing sensitivities to aromatic compounds. Mycobacterium tuberculosis expresses a single DAH7PS that is insensitive to the presence of a single amino acid yet shows extraordinary synergistic inhibition by combinations of the pathway end products Trp and Phe. The Trp+Phe-bound structure for M. tuberculosis DAH7PS, showing two separate binding sites occupied by Trp and Phe for each monomer of the tetrameric protein, was obtained by cocrystallization. Comparison of this structure with the ligand-free M. tuberculosis DAH7PS demonstrates that there is no significant change in conformation upon ligand binding, suggesting that contributions from altered dynamic properties of the enzyme may account for the allosteric inhibition. Isothermal titration calorimetry experiments demonstrate that the inhibitor binding sites are in direct communication. Molecular dynamics simulations reveal different changes in dynamic fluctuations upon single ligand binding compared to dual ligand binding. These changes account for the cross-talk between inhibitor binding sites and the active site, simultaneously potentiating both dual ligand binding and diminution of catalytic function.     

Crystal structure of pyruvate kinase from Geobacillus stearothermophilus

The pyruvate kinase (PK) from a moderate thermophile, Geobacillus stearothermophilus, is an allosteric enzyme activated by AMP and ribose 5-phosphate but not fructose 1, 6-bisphosphate (FBP), which is a common activator of PKs. It has an extra C-terminal sequence (ECTS), which contains a highly conserved phosphoenolpyruvate (PEP) binding motif, but its function and structure remain unclear. To elucidate the structural characteristics of the effector-binding site and the ECTS, the crystal structure of the C9S/C268S mutant of the enzyme was determined at 2.4 A resolution. The crystal belonged to space group P6(2)22, with unit cell parameters a, b = 145.97 A, c = 118.03 A. The enzyme was a homotetramer and its overall domain structure was similar to the previously solved structures except that the ECTS formed a new domain (C' domain). The structure of the C' domain closely resembled that of the PEP binding domain of maize pyruvate phosphate dikinase. A sulphate ion was found in a pocket in the effector-binding C domain. This site corresponds to the 6-phosphate group-binding site in yeast PK bound FBP and seems to be the effector-binding site. Through comparison of the structure of the putative effector-binding site to that of the FBP binding site of the yeast enzyme, the structural basis of the effector specificity of the G. stearothermophilus PK is discussed.     

Role of the loop containing residue 115 in the induced-fit mechanism of the bacterial cell wall biosynthetic enzyme MurA

The induced-fit mechanism in Enterobacter cloacae MurA has been investigated by kinetic studies and X-ray crystallography. The antibiotic fosfomycin, an irreversible inhibitor of MurA, induced a structural change in UDP-N-acetylglucosamine (UDPGlcNAc)-liganded enzyme with a time dependence similar to that observed for the inactivation progress. The mechanism of action of fosfomycin on MurA appeared to be of the bimolecular type, the overall rate constants of inactivation and structural change being = 104 M(-1) s(-1) and = 85 M(-1) s(-1), respectively. Fosfomycin as well as the second MurA substrate, phosphoenolpyruvate (PEP), are known to interact with the side chain of Cys115. Like wild-type MurA, the catalytically inactive single-site mutant protein Cys115Ser structurally interacted with UDPGlcNAc in a rapidly reversible reaction. However, in contrast to wild-type enzyme, binding of PEP to mutant protein induced a rate-limited, biphasic structural change. Fosfomycin did not affect the structure of the mutant protein. The crystal structure of unliganded Cys115Ser MurA at 1.9 A resolution revealed that the overall conformation of the loop comprising residues 112-121 is not influenced by the mutation. However, other than Cys115 in wild-type MurA, Ser115 exhibits two distinct side-chain conformations. A detailed view on the loop revealed the existence of an elaborate hydrogen-bonding network mainly supplied by water molecules, presumably stabilizing its conformation in the unliganded state. The comparison between the known crystal structures of MurA, together with the kinetic data obtained, suggest intermediate conformational states in the MurA reaction, in which the loop undergoes multiple structural changes upon ligand binding.     

Relation of enzyme activity to local/global stability of murine adenosine deaminase: 19F NMR studies

Adenosine deaminase (ADA, EC 3.5.4.4) is a ubiquitous (beta/alpha)8-barrel enzyme crucial for purine metabolism and normal immune competence. In this study, it was observed that loss of enzyme activity of murine ADA (mADA) precedes the global secondary and tertiary structure transition when the protein is exposed to denaturant. The structural mechanism for this phenomenon was probed using site-specific 19F NMR spectroscopy in combination with [6-19F]tryptophan labeling and inhibitor binding. There are four tryptophan residues in mADA and all are located more than 12 A from the catalytic site. The 19F NMR spectra of [6-19F]Trp-labelled mADA show that the urea-induced chemical shift change of 19F resonance of W161, one of the four tryptophan 19F nuclei, correlates with the loss of enzyme activity. The urea-induced chemical shift change of another 19F resonance of W117 correlates with the change of the apparent rate constant for the binding of transition-state analogue inhibitor deoxycoformycin to the enzyme. On the other hand, the chemical environment of the local region around W264 does not change significantly, as a consequence of perturbation by low concentrations of urea or substrate analog. The results indicate that different regions of mADA have different local stability, which controls the activity and stability of the enzyme. The results provide new insights into the relationship between the function of a protein and its conformational flexibility as well as its global stability. This study illustrates the advantage of 19F NMR spectroscopy in probing site-related conformational change information in ligand binding, enzymatic activity and protein folding.     

Role of phenylalanine B10 in plant nonsymbiotic hemoglobins

All plants contain an unusual class of hemoglobins that display bis-histidyl coordination yet are able to bind exogenous ligands such as oxygen. Structurally homologous hexacoordinate hemoglobins (hxHbs) are also found in animals (neuroglobin and cytoglobin) and some cyanobacteria, where they are thought to play a role in free radical scavenging or ligand sensing. The plant hxHbs can be distinguished from the others because they are only weakly hexcacoordinate in the ferrous state, yet no structural mechanism for regulating hexacoordination has been articulated to account for this behavior. Plant hxHbs contain a conserved Phe at position B10 (Phe(B10)), which is near the reversibly coordinated distal His(E7). We have investigated the effects of Phe(B10) mutation on kinetic and equilibrium constants for hexacoordination and exogenous ligand binding in the ferrous and ferric oxidation states. Kinetic and equilibrium constants for hexacoordination and ligand binding along with CO-FTIR spectroscopy, midpoint reduction potentials, and the crystal structures of two key mutant proteins (F40W and F40L) reveal that Phe(B10) is an important regulatory element in hexacoordination. We show that Phe at this position is the only amino acid that facilitates stable oxygen binding to the ferrous Hb and the only one that promotes ligand binding in the ferric oxidation states. This work presents a structural mechanism for regulating reversible intramolecular coordination in plant hxHbs.     

Homology modelling and molecular dynamics study of human fatty acid amide hydrolase

The three-dimensional (3D) model of the human fatty acid amide hydrolase (hFAAH) was constructed based on the crystal structure of the rat FAAH (PDB code 1MT5) in complex with a substrate using Modeller9v2 program. With the aid of molecular mechanics and molecular dynamics method, the last model was obtained and further assessed by Profile-3D, Prosa2003 and Procheck, which confirms that the refined model is reliable. Furthermore, the docking results of propofol and its structural analogue into the active site of hFAAH indicate that 2,6-di-sec-butyl phenol is a more preferred ligand than others, which is in good agreement with the experimental results. From the docking studies, we also suggest that Phe192, Ile238, Thr377, Leu380, Phe381, Phe388 and Leu404 in the hFAAH are seven important determinant residues in binding as they have strong van der Waal interactions with the ligand.     

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.     

Kinetic and expression profiling of cytosolic pyruvate kinase enzyme during seed development of Indian mustard (Brassica juncea)

In the recent years tremendous progress has been achieved in deconstructing the oil biosynthetic pathways, majority of which is in Arabidopsis. Glycolysis is fundamental to this process as it is the cardinal supplier of precursors for fatty acid metabolism. Recent reports suggest that modification and expression of pyruvate kinase (PK), a crucial regulatory enzyme involved in glycolysis is one of the plausible ways to alter seed oil content in plants. In the present study we evaluated the kinetic behavior and expression profiling of pyruvate kinase, associated with seed development in a major oilseed crop B. juncea. Developmental profiling of the enzyme showed that enzyme activity was highest during middle stage (35 DAF) of seed development which is strongly corroborated by the expression profiling of the enzyme using RT-PCR approach. Oil accumulation pattern also correlated with the enzyme expression study. Comparative activity profiling from different tissues showed seedlings to have elevated activity than other tissues. For kinetic characterization, the enzyme was partially purified by 12.3 fold using DEAE-Sephadex column and showed a narrow pH optimum of 7.0. In presence of saturated substrate concentration, the enzyme exhibited hyperbolic kinetics for both ADP and PEP with Km (Michaelis constant value) for PEP and ADP was found to be 178.5 and 96.45 μM respectively. ATP and citrate are the most significant allosteric effectors of the partially purified PK. Study on isozymes of PK resulted in a single band.     

Properties of phosphoribulokinase from Thiobacillus neapolitanus

Partially purified preparations of ribulose-5-phosphate kinase (specific activity, 50 to 125 mumoles per min per mg of protein) were employed in a series of kinetic experiments in the presence of several concentrations of H(+), Mg(2+), adenosine triphosphate (ATP), and phosphoenolpyruvate (PEP). The pH optimum of the enzyme was found to be 7.9; at this pH and above, response of the enzyme to variations in ATP concentration was hyperbolic, exhibiting a K(m) of 7 x 10(-4)m ATP. At pH values below the optimum the response to ATP was sigmoidal, as it was throughout the entire pH range in the presence of PEP at a concentration greater than 5 x 10(-4)m. In the presence of PEP the pH optimum shifted to pH 8.4. In contrast, phosphoribulokinase from spinach exhibited hyperbolic responses throughout its pH range with no inhibition caused by PEP. Thiobacillus neapolitanus phosphoribulokinase was inhibited by PEP in a sigmoidal manner; however, in the presence of suboptimal concentrations of Mg(2+) the addition of PEP caused significant stimulation of activity. It is postulated that the enzyme consists of interacting subunits with several sites on the enzyme for binding ATP and with several separate sites binding PEP. It is suggested that PEP functions as a regulator of CO(2) fixation when the organism is under conditions of unlimited concentrations of substrate and CO(2).     

Structural insight into epothilones antitumor activity based on the conformational preferences and tubulin binding modes of epothilones A and B obtained from molecular dynamics simulations

Molecular dynamics simulations were employed to analyze the conformational preferences and binding modes of epothilones A and B as a source of structural information regarding the antitumor properties of these species. Our results suggest that the conformation of free and tubulin-bound epothilones is strongly influenced by the presence of a methyl group at C12 and that epothilones A and B exploit the binding cavity in a unique and different way. The binding sites of epothilones A and B share a common region of association (Leu215, Leu217, His227, Leu228, Ala231, Phe270, Gly360, and Leu361), but lead to different ligand–residue interactions. Average interaction energies predict a larger stabilization for the epothilone B–tubulin complex, which is mainly driven by the enhancement of the electrostatic component of ligand–residue interactions compared to the epothilone A–tubulin complex. These structural and energetic results can be useful to account for the activity difference between epothilones A and B, and to design more active and potent analogs that resemble the mechanism of action of epothilones against cancer cells.     

Role of the conserved phenylalanine 181 of NADPH-cytochrome P450 oxidoreductase in FMN binding and catalytic activity.

NADPH-cytochrome P450 oxidoreductase (P450 reductase, EC 1.6.2.4) is an essential component of the P450 monooxygenase complex and binds FMN, FAD, and NADPH cofactors. Residues Tyr140 and Tyr178 are known to be involved in FMN binding. A third aromatic side chain, Phe181, is also located in the proximity of the FMN ring and is highly conserved in FMN-binding proteins, suggesting an important functional role. This role has been investigated by site-directed mutagenesis. Substitution of Phe181 with leucine or glutamine decreased the cytochrome c reductase activity of the enzyme by approximately 50%. Ferricyanide reductase activity was unaffected, indicating that the FAD domain was unperturbed. The mutant FMN domains were expressed in Escherichia coli, and the redox potentials and binding energies of their complexes with FMN were determined. The affinity for FMN was decreased approximately 50-fold in the Leu181 and Gln181 mutants. Comparison of the binding energies of the wild-type and mutant enzymes in the three redox states of FMN suggests that Phe181 stabilizes the FMN-apoprotein complex. The amide 1H and 15N resonances of the Phe181Leu FMN domain were assigned; comparison of their chemical shifts with those of the wild-type domain indicated that the effect of the substitution on FMN affinity results from perturbation of two loops which form part of the FMN binding site. The results indicate that Phe181 cooperates with Tyr140 and Tyr178 to play a major role in the binding and stability of FMN.     

Reversible ligand-induced dissociation of a tryptophan-shift mutant of phosphofructokinase from Bacillus stearothermophilus

The biophysical properties of a tryptophan-shifted mutant of phosphofructokinase from Bacillus stearothermophilus (BsPFK) have been examined. The mutant, designated W179Y/Y164W, has kinetic and thermodynamic properties similar to the wild-type enzyme. A 2-fold decrease in kcat is observed, and the mutant displays a 3-fold smaller K(0.5) for the substrate, fructose-6-phosphate (Fru-6-P), as compared to the wild-type enzyme. The dissociation constant for the inhibitor, phospho(enol)pyruvate (PEP), increases 2-fold, and the coupling parameter, Q(ay), decreases 2-fold. This suggests that while the mutant displays a slightly decreased affinity for PEP, PEP is still an effective inhibitor once bound. The new position of the tryptophan in W179Y/Y164W is approximately 6 A from the Fru-6-P portion of the active site. A 25% decrease in fluorescence intensity is observed upon Fru-6-P binding, and an 80% decrease in fluorescence intensity is observed with PEP binding. In addition, the intrinsic fluorescence polarization increases from 0.327 +/- 0.001 to 0.353 +/- 0.001 upon Fru-6-P binding, but decreases to 0.290 +/- 0.001 when PEP binds. Most notably, the presence of PEP induces dissociation of the tetramer. Dissociation of the tetramer into dimers occurs along the active site interface and can be monitored by the loss in activity or the loss in tryptophan fluorescence that is observed when the enzyme is titrated with PEP. Activity can be protected or recovered by incubating the enzyme with Fru-6-P. Recovery of activity is enzyme concentration dependent, and the rate constant for association is 6.2 +/- 0.3 M(-1) x s(-1). Ultracentrifugation experiments revealed that in the absence of PEP the mutant enzyme exists in an equilibrium between the dimer and tetramer forms with a dissociation constant of 11.8 +/- 0.5 microM, while in the presence of PEP the enzyme exists in equilibrium between the dimer and monomer forms with a dissociation constant of 7.5 +/- 0.02 microM. A 3.1 A crystal structure of the mutant enzyme suggests that the amino acid substitutions have not dramatically altered the tertiary structure of the enzyme. While it is clear that wild-type BsPFK exists as a tetramer under these same conditions, these results suggest that quaternary structural changes probably play an important role in allosteric communication.     

Phosphoenolpyruvate metabolism in Teladorsagia circumcincta: A critical junction between aerobic and anaerobic metabolism

Nematodes which have adapted to an anaerobic lifestyle in their adult stages oxidise phosphoenolpyruvate (PEP) to oxaloacetate rather than pyruvate as the final product of glycolysis. This adaptation involves selective expression of the enzyme phosphoenolpyruvate carboxykinase (PEPCK), instead of pyruvate kinase (PK). However, such adaptation is not absolute in aerobic nematode species. We have examined the activity and kinetics of PEPCK and PK in larvae (L₃) and adults of Teladorsagia circumcincta, a parasite known to exhibit oxygen uptake. Results revealed that PK and PEPCK activity existed in both L₃s and adults. The enzymes had differing affinity for nucleotide diphosphates: while both can utilise GDP, only PK utilised ADP and only PEPCK utilised IDP. In both life cycle stages, enzymes showed similar affinity for PEP. PK activity was predominant in both stages, although activity of this enzyme was lower in adults. When combined, both the activity levels and the enzyme kinetics showed that pyruvate production is probably favoured in both L₃ and adult stages of T. circumcincta and suggest that metabolism of PEP to oxaloacetate is a minor metabolic pathway in this species.