Interaction of glycolytic enzymes with F-actin is suggested to be a mechanism for compartmentation of the glycolytic pathway. Earlier work demonstrates that muscle F-actin strongly binds glycolytic enzymes, allowing for the general conclusion that "actin binds enzymes", which may be a generalized phenomenon. By taking actin from a lower form, such as yeast, which is more deviant from muscle actin than other higher animal forms, the generality of glycolytic enzyme interactions with actin and the cytoskeleton can be tested and compared with higher eukaryotes, e.g., rabbit muscle. Cosedimentation of rabbit skeletal muscle and yeast F-actin with muscle fructose-1,6-bisphosphate aldolase (aldolase) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) followed by Scatchard analysis revealed a biphasic binding, indicating high- and low-affinity domains. Muscle aldolase and GAPDH showed low-affinity for binding yeast F-actin, presumably because of fewer acidic residues at the N-terminus of yeast actin; this difference in affinity is also seen in Brownian dynamics computer simulations. Yeast GAPDH and aldolase showed low-affinity binding to yeast actin, which suggests that actin-glycolytic enzyme interactions may also occur in yeast although with lower affinity than in higher eukaryotes. The cosedimentation results were supported by viscometry results that revealed significant cross-linking at lower concentrations of rabbit muscle enzymes than yeast enzymes. Brownian dynamics simulations of yeast and muscle aldolase and GAPDH with yeast and muscle actin compared the relative association free energy. Yeast aldolase did not specifically bind to either yeast or muscle actin. Yeast GAPDH did bind to yeast actin although with a much lower affinity than when binding muscle actin. The binding of yeast enzymes to yeast actin was much less site specific and showed much lower affinities than in the case with muscle enzymes and muscle actin.

Development of new antimalarial drugs continues to be of great importance due to the resistance of the malaria parasite to currently used drugs. Glycolytic enzymes have emerged as potential targets for the development of new drugs due to the reliance of the parasite on glycolysis for energy. In this study, molecular docking was used to study the binding of some quinoline-based drugs to the glycolytic enzyme lactate dehydrogenase. The docking studies identified two potential binding sites for each ligand, one of them being the cofactor-binding site. For all ligands studied, there was the comparable binding to the cofactor-binding site as well as the secondary binding site when the cofactor was absent. All ligands showed significantly lower binding affinity than NADH for the cofactor binding site. The alternative site was the site of preference when docking was done in the presence of the cofactor. While binding to the cofactor site may support other studies suggesting potential for competitive inhibition, the fact that the binding affinities of all the ligands are significantly lower than that for NADH in this site suggests that these ligands will be ineffective competitive inhibitors. The identification of an alternative binding site with comparable affinity that is not affected by the presence of the cofactor may suggest the possibility of non-competitive inhibition that requires further exploration.

The association of glycolytic enzymes with F-actin is proposed to be one mechanism by which these enzymes are compartmentalized, and, as a result, may possibly play important roles for: regulation of the glycolytic pathway, potential substrate channeling, and increasing glycolytic flux. Historically, in vitro experiments have shown that many enzyme/actin interactions are dependent upon ionic strength. Herein, Brownian dynamics (BD) examines how ionic strength impacts the energetics of the association of F-actin with the glycolytic enzymes: lactate dehydrogenase (LDH), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), fructose-1,6-bisphosphate aldolase (aldolase), and triose phosphate isomerase (TPI). The BD simulations are steered by electrostatics calculated by Poisson-Boltzmann theory. The BD results confirm experimental observations that the degree of association diminishes as ionic strength increases but also suggest that these interactions are significant at physiological ionic strengths. Furthermore, BD agrees with experiments that muscle LDH, aldolase and GPADH interact significantly with F-actin whereas TPI does not. BD indicates similarities in binding regions for aldolase and LDH among the different species investigated. Furthermore, the residues responsible for salt bridge formation in stable complexes persist as ionic strength increases. This suggests the importance of the residues determined for these binary complexes and specificity of the interactions. That these interactions are conserved across species, and there appears to be a general trend among the enzymes, support the importance of these enzyme-F-actin interactions in creating initial complexes critical for compartmentation.

Previous Brownian dynamics (BD) simulations (Ouporov IG, Knull HR and Thomasson KA 1999. Biophys. J. 76: 17-27) of complex formation between rabbit aldolase and F-actin have identified three lysine residues (K288, K293 and K341) on aldolase and acidic residues (DEDE) at the N-terminus of actin as important to binding. BD simulations of computer models of aldolase mutants with any of these lysine residues replaced by alanine show reduced binding energy; the greatest effect of a single substitution is for K341A, and replacement of all three lysines greatly reduces binding. BD simulations of wild-type rabbit aldolase vs altered F-actin show that binding is decreased if any one of the four N-terminal acidic residues is replaced by alanine and binding is greatly reduced if three or more of the N-terminal acidic residues are replaced; none of the four actin residues appear more critical for binding than the others.

Interactions of the glycolytic enzyme, fructose‐1,6‐bisphosphate aldolase (aldolase), with F‐actin may be one mechanism for the colocalization of glycolytic enzymes. Examination of these interactions in different animal species tests this hypothesis by observing whether binding sites are conserved across species. Brownian dynamics (BD) simulations provide descriptions of such protein–protein interactions with the muscle isoforms of zebra fish and human aldolase. The results are compared with previous results obtained for rabbit muscle and yeast. The aldolase binding groove previously determined in rabbit muscle is conserved in both the human and fish muscle isoforms. The nonspecific radial free energies of interaction are similar with fish being slightly weaker than human and rabbit: human, −2.27 ± 0.05 kcal/mol; rabbit, −2.0 ± 0.04 kcal/mol; and fish, −1.5 ± 0.03 kcal/mol. BD results show a large Boltzmann population of complexes formed around the A/D and B/C grooves of aldolase with the most feasible binding mode comprising two aldolase subunits to subdomain I region of the actin subunits. These results show that the location of the important residues and binding site for fish and human aldolase is very similar to that in rabbit and that in different animals the binding site is conserved. This suggests that the binding interaction between aldolase and F‐actin is general in animal muscles and is rendered possible and energetically favorable through the conservation of this binding site. © 2006 Wiley Periodicals, Inc. Biopolymers 85: 60–71, 2007.This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com