Kinetic properties of fructose bisphosphate aldolase from Trypanosoma brucei compared to aldolase from rabbit muscle and Staphylococcus aureus

DLTKcat: deep learning-based prediction of temperature-dependent enzyme turnover rates

The enzyme turnover rate, ${k}_{cat}$, quantifies enzyme kinetics by indicating the maximum efficiency of enzyme catalysis. Despite its importance, ${k}_{cat}$ values remain scarce in databases for most organisms, primarily because of the cost of experimental measurements. To predict ${k}_{cat}$ and account for its strong temperature dependence, DLTKcat was developed in this study and demonstrated superior performance (log10-scale root mean squared error = 0.88, R-squared = 0.66) than previously published models. Through two case studies, DLTKcat showed its ability to predict the effects of protein sequence mutations and temperature changes on ${k}_{cat}$ values. Although its quantitative accuracy is not high enough yet to model the responses of cellular metabolism to temperature changes, DLTKcat has the potential to eventually become a computational tool to describe the temperature dependence of biological systems.

Defining the molecular architecture, metal dependence, and distribution of metal-dependent class II sulfofructose-1-phosphate aldolases

Sulfoquinovose (SQ, 6-deoxy-6-sulfoglucose) is a sulfosugar that is the anionic head group of plant, algal, and cyanobacterial sulfolipids: sulfoquinovosyl diacylglycerols. SQ is produced within photosynthetic tissues, forms a major terrestrial reservoir of biosulfur, and is an important species within the biogeochemical sulfur cycle. A major pathway for SQ breakdown is the sulfoglycolytic Embden-Meyerhof-Parnas pathway, which involves cleavage of the 6-carbon chain of the intermediate sulfofructose-1-phosphate (SFP) into dihydroxyacetone and sulfolactaldehyde, catalyzed by class I or II SFP aldolases. While the molecular basis of catalysis is understood for class I SFP aldolases, comparatively little is known about class II SFP aldolases. Here, we report the molecular architecture and biochemical basis of catalysis of two metal-dependent class II SFP aldolases from Hafnia paralvei and Yersinia aldovae. 3D X-ray structures of complexes with substrate SFP and product dihydroxyacetone phosphate reveal a dimer-of-dimers (tetrameric) assembly, the sulfonate-binding pocket, two metal-binding sites, and flexible loops that are implicated in catalysis. Both enzymes were metal-dependent and exhibited high KM values for SFP, consistent with their role in a unidirectional nutrient acquisition pathway. Bioinformatic analysis identified a range of sulfoglycolytic Embden-Meyerhof-Parnas gene clusters containing class I/II SFP aldolases. The class I and II SFP aldolases have mututally exclusive occurrence within Actinobacteria and Firmicutes phyla, respectively, while both classes of enzyme occur within Proteobacteria. This work emphasizes the importance of SQ as a nutrient for diverse bacterial phyla and the different chemical strategies they use to harvest carbon from this sulfosugar.

Oxidative Stress-Mediated Antibacterial Activity of the Total Flavonoid Extracted from the Agrimonia pilosa Ledeb. in Methicillin-Resistant Staphylococcusaureus (MRSA)

(1) Background: Methicillin-resistant Staphylococcus aureus (MRSA) is a zoonotic pathogen that causes endocarditis, pneumonia, and skin diseases in humans and livestock. (2) Methods: The antibacterial effect of the total flavonoid against MRSA (ATCC43300) extracted from the Agrimonia pilosa Ledeb. (A. pilosa Ledeb) was evaluated by the microdilution method. The oxidative stresses in MRSA were evaluated by the levels of intracellular hydrogen peroxide (H₂O₂), reactive oxygen species (ROS), and oxidative stress-related genes. The DNA oxidative damage was tested by the 8-hydroxy-2′-deoxyguanosine (8-OHdG) and DNA gel electrophoresis. The differentially expressed proteins were determined by the method of SDS-PAGE and NanoLC-ESI-MS/MS, while the mRNAs of differential proteins were determined by Real-Time PCR. The changes of ultra-structures in MRSA were observed by Transmission Electron Microscope (TEM). (3) Results: The minimum inhibitory concentration (MIC) of the total flavonoid against MRSA was recorded as 62.5 μg/mL. After treatment with the total flavonoid, the levels of intracellular H₂O₂ and ROS were increased and the gene expressions against oxidative stress (SodA, katA, TrxB) were decreased (p < 0.01), while the gene expression for oxidative stress (PerR) was increased (p < 0.01). The level of intracellular 8-OHdG in MRSA was increased (p < 0.01) and the DNA was damaged. The results of TEM also showed that the total flavonoid could destroy the ultra-structures in the bacteria. (4) Conclusions: The total flavonoid extracted from the A. pilosa Ledeb can induce the oxidative stress that disturbed the energy metabolism and protein synthesis in MRSA.

