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

Brownian dynamics simulations of glycolytic enzyme subsets with F-actin

Previous Brownian dynamics (BD) simulations identified specific basic residues on fructose-1,6-bisphophate aldolase (aldolase) (I. V. Ouporov et al., Biophysical Journal, 1999, Vol. 76, pp. 17-27) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (I. V. Ouporov et al., Journal of Molecular Recognition, 2001, Vol. 14, pp. 29-41) involved in binding F-actin, and suggested that the quaternary structure of the enzymes may be important. Herein, BD simulations of F-actin binding by enzyme dimers or peptides matching particular sequences of the enzyme and the intact enzyme triose phosphate isomerase (TIM) are compared. BD confirms the experimental observation that TIM has little affinity for F-actin. For aldolase, the critical residues identified by BD are found in surface grooves, formed by subunits A/D and B/C, where they face like residues of the neighboring subunit enhancing their electrostatic potentials. BD simulations between F-actin and aldolase A/D dimers give results similar to the native tetramer. Aldolase A/B dimers form complexes involving residues that are buried in the native structure and are energetically weaker; these results support the importance of quaternary structure for aldolase. GAPDH, however, placed the critical residues on the corners of the tetramer so there is no enhancement of the electrostatic potential between the subunits. Simulations using GAPDH dimers composed of either S/H or G/H subunits show reduced binding energetics compared to the tetramer, but for both dimers, the sets of residues involved in binding are similar to those found for the native tetramer. BD simulations using either aldolase or GAPDH peptides that bind F-actin experimentally show complex formation. The GAPDH peptide bound to the same F-actin domain as did the intact tetramer; however, unlike the tetramer, the aldolase peptide lacked specificity for binding a single F-actin domain.     

Interference with myosin subfragment-1 binding by site-directed mutagenesis of actin

Three N-terminal double mutants of beta-actin expressed in the yeast Saccharomyces cerevisiae have been characterized with respect to DNase-I interaction, N-terminal post-translational modification, polymerizability and myosin subfragment-1 binding. The results strongly support earlier suggestions that the acidic residues at the N-terminus of actin are part of the myosin-binding site, while they seem to be of no importance for the other aspects of actin biochemistry tested. The suitability of this expression system for production of recombinant actin in general is discussed.     

Interactions of glyceraldehyde-3-phosphate dehydrogenase with G- and F-actin predicted by Brownian dynamics

Brownian dynamics (BD) was used to simulate the binding of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) to G- and F-actin. High-resolution three-dimensional models (X-ray and homology built) of the proteins were used in the simulations. The electrostatic potential about each protein was predicted by solving the linearized Poisson-Boltzmann equation for use in BD simulations. The BD simulations resulted in complexes of GAPDH with G- or F-actin involving positively charged surface patches on GAPDH (Lyses 24, 69, 110 and 114) and negatively charged residues of the N- and C-termini (Asps 1, 25 and 363 and Glus 2, 4, 224 and 364) of actin. The actin residues all belong to subdomain 1. Although the positively charged surface patches of GAPDH are not close enough to each other to enhance their electrostatic potential, occasionally two subunits of the GAPDH tetramer may simultaneously interact with two neighboring monomers of F-actin. These results are different from those of fructose-1,6-bisphosphate aldolase, where quaternary structure directly influenced binding by two subunits combining their electrostatic potentials (see previous study, Ouporov et al., 1999, Biophys. J. 76: 17-27). Instead, GAPDH uses its quaternary structure to span the distance between two different actin subunits so that it can interact with two different actin subunits simultaneously.     

Brownian dynamics simulations of aldolase binding glyceraldehyde 3-phosphate dehydrogenase and the possibility of substrate channeling

Brownian dynamics (BD) simulations test for channeling of the substrate, glyceraldehyde 3-phosphate (GAP), as it passes between the enzymes fructose-1,6-bisphosphate aldolase (aldolase) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH). First, BD simulations determined the favorable complexes between aldolase and GAPDH; two adjacent subunits of GAPDH form salt bridges with two subunits of aldolase. These intermolecular contacts provide a strong electrostatic interaction between the enzymes. Second, BD simulates GAP moving out of the active site of the A or D aldolase subunit and entering any of the four active sites of GAPDH. The efficiency of transfer is determined as the relative number of BD trajectories that reached any active site of GAPDH. The distribution functions of the transfer time were calculated based on the duration of successful trajectories. BD simulations of the GAP binding from solution to aldolase/GAPDH complex were compared to the channeling simulations. The efficiency of transfer of GAP within an aldolase/GAPDH complex was 2 to 3% compared to 1.3% when GAP was binding to GAPDH from solution. There is a preference for GAP channeling between aldolase and GAPDH when compared to binding from solution. However, this preference is not large enough to be considered as a theoretical proof of channeling between these proteins.     

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.     

