{"id":"5a2da113-4e4e-4fec-9255-bb25f3c446af","arxiv_id":"2411.13445","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Cholesteric handedness and pitch in virus suspensions are quantitatively explained by electrostatic surface-charge chirality for stiff Y21M and by fluctuation-induced suprahelical backbone deformations for flexible M13.","lead":"This paper shows that the helical twist of liquid crystals made from rod-shaped viruses is controlled by two different microscopic mechanisms: electric charges on the virus surface for stiff viruses, and thermal bending of flexible viruses into loose helices for flexible ones. The result offers a quantitative route from nanometer-scale molecular structure to millimeter-scale helical order in chiral soft materials.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"M13 quantitative agreement relies on a fitted suprahelical pitch h=2.8Lp and an inferred handedness; without independent determination, the M13 branch of the central claim is a fit, not a prediction.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing concern: the M13 branch of the central claim rests on a free parameter (h = 2.8Lp) that is chosen to reproduce the target observable, and on an unverified handedness assumption. My independent reading confirms this is the most fragile step. The paper's Y21M electrostatic model is impressive and parameter-free in the sense that no pitch data are used to calibrate it, and its agreement with experiments at pH 5 and 8 (Fig. 3) provides real support for the electrostatic mechanism. However, the central claim is explicitly about the interplay between two mechanisms, and the M13 suprahelix model carries the full weight for the flexible-virus branch. Setting h = 2.8Lp after the fact means that the quantitative agreement of Fig. 5 is a fit, not a prediction, for the very system it is supposed to explain. The master-curve collapse is non-trivial and does add evidence that the mechanism is concentration-independent, but it does not independently fix h; any h that matches the pitch at one concentration would likely also reproduce the scaling if the model's functional form is right. The Y21M-PEG nematic result is a null prediction (pitch too large to observe) and therefore provides only weak constraint on h. Thus the central claim's quantitative force for M13 is not yet established. A direct measurement of the single-filament helical conformation would settle whether the assumed suprahelix parameters are physically correct. Given this unresolved but plausible concern, the CONDITIONAL verdict is appropriate, and my stress-test does not change it.","tokens_in":16014,"tokens_out":5579,"duration_ms":63490,"concrete_test":"Measure the full three-dimensional conformation of individual M13 filaments (e.g., by cryo-electron tomography or by tracking a fluorescent label along the contour) and determine the handedness and internal pitch h of the thermally fluctuating backbone. If the measured deformation is not a right-handed helix with h within about 30% of 2.8Lp, then the M13 suprahelix mechanism is not the quantitative origin of the measured cholesteric pitch, and the master-curve agreement is a fit rather than a prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the M13 cholesteric pitch is quantitatively explained by fluctuation-induced suprahelical backbone deformations. In Section II, the internal pitch h of the suprahelix is set such that h = 2.8Lp, i.e., chosen to match the measured M13 pitch, and the right-handedness of the deformation is inferred from the measured left-handedness of the phase via a geometric argument, not from direct observation. Consequently, the quantitative agreement in Fig. 5 is a one-parameter fit for the M13 branch; the only independent predictions are the shape of the master curve and the null result for stiff Y21M-PEG (divergent pitch). If the actual backbone deformation has a different h or handedness, the model's agreement is coincidental, and the claimed quantitative understanding of chirality transfer for semi-flexible viruses is not established. The Y21M electrostatic branch is stronger because it involves no fitted parameters targeting the pitch, but the overall 'interplay' claim depends on both branches being quantitative.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates the origin of cholesteric handedness and pitch in liquid crystals of two closely related filamentous bacteriophages, M13 and Y21M, which form cholesteric phases of opposite handedness. The authors combine experiments with two theoretical models: an atomistic electrostatic model that treats the capsid surface charges in detail, and a suprahelix model that represents thermally induced backbone deformations as a weakly curled hard helix. They report quantitative agreement for stiff Y21M using the electrostatic model with no pitch-specific fitting, and for semi-flexible M13 using the suprahelix model with internal pitch set to h = 2.8Lp, which collapses the M13 data onto a master curve. The central claim is that cholesteric self-assembly in these viruses quantitatively results from the interplay of electrostatic surface-charge chirality and fluctuation-induced backbone helicity.","tokens_in":16286,"tokens_out":3606,"duration_ms":40263,"significance":"If the claims hold, this is a substantial step toward a quantitative, bottom-up understanding of chirality transfer in colloidal liquid crystals. The Y21M electrostatic branch is particularly strong: it uses externally determined PDB structures, a standard force field, and standard protonation tools, and it reproduces the sign, magnitude, and pH-dependent unwinding of the pitch without fitted parameters. The master-curve collapse for M13 over multiple ionic strengths is also a valuable empirical result. The main weakness is that the M13 branch depends on a chosen internal pitch h = 2.8Lp and an assumed right-handedness of the backbone deformation, so the quantitative agreement in that branch is partly a fit rather than an independent prediction. Nonetheless, the paper presents a coherent, well-structured framework and will likely stimulate further work on chirality propagation in semiflexible biopolymers.","major_comments":[{"comment":"The central quantitative claim for M13 rests on the internal pitch h of the suprahelical conformation being set to h = 2.8Lp, a value chosen to match the measured M13 cholesteric pitch. This is explicitly stated in the text after Fig. 4b: \"the internal pitch h of the suprahelical conformation has been set such that h = 2.8Lp.\" Consequently, the agreement in Fig. 5 is a one-parameter fit for the M13 branch, not an independent prediction. The paper should clearly state which aspects are predictive (for example, the Y21M-PEG diverging pitch, the master-curve collapse, and the predicted scaling with contour length) and which are fitted, and should discuss how h could be determined from first principles or from independent single-filament measurements. Without such clarification, the abstract's claim that the assembly \"quantitatively results\" from the interplay is overstated.","section":"Section II, suprahelix model (Fig. 4b and Fig. 5)"},{"comment":"The right-handedness of the backbone deformation is not measured but inferred from the left-handed cholesteric phase and from geometric packing arguments (Extended Data Fig. 3). The paper itself acknowledges that \"no primary proof of such helical conformation has been reported yet.\" Because the sign of the M13 cholesteric pitch is fully determined by this assumed handedness, the M13 branch is not a complete prediction. The authors should make this assumption more prominent and propose a concrete experimental test, such as cryo-electron tomography of individual M13 filaments, that could falsify or confirm the assumed suprahelical chirality.","section":"Section II, handedness of the suprahelix (paragraph after Fig. 4)"},{"comment":"The electrostatic model uses a cutoff radius rcut = 3.5 nm for all electrostatic energies, justified by the condition kappa^{-1} << rcut, which holds at IS >= 100 mM. However, the Y21M comparison in Fig. 3 is only shown at IS = 110 mM. The sensitivity of the computed pitch to the choice of rcut is not reported, so the \"parameter-free\" characterization of the electrostatic model is not fully supported. A short test varying rcut or a statement of the observed insensitivity would strengthen the claim.","section":"Methods, numerical methods (cutoff radius rcut)"}],"minor_comments":[{"comment":"There is a typographical error \"in vitroin\" in the first paragraph of the introduction; it should read \"in vitro in.\"","section":"Introduction and Methods"},{"comment":"The text refers to the \"IIFI model\" in one place; this should be the \"1IFI model\" as used elsewhere in the paper.","section":"Methods, virus strains and capsid symmetries"},{"comment":"The Y21M electrostatic model is compared with experiment only at IS = 110 mM, while Extended Data Fig. 2 shows Y21M data at several ionic strengths. A sentence explaining why the electrostatic model is not directly compared at the other ionic strengths (for example, due to the rcut limitation) would help the reader assess the scope of the validation.","section":"Fig. 3 and Extended Data Fig. 2"},{"comment":"The paper states that error bars for computed pitches come from O(10) independent Monte-Carlo runs, but the error bars are not shown in Fig. 3 or Fig. 5. Adding representative error bars to the theoretical points would make the quantitative agreement more transparent.","section":"Methods, numerical methods"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong candidate for the journal, and the Y21M electrostatic branch alone is a solid contribution. The requested revision is not a rejection of the central mechanism, but a request for honest labeling of fitted versus predicted quantities in the M13 branch and for a concrete route to independently validate h and the backbone handedness. The authors themselves disclose the h = 2.8Lp choice, so this is a framing and evidence issue rather than a hidden error. I would be comfortable with acceptance after these revisions."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear Alex,\n\nTwo things about this one. The Y21M half is genuinely strong: a parameter-free all-atom electrostatic model reproduces the cholesteric pitch, handedness, and pH-driven unwinding of stiff Y21M rods, with no fitted parameter aimed at the pitch. That alone is a real advance. The M13 half is weaker: the suprahelix model gets its quantitative agreement by setting the internal pitch h = 2.8 Lp, chosen to match the measured pitch, and the right-handedness of the backbone deformation is inferred from the observed left-handed phase rather than measured. So the stress-test note is right: the M13 branch is a one-parameter fit, not an independent prediction.\n\nWhat's new here is the combination: two mechanisms in one framework, the master-curve collapse for M13 across ionic strengths and PEGylation, and the null result for stiff Y21M-PEG. The collapse itself does not depend on the fitted h; h only sets the vertical scale. The authors are also open about the limitation—they state the choice of h and admit there is no direct proof of the helical backbone. The Y21M model is anchored to PDB structures and a standard force field, so it has real external grounding.\n\nSoft spots: the fitted h and the assumed handedness are the main ones. I'd also note the electrostatic cutoff rcut = 3.5 nm is a simplification, though it's probably fine at the IS ≥ 100 mM conditions where the Y21M tests are done. The second-virial approximation visibly underestimates chirality at high concentration; they mention it, so it's not hidden. Minor: code is promised on GitHub but with no URL or commit hash, and I couldn't see the supplementary derivations, so full verification of the M13 model requires trusting the authors.\n\nThis is for anyone interested in cholesteric liquid crystals of biopolymers or chirality transfer generally. The Y21M result alone justifies a serious referee. The M13 branch should be treated as a plausible hypothesis with a nice master curve, not a closed case.\n\nI'd send it to review. I'd want the authors to either derive h or label it explicitly as empirical, and to discuss the handedness inference more carefully. With that, it could be an important paper.","headline":"Y21M electrostatic model is a real, parameter-free quantitative result; M13 suprahelix branch is a plausible one-parameter fit, so the paper deserves review but with the M13 caveat.","tokens_in":16776,"tokens_out":3101,"would_cite":true,"duration_ms":33118,"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":"This paper claims that chirality transfer in virus liquid crystals splits into two mechanisms—electrostatic surface charge patterns for stiff rods, thermally induced backbone coiling for flexible rods—and supports the split with…","keywords":["cholesteric liquid crystals","chirality transfer","filamentous bacteriophages","M13 virus","Y21M virus","electrostatic surface charges","suprahelical backbone fluctuations","cholesteric pitch master curve"],"falsifier":"Directly image the backbone of M13 viruses in cholesteric suspensions (for example by cryo-electron tomography or by tracking fluorescent labels along single filaments) and look for the presumed right-handed superhelix with internal pitch near 2.8 Lp ≈ 8 µm. If the thermally fluctuating backbone is straight, left-handed, or coiled at a clearly different pitch, the master-curve collapse has to be reinterpreted as a fit rather than evidence for the suprahelix mechanism.","tokens_in":15779,"feed_emoji":"🦠","tokens_out":13124,"duration_ms":132198,"temperature":0.7,"pith_summary":"Two closely related filamentous viruses, M13 and Y21M, form cholesteric liquid crystals—ordered fluids whose orientation twists in a helix—in opposite directions, and this paper claims to explain that difference quantitatively. For the stiff Y21M strain, the handedness and pitch of the twist (the distance for one full turn) are set by electrostatic interactions between the helical patterns of charge on the virus surface, described by an all-atom model built from the capsid structure. For the flexible M13 strain, those surface details are irrelevant; thermal motion bends the filament into a coherent right-handed superhelix, and the excluded-volume interactions of these coiled shapes set the twist. The authors show that all M13 pitch data collapse onto one master curve when concentration is rescaled, and that a 'suprahelix' with internal pitch fixed at 2.8 times the persistence length (a stiffness length scale) reproduces that curve. If right, this gives a general route from molecular chirality to macroscopic helical order and explains why nearly identical systems can show opposite chiral behavior.","feed_headline":"Flexible viruses twist by shape, stiff ones by surface charge","feed_subtitle":"M13's cholesteric pitch follows a universal coiling curve; Y21M's is set by capsid electrostatics.","key_machinery":"The argument runs on two machines. The first is an atomistic pair-interaction model for the capsids: starting from the deposited 1IFI and 2C0W capsid structures, every atom pair contributes screened electrostatic, van der Waals, and excluded-volume terms, and the cholesteric pitch and twist elastic constant are obtained by minimizing the second-virial free energy—a low-density expansion in pairwise interactions—of the twisted nematic state. The second is the suprahelix model, in which the thermally fluctuating backbone of a semi-flexible virus is replaced by a hard helix with radius r and internal pitch h; for a given persistence length (stiffness length scale) Lp, the internal pitch is the sole adjustable parameter, and the measured M13 behavior fixes it at h = 2.8 Lp. Both models feed the same free-energy machinery, so the two chirality mechanisms are compared on equal footing. The master-curve collapse is the key output connecting the suprahelix model to experiment.","core_discovery":"The central claim is that chirality transfer in filamentous-virus cholesterics is not dominated by a single mechanism but by two mechanisms that act at different length scales and can have opposite handedness. For the nearly rigid Y21M strain, a fully atomistic inter-particle potential—screened electrostatics, van der Waals, and steric forces summed over the roughly three million atoms of the capsid—quantitatively reproduces the measured right-handed cholesteric pitch and its dependence on pH and ionic strength, including the unwinding of the cholesteric order as the surface charge is reduced. For the semi-flexible M13 strain, the same electrostatic calculation fails, and the paper attributes the chirality instead to long-wavelength, thermally driven helical deformations of the virus backbone. Representing these fluctuations by a hard 'suprahelix' of radius r and internal pitch h, with h set to 2.8 Lp, collapses the measured inverse pitches of charged M13 and PEGylated M13-PEG at all ionic strengths onto a single master curve as a function of concentration rescaled by the isotropic binodal Ciso. The conclusion is that stiff-rod chirality is set by surface charge geometry, while flexible-rod chirality is set by entropy-driven backbone coiling, and the competition between these routes accounts for the opposite handedness of two otherwise similar viruses.","pith_inferences":["If the h = 2.8 Lp relation is universal rather than a fit, then the preferred thermal coil wavelength of any semi-flexible filament should be proportional to its persistence length; this could be tested directly by imaging individual M13 filaments under conditions that suppress or enhance thermal fluctuations.","A mutant or solvent condition that tunes M13's persistence length across the boundary between the two regimes should show a cholesteric pitch that first diverges and then changes sign—an experimentally accessible crossover the paper does not report.","The electrostatic model's success for Y21M implies a testable prediction: mutating specific charged residues on the p8 coat protein should shift the pitch quantitatively in the direction and magnitude the all-atom model computes before the experiment is done.","If both mechanisms operate simultaneously in intermediate-stiffness filaments, the observed pitch could be a sensitive probe of their relative sign and magnitude, potentially explaining why some cellulose and amyloid cholesterics show non-monotonic salt behavior."],"forward_implications":["For stiff rod-like biopolymers with known atomistic structure, the cholesteric pitch and its salt/pH response can be predicted from the capsid charge pattern alone, without invoking backbone flexibility.","For semi-flexible filaments, the detailed surface charge decoration is irrelevant to chirality: only stiffness, length, and effective diameter matter, so the pitch should follow the same master curve when concentration is rescaled by the isotropic binodal.","Because the two mechanisms can contribute with opposite handedness, mutations or conditions that change stiffness, not just surface chemistry, can flip the macroscopic handedness of a cholesteric phase.","The pitch of flexible-filament cholesterics should tighten as contour length increases, opposite to the classical prediction for rigid screws, providing an experimental signature in cellulose, amyloid, and DNA-origami systems.","Chirality studies on charged biopolymer suspensions should compare data at fixed reduced concentration C/Ciso; otherwise apparent pitch changes with ionic strength may be misread as chiral surface-charge effects rather than generic electrostatic effects on phase stability."],"supporting_citations":[{"why":"Supplies the refined atomic structure of the M13/fd capsid used as input for the electrostatic model.","marker":"[33]"},{"why":"Defines the capsid symmetries and thread angles distinguishing M13 from Y21M.","marker":"[35]"},{"why":"Frames the long-standing chirality-origin question for virus cholesterics and proposes the backbone-fluctuation idea.","marker":"[17]"},{"why":"Provides the earlier coarse-grained electrostatic/chiral model whose quantitative failure motivates the atomistic approach.","marker":"[18]"},{"why":"Supplies the fluctuation-based chiral shape framework and the free-energy minimization used for the suprahelix model.","marker":"[30]"},{"why":"Provides the measured persistence lengths and the Y21M cholesteric-pitch data used as experimental baselines.","marker":"[40]"},{"why":"Establishes that PEGylated viruses are sterically stabilized and independent of ionic strength, enabling the decoupling experiments.","marker":"[49]"},{"why":"Gives the second-virial excluded-volume theory underlying both pitch calculations.","marker":"[46]"},{"why":"Supplies the all-atom force-field parameters used to compute capsid pair interactions.","marker":"[48]"}],"fun_headline_variants":["Virus crystals: chirality from charge or coiling","Stiff rods twist by charge, flexible ones by coiling","Two mechanisms explain virus helix handedness","How viruses transfer chirality: charge vs. backbone flexibility"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The M13 result stands or falls on the assumption that thermal wiggling of the flexible virus is equivalent to one right-handed coiled shape whose coil spacing is fixed at 2.8 times the virus's stiffness length—a number chosen to match the measured pitch rather than derived from the physics; if that effective shape or its handedness is wrong, the agreement is a fit rather than a prediction.","fun_headline_variants_meta":{"raw":{"variants":["Virus crystals: chirality from charge or coiling","Stiff rods twist by charge, flexible ones by coiling","Two mechanisms explain virus helix handedness","How viruses transfer chirality: charge vs. backbone flexibility"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00021,"raw_usage":{"total_tokens":1429,"prompt_tokens":982,"completion_tokens":447,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":598,"completion_tokens_details":{"reasoning_tokens":383}},"tokens_in":598,"tokens_out":447,"duration_ms":5206,"temperature":1.0,"reasoning_tokens":383,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T16:24:26.016831+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Directly image the backbone of M13 viruses in cholesteric suspensions (for example by cryo-electron tomography or by tracking fluorescent labels along single filaments) and look for the presumed right-handed superhelix with internal pitch near 2.8 Lp ≈ 8 µm. If the thermally fluctuating backbone is straight, left-handed, or coiled at a clearly different pitch, the master-curve collapse has to be reinterpreted as a fit rather than evidence for the suprahelix mechanism.","supporting_citations":[{"cited_title":"Strano, Gerbrand Ceder, and Angela M","cited_arxiv_id":null,"evidence_quote":"Defines the capsid symmetries and thread angles distinguishing M13 from Y21M."},{"cited_title":"Smith and Valery A","cited_arxiv_id":null,"evidence_quote":"Supplies the refined atomic structure of the M13/fd capsid used as input for the electrostatic model."},{"cited_title":"Cholesteric phase in virus suspensions","cited_arxiv_id":null,"evidence_quote":"Frames the long-standing chirality-origin question for virus cholesterics and proposes the backbone-fluctuation idea."},{"cited_title":"What is the origin of chirality in the cholesteric phase of virus suspensions? Phys","cited_arxiv_id":null,"evidence_quote":"Provides the earlier coarse-grained electrostatic/chiral model whose quantitative failure motivates the atomistic approach."},{"cited_title":"Entropy-driven formation of chiral nematic phases by computer simulations","cited_arxiv_id":null,"evidence_quote":"Supplies the fluctuation-based chiral shape framework and the free-energy minimization used for the suprahelix model."},{"cited_title":"Biomimetic self-templating supramolecular struc- tures","cited_arxiv_id":null,"evidence_quote":"Provides the measured persistence lengths and the Y21M cholesteric-pitch data used as experimental baselines."},{"cited_title":"Mark, and Wilfred F","cited_arxiv_id":null,"evidence_quote":"Establishes that PEGylated viruses are sterically stabilized and independent of ionic strength, enabling the decoupling experiments."},{"cited_title":"Wensink and G","cited_arxiv_id":null,"evidence_quote":"Gives the second-virial excluded-volume theory underlying both pitch calculations."},{"cited_title":"Zhang, N","cited_arxiv_id":null,"evidence_quote":"Supplies the all-atom force-field parameters used to compute capsid pair interactions."}],"review_version":1}