Pith. sign in

REVIEW 4 major objections 5 minor 2 references

Thiol post-translational modifications modulate allosteric regulation of the OpcA-G6PDH complex through conformational gate control

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Thiol modifications on OpcA open G6PDH's active-site gate, boosting NADPH production when the cell is oxidized.

desk verdict A plausible computational extension of the OpcA-G6PDH cryo-EM story, with one genuinely new species-specific site (C398) whose modification type is inferred, not proven. read the letter →

arxiv 2507.21336 v1 pith:QFAZONV3 submitted 2025-07-28 physics.bio-ph

Hoshin Kim , Song Feng , Pavlo Bohutskyi , Xiaolu Li , Daniel Mejia-Rodriguez , Tong Zhang , Wei-Jun Qian , Margaret S. Cheung This is my paper · ORCID
classification physics.bio-ph
keywords thiolpost-translationalmodificationsOpcAG6PDHallostericregulationredoxproteomicsmoleculardynamicsoxidativepentosephosphatepathwayconformationalgate
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper seeks to explain how redox post-translational modifications on the cyanobacterial protein OpcA switch on the enzyme G6PDH, the rate-limiting step of the oxidative pentose phosphate pathway. Combining redox proteomics of Synechococcus elongatus with molecular dynamics simulations, it argues that two disulfide bonds (C162-C174 and C380-C386) plus glutathionylation of C398 near the OpcA-G6PDH interface stabilize an 'open' conformation of a molecular gate at the G6PDH active site. In the reduced state the gate is mostly closed or wide-open, which restricts substrate access or disrupts active-site hydrogen bonds; in the oxidized state three of four G6PDH subunits hold the gate open. The paper concludes that these thiol modifications act as fast allosteric switches that boost NADPH production when light is unavailable.

What carries the argument

The load-bearing object is the OpcA-G6PDH complex as modeled by AlphaFold3 and simulated with PTM-Psi parameterized MD. Within G6PDH, the central mechanism is a molecular gate: paired charged residues forming the upper (GLU241-ARG243) and lower (ASP33-ARG37) rims of the substrate channel, with ASP33-ARG243 hydrogen bonds distinguishing closed from open configurations. PTMs on OpcA act through this gate: the two disulfides and C398 glutathionylation, located at the interface with G6PDH subunits A-C, propagate long-range contacts to ARG243 and stabilize open conformations in three of four subunits. The same machinery explains why subunit D, which lacks OpcA contact, fails to adopt the open state.

What would settle it

Targeted mass spectrometry that distinguishes S-glutathionylation from sulfenylation, nitrosylation, or disulfide formation at C398 in dark-shifted S. elongatus would settle the claim; a finding that C398 is not glutathionylated, or a C398A mutant that retains the same gate-opening and activity enhancement, would falsify the paper's most specific conclusion.

Watch

Extended reading notes

Core claim

The central claim is that thiol PTMs on OpcA do not merely strengthen binding to G6PDH; they allosterically control G6PDH activity by setting the conformation of a gate near the active site. In the PTMed complex, formed by the C162-C174 and C380-C386 disulfides and glutathionylated C398, OpcA binds G6PDH about 10 kcal/mol more strongly, reduces overall flexibility, and shifts three of the four G6PDH subunits into the 'open' gate conformation that is optimal for substrate ingress and product egress. The closed and wide-open states seen in the reduced complex are less functional: closed restricts access, wide-open disrupts hydrogen bonds such as ASP197-HIS260 in the active site. Gate state is governed by hydrogen bonds between ASP33-ARG243 and GLU241-ARG243, which are present in the closed state and absent in the open state. The paper therefore identifies PTM-driven conformational gate control, not direct redox sensing by G6PDH, as the molecular mechanism of rapid metabolic switching in cyanobacteria.

Load-bearing premise

The load-bearing premise is that C398 is actually glutathionylated in S. elongatus under dark conditions; the redox proteomics data show C398 becomes oxidized but do not reveal which oxidation product it is, and the glutathionylation assignment comes from computed free energies rather than direct detection.

Editorial extensions