Binding of estrogen and progesterone-BSA conjugates to glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and the effects of the free steroids on GAPDH enzyme activity: physiological implications

In this study rat brain solubilized plasmalemma-microsomal fractions (B-P3) or cytosolic fractions were applied to P-3-BSA (progesterone linked to BSA at C-3 position) and E-6-BSA (17beta-estradiol linked to BSA at C-6 position) affinity columns. It is interesting that a 37 kDa protein was retained by both columns which was identified as glyceraldehyde-3-phosphate dehydrogenase (GAPDH) by N-terminal sequencing. The 37 kDa protein (GAPDH) was not retained by either a control BSA conjugated affinity column or a corticosterone-BSA affinity column. E-6-BSA bound to GAPDH with higher binding affinity than P-3-BSA or T-3-BSA (testosterone linked to BSA at C-3 position) affinity columns. In addition, the binding of 17beta-E-6-BSA to GAPDH was impeded by free estrogen (17beta-estradiol) completely. Binding studies of E-6-BSA and P-3-BSA to commercial GAPDH from rabbit skeletal muscle using radiolabeled ligand binding assays revealed that P-3-BSA had 10x lower GAPDH binding affinity than E-6-BSA. Next, the effects of estrogen and progesterone on GAPDH activity were studied. Rapid and significant increases in V(max) and changes in K(m) were observed by the addition of 10 nM estradiol, whereas 100 nM progesterone decreased only V(max) significantly. Testosterone, corticosterone, 17alpha-estradiol, and diethylstilbestrol did not affect the enzyme activity. The results indicate that GAPDH is a target site for 17beta-estradiol and progesterone and suggest possible roles in the regulation of cellular metabolism and synaptic remodeling in which GAPDH has been reported to be involved.     

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.     

Evolutionary conserved N-terminal region of human muscle fructose 1,6-bisphosphatase regulates its activity and the interaction with aldolase

N-terminal residues of muscle fructose 1,6-bisphosphatase (FBPase) are highly conserved among vertebrates. In this article, we present evidence that the conservation is responsible for the unique properties of the muscle FBPase isozyme: high sensitivity to AMP and Ca(2+) inhibition and the high affinity to muscle aldolase, which is a factor desensitizing muscle FBPase toward AMP and Ca(2+). The first N-terminal residue affecting the affinity of muscle FBPase to aldolase is arginine 3. On the other hand, the first residue significantly influencing the kinetics of muscle FBPase is proline 5. Truncation from 5-7 N-terminal residues of the enzyme not only decreases its affinity to aldolase but also reduces its k-(cat) and activation by Mg(2+), and desensitizes FBPase to inhibition by AMP and calcium ions. Deletion of the first 10 amino acids of muscle FBPase abolishes cooperativity of Mg(2+) activation and results in biphasic inhibition of the enzyme by AMP. Moreover, this truncation lowers affinity of muscle FBPase to aldolase about 14 times, making it resemble the liver isozyme. We suggest that the existence of highly AMP-sensitive muscle-like FBPase, activity of which is regulated by metabolite-dependent interaction with aldolase enables the precise regulation of muscle energy expenditures and might contributed to the evolutionary success of vertebrates.     

Essential functions and actin-binding surfaces of yeast cofilin revealed by systematic mutagenesis.

Cofilin stimulates actin filament turnover in vivo. The phenotypes of twenty yeast cofilin mutants generated by systematic mutagenesis were determined. Ten grew as well as the wild type and showed no cytoskeleton defects, seven were recessive-lethal and three were conditional-lethal and caused severe actin organization defects. Biochemical characterization of interactions between nine mutant yeast cofilins and yeast actin provided evidence that F-actin binding and depolymerization are essential cofilin functions. Locating the mutated residues on the yeast cofilin molecular structure allowed several important conclusions to be drawn. First, residues required for actin monomer binding are proximal to each other. Secondly, additional residues are required for interactions with actin filaments; these residues might bind an adjacent subunit in the actin filament. Thirdly, despite striking structural similarity, cofilin interacts with actin in a different manner from gelsolin segment-1. Fourthly, a previously unrecognized cofilin function or interaction is suggested by identification of spatially proximal residues important for cofilin function in vivo, but not for actin interactions in vitro. Finally, mutation of the cofilin N-terminus suggests that its sequence is conserved because of its critical role in actin interactions, not because it is sometimes a target for protein kinases.     

