Pathway and endpoint free energy calculations for cyclic nucleotide binding to HCN channels

cAMP and cGMP differentially bind to and regulate a variety of proteins, including cyclic nucleotide-gated (CNG) channels and hyperpolarization-activated cyclic nucleotide-regulated (HCN) channels. Previous site-directed mutagenesis studies have isolated two conserved residues that are critical for enabling certain channels to selectively bind cGMP relative to cAMP. However, no definitive mechanism has been identified that explains the preferential activation of other channels by cAMP. Here we apply computational binding free energy methods, including thermodynamic integration, linear interaction energy, and continuum electrostatic calculations, to gain insights into the mechanisms of cyclic nucleotide selectivity. Consistent with experimental observations, computational results for the cAMP-selective HCN channels show that the binding free energy of cAMP is lower (more favorable) than that of cGMP. Surprisingly, cAMP selectivity is not due to its preferential contacts with protein, but rather reflects the greater hydration energy of cGMP relative to cAMP, resulting in a greater energetic cost for cGMP binding.