skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on February 4, 2026

Title: Enzymes in a human cytoplasm model organize into submetabolon complexes
Enzyme–enzyme interactions are fundamental to the function of cells. Their atomistic mechanisms remain elusive mainly due to limitations of in-cell measurements. We address this challenge by atomistically modeling, for a total of ≈80 μs, a slice of the human cell cytoplasm that includes three successive enzymes along the glycolytic pathway: glyceraldehyde-3-phosphate dehydrogenase (GAPDH), phosphoglycerate kinase (PGK), and phosphoglycerate mutase (PGM). We tested the model for nonspecific protein stickiness, an artifact of current atomistic force fields in crowded environments. The simulations reveal that the human enzymes co-organize in-cell into transient submetabolon complexes, consistent with previous experimental results. Our data both reiterate known specificity between GAPDH and PGK and reveal extensive direct interactions between GAPDH and PGM. Our simulations further reveal, through force field benchmarking, the critical role of protein solvation in facilitating these enzyme–enzyme interactions. Transient interenzyme interactions with μs lifetime occur repeatedly in our simulations via specific sticky protein surface patches, with interactions often mediated by charged patch residues. Some of the residues that interact frequently with one another lie in or near the active site of the enzymes. We show that some of these patches correspond to a general mode to interact with several partners for promiscuous enzymes like GAPDH. We further show that the non-native yeast PGK is stickier than human PGK in our human cytoplasm model, supporting the idea of evolutionary pressure to reduce sticking. Our cytoplasm modeling paves the way toward capturing the atomistic dynamics of an entire enzymatic pathway in-cell.  more » « less
Award ID(s):
2243257
PAR ID:
10595355
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
122
Issue:
5
ISSN:
0027-8424
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Designing complex synthetic materials for enzyme immobilization could unlock the utility of biocatalysis in extreme environments. Inspired by biology, we investigate the use of random copolymer brushes as dynamic immobilization supports that enable supra-biological catalytic performance of immobilized enzymes. This is demonstrated by immobilizingBacillus subtilisLipase A on brushes doped with aromatic moieties, which can interact with the lipase through multiple non-covalent interactions. Incorporation of aromatic groups leads to a 50 °C increase in the optimal temperature of lipase, as well as a 50-fold enhancement in enzyme activity. Single-molecule FRET studies reveal that these supports act as biomimetic chaperones by promoting enzyme refolding and stabilizing the enzyme’s folded and catalytically active state. This effect is diminished when aromatic residues are mutated out, suggesting the importance of π-stacking and π-cation interactions for stabilization. Our results underscore how unexplored enzyme-support interactions may enable uncharted opportunities for using enzymes in industrial biotransformations. 
    more » « less
  2. Imam, Saheed (Ed.)
    ABSTRACT Upon nutrient starvation,Chlamydia trachomatisserovar L2 (CTL) shifts from its normal growth to a non-replicating form, termed persistence. It is unclear if persistence reflects an adaptive response or a lack thereof. To understand this, transcriptomics data were collected for CTL grown under nutrient-replete and nutrient-starved conditions. Applying K-means clustering on transcriptomics data revealed a global transcriptomic rewiring of CTL under stress conditions in the absence of any canonical global stress regulator. This is consistent with previous data that suggested that CTL’s stress response is due to a lack of an adaptive response mechanism. To investigate the impact of this on CTL metabolism, we reconstructed a genome-scale metabolic model of CTL (iCTL278) and contextualized it with the collected transcriptomics data. Using the metabolic bottleneck analysis on contextualized iCTL278, we observed that phosphoglycerate mutase (pgm) regulates the entry of CTL to the persistence state. Our data indicate thatpgmhas the highest thermodynamics driving force and lowest enzymatic cost. Furthermore, CRISPRi-driven knockdown ofpgmin the presence or absence of tryptophan revealed the importance of this gene in modulating persistence. Hence, this work, for the first time, introduces thermodynamics and enzyme cost as tools to gain a deeper understanding on CTL persistence. IMPORTANCEThis study uses a metabolic model to investigate factors that contribute to the persistence ofChlamydia trachomatisserovar L2 (CTL) under tryptophan and iron starvation conditions. As CTL lacks many canonical transcriptional regulators, the model was used to assess two prevailing hypotheses on persistence—that the chlamydial response to nutrient starvation represents a passive response due to the lack of regulators or that it is an active response by the bacterium. K-means clustering of stress-induced transcriptomics data revealed striking evidence in favor of the lack of adaptive (i.e., a passive) response. To find the metabolic signature of this, metabolic modeling pin-pointed pgm as a potential regulator of persistence. Thermodynamic driving force, enzyme cost, and CRISPRi knockdown of pgm supported this finding. Overall, this work introduces thermodynamic driving force and enzyme cost as a tool to understand chlamydial persistence, demonstrating how systems biology-guided CRISPRi can unravel complex bacterial phenomena. 
    more » « less
  3. Filament formation by metabolic, biosynthetic, and other enzymes has recently come into focus as a mechanism to fine-tune enzyme activity in the cell. Filamentation is key to the function of SgrAI, a sequence-specific DNA endonuclease that has served as a model system to provide some of the deepest insights into the biophysical characteristics of filamentation and its functional consequences. Structure-function analyses reveal that, in the filamentous state, SgrAI stabilizes an activated enzyme conformation that leads to accelerated DNA cleavage activity and expanded DNA sequence specificity. The latter is thought to be mediated by sequence-specific DNA structure, protein–DNA interactions, and a disorder-to-order transition in the protein, which collectively affect the relative stabilities of the inactive, non-filamentous conformation and the active, filamentous conformation of SgrAI bound to DNA. Full global kinetic modeling of the DNA cleavage pathway reveals a slow, rate-limiting, second-order association rate constant for filament assembly, and simulations of in vivo activity predict that filamentation is superior to non-filamenting mechanisms in ensuring rapid activation and sequestration of SgrAI's DNA cleavage activity on phage DNA and away from the host chromosome. In vivo studies demonstrate the critical requirement for accelerated DNA cleavage by SgrAI in its biological role to safeguard the bacterial host. Collectively, these data have advanced our understanding of how filamentation can regulate enzyme structure and function, while the experimental strategies used for SgrAI can be applied to other enzymatic systems to identify novel functional roles for filamentation. 
    more » « less
  4. Interactions between proteins lie at the heart of numerous biological processes and are essential for the proper functioning of the cell. Although the importance of hydrophobic residues in driving protein interactions is universally accepted, a characterization of protein hydrophobicity, which informs its interactions, has remained elusive. The challenge lies in capturing the collective response of the protein hydration waters to the nanoscale chemical and topographical protein patterns, which determine protein hydrophobicity. To address this challenge, here, we employ specialized molecular simulations wherein water molecules are systematically displaced from the protein hydration shell; by identifying protein regions that relinquish their waters more readily than others, we are then able to uncover the most hydrophobic protein patches. Surprisingly, such patches contain a large fraction of polar/charged atoms and have chemical compositions that are similar to the more hydrophilic protein patches. Importantly, we also find a striking correspondence between the most hydrophobic protein patches and regions that mediate protein interactions. Our work thus establishes a computational framework for characterizing the emergent hydrophobicity of amphiphilic solutes, such as proteins, which display nanoscale heterogeneity, and for uncovering their interaction interfaces. 
    more » « less
  5. Abstract Lipopolysaccharide (LPS) is essential in most Gram-negative bacteria, but mutants of several species have been isolated that can survive in its absence.Caulobacter crescentusviability in the absence of LPS is partially dependent on the anionic sphingolipid ceramide diphosphoglycerate (CPG2). Genetic analyses showed thatccna_01210, which encodes a nucleotidyltransferase, is required for CPG2 production. Using purified recombinant protein, we determined that CCNA_01210 (CpgD) is a phosphoglycerate cytidylyltransferase which uses CTP and 3-phosphoglycerate to produce CDP-glycerate, which we hypothesize is the phosphoglycerate donor for CPG2 synthesis. CpgD had optimum activity at pH 7.5-8 in the presence of magnesium. CpgD exhibited Michaelis-Menten kinetics with respect to 3-phosphoglycerate (Km,app = 10.9 ± 0.7 mM; Vmax,app = 0.72 ± 0.02 µmol/min/mg enzyme) and CTP (Km,app = 4.8 ± 0.9 mM; Vmax,app = 0.44 ± 0.03 µmol/min/mg enzyme). The characterization of this enzyme uncovers another piece of the pathway towards CPG2 synthesis. 
    more » « less