{"id":"d0bb7f65-052d-44d8-8579-85816894596a","arxiv_id":"2608.02525","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"A generalized Poisson-Boltzmann framework with self-energy cancellation enables polar solvation energy calculations for nonrigid proteins via a thermodynamic cycle correction.","lead":"Proteins change shape when they move from vacuum into water, but standard Poisson-Boltzmann solvation-energy formulas assume a fixed structure. This paper derives a correction term that lets existing solvers handle such non-rigid proteins, and validates it on a two-atom system and 70 proteins.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Accuracy claim rests on Pearson correlation, which is not an agreement metric; Table 2 suggests C-APBS may actually be closer to REG.","rationale":"The reader's verdict was CONDITIONAL with medium risk, and the reader's rationale already flagged the use of Pearson correlation rather than absolute error as a numerical validation issue. However, the reader's stated weakest_assumption concerns the theoretical assumption that partial charges and atomic radii are unchanged between vacuum and water conformations—a condition that is explicit and standard in the model. My stress-test identifies a different load-bearing concern: the empirical accuracy claim is not supported by the reported statistics, and Table 2 suggests the conclusion about C-DelPhi may be an artifact of the correlation metric. This supports the CONDITIONAL verdict (or an UNCHANGED verdict relative to the reader), since the paper should be revised to use an appropriate agreement metric before asserting which solver is more accurate. I do not find a fatal mathematical flaw in Eq. (6) itself; the self-energy cancellation is a standard renormalization, and the thermodynamic cycle is consistent with the PB electrostatic model. The main weakness is the disconnect between the strong 'accuracy' claim and the evidence presented.","tokens_in":8152,"tokens_out":19307,"duration_ms":200433,"concrete_test":"Recompute the comparison of Example 2 across all 70 proteins using an agreement metric rather than Pearson correlation: report the mean absolute deviation (MAD) and root mean square error (RMSE) of C-APBS and C-DelPhi relative to REG, along with the slope and intercept of the linear regression. Also compute Lin's concordance correlation coefficient (which penalizes both correlation and bias). If C-APBS has lower MAD/RMSE, the statement that C-DelPhi is more accurate should be revised. As a secondary check, rerun the analysis excluding the four proteins that failed to converge and report the sensitivity of the correlation values.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central derivation of Eq. (6) appears internally consistent, but the paper's empirical claim that the correction procedure is 'accurate and robust' across different PB solvers—and specifically that 'C-DelPhi is shown to be more accurate than C-APBS' (Section 5)—is not supported by the statistics reported. The evidence consists solely of Pearson correlation coefficients (0.99 for C-DelPhi vs REG, 0.83 for C-APBS vs REG). Pearson correlation measures linear association, not agreement; a systematic offset or a different slope yields r close to 1 while absolute errors remain large. In fact, the 10-protein sample in Table 2 shows C-APBS has smaller absolute deviation from REG than C-DelPhi for 6 of the 10 proteins (e.g., 1TG0: 116 vs 248 kcal/mol; 1W0N: 178 vs 344; 3LZT: 148 vs 349 kcal/mol). Thus the conclusion that C-DelPhi is more accurate may be reversed. This is load-bearing because the paper's validation is entirely internal—all methods are variations of the same PB theory, and no comparison to experiment, explicit solvent, or an independent analytical benchmark is offered. If the accuracy conclusion is an artifact of the chosen metric, the paper's claim of demonstrating accuracy across solvers is weakened, even though the theoretical formula may remain valid.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a generalization of the Poisson-Boltzmann (PB) theory to compute polar solvation energies of proteins whose structures differ between the vacuum and water states. For regularized PB models, the authors derive an expression (Eq. (6)) in which the divergent self-energy terms cancel analytically even though atom centers differ, yielding a finite sum of reaction-field, ionic, and Coulombic interaction terms. For non-regularized solvers (APBS, DelPhi), they introduce a thermodynamic cycle (Eq. (7) plus a rigid-body solvation step) and implement it as a post-processing correction. The methods are tested on a perturbed two-atom system and a set of 70 proteins, with comparisons among several PB solvers.","tokens_in":8434,"tokens_out":10658,"duration_ms":112734,"significance":"The central theoretical contribution—removing the rigidity assumption in PB solvation energy calculations—is significant and timely. The derivation is parameter-free and the proposed correction is directly applicable to widely used solvers such as APBS and DelPhi, giving the paper potential practical impact. The two-atom convergence test is a useful sanity check, and the availability of the correction code on GitHub is a strength. However, the numerical validation is entirely internal: all compared methods are PB-based, including the new REG scheme, so the 'accuracy' claims are not benchmarked against an external reference. One specific claim—that C-DelPhi is more accurate than C-APBS—is not supported by the reported statistics and may even be reversed in the 10-protein subset of Table 2.","major_comments":[{"comment":"The claim that C-DelPhi is more accurate than C-APBS is not supported by the evidence. Pearson correlation coefficients measure linear association, not agreement; a high r can coexist with large systematic offsets. In the 10-protein subset of Table 2, C-APBS is closer to REG in absolute deviation in 7 of 10 cases (e.g., 1TG0: 116 vs 248 kcal/mol; 1W0N: 178 vs 344; 3LZT: 148 vs 349 kcal/mol). Please report mean absolute errors, slopes, intercepts, or Bland-Altman statistics, and revise the conclusion accordingly.","section":"Section 5; Section 4, Figure 2 and Table 2"},{"comment":"The numerical validation demonstrates internal consistency (C-DIPB agrees with REG because they share the same diffuse-interface model; corrected sharp-interface results are bounded), but it does not establish 'accuracy' in an absolute sense because there is no external reference—all methods are PB-based and REG itself is new. The abstract's claim 'demonstrating its accuracy' is overstated. Please compare against an independent benchmark (e.g., explicit-solvent simulations or a model with an analytical solution) or temper the wording to 'internal consistency'.","section":"Section 4, Example 2; Abstract"},{"comment":"The cancellation of self-energies is stated rather than justified. The terms (1/2)e_c^2 Σ_i z_i^2/(ε_m |r_i - r_i|) are formally infinite, and subtracting two divergent series is not an algebraic operation. The paper should present this as a renormalization/regularization step (e.g., introduce a cutoff δ, show the divergent parts are position-independent, and then take the limit) or cite the established treatment in regularized PB theory. As written, the derivation may appear mathematically incomplete.","section":"Section 2, Eqs. (5)-(6)"},{"comment":"The numerical evidence for the full protein set is weakened by the exclusion of four proteins (1MC2, 2FDN, 2H5C, 4TKB) due to non-convergence, with no discussion of how this affects representativeness, and by the absence of convergence tests or error estimates for the 70-protein calculations. All protein results are reported at a single grid resolution. Please provide a grid-convergence study for a few representative proteins or otherwise quantify discretization error.","section":"Section 4, Example 2"}],"minor_comments":[{"comment":"The assumption that partial charges z_j and radii R_j are unchanged between the two conformations is critical for Eq. (6). This should be explicitly identified as a limitation, since force fields may assign charges based on protonation states that can change with conformation.","section":"Section 2"},{"comment":"The random perturbation in Eq. (9) uses uniform draws γ_{j,i}; for reproducibility, please provide a random seed or report the specific random configurations used.","section":"Section 4, Example 1"},{"comment":"The term 'valence' for the partial charge z_j is nonstandard; consider using 'partial charge' or defining 'valence' explicitly.","section":"Throughout"},{"comment":"The right half of Table 1 lists h values (0.8, 0.6, 0.4, 0.2, 0.1) in descending order, which is fine, but the row labels are not aligned with the columns; this makes the table harder to read. Consider splitting into two tables or adding explicit column headers.","section":"Section 4, Table 1"}],"recommendation":"major_revision","confidential_remarks":"The theoretical derivation is the core contribution and appears sound, but the empirical validation—particularly the relative accuracy claim for C-DelPhi vs C-APBS—needs substantial revision. The authors should also strengthen the mathematical justification for the self-energy cancellation. With these fixes, the paper could be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The theoretical core is genuinely new and looks correct. Eq. (6) removes the rigidity assumption in regularized PB by canceling the singular self-energies even when the water and vacuum structures differ. That works because the Coulomb self-energy of a point charge is position-independent, so the i=j terms in the two sums cancel exactly. The thermodynamic-cycle correction for non-regularized solvers is simple, practical, and implementable as post-processing on APBS or DelPhi. The derivation is clean, and the two-atom tests show the correction prevents the divergence seen in the uncorrected methods. The code is on GitHub. These are real contributions.\n\nThe main soft spot is the empirical accuracy claim. The paper says C-DelPhi is more accurate than C-APBS based on Pearson correlations of 0.99 vs 0.83 against REG. Pearson correlation is not an agreement metric; a systematic offset or a different slope still gives r near 1. In Table 2, the absolute deviation from REG is smaller for C-APBS in 7 of the 10 listed proteins (e.g., 1TG0: 116 vs 248 kcal/mol; 3LZT: 148 vs 349 kcal/mol). By that metric, the conclusion would reverse. The correlation values reflect a different spread, not tighter agreement. That claim needs to be reworked with proper agreement measures (mean absolute error, RMSD) and ideally a Bland-Altman-type analysis.\n\nOther weaker points: the validation is entirely internal—all methods are variations of the same PB theory, with no comparison to experiment or explicit solvent. The exclusion of four proteins because solvers fail to converge is mentioned but not discussed as a potential selection bias. There are no error bars or grid-convergence thresholds for the protein results. And Eq. (7) assumes the vacuum deformation free energy is purely Coulombic, which is consistent with the model but should be stated more prominently as a limitation.\n\nNone of this invalidates the main derivation. The formula in Eq. (6) is a genuine extension of regularized PB, and the correction scheme is practically useful. The paper deserves a serious referee. I would send it to peer review, but with the expectation of revision: fix the accuracy comparison, add convergence checks and uncertainty quantification, and soften the language from \"accurate and robust\" to something the evidence actually supports.","headline":"The analytical formula is solid and worth publishing; the empirical accuracy claim rests on a misused statistic and should be fixed.","tokens_in":8928,"tokens_out":3389,"would_cite":true,"duration_ms":37178,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N06","65N12","92C40","35Q92"],"pacs":[],"model":"deepseek-v4-flash","headline":"The polar solvation energy of a protein whose shape differs between vacuum and water can be computed in the Poisson-Boltzmann framework by analytically canceling divergent self-energies—or, in ordinary solvers, by adding a Coulombic correct","keywords":["Poisson-Boltzmann equation","polar solvation energy","self-energy cancellation","regularization","nonrigid proteins","conformational change","implicit solvent","electrostatic free energy"],"falsifier":"Compute the solvation energy for a two-atom system with identical charges and radii but a large positional shift, using the corrected non-regularized scheme and a fully regularized scheme; they should agree. A more decisive test: repeat the calculation with different partial charges in the two conformations—if the formula still returns a finite number, the cancellation logic is being misapplied, because the z_i^2 self-energies would differ.","tokens_in":8023,"feed_emoji":"🧬","tokens_out":9170,"duration_ms":87582,"temperature":0.7,"pith_summary":"This paper tries to remove the rigidity assumption that has constrained Poisson-Boltzmann (PB) solvation energy calculations: until now, the same protein structure had to be used in vacuum and water so that infinite self-energies could cancel. The authors show that even when atom centers differ between the two states, each atom's self-energy cancels exactly because it depends only on the atom's charge and not on its position. For regularized PB models, this yields a finite closed formula for the nonrigid solvation energy. For non-regularized PB solvers, they propose a thermodynamic cycle that adds a Coulombic correction term for the vacuum deformation. If correct, standard PB packages can handle nonrigid proteins via a simple post-processing step.","feed_headline":"Canceling self-energies lets PB theory handle nonrigid proteins","feed_subtitle":"A Coulombic correction gives finite solvation energies for proteins whose shape changes from vacuum to water.","key_machinery":"The central object is the regularized decomposition of the dimensionless potential, u = u_RF + G, where G is the Coulomb Green's function carrying the point-charge singularities and u_RF solves a smooth Poisson-Boltzmann equation. Cancellation of the diagonal self-energy terms in the two double sums converts a divergent expression into the finite formula (Equation 6). The companion device is a thermodynamic cycle that splits the nonrigid solvation energy into a vacuum Coulombic deformation energy plus a rigid-body PB solvation energy, so any ordinary solver needs only a post-processing correction.","core_discovery":"For a protein with atom centers r_i in water and r-hat_i in vacuum (same partial charges z_i and radii), the paper decomposes the water potential into a reaction-field part and a Coulomb Green's function, and notes the vacuum potential is exactly the Coulomb Green's function of the vacuum structure. The solvation energy then contains two double Coulomb sums. The diagonal (i=j) terms are the singular self-energies; although the centers differ, each self-energy is z_i^2/(epsilon_m |r_i-r_i|) in water and z_i^2/(epsilon_m |r-hat_i-r-hat_i|) in vacuum, both the same divergent quantity, so they cancel identically. Removing them leaves finite pairwise Coulombic terms plus a bounded reaction-field","pith_inferences":["The self-energy cancellation logic may extend to other point-charge implicit-solvent energy formulations, wherever the singular self-term is independent of the atom's location.","The thermodynamic cycle treats the conformational change as purely electrostatic in vacuum; coupling this with mechanical or entropic deformation estimates could give a more complete transfer free energy.","Because the correction is a post-processing step, it is straightforward to compute ensemble-averaged solvation energies over multiple conformations by averaging the Coulombic corrections."],"forward_implications":["Regularized PB models can now accept different vacuum and water structures without diverging, eliminating the rigidity assumption.","Any existing non-regularized PB solver can yield nonrigid solvation energies by adding the analytic Coulombic correction as post-processing.","Numerical tests on a perturbed two-atom system and 70 proteins show corrected energies stay finite and closely track the regularized scheme (correlations of 0.99 and 0.83).","The correction applies uniformly across sharp-interface and diffuse-interface PB models and across different numerical solvers."],"fun_headline_variants":["PB theory now handles nonrigid proteins via self-energy cancellation","Self-energy cancellation enables PB for proteins that change shape","Generalized PB framework for proteins with different vacuum/water structures","Coulombic correction brings PB solvation to nonrigid proteins","Nonrigid proteins get accurate PB energies with self-energy cancelation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The cancellation relies on each atom keeping the same partial charge and radius in vacuum and water; if charges or radii change between conformations, the infinite self-energies no longer cancel, and Equation (7) would ignore non-Coulombic deformation contributions.","fun_headline_variants_meta":{"raw":{"variants":["PB theory now handles nonrigid proteins via self-energy cancellation","Self-energy cancellation enables PB for proteins that change shape","Generalized PB framework for proteins with different vacuum/water structures","Coulombic correction brings PB solvation to nonrigid proteins","Nonrigid proteins get accurate PB energies with self-energy cancelation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1156,"prompt_tokens":734,"completion_tokens":422,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":478,"completion_tokens_details":{"reasoning_tokens":336}},"tokens_in":478,"tokens_out":422,"duration_ms":5484,"temperature":1.0,"reasoning_tokens":336,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T05:26:10.829121+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the solvation energy for a two-atom system with identical charges and radii but a large positional shift, using the corrected non-regularized scheme and a fully regularized scheme; they should agree. A more decisive test: repeat the calculation with different partial charges in the two conformations—if the formula still returns a finite number, the cancellation logic is being misapplied, because the z_i^2 self-energies would differ.","supporting_citations":[],"review_version":1}