Effects of sequential deletions of residues from the N- or C-terminus on the functions of epsilon subunit of the chloroplast ATP synthase

Ten truncated mutants of chloroplast ATP synthase epsilon subunit from spinach (Spinacia oleracea), which had sequentially lost 1-5 amino acid residues from the N-terminus and 6-10 residues from the C-terminus, were generated by PCR. These mutants were overexpressed in Escherichia coli, reconstituted with soluble and membrane-bound CF(1), and the ATPase activity and proton conductance of thylakoid membrane were examined. Deletions of as few as 3 amino acid residues from the N-terminus or 6 residues from the C-terminus of epsilon subunit significantly affected their ATPase-inhibitory activity in solution. Deletion of 5 residues from the N-terminus abolished its abilities to inhibit ATPase activity and to restore proton impermeability. Considering the consequence of interaction of epsilon and gamma subunit in the enzyme functions, the special interactions between the epsilon variants and the gamma subunit were detected in the yeast two-hybrid system and in vitro binding assay. In addition, the structures of these mutants were modeled through the SWISS-MODEL Protein Modeling Server. These results suggested that in chloroplast ATP synthase, both the N-terminus and C-terminus of the epsilon subunit show importance in regulation of the ATPase activity. Furthermore, the N-terminus of the epsilon subunit is more important for its interaction with gamma and some CF(o) subunits, and crucial for the blocking of proton leakage. Compared with the epsilon subunit from E. coli [Jounouchi, M., Takeyama, M., Noumi, T., Moriyama, Y., Maeda, M., and Futai, M. (1992) Arch. Biochem. Biophys. 292, 87-94; Kuki, M., Noumi, T., Maeda, M., Amemura, A., and Futai, M. (1988) J. Biol. Chem. 263, 4335-4340], the chloroplast epsilon subunit is more sensitive to N-terminal or C-terminal truncations.     

Functional consequences of N- or C-terminal deletions of the delta subunit of chloroplast ATP synthase

Mutagenesis was used to generate seven truncation mutants of the spinach (Spinacia oleracea) chloroplast ATP synthase delta subunit lacking 5, 11, 17, or 35 amino acid residues from the N-terminus or 3, 9, or 15 residues from the C-terminus. Interactions between these mutants and all other subunits of the chloroplast ATPase were investigated by a yeast two-hybrid system. The results indicate that the N-terminal deletions mainly affected interactions between the delta subunit and the other part of CF(1), but did not significantly affect interactions with the CF(0) sector. In contrast, C-terminal truncations of the delta subunit mainly affected its interaction with the CF(0) sector and caused little impairment in interactions with the other part of CF(1). The conformation of the delta subunit C-terminal domain seems to be more sensitive to the truncations, as shown by minimal expression driven by C-terminal deleted (nine residues) mutants. Further studies showed C-terminal truncations of the delta subunit greatly impaired its ability to restore cyclic photophosphorylation in NaBr vesicles, whereas N-terminal truncations had little effect on the ability of delta to plug the CF(0) channel. None of the mutants impaired ATP hydrolysis by CF(1).     

N-terminal amino acid sequence analysis of the subunits of pea photosystem I

Six 'core' subunits of pea photosystem I have been isolated and their N-terminal amino acid sequences determined by gas-phase or solid-phase sequencing. On average more than thirty residues were determined from the N-terminus of each polypeptide. This sequence analysis has revealed three polypeptides with charged N-terminal regions (21, 17 and 11 kDa subunits), one polypeptide with a predominantly hydrophobic N-terminal region (9 kDa subunit), one polypeptide which is cysteine-rich (8 kDa subunit) and one which is alanine-rich (13 kDa subunit).     

Theoretical prediction of the binding free energy for mutants of replication protein A

