A strategy for determining the equilibrium constants for heteromeric ion channels in a complex model

Abstract

Ligand-gated ion channels are oligomers containing several binding sites for the ligands. However, the signal transmission from the ligand binding site to the pore has not yet been fully elucidated for any of these channels. In heteromeric channels, the situation is even more complex than in homomeric channels. Using published data for concatamers of heteromeric cyclic nucleotide–gated channels, we show that, on theoretical grounds, multiple functional parameters of the individual subunits can be determined with high precision. The main components of our strategy are (1) the generation of a defined subunit composition by concatenating multiple subunits, (2) the construction of 16 concatameric channels, which differ in systematically permutated binding sites, (3) the determination of respectively differing concentration–activation relationships, and (4) a complex global fit analysis with corresponding intimately coupled Markovian state models. The amount of constraints in this approach is exceedingly high. Furthermore, we propose a stochastic fit analysis with a scaled unitary start vector of identical elements to avoid any bias arising from a specific start vector. Our approach enabled us to determine 23 free parameters, including 4 equilibrium constants for the closed–open isomerizations, 4 disabling factors for the mutations of the different subunits, and 15 virtual equilibrium-association constants in the context of a 4-D hypercube. From the virtual equilibrium-association constants, we could determine 32 equilibrium-association constants of the subunits at different degrees of ligand binding. Our strategy can be generalized and is therefore adaptable to other ion channels.