Computational Framework for Identifying Ion Channel Mutation-Compensating Interventions
- 1Tufts University
Abstract
We present an automated pipeline for finding therapeutic interventions in channelopathies. It starts from patch clamp recordings of the variant channel, searches over the pharmacologically accessible conductances, and says which currents must change and by how much. The intervention never touches the mutated channel: it compensates by modulating others. The feasibility of this approach is demonstrated through two computational strategies, each addressing a specific capability gap. First, we employ high-fidelity NEURON simulations combined with exploratory search algorithms to model two clinically described variants, and identify "stability windows" in which channel modulations restore the reference firing pattern. The compensating configurations are degenerate: 37 conductance triples fire the same number of action potentials at every one of the 35 injected-current levels, and the 10,780 triples reaching the best intervention value produce only six distinct firing patterns. Second, to address the computational cost of such detailed modeling, we implement a differentiable Hodgkin-Huxley model in PyTorch. Applied here to theoretical mutation screening, this approach trades granular detail for speed, treating the search as a continuous optimization problem that makes very large, non-local searches of the conductance parameter space affordable: the differentiable forward model evaluates 34,000-64,000 candidate parameter sets per second on a single consumer graphics card, against 11.9-16.9 per second for the non-differentiable pipeline on a 72-core node. The gain arises from forward-model throughput feeding a search rather than from the gradient itself, and it is what makes the 32,000,000-point survey of the solution topology reported here affordable. These computational predictions can inform practical drug development and high-throughput screening by naming the currents that must change and by how much.
Related
- Beyond individual features: Identifying computational principles in perceiving visual beauty
- Optimizing peer-led health interventions: A social network analysis approach for identifying influential fishermen in Malawi
- VARION: A Network Propagation Framework for Individual Patient Somatic Mutation Interpretation in Cancer Molecular Subtyping
- Hemodynamic-Structural Multi-Fidelity Assessment of Endovascular Flow Diversion Strategies for Cerebral Aneurysm Treatment: An Integrated Computational and Experimental Framework
- A robust mutation spectrum analysis framework for high-depth error-corrected sequencing data
BibTeX
@misc{hazan2026computational,
title = {{Computational Framework for Identifying Ion Channel Mutation-Compensating Interventions}},
author = {H. Hazan and M. Levin},
howpublished = {bioRxiv},
year = {2026},
month = sep,
doi = {10.64898/2026.09.20.753002},
url = {https://www.biorxiv.org/content/10.64898/2026.09.20.753002v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Computational Framework for Identifying Ion Channel Mutation-Compensating Interventions
AU - H. Hazan
AU - M. Levin
PY - 2026
DA - 2026/09/20
PB - bioRxiv
DO - 10.64898/2026.09.20.753002
UR - https://www.biorxiv.org/content/10.64898/2026.09.20.753002v1
AB - We present an automated pipeline for finding therapeutic interventions in channelopathies. It starts from patch clamp recordings of the variant channel, searches over the pharmacologically accessible conductances, and says which currents must change and by how much. The intervention never touches the mutated channel: it compensates by modulating others. The feasibility of this approach is demonstrated through two computational strategies, each addressing a specific capability gap. First, we employ high-fidelity NEURON simulations combined with exploratory search algorithms to model two clinically described variants, and identify "stability windows" in which channel modulations restore the reference firing pattern. The compensating configurations are degenerate: 37 conductance triples fire the same number of action potentials at every one of the 35 injected-current levels, and the 10,780 triples reaching the best intervention value produce only six distinct firing patterns. Second, to address the computational cost of such detailed modeling, we implement a differentiable Hodgkin-Huxley model in PyTorch. Applied here to theoretical mutation screening, this approach trades granular detail for speed, treating the search as a continuous optimization problem that makes very large, non-local searches of the conductance parameter space affordable: the differentiable forward model evaluates 34,000-64,000 candidate parameter sets per second on a single consumer graphics card, against 11.9-16.9 per second for the non-differentiable pipeline on a 72-core node. The gain arises from forward-model throughput feeding a search rather than from the gradient itself, and it is what makes the 32,000,000-point survey of the solution topology reported here affordable. These computational predictions can inform practical drug development and high-throughput screening by naming the currents that must change and by how much.
ER -