Amyloid-Driven Allostery.

The fields of allostery and amyloid-related pathologies, such as Parkinson's disease (PD), have been extensively explored individually, but less is known about how amyloids control allostery. Recent advancements have revealed that amyloids can drive allosteric effects in both intrinsically disordered proteins, such as alpha-synuclein (αS), and multi-domain signaling proteins, such as protein kinase A (PKA). Amyloid-driven allostery plays a central role in explaining the mechanisms of gain-of-pathological-function mutations in αS (e.g. E46K, which causes early PD onset) and loss-of-physiological-function mutations in PKA (e.g. A211D, which predisposes to tumors). This review highlights allosteric effects of disease-related mutations and how they can cause exposure of amyloidogenic regions, leading to amyloids that are either toxic or cause aberrant signaling. We also discuss multiple potential modulators of these allosteric effects, such as MgATP and kinase substrates, opening future opportunities to improve current pharmacological interventions against αS and PKA-related pathologies. Overall, we show that amyloid-driven allosteric models are useful to explain the mechanisms underlying disease-related mutations.

Conserved regulatory motifs in the juxtamembrane domain and kinase N-lobe revealed through deep mutational scanning of the MET receptor tyrosine kinase domain.

MET is a receptor tyrosine kinase (RTK) responsible for initiating signaling pathways involved in development and wound repair. MET activation relies on ligand binding to the extracellular receptor, which prompts dimerization, intracellular phosphorylation, and recruitment of associated signaling proteins. Mutations, which are predominantly observed clinically in the intracellular juxtamembrane and kinase domains, can disrupt typical MET regulatory mechanisms. Understanding how juxtamembrane variants, such as exon 14 skipping (METΔEx14), and rare kinase domain mutations can increase signaling, often leading to cancer, remains a challenge. Here, we perform a parallel deep mutational scan (DMS) of the MET intracellular kinase domain in two fusion protein backgrounds: wild-type and METΔEx14. Our comparative approach has revealed a critical hydrophobic interaction between a juxtamembrane segment and the kinase ⍺C-helix, pointing to potential differences in regulatory mechanisms between MET and other RTKs. Additionally, we have uncovered a ß5 motif that acts as a structural pivot for the kinase domain in MET and other TAM family of kinases. We also describe a number of previously unknown activating mutations, aiding the effort to annotate driver, passenger, and drug resistance mutations in the MET kinase domain.

Free energy landscape of the PI3Kα C-terminal activation.

The gene PIK3CA, encoding the catalytic subunit p110α of PI3Kα, is the second most frequently mutated gene in cancer, with the highest frequency oncogenic mutants occurring in the C-terminus of the kinase domain. The C-terminus has a dual function in regulating the kinase, playing a putative auto-inhibitory role for kinase activity and being absolutely essential for binding to the cell membrane. However, the molecular mechanisms by which these C-terminal oncogenic mutations cause PI3Kα overactivation remain unclear. To understand how a spectrum of C-terminal mutations of PI3Kα alter kinase activity compared to the WT, we perform unbiased and biased Molecular Dynamics simulations of several C-terminal mutants and report the free energy landscapes for the C-terminal "closed-to-open" transition in the WT, H1047R, G1049R, M1043L and N1068KLKR mutants. Results are consistent with HDX-MS experimental data and provide a molecular explanation why H1047R and G1049R reorient the C-terminus with a different mechanism compared to the WT and M1043L and N1068KLKR mutants. Moreover, we show that in the H1047R mutant, the cavity, where the allosteric ligands STX-478 and RLY-2608 bind, is more accessible contrary to the WT. This study provides insights into the molecular mechanisms underlying activation of oncogenic PI3Kα by C-terminal mutations and represents a valuable resource for continued efforts in the development of mutant selective inhibitors as therapeutics.

Navigating the complexity of p53-DNA binding: implications for cancer therapy.

