Ionic strength dependence of F-actin and glycolytic enzyme associations: a Brownian dynamics simulations approach.

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 on 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 GAPDH 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.

Theoretical study of interactions between muscle aldolase and F-actin: insight into different species.

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.

Glycolytic enzyme interactions with yeast and skeletal muscle F-actin.

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.

Brownian dynamics of interactions between glyceraldehyde-3-phosphate dehydrogenase (GAPDH) mutants and F-actin.

Brownian dynamics simulations of computer models of GAPDH mutants interacting with F-actin emphasized the electrostatic nature of such interactions, and confirmed the importance of four previously identified lysine residues on the GAPDH structure in these interactions. Mutants were GAPDH models in which one or more of the previously identified lysines had been replaced with alanine. Simulations showed reduced binding of these mutants to F-actin compared to wild-type GAPDH. Binding was significantly reduced by mutating the four lysines; the specific electrostatic interaction energy of the quadruple mutant was -7.3 +/- 1.0 compared to -11.4 +/- 0.5 kcal/mol for the wild enzyme. The BD simulations also reaffirmed the importance of quaternary structure for GAPDH binding F-actin.

Brownian dynamics of interactions between aldolase mutants and F-actin.

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.