{"id":"27ceab40-ce39-43d3-a53d-a4a3fdb5360c","arxiv_id":"2505.01327","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A1493 stabilizes the first codon-anticodon pair through entropy restriction and rigid wedge support, reconciling dynamic and rigid models of ribosomal decoding.","lead":"Simulations show that ribosomal nucleotide A1493 stabilizes the codon-anticodon interaction by limiting the opening motions of the first codon base and by acting as a wedge between the mRNA and tRNA backbones. This computational synthesis proposes how the ribosome tells correct tRNAs from near-correct ones and refines two competing decoding models.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 20 Å ribosome fragment with outer-atom restraints may bias the A1493 off-state PMF that underlies both the entropy-stabilization and wedge claims; a larger-shell rerun is the decisive test.","rationale":"The reader's verdict identified the 20 Å truncation as the weakest assumption, and I agree that it is the single most load-bearing issue. The reason it is load-bearing rather than merely a generic long-range-electrostatics worry is locatable: the off-state A1493 branch of the 2D PMF is not a minor correction but the exact branch used to construct Eqs. 1–2, to infer entropy restriction, and to read off the U-U distance shift that defines the wedge. If that branch is distorted by restraints, the paper's three headline mechanisms (entropy stabilization, amplification, wedge) would not necessarily survive, although the qualitative observation that the ribosome stabilizes the codon-anticodon helix might. The paper has independent value: it uses a recognized recognition-intermediate structure (5UYL), a reasonable force-field setup, and a clear 2D US design, and the authors are transparent about the small-region limitation. The proposed check is a direct computational replication at larger cutoff, which is feasible at the same class of computational cost. I would therefore keep the CONDITIONAL verdict: the concern is serious enough to require the larger-shell test or an equivalent validation, but not so severe as to reject the local structural observations without testing.","tokens_in":20138,"tokens_out":9178,"duration_ms":101176,"concrete_test":"Recompute the same 2D US (CV1 = U4-U36 opening distance, CV2 = A1493 flip angle) for the U-A and U-U systems in a 40 Å-radius fragment with positional restraints applied only beyond 30 Å, or with a stochastic boundary/elastic network on the shell, and recalculate (i) the integrated on/off U4-opening PMFs and (ii) the off-state C1' backbone distance of the U-U pair. If the on/off barrier difference changes by more than ~1 kcal/mol or the wedge distance shift (6.1→5.4 Å) moves by more than ~0.5 Å, the 20 Å/17 Å truncation is load-bearing and the entropy and wedge conclusions must be re-evaluated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claims of the paper depend on the conditional free-energy surfaces obtained by integrating the 2D umbrella-sampling PMF over A1493 \"off\" states (flip angle <100°, Eqs. 1–2; Figs. 4A and 5A/C). In the off state A1493 re-enters helix 44, and the reaction-coordinate definition includes G1494/G1491 atoms; the cost of this re-entry is therefore coupled to helix-44 bending and relaxation. But the simulation model is a 20 Å shell in which \"outer heavy atoms ... located more than 17 Å from A1493 and the first base pair were restrained to their initial positions\" (Methods, System setup). This positional restraint freezes the long-wavelength helix-44 response that would accompany A1493 flipping, so the on/off PMF difference and the U-U backbone contraction (6.1→5.4 Å, Fig. 5C) used to define the \"wedge\" may be boundary artifacts rather than intrinsic ribosome mechanics. The authors themselves acknowledge in the Discussion that restricting the simulation to a small region \"prevent[s] analysis of the competition between the overall elastic potential energy of the bent tRNA and its complete interaction with mRNA,\" and no comparison against a larger or full-ribosome system is provided. This is the most load-bearing assumption because both the entropy-stabilization argument (on vs. off A1493) and the rigid-support argument (U-U shift in the off state) pass through this truncated, restrained off-state PMF.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports molecular dynamics simulations of the ribosomal decoding center, built from a 20 Å truncated fragment of the 70S ribosome (PDB 5UYL) containing the A-site codon-anticodon helix and A1493. Using 1D and 2D umbrella sampling, the authors compute free-energy surfaces for opening of the first codon-anticodon base pair and for flipping of A1493, in cognate (U-A) and near-cognate (U-U, U-C, U-G) systems, in ribosome-free controls, with pseudouridine-modified codons, and in A-site/P-site/translocation constructs. The central claims are that A1493 stabilizes the first base pair through an entropy-restriction mechanism, that A1493 acts as a rigid wedge keeping the mRNA and tRNA backbones apart, that the ribosome amplifies the cognate versus near-cognate stability difference to about 6 kcal/mol, that A1493 flipping is non-specific, that steric complementarity with tRNA position 37 is essential, that pseudouridine stabilizes the P/A kink via a water bridge, and that the codon-anticodon helix is transiently destabilized during translocation.","tokens_in":20422,"tokens_out":3654,"duration_ms":40721,"significance":"If the claims hold, the work would provide a quantitative, structure-based reconciliation of the dynamic and rigid models of ribosomal decoding, and it would offer a molecular rationale for pseudouridine-induced mistranslation. The study has notable strengths: the use of a recognition-intermediate cryo-EM structure, explicit 2D umbrella sampling with a defined reaction-coordinate pair, multiple near-cognate systems, ribosome-free controls, comparisons against experimental structures such as 5EL4 and 4JV5, and an unusually candid Discussion section that states the model's limitations. The simulations are not fitted to experimental discrimination free energies, so the conclusions are not circular in that sense. However, several load-bearing interpretations—especially the entropy label and the magnitude of the discrimination gap—go beyond what the presented free-energy surfaces and error analysis support, and the central mechanism is computed within a heavily restrained truncated system for which no convergence or shell-size validation is provided.","major_comments":[{"comment":"The central on/off A1493 comparison is computed in a 20 Å truncated ribosome fragment in which outer heavy atoms beyond 17 Å from A1493 and the first base pair are restrained to their initial positions. When A1493 enters the 'off' state it re-enters helix 44, and the reaction coordinate explicitly includes G1494/G1491 atoms; the cost of this re-entry is therefore coupled to helix-44 bending and long-wavelength relaxation that the positional restraints freeze. The on/off free-energy difference underlying the entropy-stabilization claim (Eqs. 1-2, Fig. 4A) and the U-U backbone contraction from 11.1 to 9.1 Å with base distance 6.1 to 5.4 Å used for the wedge claim (Fig. 5C) may thus be boundary artifacts. The authors acknowledge in the Discussion that the small-region restriction prevents analysis of the overall elastic tRNA response, but no larger-shell or unrestrained comparison is provided. A decisive test would be to recompute the 2D PMF with a substantially larger shell (e.g., 30 Å or more) or with the helix-44 region left unrestrained, and to report the sensitivity of the on/off barrier and the U-U minimum shift to the restraint cutoff.","section":"Methods, System setup; Results, 'The stability contributed by A1493 depends on its flip angle'"},{"comment":"The claim that A1493 stabilizes the first base pair through an entropy mechanism is not supported by the data presented. The text states that 'A1493 reduces the entropy increase of the opening process by limiting the conformations numbers of open state,' but no enthalpy/entropy decomposition, temperature-dependent free energies, or configurational-entropy estimate is provided. The inference from reduced barrier height and constrained scatter in Fig. 3C/D is indirect: a lower barrier and a narrower open-state distribution are consistent with entropy effects, but they do not establish that the stabilization is entropic rather than enthalpic. Furthermore, Eqs. 1-2 merely re-express the same 2D free energy as conditional projections; the difference between F_on and F_off is a free-energy difference, not a decomposition into entropy and enthalpy. The 'entropy-driven stabilization' wording should be replaced by a statement about the free-energy profile unless a direct entropy calculation is added.","section":"Results, 'Mechanism of A1493 in stabilizing the codon-anticodon helix'; Summary and Discussion"},{"comment":"The comparison of Ψ-C with the unmodified U-C system is not made on equal footing. Because C36 moved instead of Ψ4, the authors write that 'the free energy of the Ψ-C pair was calculated with C36 flipping constrained.' This post-hoc constraint changes the sampled reaction coordinate for one system only, so the conclusion that Ψ-C stability reaches the level of U-A is not established without justification. Please report the unconstrained Ψ-C PMF as well, show that the constrained calculation is converged, and demonstrate that the constraint does not artificially stabilize the paired state; otherwise the claim that Ψ converts a mismatch into a cognate-like pair is not testable from the present data.","section":"Results, 'ψ enhances mRNA-tRNA pairing by pre-organizes A-form of mRNA'"},{"comment":"The ~6 kcal/mol discrimination gap is presented as 'sufficient for the ribosome to distinguish' cognate from near-cognate tRNAs, but the manuscript does not specify whether this is a barrier-height difference or a well-depth difference, does not propagate the block-analysis errors into the reported ΔΔG, and does not compare the computed difference with experimentally measured discrimination free energies. The ribosome-free U-A versus U-U difference of about 1 kcal/mol is computed in a different model (no ribosome, terminal restraints), so the 'amplification' factor is a comparison between two different free-energy surfaces rather than a single controlled variation. Please state precisely which features of the PMFs are being compared, provide error bars on the 6 kcal/mol value, and relate it quantitatively to known initial-selection and proofreading discrimination energies.","section":"Results, 'Ribosomal discrimination between cognate and near-cognate tRNA'"}],"minor_comments":[{"comment":"The figure has two panels labeled 'B': one for the aIGM/molecular-surface analysis and one for the free-energy curves of U4 opening. Please renumber the panels consistently and update the callouts in the text.","section":"Figure 5"},{"comment":"The typesetting of Eqs. 1-2 contains spacing artifacts and the subscripts 'on' and 'off' are hard to read. Please reformat these equations so that the integration domains and Boltzmann factors are unambiguous.","section":"Equations (1)-(2)"},{"comment":"Reference 43 is cited for the statement that the 5' end of RNA opens more easily, but reference 43 is the WHAM paper by Kumar et al. Please replace it with the appropriate RNA-fraying or base-pair-opening reference.","section":"Reference [43]"},{"comment":"The phrase 'pre-organizes A-form of mRNA' is grammatically awkward; consider 'pre-organizing the mRNA into the A-form' or 'pre-organization of the A-form mRNA structure.'","section":"Results, 'ψ enhances mRNA-tRNA pairing'"},{"comment":"Statistical errors are estimated by dividing each window into four blocks for WHAM analysis; this is likely to underestimate sampling error for a 2D PMF. Please state how many independent windows contribute to each region of the PMF and consider larger block counts or bootstrapping.","section":"Methods, MD simulations"}],"recommendation":"major_revision","confidential_remarks":"The paper makes strong mechanistic claims ('entropy stabilization,' 'wedge,' 'amplification') from a single truncated, restrained simulation model. The most important missing element is a shell-size/restraint-sensitivity check for the 2D PMF, since both headline mechanisms pass through the on/off A1493 states. A second missing element is any direct entropy calculation. I would not reject, because the simulations are internally consistent and well controlled in several respects, but the current wording overstates the evidence. If the authors cannot run larger systems, they should soften the mechanistic claims and present the results as conditional on the 20 Å model. I also encourage depositing the umbrella-sampling input files and topologies, since 'available from the corresponding author upon reasonable request' is weak for a computation-heavy study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things about this one: it is a genuinely new simulation result about A1493, and the main free-energy claim rests on a truncated model that has not been tested against a larger shell. Worth reading, worth refereeing, but the entropy label needs more evidence.\n\nWhat is new: the 2D umbrella sampling coupling A1493 flip angle to first base-pair opening, the wedge mechanism (A1493 prevents U-U from contracting to a two-hydrogen-bond state), the steric role of tRNA position 37 in maintaining the correct A1493 flip, and a water-bridge explanation for Ψ's stabilizing effect. The paper also offers a reasonable reconciliation of the dynamic and rigid decoding models: A1493 flipping is not tRNA-specific, but the helix is stabilized and the stability difference between cognate and near-cognate is amplified. The simulations are independent of the conclusions; no fitting to experimental discrimination free energies was done. The ribosome-free controls give useful external benchmarks.\n\nThe soft spots are real but uneven. The biggest is the 20 Å fragment with outer heavy atoms restrained. The off-state PMF, from which both the entropy-stabilization and wedge arguments are integrated, includes A1493 re-entering helix 44 and uses G1494/G1491 atoms in the reaction coordinate. Positional restraints freeze the long-wavelength helix-44 response that would accompany that re-entry. That could change the on/off difference and the U-U backbone contraction. The authors acknowledge in the Discussion that the small region prevents analysis of bent-tRNA elasticity, but they do not test the shell size. A larger-shell rerun is the decisive control.\n\nSecond, \"entropy-driven stabilization\" is asserted from the shape of the PMF, not from an enthalpy/entropy decomposition. The entropy interpretation may be right, but this paper does not demonstrate it; an enthalpic restriction of the open state would produce the same curve.\n\nThird, the Ψ-C comparison uses a post-hoc restraint on C36 flipping. That is a modeling choice, not a fatal flaw, but it needs justification. I also noticed missing references (citations 65–67 appear in text but not in the list), no numerical error bars (color strips only), and no deposited data or code. Those are fixable.\n\nBottom line: this is a solid, internally consistent simulation paper with a plausible new mechanism. It deserves a serious referee, not a desk reject. I would ask the authors for the larger-shell control and the entropy decomposition before endorsing the mechanism. If you work on translation fidelity, bring it to the group; it will generate discussion.","headline":"A1493 paper has a genuinely new 2D free-energy result and a plausible wedge/entropy mechanism, but the off-state PMF comes from a restrained 20 Å fragment, so the decisive control is a larger-shell rerun.","tokens_in":20949,"tokens_out":2656,"would_cite":true,"duration_ms":27579,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A1493 stabilizes the codon-anticodon helix by shrinking the open state's entropy and wedging mRNA and tRNA apart.","keywords":["ribosomal decoding","A1493","codon-anticodon helix","entropy stabilization","umbrella sampling","molecular dynamics","pseudouridine","translation fidelity"],"falsifier":"If a simulation of a larger ribosome model, or an experiment using atomic mutagenesis to delete or lock A1493, found that the opening free-energy barrier of the first codon-anticodon base pair does not drop when A1493 is removed and does not rise when it is locked in the fully flipped state, the claimed entropy-stabilization and wedge mechanism would be refuted.","tokens_in":19938,"feed_emoji":"🧬","tokens_out":5270,"duration_ms":48589,"temperature":0.7,"pith_summary":"This paper uses molecular dynamics simulations of the ribosomal decoding center to establish what the conserved 16S rRNA nucleotide A1493 actually does during tRNA selection. It argues that A1493 stabilizes the first codon-anticodon base pair entropically: by limiting the number of conformations the first codon base can adopt when its pair opens, it makes the closed state more probable and raises the cost of dissociation. It also claims A1493 acts as a rigid wedge between the mRNA and tRNA backbones, and that this support, together with steric complementarity involving the tRNA's 37th nucleotide, is what lets the ribosome amplify the stability difference between cognate and near-cognate tRNAs into a roughly 6 kcal/mol discrimination gap. The paper further links pseudouridine modification to water-bridged pre-organization of mRNA and describes a transient destabilization of codon-anticodon pairing during translocation.","feed_headline":"A1493 steadies codon-tRNA pairing by shrinking its open states","feed_subtitle":"MD simulations show it restricts the first codon base's open conformations and wedges the mRNA-tRNA backbones apart.","key_machinery":"The load-bearing machinery is the two-dimensional free-energy surface computed by umbrella sampling, with CV1 being the opening distance of the first codon-anticodon base pair and CV2 being the flipping angle of A1493, a dihedral defined through centers of mass of the G1494 and G1491 base pairs, the G1494 phosphate, the A1493 phosphate, and the A1493 base. This map lets the authors integrate out the A1493 angle to obtain opening profiles for the \"on\" and \"off\" states, and it reveals the wedge coupling, since in the near-cognate U-U system decreasing the A1493 angle shifts the base-pair minimum from 6.1 Å to 5.4 Å.","core_discovery":"On the paper's own terms, the central discovery is a dual mechanism for A1493. Through two-dimensional umbrella sampling that couples the opening distance of the first base pair with the flipping angle of A1493, the paper shows that a fully extrahelical A1493 restricts the accessible conformations of the open state of the first codon nucleotide, reducing the entropy gain of opening and thereby stabilizing the codon-anticodon helix. In near-cognate complexes the same nucleotide acts as a wedge: when A1493 moves toward helix 44, the mRNA and tRNA backbones approach each other and the U-U pair collapses from one hydrogen bond to two, showing that A1493 holds the backbones apart and imposes a fixed geometry. The paper also reports that A1493 flipping itself is not tRNA-specific and that the energy gap between cognate and near-cognate systems, around 6 kcal/mol, is amplified by the ribosome, with van der Waals steric complementarity, especially contact with the tRNA 37th nucleotide, more important than A-minor hydrogen bonding.","pith_inferences":["If correct, this \"conformational caging\" mechanism suggests a general fidelity strategy in RNA machines: reduce the entropy of rejected states rather than strongly stabilizing the accepted state, achieving discrimination with modest enthalpic differences.","The paper's claim that A1493 flipping is tRNA-nonspecific at the recognition intermediate implies that the specificity bottleneck lies in the coupled behavior of A1913 and helix 44; a simulation that explicitly lets A1913 insert and detach would test whether the amplification is gated by the large subunit.","The pseudouridine water-bridge mechanism predicts a measurable solvent effect: altering water activity or using heavy water should weaken pseudouridine-induced misreading, since the bridging water would be disrupted.","The translocation-destabilization result suggests that engineered ribosomes with altered U1498 stacking or modified P-site codon channels should show changed frameshift frequencies, a testable extension beyond the paper's current systems."],"forward_implications":["Ribosomal discrimination does not require A1493 to adopt different conformations for cognate versus near-cognate tRNA; its stabilizing effect is constitutive, and the distinction emerges from amplified free-energy differences.","Chemical modifications or mutations that prevent A1493 from fully flipping outward should lower the opening barrier of the first codon-anticodon base pair and increase miscoding, a prediction that atomic mutagenesis could test.","Pseudouridine at the first codon position can raise the stability of near-cognate mismatches to the level of a cognate U-A pair, providing a structural explanation for pseudouridine-induced amino acid substitutions.","Translocation transiently destabilizes the codon-anticodon helix, with stability restored only at the P site through distinct constraints such as U1498, giving a mechanistic route to frameshifting if that stabilization is lost."],"supporting_citations":[{"why":"Supplies the pre-accommodation ribosome-tRNA recognition intermediate structure used to build every A-site simulation system.","marker":"[18]"},{"why":"Defines the open-to-closed domain-closure model of decoding that the paper's discrimination argument builds on and partly supports.","marker":"[2]"},{"why":"Provides the rigid-decoding-center near-cognate structures and the claim that the first codon nucleotide is constrained against wobble, which the paper contests with its wobble U-G result.","marker":"[25]"},{"why":"Gives the prior molecular-dynamics result of cognate-specific A1492/A1493 flipping whose 3 kcal/mol difference the present study attributes to modeling the closed state rather than the recognition intermediate.","marker":"[27]"},{"why":"Documents pseudouridine-induced amino acid substitutions that the paper's pseudouridine water-bridge mechanism is built to explain.","marker":"[33]"},{"why":"Provides the post-translocation and time-resolved structures used for the P-site stability and translocation-destabilization analysis.","marker":"[15]"},{"why":"Atomic mutagenesis showing that minor-groove hydrogen-bond disruption barely affects discrimination, cited as evidence for steric complementarity over A-minor hydrogen bonding.","marker":"[28]"},{"why":"Shows that steric complementarity in the decoding center is important for tRNA selection, supporting the wedge interpretation.","marker":"[29]"}],"fun_headline_variants":["A1493 wedges mRNA-tRNA and shrinks codon open states","Entropy plus wedge: how A1493 steadies the decoding helix","A1493's steric fit and entropy effects govern codon decoding","MD shows A1493 provides rigid support while stabilizing pairing"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations assume that a 20 Å fragment of the ribosome, with its outer atoms held fixed, reproduces the thermodynamics of the intact decoding center; if long-range electrostatic, allosteric, or mechanical coupling from the rest of the ribosome changes these free-energy surfaces, the entropy and wedge mechanisms could be artifacts of the truncation.","fun_headline_variants_meta":{"raw":{"variants":["A1493 wedges mRNA-tRNA and shrinks codon open states","Entropy plus wedge: how A1493 steadies the decoding helix","A1493's steric fit and entropy effects govern codon decoding","MD shows A1493 provides rigid support while stabilizing pairing"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001074,"raw_usage":{"total_tokens":4500,"prompt_tokens":953,"completion_tokens":3547,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":3472}},"tokens_in":569,"tokens_out":3547,"duration_ms":23457,"temperature":1.0,"reasoning_tokens":3472,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T04:20:25.551732+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"If a simulation of a larger ribosome model, or an experiment using atomic mutagenesis to delete or lock A1493, found that the opening free-energy barrier of the first codon-anticodon base pair does not drop when A1493 is removed and does not rise when it is locked in the fully flipped state, the claimed entropy-stabilization and wedge mechanism would be refuted.","supporting_citations":[{"cited_title":"(2017) Ensemble cryo-EM elucidates the mechanism of translation fidelity","cited_arxiv_id":null,"evidence_quote":"Supplies the pre-accommodation ribosome-tRNA recognition intermediate structure used to build every A-site simulation system."},{"cited_title":"and Ramakrishnan,V","cited_arxiv_id":null,"evidence_quote":"Defines the open-to-closed domain-closure model of decoding that the paper's discrimination argument builds on and partly supports."},{"cited_title":"and Yusupova,G","cited_arxiv_id":null,"evidence_quote":"Provides the rigid-decoding-center near-cognate structures and the claim that the first codon nucleotide is constrained against wobble, which the paper contests with its wobble U-G result."},{"cited_title":"and Brooks,C.L","cited_arxiv_id":null,"evidence_quote":"Gives the prior molecular-dynamics result of cognate-specific A1492/A1493 flipping whose 3 kcal/mol difference the present study attributes to modeling the closed state rather than the recognition intermediate."},{"cited_title":"and Koutmou,K.S","cited_arxiv_id":null,"evidence_quote":"Documents pseudouridine-induced amino acid substitutions that the paper's pseudouridine water-bridge mechanism is built to explain."},{"cited_title":"and Korostelev,A.A","cited_arxiv_id":null,"evidence_quote":"Provides the post-translocation and time-resolved structures used for the P-site stability and translocation-destabilization analysis."},{"cited_title":"and Erlacher,M","cited_arxiv_id":null,"evidence_quote":"Atomic mutagenesis showing that minor-groove hydrogen-bond disruption barely affects discrimination, cited as evidence for steric complementarity over A-minor hydrogen bonding."},{"cited_title":"and Joseph,S","cited_arxiv_id":null,"evidence_quote":"Shows that steric complementarity in the decoding center is important for tRNA selection, supporting the wedge interpretation."}],"review_version":1}