Abstract: The tumor suppressor protein p53, a transcription factor playing a key role in cancer prevention, interacts with DNA as its primary means of determining cell fate in the event of DNA damage. When it becomes mutated, it opens damaged cells to the possibility of reproducing unchecked, which can lead to formation of cancerous tumors. Despite its critical role, therapies at the molecular level to restore p53 native function remain elusive, due to its complex nature. Nevertheless, considerable information has been amassed, and new means of investigating the problem have become available. Objectives: We consider structural, biophysical, and bioinformatic insights and their implications for the role of direct and indirect readout and how they contribute to binding site recognition, particularly those of low consensus. We then pivot to consider advances in computational approaches to drug discovery. Materials and methods: We have conducted a review of recent literature pertinent to the p53 protein. Results: Considerable literature corroborates the idea that p53 is a complex allosteric protein that discriminates its binding sites not only via consensus sequence through direct H-bond contacts, but also a complex combination of factors involving the flexibility of the binding site. New computational methods have emerged capable of capturing such information, which can then be utilized as input to machine learning algorithms towards the goal of more intelligent and efficient de novo allosteric drug design. Conclusions: Recent improvements in machine learning coupled with graph theory and sector analysis hold promise for advances to more intelligently design allosteric effectors that may be able to restore native p53-DNA binding activity to mutant proteins. Clinical relevance: The ideas brought to light by this review constitute a significant advance that can be applied to ongoing biophysical studies of drugs for p53, paving the way for the continued development of new methodologies for allosteric drugs. Our discoveries hold promise to provide molecular therapeutics which restore p53 native activity, thereby offering new insights for cancer therapies. Graphical Abstract: Structural representation of the p53 DBD (PDBID 1TUP). DNA consensus sequence is shown in gray, and the protein is shown in blue. Red beads indicate hotspot residue mutations, green beads represent DNA interacting residues, and yellow beads represent both.

Reversible Control of Native GluN2B-Containing NMDA Receptors with Visible Light.

NMDA receptors (NMDARs) are glutamate-gated ion channels playing a central role in synaptic transmission and plasticity. NMDAR dysregulation is linked to various neuropsychiatric disorders. This is particularly true for GluN2B-containing NMDARs (GluN2B-NMDARs), which have major pro-cognitive, but also pro-excitotoxic roles, although their exact involvement in these processes remains debated. Traditional GluN2B-selective antagonists suffer from slow and irreversible effects, limiting their use in native tissues. We therefore developed OptoNAM-3, a photoswitchable negative allosteric modulator selective for GluN2B-NMDARs. OptoNAM-3 provided light-induced reversible inhibition of GluN2B-NMDAR activity with precise temporal control both in vitro and in vivo on the behavior of freely moving Xenopus tadpoles. When bound to GluN2B-NMDARs, OptoNAM-3 displayed remarkable red-shifting of its photoswitching properties allowing the use of blue light instead of UV light to turn-off its activity, which we attributed to geometric constraints imposed by the binding site onto the azobenzene moiety of the ligand. This study therefore highlights the importance of the binding site in shaping the photochemical properties of azobenzene-based photoswitches. In addition, by enabling selective, fast, and reversible photocontrol of native GluN2B-NMDARs with in vivo compatible photochemical properties (visible light), OptoNAM-3 should be a useful tool for the investigation of the GluN2B-NMDAR physiology in native tissues.

A platform for the early selection of non-competitive antibody-fragments from yeast surface display libraries.