The replication protein A (RPA) is a heterotrimeric (70, 32, and 14 kDa subunits), single stranded DNA (ssDNA) binding protein required for pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Studies based on deletions and mutations have identified the high affinity ssDNA binding domains in the 70 kDa subunit of RPA, regions A and B. Individually, the domain A and B have a low affinity for ssDNA, while tandems composed of AA, AB, BB, and BA sequences bind the ssDNA with moderate to high affinity. Single and double point mutations on polar residues in the binding domains leads to a reduction in affinity of RPA for ssDNA, in particular when two hydrophilic residues are involved. In view of these results, we performed a study based on molecular dynamics simulation aimed to reproduce the experimental change in binding free energy, ΔΔG, of RPA70 mutants to further elucidate the nature of the protein-ssDNA interaction. The MM-PB(GB)SA methods implemented in Amber10 and the code FoldX were used to estimate the binding free energy. The theoretical and experimental ΔΔG values correlate better when the results are obtained by MM-PBSA calculated on individual trajectories for each mutant. In these conditions, the correlation coefficient between experimental and theoretical ΔΔG reaches a value of 0.95 despite the overestimation of the energy change by one order of magnitude. The decomposition of the MM-GBSA energy per residue allows us to correlate the change of the affinity with the residue polarity and energy contribution to the binding. The method revealed reliable predictions of the change in the affinity in function of mutations, and can be used to identify new mutants with distinct binding properties.     

Identification of intracellular target proteins of the calcium-signaling protein S100A12.

In this report, we have focused our attention on identifying intracellular mammalian proteins that bind S100A12 in a Ca2+-dependent manner. Using S100A12 affinity chromatography, we have identified cytosolic NADP+-dependent isocitrate dehydrogenase (IDH), fructose-1,6-bisphosphate aldolase A (aldolase), glyceraldehyde-3-phosphate dehydrogenese (GAPDH), annexin V, S100A9, and S100A12 itself as S100A12-binding proteins. Immunoprecipitation experiments indicated the formation of stable complexes between S100A12 and IDH, aldolase, GAPDH, annexin V and S100A9 in vivo. Surface plasmon resonance analysis showed that the binding to S100A12, of S100A12, S100A9 and annexin V, was strictly Ca2+-dependent, whereas that of GAPDH and IDH was only weakly Ca2+-dependent. To localize the site of S100A12 interaction, we examined the binding of a series of C-terminal truncation mutants to the S100A12-immobilized sensor chip. The results indicated that the S100A12-binding site on S100A12 itself is located at the C-terminus (residues 87-92). However, cross-linking experiments with the truncation mutants indicated that residues 87-92 were not essential for S100A12 dimerization. Thus, the interaction between S100A12 and S100A9 or immobilized S100A12 should not be viewed as a typical S100 homo- or heterodimerization model. Ca2+-dependent affinity chromatography revealed that C-terminal residues 75-92 are not necessary for the interaction of S100A12 with IDH, aldolase, GAPDH and annexin V. To analyze the functional properties of S100A12, we studied its action in protein folding reactions in vitro. The thermal aggregation of IDH or GAPDH was facilitated by S100A12 in the absence of Ca2+, whereas in the presence of Ca2+ the protein suppressed the aggregation of aldolase to less than 50%. These results suggest that S100A12 may have a chaperone/antichaperone-like function which is Ca2+-dependent.     

Effect of the substitution of muscle actin-specific subdomain 1 and 2 residues in yeast actin on actin function

Muscle and yeast actins display distinct behavioral characteristics. To better understand the allosteric interactions that regulate actin function, we created a muscle/yeast hybrid actin containing a muscle-specific outer domain (subdomains 1 and 2) and a yeast inner domain (subdomains 3 and 4). Actin with muscle subdomain 1 and the two yeast N-terminal negative charges supported viability. The four negative charge muscle N terminus in a muscle subdomain 1 background caused death, but in the same background actin with three N-terminal acidic residues (3Ac/Sub1) led to sick but viable cells. Addition of three muscle subdomain 2 residues (3Ac/Sub12) produced no further deleterious effects. These hybrid actins caused depolarized cytoskeletons, abnormal vacuoles, and mitochondrial and endocytosis defects. 3Ac/Sub1 G-actin exchanged bound epsilonATP more slowly than wild type actin, and the exchange rate for 3Ac/Sub12 was even slower, similar to that for muscle actin. The mutant actins polymerized faster and produced less stable and shorter filaments than yeast actin, the opposite of that expected for muscle actin. Unlike wild type actin, in the absence of unbound ATP, polymerization led to ADP-F-actin, which rapidly depolymerized. Like yeast actin, the hybrid actins activated muscle myosin S1 ATPase activity only about one-eighth as well as muscle actin, despite having essentially a muscle actin-specific myosin-binding site. Finally, the hybrid actins behaved abnormally in a yeast Arp2/3-dependent polymerization assay. Our results demonstrate a unique sensitivity of yeast to actin N-terminal negative charge density. They also provide insight into the role of each domain in the control of the various functions of actin.