Brownian dynamics of interactions between aldolase mutants and F-actin

Therapeutic Targeting of Aldolase A Interactions Inhibits Lung Cancer Metastasis and Prolongs Survival

Cancer metabolic reprogramming promotes tumorigenesis and metastasis; however, the underlying molecular mechanisms are still being uncovered. In this study, we show that the glycolytic enzyme aldolase A (ALDOA) is a key enzyme involved in lung cancer metabolic reprogramming and metastasis. Overexpression of ALDOA increased migration and invasion of lung cancer cell lines in vitro and formation of metastatic lung cancer foci in vivo. ALDOA promoted metastasis independent of its enzymatic activity. Immunoprecipitation and proteomic analyses revealed γ-actin binds to ALDOA; blocking this interaction using specific peptides decreased metastasis both in vitro and in vivo. Screening of clinically available drugs based on the crystal structure of ALDOA identified raltegravir, an antiretroviral agent that targets HIV integrase, as a pharmacologic inhibitor of ALDOA-γ-actin binding that produced antimetastatic and survival benefits in a xenograft model with no significant toxicity. In summary, ALDOA promotes lung cancer metastasis by interacting with γ-actin. Targeting this interaction provides a new therapeutic strategy to treat lung cancer metastasis. SIGNIFICANCE: This study demonstrates the role of aldolase A and its interaction with γ-actin in the metastasis of non-small lung cancer and that blocking this interaction could be an effective cancer treatment.

Diffusion coefficients of endogenous cytosolic proteins from rabbit skinned muscle fibers

Efflux time courses of endogenous cytosolic proteins were obtained from rabbit psoas muscle fibers skinned in oil and transferred to physiological salt solution. Proteins were separated by gel electrophoresis and compared to load-matched standards for quantitative analysis. A radial diffusion model incorporating the dissociation and dissipation of supramolecular complexes accounts for an initial lag and subsequent efflux of glycolytic and glycogenolytic enzymes. The model includes terms representing protein crowding, myofilament lattice hindrance, and binding to the cytomatrix. Optimization algorithms returned estimates of the apparent diffusion coefficients, D(r,t), that were very low at the onset of diffusion (∼10(-10) cm(2) s(-1)) but increased with time as cytosolic protein density, which was initially high, decreased. D(r,t) at later times ranged from 2.11 × 10(-7) cm(2) s(-1) (parvalbumin) to 0.20 × 10(-7) cm(2) s(-1) (phosphofructose kinase), values that are 3.6- to 12.3-fold lower than those predicted in bulk water. The low initial values are consistent with the presence of complexes in situ; the higher later values are consistent with molecular sieving and transient binding of dissociated proteins. Channeling of metabolic intermediates via enzyme complexes may enhance production of adenosine triphosphate at rates beyond that possible with randomly and/or sparsely distributed enzymes, thereby matching supply with demand.

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.

Molecular modeling of the cytoskeleton

Molecular modeling techniques have truly come of age in recent decades, and here we cover several of the most commonly used techniques, namely molecular dynamics, Brownian dynamics, and molecular docking. In each case, we explain the physical basis and limitations of the various techniques and then illustrate their application to various problems related to the cytoskeleton. This set of studies covers a relatively wide range of examples and is comprehensive enough to clearly see how these techniques could be applied to other systems. Finally, we cover several related methodologies that expand on these basic techniques to allow for more detailed and specific simulation and analysis.

Cofactor CLIM2 promotes the repressive action of LIM homeodomain transcription factor Lhx2 in the expression of porcine pituitary glycoprotein hormone alpha subunit gene

We have cloned a porcine orthologue of cofactor CLIM2 (Ldb1/NLI) from the porcine pituitary cDNA library by protein-protein interaction with the Yeast Two-Hybrid System using porcine Lhx2 as a bait protein. Porcine CLIM2 shows a high identity (99%) in the dimerization domain, nuclear localization signal and LIM binding domain with those of man and mouse. The expression of CLIM2 gene in the anterior pituitary lobe was detected during the porcine fetal and postnatal period by RT-PCR analysis, suggesting that this protein is constitutively expressing and plays a basic role in the anterior pituitary. Transfection assay to the pituitary tumor derived LbetaT2 cells, and the Chinese hamster ovary cells demonstrated that CLIM2 acts as a corepressor of the porcine Lhx2 function. Interestingly, CLIM2 alone apparently repressed the high level of alphaGSU gene expression in LbetaT2 cells. These data suggest that CLIM2 is a basic factor in the pituitary development and function, and plays the role of repressor to modify the function of Lhx2 on the alphaGSU gene expression.

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.

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