In this work, we report the development of a platform for the early selection of non-competitive antibody-fragments against cell surface receptors that do not compete for binding of their natural ligand. For the isolation of such subtype of blocking antibody-fragments, we applied special fluorescence-activated cell sorting strategies for antibody fragments isolation from yeast surface display libraries. Given that most of the monoclonal antibodies approved on the market are blocking ligand-receptor interactions often leading to resistance and/or side effects, targeting allosteric sites represents a promising mechanism of action to open new avenues for treatment. To directly identify these antibody-fragments during library screening, we employed immune libraries targeting the epidermal growth factor receptor as proof of concept. Incorporating a labeled orthosteric ligand during library sorting enables the early selection of non-competitive binders and introduces an additional criterion to refine the selection of candidates exhibiting noteworthy properties. Furthermore, after sequencing, more candidates were identified compared to classical sorting based solely on target binding. Hence, this platform can significantly improve the drug discovery process by the early selection of more candidates with desired properties.

SecM leader peptide as an allosteric translation inhibitor: a molecular dynamics study.

The SecM leader peptide regulates translation of the SecA protein, being a part of the Sec translocase, that reversibly arrests the ribosome. In the present study the structure of the SecM complex with the E. coli A/A,P/P-ribosome was obtained by means of docking and molecular dynamics simulation methods. It has been established that binding of the SecM leader peptide in the nascent peptide exit tunnel leads to a turn of the aminoacylating proline residue away from the C-terminal SecM glycine residue, which is adverse to the peptidyltransferase reaction. Besides, the SecM binding leads to a disturbance of the A-tRNA contacts with the tip of the H38 helix of the 23S rRNA (the A-site finger, ASF) and ribosomal protein uL16. Allosteric interrelation between these events has been proved by a construction of networks of concerted changes in non-covalent interactions throughout the whole ribosome, whereupon the A1614 and A751 residues of the 23S rRNA in the exit tunnel that formed stacking interactions with the SecM residues during the MD simulations, were found to be the principal triggers, inducing crucial alterations in the A-tRNA binding. The allosteric signal from the SecM peptide to the ASF, according to our model, is transmitted through ribosomal protein uL22, and there is reason to believe that this sensor is used not only by the SecM leader peptide, but also by other peptides that cause translation arrest.

Dynamic and structural insights into allosteric regulation on MKP5 a dual-specificity phosphatase.

Dual-specificity mitogen-activated protein kinase (MAPK) phosphatases (MKPs) directly dephosphorylate and inactivate the MAPKs. Although the catalytic mechanism of dephosphorylation of the MAPKs by the MKPs is established, a complete molecular picture of the regulatory interplay between the MAPKs and MKPs still remains to be fully explored. Here, we sought to define the molecular mechanism of MKP5 regulation through an allosteric site within its catalytic domain. We demonstrate using crystallographic and NMR spectroscopy approaches that residue Y435 is required to maintain the structural integrity of the allosteric pocket. Along with molecular dynamics simulations, these data provide insight into how changes in the allosteric pocket propagate conformational flexibility in the surrounding loops to reorganize catalytically crucial residues in the active site. Furthermore, Y435 contributes to the interaction with p38 MAPK and JNK, thereby promoting dephosphorylation. Collectively, these results highlight the role of Y435 in the allosteric site as a novel mode of MKP5 regulation by p38 MAPK and JNK.

Dynamic allostery drives autocrine and paracrine TGF-ß signaling.

TGF-ß, essential for development and immunity, is expressed as a latent complex (L-TGF-ß) non-covalently associated with its prodomain and presented on immune cell surfaces by covalent association with GARP. Binding to integrin αvß8 activates L-TGF-ß1/GARP. The dogma is that mature TGF-ß must physically dissociate from L-TGF-ß1 for signaling to occur. Our previous studies discovered that αvß8-mediated TGF-ß autocrine signaling can occur without TGF-ß1 release from its latent form. Here, we show that mice engineered to express TGF-ß1 that cannot release from L-TGF-ß1 survive without early lethal tissue inflammation, unlike those with TGF-ß1 deficiency. Combining cryogenic electron microscopy with cell-based assays, we reveal a dynamic allosteric mechanism of autocrine TGF-ß1 signaling without release where αvß8 binding redistributes the intrinsic flexibility of L-TGF-ß1 to expose TGF-ß1 to its receptors. Dynamic allostery explains the TGF-ß3 latency/activation mechanism and why TGF-ß3 functions distinctly from TGF-ß1, suggesting that it broadly applies to other flexible cell surface receptor/ligand systems.

