Current-Voltage Characteristics of K+ Channels Estimated by MD Simulations and Markov State Models.
- 1Department of Electrical, Electronic, and Information Engineering "Guglielmo Marconi", Alma Mater Studiorum─University of Bologna, via dell'Università 50, Cesena (FC)47521, Italy.
- 2Department of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, Perugia06123, Italy.
Abstract
Molecular dynamics (MD) simulations have yielded important insights into ion conduction in potassium channels, but quantitative comparison with electrophysiological experiments remains challenging. Due to their high computational cost, MD simulations are typically performed at membrane potentials well above physiological values, and at only a limited number of voltages. Since current-voltage relationships are not necessarily linear, this limits direct comparison between simulations and experiments. Here, we introduce a method to estimate the current-voltage characteristics of ion channels from Markov state models (MSMs) constructed from MD simulations performed at only few membrane potentials. Time-discrete MSMs of ion conduction are converted into continuous-time rate matrices, whose voltage dependence is modeled using a rate theory formulation with free energy barriers depending on membrane potential. This approach enables the prediction of channel currents over a wide voltage range without additional simulations. We validated the method using MD simulations of the potassium channels KcsA and MthK. In both cases, the currents predicted at low membrane potentials are in good agreement with those obtained directly from MD simulations, demonstrating the robustness and efficiency of the approach.
Related
- Periodic paralysis and voltage-gated ion channels
- The extracellular turret as a structural hub for voltage sensor-pore coupling in HCN channels
- Voltage versus concentration in ion channels: Physical signatures reveal rate-limiting steps
- Revealing intermediate states of voltage-gated potassium channels with physics-informed deep learning: application to Kv7.1
- Environmental monitoring of static magnetic fields around operating high-voltage direct-current transmission facilities in Japan
BibTeX
@article{furini2026currentvoltage,
title = {{Current-Voltage Characteristics of K+ Channels Estimated by MD Simulations and Markov State Models.}},
author = {Simone Furini and Luigi Catacuzzeno},
journal = {Journal of chemical information and modeling},
year = {2026},
month = sep,
doi = {10.1021/acs.jcim.6c01258},
url = {https://europepmc.org/article/MED/42906029},
note = {CC BY},
}
RIS
TY - JOUR
TI - Current-Voltage Characteristics of K+ Channels Estimated by MD Simulations and Markov State Models.
AU - Simone Furini
AU - Luigi Catacuzzeno
PY - 2026
DA - 2026/09/01
JO - Journal of chemical information and modeling
DO - 10.1021/acs.jcim.6c01258
UR - https://europepmc.org/article/MED/42906029
AB - Molecular dynamics (MD) simulations have yielded important insights into ion conduction in potassium channels, but quantitative comparison with electrophysiological experiments remains challenging. Due to their high computational cost, MD simulations are typically performed at membrane potentials well above physiological values, and at only a limited number of voltages. Since current-voltage relationships are not necessarily linear, this limits direct comparison between simulations and experiments. Here, we introduce a method to estimate the current-voltage characteristics of ion channels from Markov state models (MSMs) constructed from MD simulations performed at only few membrane potentials. Time-discrete MSMs of ion conduction are converted into continuous-time rate matrices, whose voltage dependence is modeled using a rate theory formulation with free energy barriers depending on membrane potential. This approach enables the prediction of channel currents over a wide voltage range without additional simulations. We validated the method using MD simulations of the potassium channels KcsA and MthK. In both cases, the currents predicted at low membrane potentials are in good agreement with those obtained directly from MD simulations, demonstrating the robustness and efficiency of the approach.
SN - 1549-9596
SN - 1549-960X
ER -