Substitutions at a rheostat position in human aldolase A cause a shift in the conformational population

Some protein positions play special roles in determining the magnitude of protein function: at such "rheostat" positions, varied amino acid substitutions give rise to a continuum of functional outcomes, from wild type (or enhanced), to intermediate, to loss of function. This observed range raises interesting questions about the biophysical bases by which changes at single positions have such varied outcomes. Here, we assessed variants at position 98 in human aldolase A ("I98X"). Despite being ~17 Å from the active site and far from subunit interfaces, substitutions at position 98 have rheostatic contributions to the apparent cooperativity (n_H ) associated with fructose-1,6-bisphosphate substrate binding and moderately affected binding affinity. Next, we crystallized representative I98X variants to assess structural consequences. Residues smaller than the native isoleucine (cysteine and serine) were readily accommodated, and the larger phenylalanine caused only a slight separation of the two parallel helixes. However, the diffraction quality was reduced for I98F, and further reduced for I98Y. Intriguingly, the resolutions of the I98X structures correlated with their n_H values. We propose that substitution effects on both n_H and crystal lattice disruption arise from changes in the population of aldolase A conformations in solution. In combination with results computed for rheostat positions in other proteins, the results from this study suggest that rheostat positions accommodate a wide range of side chains and that structural consequences manifest as shifted ensemble populations and/or dynamics changes.

Enhanced Diffusion of Passive Tracers in Active Enzyme Solutions

Colloidal suspensions containing microscopic swimmers have been the focus of recent studies aimed at understanding the principles of energy transfer in fluidic media at low Reynolds number conditions. Going down in scale, active enzymes have been shown to be force-generating, nonequilibrium systems, thus offering opportunity to examine energy transfer at the ultralow Reynolds number regime. By monitoring the change of diffusion of inert tracers dispersed in active enzyme solutions, we demonstrate that the nature of energy transfer in these systems is similar to that reported for larger microscopic active systems, despite the large differences in scale, modes of energy transduction, and propulsion. Additionally, even an enzyme that catalyzes an endothermic reaction behaves analogously, suggesting that heat generation is not the primary factor for the observed enhanced tracer diffusion. Our results provide new insights into the mechanism of energy transfer at the molecular level.

qPIPSA: relating enzymatic kinetic parameters and interaction fields

Background

The simulation of metabolic networks in quantitative systems biology requires the assignment of enzymatic kinetic parameters. Experimentally determined values are often not available and therefore computational methods to estimate these parameters are needed. It is possible to use the three-dimensional structure of an enzyme to perform simulations of a reaction and derive kinetic parameters. However, this is computationally demanding and requires detailed knowledge of the enzyme mechanism. We have therefore sought to develop a general, simple and computationally efficient procedure to relate protein structural information to enzymatic kinetic parameters that allows consistency between the kinetic and structural information to be checked and estimation of kinetic constants for structurally and mechanistically similar enzymes.

Results

We describe qPIPSA: quantitative Protein Interaction Property Similarity Analysis. In this analysis, molecular interaction fields, for example, electrostatic potentials, are computed from the enzyme structures. Differences in molecular interaction fields between enzymes are then related to the ratios of their kinetic parameters. This procedure can be used to estimate unknown kinetic parameters when enzyme structural information is available and kinetic parameters have been measured for related enzymes or were obtained under different conditions. The detailed interaction of the enzyme with substrate or cofactors is not modeled and is assumed to be similar for all the proteins compared. The protein structure modeling protocol employed ensures that differences between models reflect genuine differences between the protein sequences, rather than random fluctuations in protein structure.