If this is right

  • Oxidized OpcA enhances G6PDH activity by shifting three of four subunits into the open-gate conformation, so NADPH production can be ramped up on a sub-second-to-minute timescale without new protein synthesis.
  • The conserved gate residues imply the same open/closed/wide-open logic may govern G6PDH activity in other organisms, including humans, even where OpcA is absent.
  • C398 glutathionylation adds a third, species-specific redox switch beyond the two conserved disulfides, allowing S. elongatus to fine-tune its dark response.
  • Because gate conformation rather than active-site chemistry is the controlled variable, the study predicts that mutations at ASP33, ARG243, GLU241, or ARG243 will alter G6PDH activity by changing gate stability.
  • The roughly 10 kcal/mol stronger OpcA-G6PDH interaction upon PTM formation provides a quantitative anchor for why oxidized OpcA remains bound and activates the complex.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If targeted mass spectrometry confirms S-glutathionylation at C398, the C398 switch could be engineered independently of the disulfides, giving synthetic biologists a single-residue dial for NADPH output in cyanobacterial strains.
  • The gate-distance classification (closed below 0.4 nm, open 0.4-0.7 nm, wide-open above 0.7 nm) could be reused as a screening metric to predict whether other PTMs, drugs, or mutations activate or suppress G6PDH.
  • The paper's reliance on AlphaFold3 for the complex geometry leaves open the possibility that alternative OpcA-G6PDH stoichiometries or docking modes produce different gate outcomes; testing this would require experimental structures of S. elongatus complexes.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The manuscript combines redox proteomics of Synechococcus elongatus PCC 7942 with PTM-Psi-based molecular dynamics simulations to propose that thiol post-translational modifications of OpcA—specifically the intramolecular disulfides C162-C174 and C380-C386 plus S-glutathionylation of C398—allosterically regulate the G6PDH tetramer by stabilizing an "open" active-site gate conformation in three of four subunits. The authors classify gate distances into closed/open/wide-open categories, report that PTMs increase the OpcA-G6PDH interaction energy by roughly 10 kcal/mol, and argue that open gates facilitate substrate ingress and product egress, thereby enhancing G6PDH catalytic activity under oxidizing conditions. The central mechanistic claim is that PTM-induced conformational changes at the protein-protein interface propagate to the enzyme active site through long-range contacts to specific gate residues.

Significance. If correct, the paper would provide a molecular-level mechanism for redox regulation of the oxidative pentose phosphate pathway in cyanobacteria, with a species-specific C398 glutathionylation site as a novel regulatory element beyond the conserved disulfides. The work integrates experimental redox proteomics with atomistic simulation and uses QM-validated force field parameters for non-standard cysteine modifications, which is a clear strength. The gate hypothesis is falsifiable and offers concrete targets for mutagenesis. However, the central claim that PTMs enhance catalytic activity is not directly measured; it is inferred from gate geometry and hydrogen-bonding patterns, and the physiological identity of the C398 modification is not established by the proteomics data. These issues limit the certainty with which the conclusions can be accepted.

