{"id":"b90166f8-102e-4ee0-a571-a3ae2f87f0a6","arxiv_id":"2505.07273","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Numerical simulations suggest that strong clustering of CaV1.3 calcium channels can synchronize multivesicular release and increase monophasic EPSC amplitude in inner hair cells.","lead":"This paper simulates the full chain of events in an auditory hair cell synapse, from calcium channel opening to postsynaptic current, and asks why some signals are unusually large. It argues that tight clustering of calcium channels could make multiple vesicles release together, creating large signals, though the simulation only partially matches experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Simulation does not reproduce the required monophasic/multiphasic amplitude ordering; the strong-clustering conclusion extrapolates beyond the computed parameter range.","rationale":"I focused on the amplitude-ordering failure because it is internal to the paper and directly bears on the headline claim. The reader's chosen weakest assumption—that tight CaV1.3 clusters exist in vivo—is real but less decisive: the paper's conclusion is a 'can' claim, and a simulation may legitimately posit a structure to test a mechanism. What is not legitimate is to conclude 'can' when the simulated strong-clustering case still gives monophasic amplitudes below multiphasic ones, contrary to the experimental fact the mechanism is supposed to explain. That mismatch is acknowledged in the authors' own Fig. 2c and in the Section II text, so the stress-test concern is not about consensus; it is an internal inconsistency between the plotted results and the stated conclusion. I still keep the reader's CONDITIONAL verdict rather than moving to REJECT because the missing evidence is a well-defined, feasible computation: a parameter sweep over clustering strength would show whether the claimed crossover exists. If that sweep fails, the paper's main conclusion would be unsupported. I also note the unresolved eta correction and missing code, but those are secondary to the fact that the central effect is not yet demonstrated.","tokens_in":9947,"tokens_out":7407,"duration_ms":71143,"concrete_test":"Run a systematic sweep of the strong-clustering geometry from Fig. 3 (e.g., inter-channel spacings of 25, 15, 10, 5, 2, and 1 nm, preserving the 80-channel/4x20 layout and all other parameters) and classify EPSCs into monophasic/multiphasic by the same amplitude-charge logic used in Fig. 2c. Plot mean monophasic versus mean multiphasic amplitude as a function of spacing. If the monophasic mean never exceeds the multiphasic mean for any spacing, the central claim is falsified; if it does, the paper needs to report the crossover spacing and compare it with experimentally measured CaV1.3 nearest-neighbor distances (e.g., dSTORM) to establish plausibility.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At the center of the paper is the claim (Conclusions) that strong clustering of CaV1.3 channels enables multivesicular release to produce large-amplitude monophasic EPSCs. The actual simulation results do not show this. The experimental pattern the paper must explain is stated in Fig. 2c: multiphasic EPSC amplitudes are lower than monophasic EPSC amplitudes. The weak-clustering simulation (Fig. 2 right) shows the opposite ordering, and the authors' caption says so: EPSCs with fewer vesicles are monophasic, those with more vesicles are multiphasic. In the strong-clustering run (Fig. 3), monophasic amplitudes increase, but the authors concede: 'Although the amplitude of monophasic EPSCs does not yet exceed that of multiphasic EPSCs in the current simulation results' (Section II, after Fig. 3d). Therefore the plotted results fail the very comparison that motivated the mechanism. The paper's conclusion rests on the hope that a stronger, un-simulated clustering regime will reverse the ordering. That is an extrapolation, not a demonstration. Moreover, the 'strong' cluster geometry is not parameterized anywhere (no inter-channel spacing or packing numbers are given), so the extrapolated regime is undefined. This is an internal gap: the claimed capability is not exhibited by the simulations as run.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript presents a multiscale computational model of synaptic vesicle fusion at inner hair cell ribbon synapses, coupling CaV1.3 channel gating, buffered Ca2+ diffusion, a five-Ca2+-binding synaptotagmin-like sensor, vesicle fusion pore expansion, glutamate release, and a 16-state AMPA receptor model. The authors compare simulated EPSC amplitude versus charge distributions against experimental data from Chapochnikov et al. (Ref. 11) under two hypothesized spatial organizations of calcium channels: weak and strong clustering. They argue that strong clustering synchronizes multivesicular release and thereby explains large monophasic EPSCs, favoring multivesicular release over univesicular release as the origin of large EPSCs.","tokens_in":10102,"tokens_out":3486,"duration_ms":33775,"significance":"If the central claim could be established, the paper would offer a unified mechanistic account of two puzzling features of inner hair cell EPSCs—the coexistence of monophasic and multiphasic waveforms and the unusually large amplitudes of monophasic events—without invoking oversized vesicles. The authors are to be credited for building an unusually complete forward simulation chain from Ca2+ microdomains to postsynaptic currents and for comparing against quantitative experimental distributions (Fig. 2c). However, the key mechanistic prediction is not actually exhibited by the simulations as run: the strong-clustering results still show monophasic amplitudes below multiphasic ones. The claim therefore currently rests on an extrapolation beyond the parameter range that the authors simulated, which is a load-bearing gap.","major_comments":[{"comment":"The experimental target, as stated in the text, is that multiphasic EPSC amplitudes are lower than monophasic amplitudes. The weak-clustering simulation (Fig. 2, right) shows the opposite ordering, and the strong-clustering run (Fig. 3d) does not reverse it; the authors concede that \"the amplitude of monophasic EPSCs does not yet exceed that of multiphasic EPSCs in the current simulation results.\" The concluding claim that strong clustering can produce large-amplitude monophasic EPSCs is therefore an extrapolation rather than a demonstrated result. Please run simulations in a stronger clustering regime and show that the ordering flips, or soften the conclusion to match the data.","section":"Section II, Figs. 2c and 3d"},{"comment":"The \"strong clustering\" condition is not parameterized. The text states only that weak clustering corresponds to a 25 nm inter-channel distance, but no inter-channel spacing, number of channels per cluster, or cluster geometry is given for the strong-clustering case. Without these values, the reader cannot assess how extreme the hypothesized clustering is, nor reproduce the simulations. Please state the geometric parameters for both conditions.","section":"Section III.2 and Fig. 3a"},{"comment":"The analytic Ca2+ solution introduces an ad hoc correction factor eta = 0.19 sqrt(r/r0 - 1) + 1, inserted to account for the ribbon ceiling. The manuscript says this reduces the difference between analytic and numerical solutions, but no derivation or sensitivity analysis is provided. Because this factor scales the entire Ca2+ transient, it can affect vesicle fusion timing and hence EPSC shape and amplitude. Please either derive eta from the boundary conditions or evaluate the sensitivity of the main results to its functional form.","section":"Section III.2, Eqs. (10)-(11)"}],"minor_comments":[{"comment":"There are several typos: \"potasium\" and \"pottasium\" should be \"potassium,\" and \"subsequant\" should be \"subsequent.\"","section":"Section III, first paragraph"},{"comment":"The resting calcium concentration is listed as 5 x 10^5 mM, which is physically impossible; it should presumably be 5 x 10^-5 mM (50 nM). Also, the units for dCa2+ are written as s/m^2 but should be m^2/s.","section":"Table 2"},{"comment":"The phrase \"strong calcium channeling of the calcium channel clusters\" is awkward and presumably should be \"strong clustering of calcium channels.\"","section":"Abstract"},{"comment":"The caption for Fig. 2c does not clearly identify which panel is experimental and which is simulated; please label the panels explicitly in the figure or caption.","section":"Section II and Fig. 2"},{"comment":"No code or data availability statement is provided. For a simulation paper with many parameters, making the code available would substantially aid reproducibility and trust in the results.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a valuable and unusually complete simulation framework, but its central claim is not yet supported by the plotted results. The authors' own statement that monophasic amplitudes do not yet exceed multiphasic amplitudes in the strong-clustering case is a significant gap. I see no reason to reject outright, because the gap could be closed by additional simulations in a well-parameterized stronger-clustering regime, but the manuscript should not be accepted in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nWhat's actually new here is the end-to-end simulation chain: stochastic CaV1.3 gating, buffered calcium diffusion, sensor binding, vesicle fusion pore expansion, glutamate release, and a 16-state AMPA model, all run together to ask whether clustered channels synchronize multivesicular release and make large monophasic EPSCs. That integration is a genuine contribution, and the authors are transparent about mismatches with experiment. They also disclose that their analytic calcium solution drifts from the numerical one and introduce an ad hoc correction η, which is a patch but at least an acknowledged one.\n\nThe central problem is the one the stress-test identifies. The experimental pattern is that monophasic EPSCs are larger than multiphasic EPSCs. The weak-clustering simulation gives the opposite ordering, and in the strong-clustering run the authors concede that monophasic amplitudes still do not exceed multiphasic amplitudes. So the plotted results fail the very comparison that motivated the mechanism. The conclusion rests on the hope that an even stronger, un-simulated clustering regime will reverse the ordering. That is extrapolation, not demonstration. The strong-cluster geometry is also not parameterized: no channel spacing or packing density is given for the 'strong' condition, so the extrapolated regime is undefined. This is a load-bearing gap, not a cosmetic one.\n\nOther soft spots, in proportion. The clustering hypothesis leans on depletion-interaction theory, but the paper presents no experimental evidence that CaV1.3 channels form tight clusters in vivo; the authors state it as a working hypothesis. The sensor-ion kinetics are inherited from a model fitted to hair-cell exocytosis, so the forward-model claim is weakened by that circularity, though not fatally. No code or data are provided for replication, and the parameter tables contain apparent unit problems, e.g. [Ca2+]r listed as 5×10^5 mM. These are fixable but should not be ignored.\n\nWho is this for? A modeling-oriented synaptic physiologist who wants a scaffold for testing whether channel clustering can drive coherent release. The integrated model is useful even if the headline claim is not yet supported. With the amplitude discrepancy resolved, a defined strong-clustering geometry, and released code, this could become a solid paper. As it stands, I would not cite it for the conclusion, but I would not desk-reject it either.\n\nRecommendation: send to peer review, but tell the authors the central claim needs either a simulated regime that actually produces the experimental ordering or a substantially toned-down conclusion.","headline":"The integrated simulation is a real piece of work, but the paper's central claim about channel clustering relies on simulations that do not yet reproduce the amplitude ordering it sets out to explain.","tokens_in":10701,"tokens_out":1824,"would_cite":false,"duration_ms":18502,"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":"The paper argues that strong clustering of CaV1.3 calcium channels synchronizes multivesicular release, producing the large monophasic EPSCs observed in auditory hair cells.","keywords":["hair cell ribbon synapse","multivesicular release","calcium channel clustering","CaV1.3","EPSC waveforms","synaptic vesicle fusion simulation","cochlear hearing"],"falsifier":"Super-resolution or electron microscopy of inner hair cell active zones showing that CaV1.3 channels are not packed into sub-diffraction clusters would falsify the proposed mechanism; alternatively, an experiment that disrupts the presumed clustering forces without changing the large monophasic EPSC amplitude distribution would contradict the claim.","tokens_in":9637,"feed_emoji":"⚡","tokens_out":5522,"duration_ms":55877,"temperature":0.7,"pith_summary":"This paper argues that the large, single-peaked electrical signals recorded at inner hair cell synapses come from several vesicles fusing at nearly the same time, not from one giant vesicle. The synchronizing agent, the authors propose, is strong local clustering of CaV1.3 calcium channels: a tight cluster admits a concentrated calcium cloud that binds the vesicle sensor and triggers coherent multivesicular release. The claim matters because it would settle a long-standing debate about the origin and shape of large EPSCs, and it would connect the physics of channel-channel interactions to hearing. The full simulation chain, from channel gating to AMPA receptor kinetics, reproduces the observed amplitude-over-charge pattern and shows that stronger clustering raises the amplitude of monophasic EPSCs.","feed_headline":"Clustered calcium channels may explain giant auditory signals","feed_subtitle":"A full simulation links strong CaV1.3 clustering to synchronized multivesicular release, producing large monophasic EPSCs.","key_machinery":"The load-bearing element is the spatial arrangement of CaV1.3 calcium channels, specifically tight sub-clusters within the active zone. The model couples a three-state Markov gating model of the channels, an analytic linearized-buffer solution for calcium diffusion from point sources, a five-calcium-ion sensor that triggers vesicle fusion, a fusion-pore expansion equation, and a 16-state AMPA receptor kinetic scheme. Strong clustering concentrates the calcium nanodomain so that several nearby vesicles reach the five-$\\mathrm{Ca}^{2+}$ threshold within a narrow time window, generating a large, monophasic postsynaptic current.","core_discovery":"The paper's central discovery claim is that multivesicular release, driven by clustered CaV1.3 channels, can account for the large-amplitude monophasic EPSCs observed in inner hair cells, whereas univesicular release cannot plausibly explain the multiphasic EPSC shapes. In the simulations, weak clustering reproduces the experimentally known distributions of EPSC amplitude and charge, while strong clustering, modeled as tightly packed channels, produces higher local calcium concentrations and enhances monophasic EPSC amplitude. The authors conclude that a large EPSC is the summed postsynaptic response of several vesicles released nearly simultaneously, and that strong channel clustering is the mechanism that makes the release coherent. They report one residual discrepancy: in the current simulations multiphasic EPSC amplitudes still exceed monophasic ones, opposite to the experimental ordering, and they frame the result as a trend that stronger clustering would complete.","pith_inferences":["If CaV1.3 clustering is a tunable property, interventions that alter channel mobility or the membrane environment could shift the balance between monophasic and multiphasic EPSCs, effectively changing temporal coding in the auditory nerve.","The same clustering-synchronization logic might apply to other ribbon synapses, such as retinal bipolar cells or vestibular hair cells, where large EPSCs have been reported; a simulation using their channel densities would be a direct test.","A specific prediction is that the amplitude distribution of monophasic EPSCs should correlate with the measured degree of CaV1.3 sub-clustering across developmental stages or experimental manipulations.","The residual sign discrepancy between simulated and experimental amplitude ordering suggests that, in addition to presynaptic clustering, a postsynaptic ingredient such as receptor saturation or desensitization shaped by the glutamate transient may be needed to complete the explanation."],"forward_implications":["If the claim holds, large EPSCs in inner hair cells should be counted as compound events, so quantal analysis of these synapses would need to treat multivesicular sums as the elementary signal.","Monophasic and multiphasic EPSCs become two ends of a synchronization continuum, with the degree of calcium-channel clustering controlling the degree of vesicle synchrony.","The model's parameter sweep shows that fusion of three vesicles yields about 250 pA, matching the experimentally observed average EPSC amplitude and making multivesicular sums quantitatively consistent with measured amplitudes.","Univesicular release would require vesicles larger than those observed and cannot naturally generate multiphasic waveforms, so the simulation shifts the burden of explanation onto multivesicular release.","The apparent contradiction between the roughly 700 vesicles/s sustained release rate and the 200-400 spikes/s afferent firing rates becomes a well-posed question about how few vesicles need to fuse coherently to produce a large EPSC."],"supporting_citations":[{"why":"Defines the two release modes and the monophasic and multiphasic EPSC classification that the paper seeks to explain.","marker":"[10]"},{"why":"Supplies the experimental EPSC amplitude and charge distributions and the dynamic fusion pore (univesicular release) hypothesis that the paper argues against.","marker":"[11]"},{"why":"Provides CaV1.3 channel cluster parameters, endogenous buffer species and concentrations, and the gating-rate formalism used in the simulation.","marker":"[13]"},{"why":"Supplies the depletion-interaction theory that motivates the hypothesis that calcium channels form strong local clusters.","marker":"[16]"},{"why":"Provides the sustained vesicle release rate of about 700 vesicles/s, used for vesicle replenishment and for the rate-based argument in the introduction.","marker":"[20]"},{"why":"Supplies the 200-400 spikes/s afferent fiber firing rates that motivate the question of whether fewer than six vesicles can generate a strong EPSC.","marker":"[21]"},{"why":"Provides the empirical activation time constant formula for CaV1.3 gating, including the resting-voltage value the simulation needs.","marker":"[27]"},{"why":"Supplies the linearized buffered calcium diffusion solution used to model calcium nanodomains near open channels.","marker":"[29]"},{"why":"Provides the five-calcium-ion binding kinetic model for the calcium sensor that triggers vesicle fusion.","marker":"[31]"},{"why":"Supplies the charged-vesicle discharge model used for glutamate and ion exchange through the expanding fusion pore.","marker":"[32]"}],"fun_headline_variants":["Calcium clusters sync vesicle release for giant auditory EPSCs","Strong CaV1.3 clustering yields coherent multivesicular release","Giant EPSCs from calcium channel clustering in hair cells","Multivesicular release explained by strong calcium cluster","Clustered calcium channels sync fusion for large EPSCs"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire mechanism depends on CaV1.3 channels actually forming tight local clusters at inner hair cell active zones, which the paper assumes from depletion-interaction theory rather than demonstrates experimentally; if the channels are spread uniformly, the proposed synchronization does not occur.","fun_headline_variants_meta":{"raw":{"variants":["Calcium clusters sync vesicle release for giant auditory EPSCs","Strong CaV1.3 clustering yields coherent multivesicular release","Giant EPSCs from calcium channel clustering in hair cells","Multivesicular release explained by strong calcium cluster","Clustered calcium channels sync fusion for large EPSCs"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000691,"raw_usage":{"total_tokens":3095,"prompt_tokens":878,"completion_tokens":2217,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":494,"completion_tokens_details":{"reasoning_tokens":2145}},"tokens_in":494,"tokens_out":2217,"duration_ms":14926,"temperature":1.0,"reasoning_tokens":2145,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:20:30.251463+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Super-resolution or electron microscopy of inner hair cell active zones showing that CaV1.3 channels are not packed into sub-diffraction clusters would falsify the proposed mechanism; alternatively, an experiment that disrupts the presumed clustering forces without changing the large monophasic EPSC amplitude distribution would contradict the claim.","supporting_citations":[{"cited_title":"”Two modes of release shape the postsynaptic response at the inner hair cell ribbon synapse.” Journal of Neuroscience 30.12 (2010): 4210-4220","cited_arxiv_id":null,"evidence_quote":"Defines the two release modes and the monophasic and multiphasic EPSC classification that the paper seeks to explain."},{"cited_title":"”Uniquantal release through a dynamic fusion pore is a candidate mechanism of hair cell exocytosis.” Neuron 83.6 (2014): 1389-1403","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental EPSC amplitude and charge distributions and the dynamic fusion pore (univesicular release) hypothesis that the paper argues against."},{"cited_title":"”Developmental refinement of hair cell synapses tightens the coupling of Ca2+ influx to exocytosis.” The EMBO journal 33.3 (2014): 247-264","cited_arxiv_id":null,"evidence_quote":"Provides CaV1.3 channel cluster parameters, endogenous buffer species and concentrations, and the gating-rate formalism used in the simulation."},{"cited_title":"”Polymer brush-induced depletion interactions and clustering of membrane proteins.” The Journal of Chemical Physics 154.21 (2021): 214901","cited_arxiv_id":null,"evidence_quote":"Supplies the depletion-interaction theory that motivates the hypothesis that calcium channels form strong local clusters."},{"cited_title":"”Hearing requires otoferlin-dependent efficient replenishment of synaptic vesicles in hair cells.” Nature neuroscience 13.7 (2010): 869","cited_arxiv_id":null,"evidence_quote":"Provides the sustained vesicle release rate of about 700 vesicles/s, used for vesicle replenishment and for the rate-based argument in the introduction."},{"cited_title":"Charles Liberman","cited_arxiv_id":null,"evidence_quote":"Supplies the 200-400 spikes/s afferent fiber firing rates that motivate the question of whether fewer than six vesicles can generate a strong EPSC."},{"cited_title":"”Biophysical properties of CaV1","cited_arxiv_id":null,"evidence_quote":"Provides the empirical activation time constant formula for CaV1.3 gating, including the resting-voltage value the simulation needs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the linearized buffered calcium diffusion solution used to model calcium nanodomains near open channels."},{"cited_title":"”Calcium dependence of exocytosis and endocytosis at the cochlear inner hair cell afferent synapse.” Neuron 29.3 (2001): 681-690","cited_arxiv_id":null,"evidence_quote":"Provides the five-calcium-ion binding kinetic model for the calcium sensor that triggers vesicle fusion."},{"cited_title":"”A mechanism for discharge of charged excitatory neurotransmitter.” Bio- 11 physical journal 72.2 (1997): 507-521","cited_arxiv_id":null,"evidence_quote":"Supplies the charged-vesicle discharge model used for glutamate and ion exchange through the expanding fusion pore."}],"review_version":1}