Conclusion

Provided that the experimental conditions and the protein structural models refer to the same protein state or conformation, correlations between interaction fields and kinetic parameters can be established for sets of related enzymes. Outliers may arise due to variation in the importance of different contributions to the kinetic parameters, such as protein stability and conformational changes. The qPIPSA approach can assist in the validation as well as estimation of kinetic parameters, and provide insights into enzyme mechanism.

Simulated complex dynamics of glycolysis in the protozoan parasite Trypanosoma brucei

Glycolysis in Trypanosoma brucei was modeled using a reaction transport simulator and tested for possible complex dynamics. The glycolytic model is multi-compartmentalized and accounts for the exchange of metabolites between the glycosomes, cytosol, mitochondrion and the host medium. The model is used to examine the effects of a range of culture medium concentrations of oxygen on the glycolysis of T. brucei. Our results are in good agreement with steady-state experiments. We also find that under aerobic conditions, increasing the activity of glycerol-3-phosphate dehydrogenase induces complex dynamics in the system. We report the presence of three distinct types of these dynamics. Varying the oxygen concentration in the medium can induce the transition between these dynamics.

Can yeast glycolysis be understood in terms of in vitro kinetics of the constituent enzymes? Testing biochemistry

This paper examines whether the in vivo behavior of yeast glycolysis can be understood in terms of the in vitro kinetic properties of the constituent enzymes. In nongrowing, anaerobic, compressed Saccharomyces cerevisiae the values of the kinetic parameters of most glycolytic enzymes were determined. For the other enzymes appropriate literature values were collected. By inserting these values into a kinetic model for glycolysis, fluxes and metabolites were calculated. Under the same conditions fluxes and metabolite levels were measured. In our first model, branch reactions were ignored. This model failed to reach the stable steady state that was observed in the experimental flux measurements. Introduction of branches towards trehalose, glycogen, glycerol and succinate did allow such a steady state. The predictions of this branched model were compared with the empirical behavior. Half of the enzymes matched their predicted flux in vivo within a factor of 2. For the other enzymes it was calculated what deviation between in vivo and in vitro kinetic characteristics could explain the discrepancy between in vitro rate and in vivo flux.

The length of the combined 3' untranslated region and poly(A) tail does not control rates of glyceraldehyde-3-phosphate dehydrogenase mRNA translation in three species of parasitic protists

Experimental observations suggested that the length of the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mRNA 3' end has a role in regulating rates of translation in the parasitic protists Trypanosoma brucei, Leishmania donovani, and Trichomonas vaginalis. Using a PCR assay for poly(A) tail length, we measured the size of the RNA 3' end under different growth conditions in all three species. Our results showed that the combined 3' untranslated region and poly(A) tail of GAPDH mRNA do not vary with different rates of translation.

Regulation and adaptation of glucose metabolism of the parasitic protist Leishmania donovani at the enzyme and mRNA levels

Adaptation of the glucose metabolism of Leishmania donovani promastigotes (insect stage) was investigated by simultaneously measuring metabolic rates, enzyme activities, message levels, and cellular parameters under various conditions. Chemostats were used to adapt cells to different growth rates with growth rate-limiting or excess glucose concentrations. L. donovani catabolized glucose to CO(2), succinate, acetate, and pyruvate in ratios that depended on growth rate and glucose availability. Rates of glucose consumption were a linear function of growth rate and were twice as high in excess glucose-grown cells as in glucose-limited organisms. The major end product was CO(2), but organic end products were also formed in ratios that varied strongly with growth conditions. The specific activities of the 14 metabolic enzymes measured varied by factors of 3 to 17. Two groups of enzymes adapted specific activities in parallel, but there was no correlation between the groups. The activities of only one group correlated with specific rates of glucose metabolism. Total RNA content per cellular protein varied by a factor of 6 and showed a linear relationship with the rate of glucose consumption. There was no correlation between steady-state message levels and activities of the corresponding enzymes, suggesting regulation at the posttranscriptional level. A comparison of the adaptation of energy metabolism in L. donovani and other species suggests that the energy metabolism of L. donovani is inefficient but is well suited to the environmental challenges that it encounters during residence in the sandfly, its insect vector.

Molecular evolution of amphioxus fructose-1,6-bisphosphate aldolase