Insights into Ligand-Mediated Activation of an Oligomeric Ring-Shaped Gene-Regulatory Protein from Solution- and Solid-State NMR.

The 91 kDa oligomeric ring-shaped ligand binding protein TRAP (trp RNA binding attenuation protein) regulates the expression of a series of genes involved in tryptophan (Trp) biosynthesis in bacilli. When cellular Trp levels rise, the free amino acid binds to sites buried in the interfaces between each of the 11 (or 12, depending on the species) protomers in the ring. Crystal structures of Trp-bound TRAP show the Trp ligands are sequestered from solvent by a pair of loops from adjacent protomers that bury the bound ligand via polar contacts to several threonine residues. Binding of the Trp ligands occurs cooperatively, such that successive binding events occur with higher apparent affinity but the structural basis for this cooperativity is poorly understood. We used solution methyl-TROSY NMR relaxation experiments focused on threonine and isoleucine sidechains, as well as magic angle spinning solid-state NMR 13C-13C and 15N-13C chemical shift correlation spectra on uniformly labeled samples recorded at 800 and 1200 MHz, to characterize the structure and dynamics of the protein. Methyl 13C relaxation dispersion experiments on ligand-free apo TRAP revealed concerted exchange dynamics on the µs-ms time scale, consistent with transient sampling of conformations that could allow ligand binding. Cross-correlated relaxation experiments revealed widespread disorder on fast timescales. Chemical shifts for methyl-bearing side chains in apo- and Trp-bound TRAP revealed subtle changes in the distribution of sampled sidechain rotameric states. These observations reveal a pathway and mechanism for induced conformational changes to generate homotropic Trp-Trp binding cooperativity.

Discovery of potent allosteric antibodies inhibiting EGFR.

In this work, we report the discovery of potent anti-epidermal growth factor receptor (EGFR) allosteric heavy-chain antibodies by combining camelid immunization and fluorescence-activated cell sorting (FACS). After immunization and yeast surface display library construction, allosteric clones were obtained by introducing the labeled EGF Fc fusion protein as an additional criterion for FACS. This sorting method enabled the identification of 11 heavy-chain antibodies that did not compete with the orthosteric ligand EGF for the binding to EGFR. These antibodies bind to a triple-negative breast cancer cell line expressing EGFR with affinities in the picomolar to nanomolar range. Those camelid-derived antibodies also exhibit interesting properties by modulating EGFR affinity for EGF. Moreover, they are also able to inhibit EGF-induced downstream signaling pathways. In particular, we identified one clone that is more potent than the approved blocking antibody cetuximab in inhibiting both PI3K/AKT and MAPK/ERK pathways. Our results suggest that allosteric antibodies may be potential new modalities for therapeutics.

Cooperative Folding as a Molecular Switch in an Evolved Antibody Binder.

Designing proteins with tunable activities from easily accessible external cues remains a biotechnological challenge. Here, we set out to create a small antibody-binding domain equipped with a molecular switch inspired by the allosteric response to calcium seen in naturally derived proteins like calmodulin. We have focused on one of the three domains of Protein G that show inherent affinity to antibodies. By combining a semi-rational protein design with directed evolution, we engineered novel variants containing a calcium-binding loop rendering the inherent antibody affinity calcium-dependent. The evolved variants resulted from a designed selection strategy subjecting them to negative and positive selection pressures focused on conditional antibody-binding. Hence, these variants contained molecular "on/off" switches, controlling the target affinity towards antibody fragments simply by the presence or absence of calcium. From NMR spectroscopy we found that the molecular mechanism underlying the evolved switching behavior was a coupled calcium-binding and folding event where the target binding surface was intact and functional only in the presence of bound calcium. Notably, it was observed that the response to the employed selection pressures gave rise to the evolution of a cooperative folding mechanism. This observation illustrates why the cooperative folding reaction is an effective solution seen repeatedly in the natural evolution of fine-tuned macromolecular recognition. Engineering binding moieties to confer conditional target interaction has great potential due to the exquisite interaction control that is tunable to application requirements. Improved understanding of the molecular mechanisms behind regulated interactions is crucial to unlock how to engineer switchable proteins useful in a variety of biotechnological applications.

