Abstract
Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are opened in an allosteric manner by membrane hyperpolarization and cyclic nucleotides such as cAMP. Because of conflicting reports from experimental studies on whether cAMP binding to the four available binding sites in the channel tetramer operates cooperatively in gating, we employ here a computational approach as a promising route to examine ligand-induced conformational changes after binding to individual sites. By combining an elastic network model (ENM) with linear response theory (LRT) for modeling the apo-holo transition of the cyclic nucleotide-binding domain (CNBD) in HCN channels, we observe a distinct pattern of cooperativity matching the "positive-negative-positive" cooperativity reported from functional studies. This cooperativity pattern is highly conserved among HCN subtypes (HCN4, HCN1), but only to a lesser extent visible in structurally related channels, which are only gated by voltage (KAT1) or cyclic nucleotides (TAX4). This suggests an inherent cooperativity between subunits in HCN channels as part of a ligand-triggered gating mechanism in these channels.