major comments (4)
  1. [Results ("Redox proteomics data...") and Methods ("MD simulations of the Reduced and PTM-ed OpcA-G6PDH binary…] C398 glutathionylation is load-bearing but not experimentally established. The redox proteomics data in Figure 3 show increased oxidation of C398 under dark conditions, but the RAC-based workflow enriches reversibly oxidized cysteines generically and does not identify the chemical modification. The assignment of S-glutathionylation rests solely on the PTM-Psi RBFE calculation (Results: "glutathionylation at C398 is the most energetically favorable PTM"). This thermodynamic preference in a monomer does not prove that glutathionylation occurs in vivo, where kinetics, glutathione redox state, and enzyme accessibility matter. Because the PTMed complex simulated in Figure 4 and all subsequent conclusions include C398 glutathionylation, the central claim is conditional on this unverified modification. Please provide direct mass-spectrometric evidence for S-glutathionylation, or run control simulations with C398 reduced or alternatively modified, and soften the abstract and conclusions accordingly.
  2. [Results ("Thiol PTMs of OpcA in proximity..."), Figure 4, and Methods ("MD simulations of the Reduced and PTM-ed...")] The gate population statistics in Figure 4 are not adequately supported. The thresholds 0.4 nm and 0.7 nm are taken from cryo-EM structures of different cyanobacteria and applied as hard cutoffs, but no uncertainties are reported. The simulations are 100 ns and the number of replicas for the two complexes is not stated in the Methods; the text says "The aforementioned procedures were applied" without specifying length or replicas. The statement that "three out of four G6PDH subunits showed open conformations" depends on these thresholds and on limited sampling. Please report the full distance distributions with block-error estimates, state the simulation length and replica count for the complexes, and show that the conclusions are robust to reasonable variations of the open/closed thresholds.
  3. [Abstract and Conclusions] The abstract's claim that "these PTMs promote rapid metabolic switching by enhancing G6PDH catalytic activity when OpcA is oxidized" is not directly supported by the data. No enzymatic activity, kcat, or Km is measured; activity is inferred from gate geometry and hydrogen-bond occupancy (Figures 4 and 5E). While the gate mechanism is plausible and cited literature supports gate-activity relationships, the inference is indirect. Please either temper the claim to "likely enhance" or "are consistent with enhanced activity," or provide a quantitative computational estimate (e.g., substrate binding free energy or barrier height in open versus closed states).
  4. [Results and Figure S4] The ~10 kcal/mol increase in OpcA-G6PDH interaction energy (quoted in Results and Conclusions) is reported without uncertainty and is obtained from gmx_MMPBSA, an end-state method known to be sensitive to system setup. Please provide error bars (e.g., over frames or replicas) and, if possible, corroborate with an alchemical or PMF-based binding free energy difference.
minor comments (5)
  1. [Figure 1B] The color bar in Figure 1B is labeled "RMSD (nm)" while the analysis and caption describe RMSF; please relabel the color bar as RMSF.
  2. [Methods ("MD simulations of the Reduced and PTM-ed...")] The simulation length and number of replicas are specified for the OpcA monomer MD (100 ns, three replicas) but are not stated for the OpcA-G6PDH complex simulations; please add these details.
  3. [Figure 3 and Methods ("Computational Analyses")] The right panel of Figure 3 reports RBFE values without error bars, and the definition of the L1 norm is not described in the text; please define the L1 norm in the caption or Methods and include uncertainties for the free energy estimates.
  4. [Results ("Redox proteomics data...")] The statement that "glutathionylation at C398 is the most energetically favorable PTM compared to other potential thiol PTM sites during light/dark perturbations" overstates the RBFE, which is a relative free energy in the monomer, not a rate or an experimentally measured occupancy under light/dark conditions; please revise to "computed to be the most energetically favorable."
  5. [Figure S3] The multiple sequence alignment in Figure S3 compares only two OpcA sequences; the claim that C398 is "not conserved in OpcA proteins from other species" would be better supported by a broader alignment across more cyanobacterial species.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: gate populations are emergent MD observables and the only self-citations are to a QM-validated toolkit.

full rationale

The derivation chain is self-contained. The central claim—that C162–C174 and C380–C386 disulfides plus C398 glutathionylation stabilize open G6PDH gates—rests on MD simulations of two constructed complexes, and the gate populations are emergent observables, not imposed by the simulation setup. The two disulfide sites are taken from an external cryo-EM study (Doello et al. 2024) and from the authors' redox proteomics, which independently show increased oxidation of C162, C174, C380, C386, and C398 under dark conditions; the assignment of C398 as glutathionylated is a thermodynamic RBFE prediction from PTM-Psi, not a fit to the gate data. The PTM-Psi toolkit is cited from the authors' prior work, but the force-field parameters are QM-validated and the RBFE protocol is described in the Methods, so the self-citation does not smuggle in the target result. No equation equates a predicted quantity with an input by construction, no parameter is fitted to the quantity later called a prediction, and no uniqueness theorem from the authors is invoked to force the choice of PTM state. Areas of concern, such as the absence of direct mass-spectrometric detection of S-glutathionylation at C398 and the lack of a control simulation with an alternative C398 modification, are physiological-assumption and completeness issues, not circularity.

Assumptions & free parameters 1 free parameters · 4 assumptions · 0 invented entities

The central claims rest on assumptions about model accuracy, sampling adequacy, the gate-activity link, and the specific identity of the C398 modification. The gate itself is a structural feature inferred from MD, not a newly invented physical entity.

free parameters (1)
  • Gate classification distance thresholds = 0.4 nm and 0.7 nm
    Used to bin gate conformations into closed, open, and wide-open states; chosen from cryo-EM distances of G6PDH complexes, not from a first-principles rule. The population comparison in Figure 4 depends on these thresholds.
assumptions (4)
  • domain assumption AlphaFold3-predicted OpcA-G6PDH complex is sufficiently accurate for allosteric analysis
    The complex was not experimentally determined for S. elongatus; the model is validated by TM-score against a Synechocystis cryo-EM structure, but subtle interface geometry relevant to allostery may differ.
  • domain assumption 100 ns MD simulations sample the relevant gate conformational states
    Gate opening and closing may occur on longer timescales; no convergence analysis or replica exchange is reported. Gate population distributions from the trajectories underpin the main claim.
  • domain assumption Open gate conformation correlates with higher G6PDH catalytic activity
    The paper infers enhanced activity from open gates based on prior enzyme gate literature, but no G6PDH activity assay is performed in this study. Wide-open states are treated as detrimental based on loss of ASP197-HIS260 hydrogen bonds.
  • domain assumption C398 is glutathionylated under dark conditions
    Redox proteomics detects oxidation at C398 but does not identify the modification type. The selection of glutathionylation is based on computed RBFE values, not direct mass spectrometry evidence.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Thiol post-translational modifications modulate allosteric regulation of the OpcA-G6PDH complex through conformational gate control." pith.science (2026). https://pith.science/paper/QFAZONV3

@misc{pith2026250721336,
  author       = {Pith},
  title        = {Pith review of: Thiol post-translational modifications modulate allosteric regulation of the OpcA-G6PDH complex through conformational gate control},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/QFAZONV3}},
  note         = {Machine review of arXiv:2507.21336}
}
read the original abstract

