{"id":"7ffdd4a1-5059-4c51-b616-a181bba55b22","arxiv_id":"2507.21336","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Thiol modifications on OpcA, including a S. elongatus-specific glutathionylation at C398, are predicted to stabilize open gate conformations in three of four G6PDH subunits and enhance allosteric activation.","lead":"This study combines computer simulations and redox proteomics to show how chemical modifications on the protein OpcA change the shape of the G6PDH enzyme complex that produces a key energy-carrying molecule in cyanobacteria. The work identifies a specific gate mechanism that may explain how cells switch metabolism on and off within seconds.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"C398 glutathionylation is the linchpin of the paper's novel claim, but it rests on computed RBFE rather than direct MS detection; if C398 is not S-glutathionylated, the simulated PTMed complex may not represent physiology.","rationale":"The reader's weakest_assumption correctly identifies the unsupported C398 glutathionylation as the load-bearing fragility. My review of the full text confirms that the redox proteomics only reports cysteine oxidation without modification-type resolution, and the methods (Supporting Information) describe TI-based RBFE as the basis for the glutathionylation assignment. This is a critical gap because the PTMed complex simulation uses C398-glutathionylation as one of three fixed modifications; the gate distributions are conditional on that state. The paper gives no direct MS/MS evidence of the GSH adduct, no in vitro assay, and no control with C398 unmodified. I agree that this warrants a CONDITIONAL verdict: the model is plausible and internally consistent, but the physiological relevance of the most novel modification is unverified. The paper does have independent support: the AlphaFold3 complex is checked against cryo-EM structures (TM-scores in Fig. S2), the disulfide pairs match a prior cryo-EM study, and the redox proteomics is experimental. These reduce risk but do not remove the C398-specific gap. A targeted MS/MS re-analysis or a comparative simulation would settle the question. Therefore I recommend no change to the reader's verdict.","tokens_in":19769,"tokens_out":8156,"duration_ms":95996,"concrete_test":"Re-analyze the raw LC-MS/MS data from the dark-treated S. elongatus samples for the C398-containing peptide with a +305.07 Da (glutathionyl) mass shift, including diagnostic fragment ions (e.g., y or b ions with the added GSH moiety). If S-glutathionylation is not detectable at C398, then simulate the reduced-C398 control (or alternative oxidation states) and test whether the two disulfides alone still produce the open-gate distribution in three of four G6PDH subunits. If C398 glutathionylation is confirmed, the current conditional verdict can be upgraded; if not, the central claim must be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim—that thiol PTMs (C162-C174 and C380-C386 disulfides plus C398 glutathionylation) allosterically gate G6PDH—depends on the physiological reality of these modifications. The redox proteomics data (Results, Figure 3) demonstrate increased oxidation of C398 under dark conditions but do not identify the chemical nature of the modification; the RAC-based workflow enriches reversibly oxidized cysteines generically. The assignment of C398 as S-glutathionylated comes solely from computed relative binding free energies (RBFE) in the PTM-Psi toolkit, comparing glutathionylation against sulfenylation and nitrosylation. This thermodynamic preference in a monomeric model does not establish that glutathionylation occurs in vivo, where kinetics, glutathione redox state, and enzyme accessibility matter. If C398 is actually present as a sulfenic acid, S-nitrosothiol, or intramolecular disulfide, the simulated 'PTMed' complex contains a modification that does not exist under the measured conditions, and the gate-opening statistics (Figure 4B) may be an artifact of that hypothetical modification. The paper does not include a control simulation with C398 reduced or with an alternative modification, so the specific contribution of C398 glutathionylation to the three-of-four open gates is untested. The novelty of the paper hinges on this species-specific site, making the unverified modification type the single most load-bearing assumption.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":20059,"tokens_out":4302,"duration_ms":46762,"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":[{"comment":"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.","section":"Results (\"Redox proteomics data...\") and Methods (\"MD simulations of the Reduced and PTM-ed OpcA-G6PDH binary…"},{"comment":"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.","section":"Results (\"Thiol PTMs of OpcA in proximity...\"), Figure 4, and Methods (\"MD simulations of the Reduced and PTM-ed...\")"},{"comment":"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).","section":"Abstract and Conclusions"},{"comment":"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.","section":"Results and Figure S4"}],"minor_comments":[{"comment":"The color bar in Figure 1B is labeled \"RMSD (nm)\" while the analysis and caption describe RMSF; please relabel the color bar as RMSF.","section":"Figure 1B"},{"comment":"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.","section":"Methods (\"MD simulations of the Reduced and PTM-ed...\")"},{"comment":"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.","section":"Figure 3 and Methods (\"Computational Analyses\")"},{"comment":"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.\"","section":"Results (\"Redox proteomics data...\")"},{"comment":"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.","section":"Figure S3"}],"recommendation":"major_revision","confidential_remarks":"The stress-test concern about C398 glutathionylation is legitimate and is a primary basis for the major revision. The paper is otherwise a useful integration of redox proteomics with MD simulation, but the unverified modification identity and the indirect activity inference need to be addressed before publication. The authors should also consider whether the journal's readership will accept an AlphaFold3-derived complex as the sole structural basis for the S. elongatus system without experimental validation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nThis is worth a look if you work on redox regulation of cyanobacterial metabolism. The paper does something simple and useful: it takes the known OpcA-G6PDH disulfide switches and asks how they change the dynamics of the complex. The genuinely new bits are the ASP33-ARG243/GLU241-ARG243 gate hydrogen-bond mechanism and the S. elongatus-specific C398 glutathionylation. The gate mechanism is a reasonable hypothesis grounded in the existing cryo-EM structures, and the redox proteomics data showing dark-induced oxidation at C162, C174, C380, C386, and C398 is real, usable evidence. The PTM-Psi toolkit is their own, but it is QM-validated rather than fitted to this system, so the circularity worry is minor.\n\nThe soft spots are real but not fatal. First, the C398 modification type is not experimentally identified. The redox proteomics enriches reversibly oxidized cysteines generically; the assignment as S-glutathionylation comes solely from RBFE calculations. That is load-bearing for the most novel claim, and a control simulation with reduced C398 or an alternative modification would have made it much stronger. Second, the activity enhancement is inferred from gate geometry, not measured. The gate classification thresholds (0.4 and 0.7 nm) are hand-chosen, and the 100 ns simulations lack reported uncertainties. That weakens confidence in quantitative claims, e.g., the 10 kcal/mol binding difference. Third, no code or data are shipped, which makes it hard to check the gate statistics.\n\nNone of this kills the paper. The central idea—that oxidized OpcA stabilizes open active-site gates in three of four G6PDH subunits—is consistent with the cryo-EM evidence and with the proteomics, and the gate mechanism is specific enough to test experimentally (e.g., by mutating ASP33 or ARG243). I'd send it to a serious referee. The referee should push for MS-based confirmation of C398 modification, error bars or longer simulations, and ideally an activity assay. But the paper deserves engagement, not a desk reject.\n\nRecommendation: accept for peer review, likely major revisions.","headline":"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.","tokens_in":20609,"tokens_out":1953,"would_cite":true,"duration_ms":21439,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Thiol modifications on OpcA open G6PDH's active-site gate, boosting NADPH production when the cell is oxidized.","keywords":["thiol post-translational modifications","OpcA","G6PDH","allosteric regulation","redox proteomics","molecular dynamics","oxidative pentose phosphate pathway","conformational gate"],"falsifier":"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.","tokens_in":19582,"feed_emoji":"⚡","tokens_out":6913,"duration_ms":69529,"temperature":0.7,"pith_summary":"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.","feed_headline":"Oxidized OpcA opens G6PDH's active-site gate, boosting NADPH output","feed_subtitle":"Redox proteomics plus simulations show the modified OpcA holds three of four G6PDH gates open for NADPH production.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the cryo-EM structure of the cyanobacterial OpcA-G6PDH complex and identifies the two intramolecular disulfides that allosterically activate G6PDH.","marker":"(Doello et al. 2024)"},{"why":"Provides the PTM-Psi toolkit and QM-parameterized force fields used to model sulfenylated, nitrosylated, and glutathionylated cysteines and to compute relative binding free energies.","marker":"(Mejia-Rodriguez et al. 2023)"},{"why":"AlphaFold3 supplies the predicted OpcA-G6PDH complex structure on which the MD simulations are built.","marker":"(Abramson et al. 2024)"},{"why":"Establishes the mass-spectrometry approach for detecting multiple types of protein thiol modifications, underlying the interpretation of the redox proteomics data.","marker":"(Li et al. 2021)"},{"why":"Defines the S-glutathionylation proteome methods used to contextualize C398 as a potential glutathionylation site.","marker":"(Li et al. 2022)"},{"why":"Provides the resin-assisted capture enrichment protocol that produced the cysteine oxidation data for OpcA.","marker":"(Guo et al. 2014)"},{"why":"Shows thioredoxin regulates G6PDH activity through OpcA redox state in Anabaena and identifies the conserved cysteine pair critical for this regulation.","marker":"(Mihara et al. 2018)"},{"why":"Provides the diel transcriptome data showing opcA and zwf are strongly differentially expressed between dusk and noon, supporting the slow abundance-based regulatory layer.","marker":"(Gilliam et al. 2025)"}],"fun_headline_variants":["Thiol PTMs on OpcA pry open G6PDH's gate for NADPH","OpcA redox modifications swing G6PDH gate open","Redox-modified OpcA holds G6PDH gate open for NADPH output","OpcA thiol PTMs flip G6PDH gate to open state"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Thiol PTMs on OpcA pry open G6PDH's gate for NADPH","OpcA redox modifications swing G6PDH gate open","Redox-modified OpcA holds G6PDH gate open for NADPH output","OpcA thiol PTMs flip G6PDH gate to open state"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000707,"raw_usage":{"total_tokens":3274,"prompt_tokens":1123,"completion_tokens":2151,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":2064}},"tokens_in":739,"tokens_out":2151,"duration_ms":19531,"temperature":1.0,"reasoning_tokens":2064,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T12:52:26.808338+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[],"review_version":1}