GPCR Signaling: A Study of the Interplay Between Structure, Energy, and Function.

G protein-coupled receptors (GPCRs) exemplify sophisticated allosteric communication, transducing extracellular signals through ligand-induced structural rearrangements that resonate through the molecular scaffold. Despite extensive study, the biophysical underpinnings of how conformational changes spread remain unclear. This work employs a novel physics-based framework to characterize the role of energy dissipation in directing intramolecular signaling pathways. By modeling each residue as a network of coupled oscillators, we generate a localization landscape depicting the vibrational energy distribution throughout the protein scaffold. Quantifying directional energy flux between residues reveals distinct pathways for energy and information transfer, illuminating sequences of allosteric communication. Our analysis of CB1 and CCR5 crystal structures unveils an anisotropic pattern of energy dissipation aligning with key functional dynamics, such as activation-related conformational changes. These anisotropic patterns of vibrational energy flow constitute pre-configured channels for allosteric signaling. Elucidating the relationship between structural topology and energy dissipation patterns provides key insights into the thermodynamic drivers of conformational signaling. This methodology significantly advances our mechanistic understanding of allostery in GPCRs and presents a broadly applicable approach for rationally dissecting allosteric communication pathways, with potential implications for structure-based drug design targeting these critical receptors.

Subtype-specific conformational landscape of NMDA receptor gating.

N-methyl-D-aspartate receptors are ionotropic glutamate receptors that mediate synaptic transmission and plasticity. Variable GluN2 subunits in diheterotetrameric receptors with identical GluN1 subunits set very different functional properties. To understand this diversity, we use single-molecule fluorescence resonance energy transfer (smFRET) to measure the conformations of the ligand binding domain and modulatory amino-terminal domain of the common GluN1 subunit in receptors with different GluN2 subunits. Our results demonstrate a strong influence of the GluN2 subunits on GluN1 rearrangements, both in non-agonized and partially agonized activation intermediates, which have been elusive to structural analysis, and in the fully liganded state. Chimeric analysis reveals structural determinants that contribute to these subtype differences. Our study provides a framework for understanding the conformational landscape that supports highly divergent levels of activity, desensitization, and agonist potency in receptors with different GluN2s and could open avenues for the development of subtype-specific modulators.

Modulation of circulating free testosterone fraction by testosterone, dihydrotestosterone, and estradiol during testosterone replacement therapy.