The cDNA for amphioxus fructose-1,6-bisphosphate (FBP)-aldolase was isolated and its nucleotide sequence was determined. In the cDNA, there existed a probable open reading frame comprising 1080 bp; hence, 359 amino acid residues were deduced. The amino acid sequence indicates the deletion of 4 residues from N-terminus, in comparison with the sequence of FBP-aldolase isozymes from other sources. There was only one FBP-aldolase gene, and one enzyme species corresponding, in the amphioxus; this is the first report of the existence of a single FBP-aldolase species in animals. Enzymatic studies of both native and the recombinant FBP-aldolase suggest that the amphioxus enzyme belongs to an ancestral class I type which is not discovered among vertebrate aldolase isozymes.

Partial purification and characterization of a murine malaria parasite, Plasmodium berghei specific aldolase

Cell-free P. berghei contains 26.1 times more aldolase activity as compared to normal mouse erythrocytes. Subcellular fractionation of cell-free parasite showed maximum enzyme activity in the soluble fraction. The parasite enzyme was active in a narrow pH range of 7.8-8.0. Of the enzyme activity 90% was lost within 2 weeks at 4 degrees C. Slight inhibition was observed with specific inhibitors ATP, pyrophosphate (PPi) and PEP. The F1, 6DP Km was 0.025 mM.

Lamprey fructose-1,6-bisphosphate aldolase: characterization of the muscle-type and non-muscle-type isozymes

To study evolutionary aspects of fructose-1,6-bisphosphate (Fru-1,6-P2) aldolase during deuterostomian evolution, we have purified and characterized aldolases from the muscle and liver of lamprey (Entosphenus japonicus). Aldolase from the skeletal muscle and liver was identified to be the muscle-type isozyme and the non-muscle-type isozyme that was encoded by cDNAs M8 and L3, respectively, as described previously (Zhang, R., Yatsuki, H., Kusakabe, T., Iwabe, Miyata, T., Imai, T., Yoshida, M., and Hori, K., J. Biochem. 117, 545-553, 1995). The muscle-type isozyme has properties similar to vertebrate aldolase A, while the non-muscle-type isozyme shows a similarity to bacterial class I aldolase and vertebrate aldolase C but not to aldolase B, the liver-type aldolase, in terms of kinetic parameters: the Kcat values toward Fru-1,6-P2 and Fru-1-P, the Fru-1,6-P2/Fru-1-P activity ratio, and the Km values toward these substrates. The two enzymes have tetrameric forms with a molecular mass of approximately 160,000 and have similar pH optimum. The muscle-type and non-muscle-type isozymes from the tissues show different electrophoretic mobility; the muscle-type isozyme moves much faster than the non-muscle-type isozyme toward anodic side. The recombinant muscle-type and non-muscle-type aldolases gave similar characteristics as those from the tissues. The results presented in this paper, together with the data presented in the previous paper, strongly suggest that in lamprey it is possible to have two types of aldolase isozymes rather than one or three isozymes.

Caenorhabditis elegans has two isozymic forms, CE-1 and CE-2, of fructose-1,6-bisphosphate aldolase which are encoded by different genes

Two distinct types of cDNAs for fructose-1,6-bisphosphate (FBP) aldolase, Ce-1 and Ce-2, have been isolated from nematode Caenorhabditis elegans, and the respective recombinant aldolase isozymes, CE-1 and CE-2, have been purified and characterized. The Ce-1 and Ce-2 are 1282 and 1248 bp in total length, respectively, and both have an open reading frame of 1098 bp, which encodes 366 amino acid residues. The entire amino acid sequences deduced from Ce-1 and Ce-2 show a high degree of identity to one another and to those of vertebrate and invertebrate aldolases. The highest sequence diversity was found in the carboxyl-terminal region that corresponds to one of the isozyme group-specific sequences of vertebrate aldolase isozymes that play a role in determining isozyme-specific functions. Southern blot analysis suggests that CE-1 and CE-2 are encoded by different genes. Concerning general or kinetic properties, CE-2 is quite different from CE-1. CE-1 exhibits unique characteristics which are not identical to any aldolase isozymes previously reported, whereas CE-2 is similar to vertebrate aldolase C. These results suggest that CE-2 might preserve the properties of a progenitor aldolase with a moderate preference for FBP over fructose 1-phosphate (F1P) as a substrate, whereas CE-1 evolved to act as an intrinsic enzyme that exhibits a much broader substrate specificity than dose CE-2.