Cyanobacteria require ultra-fast metabolic switching to maintain reducing power balance during environmental fluctuations. Glucose-6-phosphate dehydrogenase (G6PDH), catalyzing the rate-limiting step of the oxidative pentose phosphate pathway (OPPP), provides essential NADPH and metabolic intermediates for biosynthetic processes and redox homeostasis. In cyanobacteria, the unique redox-sensitive protein OpcA acts as a metabolic switch for G6PDH, enabling rapid adjustment of reducing power generation from glycogen catabolism and resulting in precise regulation of carbon flux between anabolic and catabolic pathways. While the redox-sensitive cysteine structures of OpcA are known to regulate G6PDH, the detailed mechanisms of how redox post-translational modifications (PTMs) influence OpcA's allosteric effects on G6PDH structures and function remain elusive. To investigate this mechanism, we utilized computational modeling combined with experimental redox proteomics using Synechococcus elongatus PCC 7942 as a model system. Redox proteomics captured modified cysteine residues under light/dark or circadian shifts. Computational simulation revealed that thiol PTMs near the OpcA-G6PDH interface are crucial to allosteric regulation of regions affecting the G6PDH activity, including a potential gate region for substrate ingress and product egress, as well as critical hydrogen bond networks within the active site. These PTMs promote rapid metabolic switching by enhancing G6PDH catalytic activity when OpcA is oxidized. This study provides evidence for novel molecular mechanisms that elucidate the importance of thiol PTMs of OpcA in modulating G6PDH structure and function in an allosteric manner, demonstrating how PTM-level regulation provides a critical control mechanism that enables cyanobacteria to rapidly adapt to environmental fluctuations through precise metabolic fine-tuning.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    E., Jalalypour, F ., Jordan, J., Kutzner, C., Lemkul, J

    Abraham M, Alekseenko, A., Basov, V ., Bergh, C., Briand, E., Brown, A., Doijade, M., Fiorin, G., Fleischmann, S., Gorelov, S., Gouaillardet, G., Grey, A., Irrgang, M. E., Jalalypour, F ., Jordan, J., Kutzner, C., Lemkul, J. A., Lundborg, M., Merz, P ., Lindahl, E. 2024. Gromacs 2024.3 manual (2024.3). Abramson J, Adler J, Dunger J, Evans R, Green T, Prit...

    work page 2024

  2. [2]

    AJ, Bambrick J et al. 2024. Accurate structure prediction of biomolecular interactions with alphafold 3. Nature. 630(8016):493-500. Bennett CH. 1976. Efficient estimation of free -energy differences from monte -carlo data. J Comput Phys. 22(2):245-268. Bernetti M, Bussi G. 2020. Pressure control using stochastic cell rescaling. The Journal of Chemical Phy...

    work page 2024

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.