BACKGROUND: Testosterone, estradiol, and dihydrotestosterone share common ligand binding sites on sex hormone binding globulin and albumin. It is unknown whether and how changes in testosterone, dihydrotestosterone, and estradiol concentrations during testosterone replacement therapy affect free testosterone fraction. OBJECTIVE: To determine the effect of changes in testosterone, dihydrotestosterone, and estradiol concentrations on free testosterone fraction during testosterone replacement therapy of men with hypogonadism. METHODS: Using data from the Testosterone Trials, we assessed the association of changes in total testosterone, estradiol, and dihydrotestosterone concentrations over 12 months of testosterone replacement therapy with changes in free testosterone fraction, measured using equilibrium dialysis. We used random forests to evaluate the associations of predicted mean changes in free testosterone fraction with changes in circulating concentrations of each hormone at low, mean, or high change in the other two hormones. RESULTS: Testosterone replacement therapy not only increased total testosterone, dihydrotestosterone, estradiol, and free testosterone concentrations, but also the percent free testosterone, even though sex hormone binding globulin levels did not change. The predicted changes in free testosterone fraction during testosterone replacement therapy exhibited a non-linear relationship with changes in each of total testosterone, dihydrotestosterone, and estradiol concentrations. Greater increases in testosterone, dihydrotestosterone, and estradiol levels during testosterone replacement therapy were each associated with higher model-predicted percent free testosterone. Substantially smaller changes in molar concentrations of estradiol and dihydrotestosterone had a greater effect on percent free testosterone than those in testosterone. CONCLUSION: During testosterone replacement therapy of men with hypogonadism, changes in testosterone, dihydrotestosterone, and estradiol concentrations each altered percent free testosterone non-linearly. Small changes in estradiol concentrations exerted much larger effect on the free testosterone fraction than testosterone and dihydrotestosterone, suggesting complex interactions of the three hormones with the binding proteins. Assessment of changes in free testosterone during testosterone replacement therapy should include consideration of changes in all three hormones.

Revealing the principles of inter- and intra-domain regulation in a signaling enzyme via scanning mutagenesis.

Multi-domain enzymes can be regulated by both inter-domain interactions and structural features intrinsic to the catalytic domain. The tyrosine phosphatase SHP2 is a quintessential example of a multi-domain protein that is regulated by inter-domain interactions. This enzyme has a protein tyrosine phosphatase (PTP) domain and two phosphotyrosine-recognition domains (N-SH2 and C-SH2) that regulate phosphatase activity through autoinhibitory interactions. SHP2 is canonically activated by phosphoprotein binding to the SH2 domains, which causes large inter-domain rearrangements, but autoinhibition can also be disrupted by disease-associated mutations. Many details of the SHP2 activation mechanism are still unclear, the physiologically-relevant active conformations remain elusive, and hundreds of human variants of SHP2 have not been functionally characterized. Here, we perform deep mutational scanning on both full-length SHP2 and its isolated PTP domain to examine mutational effects on inter-domain regulation and catalytic activity. Our experiments provide a comprehensive map of SHP2 mutational sensitivity, both in the presence and absence of inter-domain regulation. Coupled with molecular dynamics simulations, our investigation reveals novel structural features that govern the stability of the autoinhibited and active states of SHP2. Our analysis also identifies key residues beyond the SHP2 active site that control PTP domain dynamics and intrinsic catalytic activity. This work expands our understanding of SHP2 regulation and provides new insights into SHP2 pathogenicity.

Taste Sensor Assessment of Bitterness in Medicines: Overview and Recent Topics.

In recent decades, taste sensors have been increasingly utilized to assess the taste of oral medicines, particularly focusing on bitterness, a major obstacle to patient acceptance and adherence. This objective and safe method holds promise for enhancing the development of patient-friendly medicines in pharmaceutical companies. This review article introduces its application in measuring the intensity of bitterness in medicine, confirming the achievement of taste masking, distinguishing taste differences between branded and generic medicines, and identifying substances to suppress bitterness in target medicines. Another application of the sensor is to predict a significant increase in bitterness when medicine is taken with certain foods/beverages or concomitant medication. Additionally, to verify the sensor's predictability, a significant correlation has been demonstrated between the output of a bitter-sensitive sensor designed for drug bitterness (BT0) and the bitterness responses of the human taste receptor hT2R14 from BitterDB (huji.ac.il). As a recent advancement, a novel taste sensor equipped with lipid/polymer membranes modified by 3-Br-2,6-dihydroxybenzoic acid (2,6-DHBA), based on the concept of allostery, is introduced. This sensor successfully predicts the bitterness of non-charged pharmaceuticals with xanthine skeletons, such as caffeine or related compounds. Finally, the future prospects of taste sensors are discussed.