Glycolysis in bloodstream form Trypanosoma brucei can be understood in terms of the kinetics of the glycolytic enzymes

In trypanosomes the first part of glycolysis takes place in specialized microbodies, the glycosomes. Most glycolytic enzymes of Trypanosoma brucei have been purified and characterized kinetically. In this paper a mathematical model of glycolysis in the bloodstream form of this organism is developed on the basis of all available kinetic data. The fluxes and the cytosolic metabolite concentrations as predicted by the model were in accordance with available data as measured in non-growing trypanosomes, both under aerobic and under anaerobic conditions. The model also reproduced the inhibition of anaerobic glycolysis by glycerol, although the amount of glycerol needed to inhibit glycolysis completely was lower than experimentally determined. At low extracellular glucose concentrations the intracellular glucose concentration remained very low, and only at 5 mM of extracellular glucose, free glucose started to accumulate intracellularly, in close agreement with experimental observations. This biphasic relation could be related to the large difference between the affinities of the glucose transporter and hexokinase for intracellular glucose. The calculated intraglycosomal metabolite concentrations demonstrated that enzymes that have been shown to be near-equilibrium in the cytosol must work far from equilibrium in the glycosome in order to maintain the high glycolytic flux in the latter.

Regulation and control of compartmentalized glycolysis in bloodstream form Trypanosoma brucei

Unlike other eukaryotic cells, trypanosomes possess a compartmentalized glycolytic pathway. The conversion of glucose into 3-phosphoglycerate takes place in specialized peroxisomes, called glycosomes. Further conversion of this intermediate into pyruvate occurs in the cytosol. Due to this compartmentation, many regulatory mechanisms operating in other cell types cannot work in trypanosomes. This is reflected by the insensitivity of the glycosomal enzymes to compounds that act as activity regulators in other cell types. Several speculations have been raised about the function of compartmentation of glycolysis in trypanosomes. We calculate that even in a noncompartmentalized trypanosome the flux through glycolysis should not be limited by diffusion. Therefore, the sequestration of glycolytic enzymes in an organelle may not serve to overcome a diffusion limitation. We also search the available data for a possible relation between compartmentation and the distribution of control of the glycolytic flux among the glycolytic enzymes. Under physiological conditions, the rate of glycolytic ATP production in the bloodstream form of the parasite is possibly controlled by the oxygen tension, but not by the glucose concentration. Within the framework of Metabolic Control Analysis, we discuss evidence that glucose transport, although it does not qualify as the sole rate-limiting step, does have a high flux control coefficient. This, however, does not distinguish trypanosomes from other eukaryotic cell types without glycosomes.

Glucose metabolism in Escherichia coli and the effect of increased amount of aldolase

We present a comparative study of Escherichia coli with normal and increased amounts of fructose-1,6-bisphosphate aldolase. Most experiments employed a resting cell system involving a high cell density (so as to obtain the soluble pool by direct extraction) and anaerobic incubation in the presence of chloramphenicol. Glucose use is linear with time with a rate ca. half of that in growth, fermentation is almost quantitative, and metabolite concentrations reach a quasi steady state. Increased amount of aldolase had little effect on glucose flux; fructose-1,6-P2 concentration decreased by ca. one-third, and the extent of equilibration of its two halves, measured by a dismutation procedure on samples taken during metabolism of [6-14C]glucose, increased from 0.33 [(cpm in C1-3)/(cpm in C1-6)] to 0.43. Using the simplest model, that increased amount of aldolase does not perturb net flux or later metabolites, together with the steady-state rate equations for aldolase and triose-P isomerase, we show that the results with resting cells fit with the extra enzyme being fully active, and do not necessitate special assumptions concerning a glycolytic complex, metabolite compartmentation, or secondary mechanisms assuring high metabolite concentration. However, the fit does require that the measured Vmax values substantially underestimate the actual ones. Calculation also shows that the forms of the predicted curves--and hence the fit with experimental data--of fructose-1,6-P2 concentration and labeling as a function of the amount of aldolase are highly dependent on glyceraldehyde-3-P concentration but independent of the kinetic parameters of aldolase.