Varying Selection Pressure for a Na+ Sensing Site in Epithelial Na+ Channel Subunits Reflect Divergent Roles in Na+ Homeostasis.

The epithelial Na+ channel (ENaC) emerged early in vertebrates and has played a role in Na+ and fluid homeostasis throughout vertebrate evolution. We previously showed that proteolytic activation of the channel evolved at the water-to-land transition of vertebrates. Sensitivity to extracellular Na+, known as Na+ self-inhibition, reduces ENaC function when Na+ concentrations are high and is a distinctive feature of the channel. A fourth ENaC subunit, δ, emerged in jawed fishes from an α subunit gene duplication. Here, we analyzed 849 α and δ subunit sequences and found that a key Asp in a postulated Na+ binding site was nearly always present in the α subunit, but frequently lost in the δ subunit (e.g. human). Analysis of site evolution and codon substitution rates provide evidence that the ancestral α subunit had the site and that purifying selection for the site relaxed in the δ subunit after its divergence from the α subunit, coinciding with a loss of δ subunit expression in renal tissues. We also show that the proposed Na+ binding site in the α subunit is a bona fide site by conferring novel function to channels comprising human δ subunits. Together, our findings provide evidence that ENaC Na+ self-inhibition improves fitness through its role in Na+ homeostasis in vertebrates.

Allostery at a Protein-Protein Interface Harboring an Intermolecular Motional Network.

Motional properties of proteins govern recognition, catalysis, and regulation. The dynamics of tightly interacting residues can form intramolecular dynamic networks, dependencies fine-tuned by evolution to optimize a plethora of functional aspects. The constructive interaction of residues from different proteins to assemble intermolecular dynamic networks is a similarly likely case but has escaped thorough experimental assessment due to interfering association/dissociation dynamics. Here, we use fast-MAS solid-state 15N R1ρ NMR relaxation dispersion aided by molecular-dynamics simulations to mechanistically assess the hierarchy of individual µs timescale motions arising from a crystal-crystal contact, in the absence of translational motion. In contrast to the monomer, where particular mutations entail isolated perturbations, specific intermolecular interactions couple the motional properties between distant residues in the same protein. The mechanistic insights obtained from this conceptual work may improve our understanding on how intramolecular allostery can be tuned by intermolecular interactions via assembly of dynamic networks from previously isolated elements.

Deciphering the interaction between PKM2 and the built-in thermodynamic properties of the glycolytic pathway in cancer cells.

Most cancer cells exhibit high glycolysis rates under conditions of abundant oxygen. Maintaining a stable glycolytic rate is critical for cancer cell growth as it ensures sufficient conversion of glucose carbons to energy, biosynthesis, and redox balance. Here we deciphered the interaction between PKM2 and the thermodynamic properties of the glycolytic pathway. Knocking down or knocking out PKM2 induced a thermodynamic equilibration in the glycolytic pathway, characterized by the reciprocal changes of the Gibbs free energy (ΔG) of the reactions catalyzed by PFK1 and PK, leading to a less exergonic PFK1-catalyzed reaction and a more exergonic PK-catalyzed reaction. The changes in the ΔGs of the two reactions cause the accumulation of intermediates, including the substrate PEP (the substrate of PK), in the segment between PFK1 and PK. The increased concentration of PEP in turn increased PK activity in the glycolytic pathway. Thus, the interaction between PKM2 and the thermodynamic properties of the glycolytic pathway maintains the reciprocal relationship between PK concentration and its substrate PEP concentration, by which, PK activity in the glycolytic pathway can be stabilized and effectively counteracts the effect of PKM2 KD or KO on glycolytic rate. In line with our previous reports, this study further validates the roles of the thermodynamics of the glycolytic pathway in stabilizing glycolysis in cancer cells. Deciphering the interaction between glycolytic enzymes and the thermodynamics of the glycolytic pathway will promote a better understanding of the flux control of glycolysis in cancer cells.