Pyruvate kinase from Trichomonas vaginalis, an allosteric enzyme stimulated by ribose 5-phosphate and glycerate 3-phosphate

Trichomonas vaginalis pyruvate kinase was purified over 1750 fold to a specific activity greater than 100 mumol min-1 (mg protein)-1. The enzyme is a tetramer of M(r) 266,000, consisting of subunits of M(r) 53,000 and 56,000 in equivalent amounts. Its activity was dependent on the presence of magnesium but was not stimulated by potassium or ammonium. The enzyme exhibited positive cooperativity towards phosphoenolpyruvate and was inhibited by inorganic phosphate, which increased the sigmoidicity of the saturation curve for phosphoenolpyruvate without affecting maximal activity. It was heterotropically stimulated by ribose 5-phosphate and glycerate 3-phosphate, not previously known to act on eukaryotic pyruvate kinases, but was unaffected by known effectors of most pyruvate kinases, including fructose 1,6-bisphosphate and fructose 2,6-bisphosphate.

Inhibition of triosephosphate isomerase from Trypanosoma brucei with cyclic hexapeptides

Two series of oligopeptides have been synthesized. Their effects on the activity of purified triosephosphate isomerase from Trypanosoma brucei and various other organisms have been studied. Using detailed three-dimensional structure information, the first series consisted of both cyclic and linear hydrophilic peptides that were designed to mimic the beta turns of the subunit interface loops of the trypanosome triosephosphate isomerase dimer. None of these exerted any inhibitory effect. The second series consisted of more hydrophobic cyclic peptides, originally designed to inhibit a hepatic transport system. Several of these were very effective in inhibiting the trypanosome triosephosphate isomerase, but not the homologous enzymes from rabbit, dog, yeast or Escherichia coli. The most active peptide, cyclo[-Trp-Phe-D-Pro-Phe-Phe-Lys(Z)-], exerted 50% inhibitory activity at a concentration of 3 microM. The nature of the inhibitory action of one of these compounds cyclo[-Trp-Tyr(OSO3Na)-D-Pro-Phe-Thr(OSO3Na)-Lys(Z)-] was studied in more detail. Its inhibition was noncompetitive and reversible and more than one peptide was able to bind/active site.

Chemical modification of fructose bisphosphate aldolase from Trypanosoma brucei compared to aldolase from rabbit muscle and Staphylococcus aureus

Chemical modifications of Class I aldolases from Trypanosoma brucei, rabbit muscle and Staphylococcus aureus with carboxypeptidase A, glyceraldehyde 3-phosphate and cysteine-specific reagents revealed the following differences between the three homologous enzymes. Aldolase from S. aureus was not affected by any of these reagents. Carboxypeptidase-A treatment of rabbit-muscle and T. brucei aldolase inhibited the activity of both enzymes towards fructose-1,6-bisphosphate (Fru(1,6)P2), while the activity towards fructose-1-phosphate (Fru-1-P) was affected only in the case of the trypanosomal enzyme. Moreover carboxypeptidase-A treatment reduced the turnover numbers of these two aldolases for both Fru(1,6)P2 and Fru-1-P to a similar level. Glyceraldehyde 3-phosphate, in the absence of dihydroxyacetone phosphate, also inactivated aldolases from rabbit muscle and T. brucei with second order rate constants of 1054 and 254 min-1 M-1, respectively. Using 5,5'-dithiobis-(2-nitrobenzoic acid) with rabbit-muscle aldolase, a total of 4 thiol groups could be titrated per subunit, resulting in a total inactivation. The presence of substrate completely protected the enzyme from inactivation. Methyl methanethiosulfonate also reacted with four cysteine residues, but this led to very little inactivation. This indicates that the inactivation by modification with DTNB is due to conformational changes in the enzyme. In T. brucei aldolase only one thiol group could be titrated with methyl methanesulfonate and there was no loss of activity. With 5,5'-dithiobis-(2-nitrobenzoic acid) five cysteines were titrated with an immediate and complete loss of activity.(ABSTRACT TRUNCATED AT 250 WORDS)