# A synthetic peptide that prevents cAMP regulation in mammalian hyperpolarization-activated cyclic nucleotide-gated (HCN) channels

A. Saponaro, F. Cantini, A. Porro, A. Bucchi, D. Difrancesco, V. Maione, C. Donadoni, B. Introini, Pietro Mesirca, Matteo E. Mangoni, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

A. Saponaro, F. Cantini, A. Porro, A. Bucchi, D. Difrancesco, et al. A synthetic peptide that prevents cAMP regulation in mammalian hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. eLife, 2018, 7, 10.7554/eLife.35753. hal-02073210

HAL Id: hal-02073210

https://hal.umontpellier.fr/hal-02073210

Submitted on 7 Jun 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# A synthetic peptide that prevents cAMP regulation in mammalian hyperpolarization-activated cyclic nucleotide-gated (HCN) channels

Andrea Saponaro¹, Francesca Cantini²,³, Alessandro Porro¹, Annalisa Bucchi¹, Dario DiFrancesco¹, Vincenzo Maione⁴, Chiara Donadoni¹, Bianca Introini¹, Pietro Mesirca⁵,⁶, Matteo E Mangoni⁵,⁶, Gerhard Thiel⁷, Lucia Banci²,³,⁴,⁸, Bina Santoro⁹, Anna Moroni¹,¹⁰*

¹Department of Biosciences, University of Milan, Milan, Italy; ²Department of Chemistry, University of Florence, Florence, Italy; ³Magnetic Resonance Center, University of Florence, Florence, Italy; ⁴Interuniversity Consortium for Magnetic Resonance of Metalloproteins, Sesto Fiorentino, Italy; ⁵Institut de Génomique Fonctionnelle, CNRS, INSERM F-34094, Université de Montpellier, Montpellier, France; ⁶Laboratory of Excellence Ion Channels Science and Therapeutics, Valbonne, France; ⁷Department of Biology, TU-Darmstadt, Darmstadt, Germany; ⁸Institute of Neurosciences, Consiglio Nazionale delle Ricerche, Florence, Italy; ⁹Department of Neuroscience, Columbia University, New York, United States; ¹⁰Institute of Biophysics, Consiglio Nazionale delle Ricerche, Milan, Italy

**Abstract** Binding of TRIP8b to the cyclic nucleotide binding domain (CNBD) of mammalian hyperpolarization-activated cyclic nucleotide-gated (HCN) channels prevents their regulation by cAMP. Since TRIP8b is expressed exclusively in the brain, we envisage that it can be used for orthogonal control of HCN channels beyond the central nervous system. To this end, we have identified by rational design a 40-aa long peptide (TRIP8b_nano) that recapitulates affinity and gating effects of TRIP8b in HCN isoforms (hHCN1, mHCN2, rbHCN4) and in the cardiac current I<sub>f</sub> in rabbit and mouse sinoatrial node cardiomyocytes. Guided by an NMR-derived structural model that identifies the key molecular interactions between TRIP8b_nano and the HCN CNBD, we further designed a cell-penetrating peptide (TAT-TRIP8b_nano) which successfully prevented β-adrenergic activation of mouse I<sub>f</sub> leaving the stimulation of the L-type calcium current (I<sub>CaL</sub>) unaffected. TRIP8b_nano represents a novel approach to selectively control HCN activation, which yields the promise of a more targeted pharmacology compared to pore blockers.

DOI: https://doi.org/10.7554/eLife.35753.001

## Introduction

Hyperpolarization-activated cyclic nucleotide-gated (HCN1-4) channels are the molecular correlate of the I<sub>f</sub>/I<sub>h</sub> current, which plays a key role in controlling several higher order physiological functions, including dendritic integration and intrinsic rhythmicity both in cardiac and neuronal cells (Robinson and Siegelbaum, 2003). Unique among the voltage-gated ion channel superfamily, HCN channels are modulated by the direct binding of cAMP to their cyclic nucleotide binding domain (CNBD). Binding of the cyclic nucleotide increases the channel open probability upon

<footnote>*For correspondence: anna.moroni@unimi.it</footnote>

<footnote>**Competing interests:** The authors declare that no competing interests exist.</footnote>

<footnote>**Funding:** See page 19</footnote>

<footnote>**Received:** 07 February 2018</footnote>

<footnote>**Accepted:** 14 June 2018</footnote>

<footnote>**Published:** 20 June 2018</footnote>

<footnote>**Reviewing editor:** Volker Dötsch, J.W. Goethe-University, Germany</footnote>

<footnote>© Copyright Saponaro et al. This article is distributed under the terms of the [Creative Commons Attribution License](https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
hyperpolarization via conformational changes in the CNBD that are propagated to the pore through the C-linker domain (DiFrancesco and Tortora, 1991; Wainger et al., 2001; Zagotta et al., 2003).

In addition to cAMP, HCN channels are regulated by TRIP8b, a brain-specific auxiliary (β) subunit, which modulates two independent features of the channel, namely trafficking and gating (Santoro et al., 2009; Zolles et al., 2009). For this dual regulation, TRIP8b binds HCN channels through two distinct sites: via the tetratricopeptide repeat (TPR) domain, which interacts with the last three amino acids (SNL) of HCN channels and regulates their trafficking; and via the TRIP8b<sub>core</sub> domain, which interacts with the CNBD and antagonizes the effect of cAMP on the voltage dependency of the channel (Santoro et al., 2011; Han et al., 2011; Hu et al., 2013).

Here, we focus our attention on the specific action of TRIP8b in preventing cAMP regulation of HCN channels. Given the brain-specific localization of TRIP8b, we posit that a TRIP8b-derived peptide drug, able to reproduce the effect of the full length protein on HCN channel gating, can be developed for orthogonal selective regulation of HCN in cells/tissues in which TRIP8b is not expressed. cAMP-dependent modulation of HCN channels underlies distinct roles of cAMP in heart rate regulation (DiFrancesco, 1993) and development of peripheral neuropathic pain (Emery et al., 2012; Herrmann et al., 2017), which can be dissected by using a TRIP8b-based tool. In this regard, peptide-based drugs (2–50 aa long) are emerging as a fascinating application area as they open new therapeutic possibilities with an advantage over small molecules in terms of specificity and affinity for the target (Fosgerau and Hoffmann, 2015; Henninot et al., 2018). To this end, we searched for the minimal peptide that binds to the CNBD and recapitulates the gating effect of full length TRIP8b in three HCN isoforms (HCN1, HCN2 and HCN4) and in the native $I_f$ current. In previous studies, we identified the core portion of TRIP8b (TRIP8b<sub>core</sub>, 80 aa long) that interacts with the HCN CNBD and prevents cAMP modulation in full length channels (Santoro et al., 2011; Hu et al., 2013; Saponaro et al., 2014). A recent paper (Lyman et al., 2017) reported an even shorter binding sequence of TRIP8b (37 aa). However, this peptide, which was identified by progressive truncation of TRIP8b<sub>core</sub>, failed to reproduce the binding affinity of the starting construct. Moreover, evidence for activity of this peptide on HCN currents is lacking. In the present study, we adopted a structure-driven rational design approach to engineer a 40-aa long peptide, TRIP8b<sub>nano</sub>, that efficiently prevents cAMP regulation of HCN channels. The rational design of this peptide, based on secondary structure predictions and on NMR data of TRIP8b<sub>core</sub>, was supported by an NMR-based 3D model structure of the complex formed by the TRIP8b<sub>nano</sub> peptide and CNBD of the human HCN2 channel isoform. This structural information identifies crucial interactions between the two partners and explains both direct (Han et al., 2011; DeBerg et al., 2015; Bankston et al., 2017) and indirect (allosteric) (Hu et al., 2013; Saponaro et al., 2014) modes of competition between TRIP8b and cAMP for binding to the CNBD. The evidence that TRIP8b<sub>nano</sub> establishes all relevant interactions with the CNBD is reflected by the finding that, contrary to shorter core sequences (Lyman et al., 2017), it binds to the isolated CNBD with the same affinity as TRIP8b<sub>core</sub> and acts with even higher efficacy than TRIP8b<sub>core</sub> in preventing cAMP modulation of full length HCN channels (Hu et al., 2013). In pacemaker myocytes of the sino-atrial node (SAN), TRIP8b<sub>nano</sub> equally prevented cAMP stimulation of native f-channels leading to a 30% reduction in spontaneous firing rate.

To develop TRIP8b<sub>nano</sub> as a membrane permeable drug, we linked it with the positively charged TAT sequence (Herce et al., 2014). TAT-TRIP8b<sub>nano</sub> was tested in SAN pacemaker myocytes where its addition to the extracellular buffer prevented adrenergic stimulation of the $I_f$ current leaving the activation of the L-type calcium current ($I_{CaL}$) unaffected. Our study opens the possibility of selective in vivo control of the cAMP-dependent facilitation of HCN channel opening, by local supply of TAT-TRIP8b<sub>nano</sub> peptide.

## Results

We have previously shown that TRIP8b<sub>core</sub> (residues 223–303 of mouse TRIP8b splice variant 1a4, hereafter TRIP8b) interacts with two elements of the isolated CNBD protein fragment from HCN channels (residues 521–672 of human HCN2, hereafter CNBD): the C-helix and the N-bundle loop, a sequence connecting helix E' of the C-linker with helix A of the CNBD (Saponaro et al., 2014). Biochemical assays confirmed that each of these two elements, that is, the N-bundle loop and C-helix, is necessary but not sufficient for binding (Saponaro et al., 2014).
To understand the interaction in atomic detail, we used solution NMR spectroscopy to characterize the structural properties of the CNBD - TRIP8b<sub>core</sub> complex. However, the NMR spectra of TRIP8b<sub>core</sub> showed very few signals. In order to improve the quality of the NMR spectra, we reduced the length of the TRIP8b fragment by progressively removing residues at the N- and C-termini with no predicted secondary structure. The truncated peptides were then tested for CNBD binding activity by isothermal titration calorimetry (ITC). We thus identified a 40-aa peptide (TRIP8b<sub>nano</sub>, comprising residues 235-275 of TRIP8b, Figure 1A) with a binding K<sub>D</sub> of 1.5 ± 0.1 µM, a value similar to the K<sub>D</sub> of 1.2 ± 0.1 µM obtained with TRIP8b<sub>core</sub> (Figure 1B). TRIP8b<sub>nano</sub> was therefore employed for all subsequent NMR experiments, resulting in a remarkable improvement in the spectral quality and sample stability.

## Structural characterization of TRIP8b<sub>nano</sub> bound to CNBD

The comparison of the ¹H-¹⁵N HSQC spectra of TRIP8b<sub>nano</sub> with and without CNBD bound shows that the peptide folds upon interaction with the CNBD. Thus, the ¹H-¹⁵N HSQC spectrum of TRIP8b-nano without CNBD shows a limited ¹H resonance dispersion, characteristic of unstructured proteins (Dyson and Wright, 2004), while a larger number of well-dispersed amide signals appear in the spectrum of the CNBD-bound form (Figure 1C). Importantly, we were now able to assign the backbone chemical shift resonances of TRIP8b<sub>nano</sub> bound to the CNBD. The φ and ψ dihedral angles obtained from the NMR assignment indicate that the peptide displays two α-helices (stretch L<sub>238</sub>-E<sub>250</sub> named helix N and stretch T<sub>253</sub>-R<sub>269</sub> named helix C) when bound to CNBD. The helices are separated by two amino acids; three and six residues at the N- and C- termini, respectively, are unstructured (Figure 1D).

## Structural characterization of CNBD bound to TRIP8b<sub>nano</sub>

NMR-analysis of the CNBD fragment bound to TRIP8b<sub>nano</sub> revealed that the interaction with the peptide does not affect the overall fold of the protein. Thus, the CNBD adopts the typical fold of the cAMP-free state, in line with previous evidence that this is the form bound by TRIP8b (Saponaro et al., 2014; DeBerg et al., 2015). More specifically, the secondary structure elements of the cAMP-free CNBD are all conserved in the TRIP8b<sub>nano</sub>-bound CNBD (Figure 2). This finding generally agrees with a previous double electron-electron resonance (DEER) analysis of the CNBD - TRIP8b interaction, which showed that TRIP8b binds to a conformation largely similar to the cAMP-free state (DeBerg et al., 2015). Despite the overall agreement with the DEER study, the NMR data also reveal a new and unexpected feature of TRIP8b binding to the CNBD. Indeed, our results show that TRIP8b<sub>nano</sub> binding to the CNBD induces a well-defined secondary structure of the distal region of the C-helix (Figure 2). This means that the distal region of the C-helix (residues 657-662), which is unstructured in the free form of the CNBD (Saponaro et al., 2014; Lee and MacKinnon, 2017), extends into a helical structure upon ligand binding irrespectively of whether the ligand is cAMP (Puljung and Zagotta, 2013; Saponaro et al., 2014; Lee and MacKinnon, 2017) or TRIP8b (Figure 2). In contrast, and very differently from cAMP, which directly contacts the P-helix in the Phosphate Binding cassette (PBC) and causes its folding (Saponaro et al., 2014; Lee and MacKinnon, 2017), the NMR data show that TRIP8b<sub>nano</sub> binding to the CNBD does not induce P-helix formation (Figure 2).

## Modeling the CNBD - TRIP8b<sub>nano</sub> complex

Despite the significant improvement in sample stability and NMR spectra quality achieved upon TRIP8b<sub>nano</sub> binding, we were still unable to assign the side chains of both proteins in the complex and thus could not solve the solution structure of the complex by the canonical NMR procedure. We therefore built a model (Figure 3) of the CNBD - TRIP8b<sub>nano</sub> complex by docking the two NMR-derived structures described above using the Haddock program (a detailed description of how the respective structures were generated is provided in Materials and methods and Figure 3—source data 1).

In order to define the active residues (ambiguous interaction restraints) on the CNBD we used the chemical shift perturbation values as described in Figure 3—figure supplement 1. For TRIP8b-nano, we defined as active a stretch of residues, E<sub>239</sub>-E<sub>243</sub>, previously identified as critical for the interaction (Santoro et al., 2011). Output clusters of this first molecular docking calculation (settings can
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Functional and structural characterization of TRIP8b<sub>nano</sub>. (A) Primary sequence of TRIP8b<sub>nano</sub>. Amino acid numbering refers to full length mouse TRIP8b (1a4). (B) Binding of TRIP8b<sub>core</sub> and TRIP8b<sub>nano</sub> to purified His<sub>6</sub>-MBP-CNBD measured by Isothermal titration calorimetry (ITC). Upper panel, heat changes (μcal/sec) during successive injections of 8 μL of the corresponding TRIP8b peptide (200 μM) into the chamber containing His<sub>6</sub>-MBP-CNBD (20 Figure 1 continued on next page</figcaption></figure>
Figure 1 continued

µM). Lower panel, binding curve obtained from data displayed in the upper panel. The peaks were integrated, normalized to TRIP8b peptide concentration, and plotted against the molar ratio (TRIP8b peptide/His<sub>6</sub>-MBP-CNBD). Solid line represents a nonlinear least-squares fit to a single-site binding model, yielding, in the present examples, a K<sub>D</sub> = 1.2 ± 0.1 µM for TRIP8b<sub>core</sub> and K<sub>D</sub> = 1.4 ± 0.1 µM for TRIP8b<sub>nano</sub>. (C) Evidence for TRIP8b<sub>nano</sub> folding upon CNBD binding based on the superimposition of the [¹H, ¹⁵N] heteronuclear single quantum coherence (HSQC) NMR spectrum of CNBD-free TRIP8b<sub>nano</sub> (black) and CNBD-bound TRIP8b<sub>nano</sub> (red). The latter experiment was performed at the molar ratio ([CNBD]/[TRIP8b<sub>nano</sub>])=3. The backbone amide (HN) signals of the residues of CNBD-bound TRIP8b<sub>nano</sub> are labelled in red. (D) (Top) Ribbon representation of the 10 lowest energy conformers of TRIP8b<sub>nano</sub> bound to CNBD used for in silico modelling of CNBD-TRIP8b<sub>nano</sub> complex. The unfolded regions at the N- and C-termini of the construct (residues 235–237 and 270–275) are omitted for clarity. (Bottom) Chemical Shift Index (CSI, calculated using TALOS+) plotted as a function of the residue number of TRIP8b<sub>nano</sub> bound to CNBD. Positive values represent helical propensity.

DOI: https://doi.org/10.7554/eLife.35753.002

be found in Materials and methods) were further screened for TRIP8b_nano orientations in agreement with a previous DEER analysis, which placed TRIP8b residue A<sub>248</sub> closer to the proximal portion and TRIP8b residue A<sub>261</sub> closer to the distal portion of the CNBD C-helix (DeBerg et al., 2015). Remarkably, in all clusters selected in this way, residues E<sub>264</sub> or E<sub>265</sub> in TRIP8b were found to interact with residues K<sub>665</sub> or K<sub>666</sub> of the CNBD (Figure 3—figure supplement 2). This finding was notable, because we previously identified K<sub>665</sub>/K<sub>666</sub> as being critical for TRIP8b interaction in a biochemical binding assay (Saponaro et al., 2014). We thus proceeded to individually mutate each of these four positions, and test their effect on binding affinity through ITC. As expected, reverse charge mutations K<sub>665</sub>E or K<sub>666</sub>E (CNBD) as well as E<sub>264</sub>K or E<sub>265</sub>K (TRIP8b_nano) each strongly reduced the CNBD/TRIP8b_nano binding affinity (Figure 3—figure supplement 3).

Based on these observations, we performed a second molecular docking calculation, including E<sub>264</sub> and E<sub>265</sub> as additional active residues for TRIP8b<sub>nano</sub>. This procedure resulted in the model shown in Figure 3, which represents the top-ranking cluster for energetic and scoring function (Figure 3—source data 2) and was fully validated by mutagenesis analysis as described below. Scrutiny of the model shows that TRIP8b<sub>nano</sub> binds to both the C-helix and the N-bundle loop (Figure 3A). Binding to the C-helix is mainly guided by electrostatic interactions between the negative charges on TRIP8b<sub>nano</sub>, and the positive charges on the CNBD (Figure 3A). As shown in Figure 3B, the model highlights a double saline bridge (K<sub>665</sub> and K<sub>666</sub> of CNBD with E<sub>265</sub> and E<sub>264</sub> of TRIP8b<sub>nano</sub>) in line with the ITC results described above (Figure 3—figure supplement 3). Of note, the contribution of residue R<sub>662</sub> to the binding is also consistent with previous experiments showing residual TRIP8b interaction in a CNBD deletion mutant ending at position 663 (Saponaro et al., 2014). Our modeling data suggest that, upon folding of the distal portion of the C-helix, the side chains of residues R<sub>662</sub> and K<sub>665</sub> face to the inside when contacting cAMP, but face to the outside when binding TRIP8b (Figure 3—figure supplement 4). This indicates that cAMP and TRIP8b directly compete for the binding to the distal region of C-helix.

In addition to clarifying the role of residues in the distal portion of the CNBD C-helix, the model also highlights a second important cluster of electrostatic interactions, with R<sub>650</sub> in the proximal portion of the CNBD C-helix contacting E<sub>240</sub> and E<sub>241</sub> in helix N of TRIP8b<sub>nano</sub> (Figure 3C). To confirm the contribution of these residues, we reversed charges and tested each residue mutation for binding in ITC. The results in Figure 3—figure supplement 3 show that R<sub>650</sub>E caused a more than sixfold reduction in binding affinity for TRIP8b<sub>nano</sub>, with smaller but significant effects seen also for E<sub>240</sub>R and E<sub>241</sub>R.

A third important contact highlighted by the model is the interaction between N<sub>547</sub> in the N-bundle loop of the CNBD and D<sub>252</sub> in the link between helix N and helix C of TRIP8b<sub>nano</sub> (Figure 3D). We tested this potential interaction by disrupting the expected hydrogen bond between N<sub>547</sub> and the carboxyl group of the negative residue (D<sub>252</sub>) in TRIP8b<sub>nano</sub>. The asparagine in CNBD was mutated into aspartate (N<sub>547</sub>D) to generate an electrostatic repulsion for D<sub>252</sub>, and the carboxyl group in D<sub>252</sub> of TRIP8b<sub>nano</sub> was removed by mutation into asparagine (D<sub>252</sub>N). As predicted, N<sub>547</sub>D greatly reduced binding to TRIP8b in ITC assays (Figure 3—figure supplement 3), with a smaller but significant effect observed also for D<sub>252</sub>N (Figure 3—figure supplement 3). These results
<figure><img src="image_5.png" /><figcaption>**Figure 2.** NMR structure of CNBD bound to TRIP8b_nano. (A) (Top) comparison of secondary structure elements of cAMP-free CNBD (*Saponaro et al., 2014*), cAMP-bound CNBD (*Zagotta et al., 2003*) and cAMP-free CNBD bound to TRIP8b_nano (this study). Secondary structure elements are indicated by arrows (β-strands) and cylinders (α-helices) and labeled. The loop between β<sub>6</sub> and β<sub>7</sub> constitutes the Phosphate Binding Cassette (PBC). The elements that fold upon binding of cAMP and TRIP8b_nano are shown in red. (Bottom) Chemical Shift Index (CSI, calculated using TALOS+) plotted as a function of the residue number of CNBD bound to TRIP8b_nano. Positive values represent helical propensity, while negative values represent strands. (B) Ribbon representation of the 10 lowest energy conformers of CNBD bound to TRIP8b_nano used for in silico modeling of CNBD-TRIP8b_nano complex. Secondary structure elements are coloured in light gray and labeled. Loop regions are colored in dark gray. The distal region of the C-helix (residues 657–662), which is unfolded in the free form of the CNBD (*Saponaro et al., 2014*) and folds upon TRIP8b_nano binding, is coloured in red. The unfolded regions at the N- and C-termini of the construct (residues 521–532 and 663–672 respectively) are omitted for clarity.</figcaption></figure>

DOI: https://doi.org/10.7554/eLife.35753.003
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Structural model of CNBD – TRIP8b_nano complex. (A) Ribbon representation of the complex where CNBD is in gray and TRIP8b_nano is in orange. Helix N (N) and helix C (C) of TRIP8b_nano are labeled. C-helix of CNBD (αC) is labeled, while N-bundle loop is colored in yellow. Positively charged residues of C-helix CNBD (blue) and negatively charged residues of TRIP8b_nano (red) involved in salt bridges are shown as sticks and labeled. N547 of the N-bundle loop (yellow) and D252 of TRIP8b_nano (orange) are shown as sticks and labeled. (B) Close view of K665 and K666 of CNBD that interact, respectively, with E265 and E264 of TRIP8b_nano. C-helix (αC) of CNBD, and Helix C (C) of TRIP8b_nano are labeled. (C) Close view of R650 of CNBD that is positioned between E240 and E241 of TRIP8b_nano. C-helix (αC) of CNBD, and Helix N (N) of TRIP8b_nano are labeled. (D) Close view of N547 of N-bundle loop that forms a hydrogen bond (red dashed line) with D252 of TRIP8b_nano. Helix E' (αE') and C-helix (αC) of CNBD, and Helix N (N) of TRIP8b_nano are labeled.</figcaption></figure>

DOI: https://doi.org/10.7554/eLife.35753.004

The following source data and figure supplements are available for figure 3:

**Source data 1.** Acquisition parameters for NMR experiments performed on cAMP-free human HCN2 CNBD in complex with TRIP8b_nano and vice-versa.

DOI: https://doi.org/10.7554/eLife.35753.011

**Source data 2.** Docking calculation.

DOI: https://doi.org/10.7554/eLife.35753.012

**Figure supplement 1.** CNBD residues involved in TRIP8b_nano binding.

DOI: https://doi.org/10.7554/eLife.35753.005

**Figure supplement 2.** Representative families of clusters obtained from the first docking calculation.

Figure 3 continued on next page
Figure 3 continued

DOI: https://doi.org/10.7554/eLife.35753.006

**Figure supplement 3.** Biochemical validation of CNBD – TRIP8b<sub>nano</sub> complex.

DOI: https://doi.org/10.7554/eLife.35753.007

**Figure supplement 4.** Different orientation of R<sub>662</sub> and K<sub>665</sub> in the cAMP-bound and TRIP8b<sub>nano</sub>-bound conformation of the CNBD.

DOI: https://doi.org/10.7554/eLife.35753.008

**Figure supplement 5.** Mutation N<sub>520</sub>D affects cAMP affinity in full-length HCN2 channel.

DOI: https://doi.org/10.7554/eLife.35753.009

**Figure supplement 6.** Structural characterization of N<sub>547</sub>D CNBD protein.

DOI: https://doi.org/10.7554/eLife.35753.010

confirm and extend our previous finding that the N-bundle loop contributes in a substantial manner to the binding of TRIP8b (Saponaro et al., 2014). To understand the mechanism for the allosteric effect of TRIP8b on cAMP binding, which has been postulated on the basis of electrophysiological and structural data (Hu et al., 2013; Saponaro et al., 2014), we further tested by ITC the affinity of the N<sub>547</sub>D CNBD mutant for cAMP. Somewhat surprisingly, the affinity of the mutant is much lower than that of the wt (N<sub>547</sub>D K<sub>D</sub> = 5.5 ± 0.4 µM (n = 3) vs. wt K<sub>D</sub> = 1.4 ± 0.1 µM (n = 3)). Moreover, we measured a reduced sensitivity to cAMP also in patch clamp experiments where addition of 5 µM cAMP caused a right shift in the V<sub>1/2</sub> of the mutant HCN2 channel of only 5 mV while the wt channel shifted by 12 mV (Figure 3—figure supplement 5). To exclude that the N<sub>547</sub>D mutation affects the overall structure of the CNBD, we performed the NMR (¹H-¹⁵N HSQC spectrum) analysis of the N<sub>547</sub>D CNBD mutant. Our data show that the protein is appropriately folded (Figure 3—figure supplement 6). In conclusion, since the N-bundle loop does not directly contact any of the residues of the cAMP binding pocket, these findings underscore a previously unaddressed role of the N-bundle loop in allosterically modulating cAMP binding to the CNBD (see Discussion).

## TRIP8b<sub>nano</sub> as a tool for the direct regulation of native HCN currents

Next, we asked whether the relatively short TRIP8b<sub>nano</sub> could be used to block cAMP-dependent modulation of HCN channels by delivering the peptide to full length channels. To this end, we dia-lyzed TRIP8b<sub>nano</sub> into the cytosol of HEK 293 T cells transfected either with HCN1, HCN2, or HCN4 channels. The peptide was added (10 µM) in the recording pipette together with a non-saturating concentration of cAMP (5 µM for HCN2, 1 µM for HCN4) expected to induce a ~ 10 mV rightward shift in the half-activation potential (V<sub>1/2</sub>) of the channels (Figure 4). No cAMP was added in the case of HCN1, because, in HEK 293 T cells, this isoform is already fully shifted to the right by the endoge-nous cAMP and does not respond further (Figure 4—figure supplement 1). Indeed, it is possible to induce a ~10 mV left shift in HCN1 V<sub>1/2</sub> by introducing the mutation R<sub>549</sub>E that prevents cAMP bind-ing to the CNBD (Figure 4—figure supplement 1).

**Figure 4A-C** show representative current traces recorded at four given voltages, in control solution, cAMP, and cAMP +10 µM TRIP8b<sub>nano</sub> (HCN4 and HCN2) or +10 µM TRIP8b<sub>nano</sub> only (HCN1) in the patch pipette. Already from a visual comparison of the most positive voltage at which the current appears measurable, it is evident that TRIP8b<sub>nano</sub> counteracts the activ-ating effect of cAMP on the voltage-dependent gating. In the case of HCN1, the effect of TRIP8b- nano can be observed without added cAMP for the aforementioned reasons. **Figure 4D-F** show the mean channel activation curves obtained from the above and other experiments. Fitting the Boltz-mann equation to the data (solid and dashed lines of Figure 4D-F, see Materials and methods for equation) yielded the half-activation potential values (V<sub>1/2</sub>) plotted in **Figure 4G**. The addition of TRIP8b<sub>nano</sub> prevents the cAMP-induced right shift of about 13 mV in HCN4 (V<sub>1/2</sub> = -102.8, -89.2, -102.1 mV, for control, cAMP and cAMP + TRIP8b<sub>nano</sub>, respectively), and of about 11 mV in HCN2 (V<sub>1/2</sub> = -93.7, -83.5, -94.5 mV, for control, cAMP and cAMP + TRIP8b<sub>nano</sub>, respectively). In HCN1, TRIP8b<sub>nano</sub> induced a left shift in V<sub>1/2</sub> of about 10 mV (from - 72.8 to -82 mV) which is comparable to that induced by the R<sub>549</sub>E mutation (from -72.7 to -80.4 mV) (Figure 4—figure supplement 1).

**Figure 4G** also shows the result of a control experiment performed on HCN4 where 10 µM TRIP8b<sub>nano</sub> was added to the extracellular medium (TRIP8b<sub>nano</sub> bath) in order to test if the peptide was able to cross the cell membrane (current traces and activation curves not shown). The V<sub>1/2</sub> value, which is similar to that of cAMP alone (-91.7 vs. -89.2 mV), confirmed that TRIP8b<sub>nano</sub> peptide
<figure><img src="image_7.png" /><figcaption>HCN4
+ cAMP</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>HCN4
+ cAMP
+ TRIP8b<sub>nano</sub></figcaption></figure>

<figure><img src="image_9.png" /><figcaption>B
HCN2</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>HCN2
+ cAMP</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>HCN2
+ cAMP
+ TRIP8b<sub>nano</sub></figcaption></figure>

<figure><img src="image_12.png" /><figcaption>C
HCN1</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>HCN1
+ TRIP8b<sub>nano</sub></figcaption></figure>

<figure><img src="image_14.png" /><figcaption>F

Voltage (mV)</figcaption></figure>

<figure><img src="image_15.png" /></figure>

<figure><img src="image_16.png" /></figure>

<figure><img src="image_17.png" /></figure>

<figure><img src="image_18.png" /></figure>

**Figure 4.** TRIP8b<sub>nano</sub> abolishes cAMP effect on HCN channel gating. (A-C) Representative whole-cell HCN4, HCN2 and HCN1 currents recorded, at the indicated voltages, with control solution or with cAMP (1 µM for HCN4 and 5 µM for HCN2) or with cAMP + 10 µM TRIP8b<sub>nano</sub> in the patch pipette (for HCN1, 10 µM TRIP8b<sub>nano</sub> only was added). (D-F) Mean activation curves measured from HCN4, HCN2 and HCN1 in control solution (open circles), cAMP (filled circles), cAMP +TRIP8b<sub>nano</sub>, or TRIP8b<sub>nano</sub> only in the case of HCN1 (filled triangles). Solid and dashed lines indicate Boltzmann fitting to the data (see Materials and methods). (G) Half activation potential (V<sub>1/2</sub>) values of HCN4 (blue), HCN2 (black) HCN1 (red) in control solution (open circle), cAMP (filled circle) and cAMP +TRIP8b<sub>nano</sub>, or TRIP8b<sub>nano</sub> only in the case of HCN1 (filled triangle). HCN4, control = -102.8 ± 0.3 mV; HCN4 +1 µM cAMP = -89.2 ± 0.6 mV; HCN4 +1 µM cAMP +10 µM TRIP8b<sub>nano</sub> = -102.1 ± 0.6 mV, HCN4 +1 µM cAMP in the patch pipette +10 µM TRIP8b<sub>nano</sub> in the bath solution = -91.7 ± 0.3 mV; HCN2, control = -93.7 ± 0.3 mV; HCN2 +5 µM cAMP = -83.5 ± 0.3 mV; HCN2 +5 µM cAMP +10 µM TRIP8b<sub>nano</sub> = -94.5 ± 0.6 mV; HCN1, control = -72.8 ± 0.2 mV; HCN1 +10 µM TRIP8b<sub>nano</sub> = -82 ± 0.5 mV. Data are presented as mean ± SEM.
Number of cells (N) ≥ 11. There is no significant difference between the controls and the addition of TRIP8b<sub>nano</sub> with (HCN4, HCN2) or without (HCN1) cAMP in the pipette. No significant difference was observed following the addition of TRIP8b<sub>nano</sub> in the bath. Statistical analysis performed with two-way ANOVA, followed by post-hoc Tukey test (p<0.001).

DOI: https://doi.org/10.7554/eLife.35753.013

The following figure supplements are available for figure 4:

**Figure supplement 1.** Comparison of half activation potentials (V<sub>1/2</sub>) of HCN1 WT and HCN1 R<sub>549</sub>E mutant (Chen et al., 2001).

DOI: https://doi.org/10.7554/eLife.35753.014

**Figure supplement 2.** Voltage-dependency of activation time constant (τ<sub>on</sub>) of HCN4 channels in control solution (open circles), 1 µM cAMP (filled circles), 1 µM cAMP + 10 µM TRIP8b<sub>nano</sub> (filled triangles).

DOI: https://doi.org/10.7554/eLife.35753.015
affects channel gating only if added to the intracellular solution presumably because it does not dif-
fuse through the cell membrane.

It is worth noting that TRIP8b_nano prevents other related effects of cAMP activation in HCN channels, such as the acceleration of activation kinetics (Wainger et al., 2001) and, for HCN2 only, the increase in maximal current (Chen et al., 2007; Hu et al., 2013). For example, the activation kinetics (τon) of HCN4 measured at –120 mV was: control = 2 ± 0.2 s, 1 µM cAMP = 1.2 ± 0.1 s, 1 µM cAMP + 10 µM TRIP8b_nano = 2 ± 0.3 s (Figure 4—figure supplement 2). Moreover, Figure 4B clearly shows that TRIP8b_nano fully prevented the increase in maximal current in HCN2.

Based on these results, we reckoned the peptide may be employed as a regulatory tool for native $I_f/I_h$ currents. As proof of principle, we tested whether TRIP8b_nano can modulate the frequency of action potential firing in SAN myocytes. In these cells, $I_f$ is key contributor of the diastolic depolarization phase of the pacemaker action potential cycle. Moreover, the autonomic nervous system modulates the frequency of action potential firing by changing intracellular cAMP levels, which in turn acts on f-HCN channel open probability (DiFrancesco, 1993). We thus recorded the native $I_f$ current in acutely isolated rabbit SAN myocytes with and without 10 µM TRIP8b_nano in the pipette solution (Figure 5A). Figure 5B shows that the averaged $I_f$ activation curve measured in presence of TRIP8b_nano is significantly shifted to hyperpolarized voltages compared to the control. This indicates that the peptide is displacing the binding of endogenous cAMP to native HCN channels. Moreover, when the experiment was repeated in the presence of 1 µM cAMP, TRIP8b_nano prevented the typical cAMP-dependent potentiation of the native $I_f$ current (Figure 5B). In light of these results, we tested whether TRIP8b_nano is also able to modulate cardiac automaticity by antagonizing basal cAMP. The data in Figure 5C show that TRIP8b_nano indeed significantly decreased the rate of action potential firing in single SAN cells. Strikingly, the observed 30% decrease in action potential rate corresponds to the effect induced by physiological concentrations of acetylcholine (DiFrancesco et al., 1989).

To conclusively prove that the inhibition of the native $I_f$ current was specifically due to TRIP8b_nano rather than caused by the dilution of the cellular content following whole cell configuration, we created a TAT version of TRIP8b_nano (hereafter TAT-TRIP8b_nano). Indeed, the TAT sequence allows the entry of biomolecules into a cell via endocytosis and/or direct translocation across the plasma membrane, thus leaving the cytosolic content unaltered (Guidotti et al., 2017).

We therefore tested whether both TRIP8b_nano and TAT-TRIP8b_nano were able to selectively inhibit the β-adrenergic stimulation of $I_f$ current, while leaving the potentiation of L-type Ca²⁺ current ($I_{Ca,L}$) unaltered. To this end, we recorded either the native $I_f$ or $I_{Ca,L}$ current from cardiomyocytes acutely isolated from mouse sinoatrial node (SAN) in the presence and in the absence of 10 µM TRIP8b_nano or TAT-TRIP8b_nano, before and after stimulation with 100 nM isoproterenol (ISO), a β-adrenergic receptor agonist (Figure 6). Strikingly, TRIP8b_nano prevented the isoproterenol-induced increase of $I_f$ current density, both when the peptide was added in the recording pipette solution (Figure 6A and B), and when it was used in the TAT version added to the bath (Figure 6A and C). The specificity of TRIP8b_nano for $I_f$ current was confirmed by the absent inhibition of basal $I_{Ca,L}$ (Figure 6D). In addition, we failed to record a significant difference in the isoproterenol-stimulated increase of the $I_{Ca,L}$ current density between the control condition and 10 µM TRIP8b_nano (Figure 6E) or TAT-TRIP8b_nano (Figure 6F) conditions. To test whether the TAT-TRIP8b_nano effect described above was exclusively due to TRIP8b_nano peptide, we repeated the experiments with a scrambled version of the peptide (TAT- (SCRAMBLED) TRIP8b_nano) to exclude that the effect could be due to the TAT sequence (Figure 6—figure supplement 1). We failed to observe a significant reduction in the responsiveness of $I_f$ to isoproterenol in the presence of TAT- (SCRAMBLED) TRIP8b_nano confirming that prevention of cAMP induced f- current stimulation was specific of the TRIP8b_nano sequence.

Discussion

TRIP8b-CNBD complex

In this study, we have identified the minimal binding peptide that reproduces the effects of TRIP8b
on HCN channel gating. The peptide is 40 aa long and binds the HCN CNBD with high affinity
($K_D$ = 1.4 µM). By solving the NMR structures of TRIP8b_nano and HCN CNBD in the bound form, we
generated a structural model of their complex. The model provides detailed information on this pro-
tein-protein interaction at the atomic level with implications on their physiological function. The data
<figure><img src="image_19.png" /><figcaption>**Figure 5.** Effects of TRIP8b_nano on voltage-dependent activation of I<sub>f</sub> and spontaneous rate in rabbit sinoatrial node (SAN) myocytes. (A) Representative whole-cell I<sub>f</sub> currents recorded, at the indicated voltages, in the control solution and in the presence of 10 µM TRIP8b_nano, without (top) and with 1 µM cAMP in the pipette (bottom). (B) Mean I<sub>f</sub> activation curves were measured using a two-step protocol (see Materials and methods) in control (filled circles) or in the presence of: 1 µM cAMP (open circles); 10 µM TRIP8b_nano (filled squares); 1 µM cAMP +10 µM TRIP8b_nano (open squares). Ligands were added in the patch pipette. Half activation potential (V<sub>1/2</sub>) of I<sub>f</sub> activation curves measured in control = -64.1 ± 0.4 mV or in the presence of: 1 µM cAMP = -59.9 ± 0.4 mV; 10 µM TRIP8b_nano = -67.7 ± 0.4 mV; 1 µM cAMP +10 µM TRIP8b_nano = -67.6 ± 0.7 mV. Data are presented as mean ± SEM. Number of cells (N) was ≥15. V<sub>1/2</sub> values are significantly different between each other's whit the exception of V<sub>1/2</sub> obtained in the presence TRIP8b_nano and cAMP +TRIP8b_nano. Statistical analysis performed with two-way ANOVA, followed by post-hoc Bonferroni test (*p<0.05) (C) (Left) Representative recordings of single SAN cell spontaneous activity in control and in the presence of 10 µM TRIP8b_nano. (Right) Mean spontaneous rate (Hz) recorded in control solution = 3.65 ± 0.29 Hz and in the presence of 10 µM TRIP8b_nano added to the pipette = 2.69 ± 0.27 Hz. Data are presented as mean ± SEM. Number of cells (N) was ≥7. Statistical analysis performed with t test (*p<0.05).</figcaption></figure>

DOI: https://doi.org/10.7554/eLife.35753.016

show that the minimal binding unit of TRIP8b, TRIP8b_nano, folds in two helices upon binding and sug-
gest that this region is intrinsically disordered when it is not bound. The model structurally validates
previous indirect evidence, which suggested that TRIP8b binds to two discrete elements of the
CNBD: the N-bundle loop and the C-helix (Saponaro et al., 2014). The complex forms by electro-
static interactions, which are spread throughout the contact surface. As a consequence of the inter-
action with TRIP8b_nano, the C-helix of CNBD increases in length, a behavior previously observed in
the case of cAMP binding (Puljung and Zagotta, 2013). This portion of C-helix includes the two resi-
dues R<sub>662</sub> and K<sub>665</sub> engaged in salt bridge formation with respectively E<sub>250</sub>/D<sub>257</sub> and E<sub>264</sub> of TRIP8b-
nano. It is important to note that these two cationic residues are also involved in cAMP binding
(Zagotta et al., 2003; Zhou and Siegelbaum, 2007; Lolicato et al., 2011). The finding that TRIP8b
and cAMP share the same binding sites on the C-helix provides a solid molecular explanation for
<figure><img src="image_20.png" /><figcaption>**Figure 6.** Effect of TRIP8b<sub>nano</sub> and TAT-TRIP8b<sub>nano</sub> on I<sub>f</sub> and I<sub>Ca,L</sub> in mouse sinoatrial node (SAN) myocytes. (A) Representative examples of I<sub>f</sub> recordings at −95 mV in control conditions (top), in 10 µM TRIP8b<sub>nano</sub> dialyzed cell (middle) and in cells perfused with 10 µM TAT-TRIP8b<sub>nano</sub> (bottom), before (black trace) and after (red trace) application of ISO 100 nM. The voltage-clamp protocol used for recordings is shown above current traces. (B) Mean normalized I<sub>f</sub> current density recorded at −95 mV in absence and in presence of 10 µM TRIP8b<sub>nano</sub> in the patch pipette, before (open bars) and after *Figure 6 continued on next page*</figcaption></figure>
Figure 6 continued

(filled bars) 100 nM ISO perfusion. Data are presented as mean ± SEM. Number of cells (N) ≥ 6. Statistical analysis performed with two-way ANOVA test, followed by Sidak multiple comparisons test (\*\*p<0.01). (C) Mean normalised $I_f$ current density recorded at –95 mV in control solution or in the solution containing 10 µM TAT-TRIP8b_nano, in absence (open bars) and in the presence (filled bars) of 100 nM ISO. Data are presented as mean ± SEM. Number of cells (N) ≥ 8. Statistical analysis performed with two-way ANOVA test, followed by Sidak multiple comparisons test (\*\*\*p<0.0001). (D) Representative examples of $I_{Ca,L}$ recordings at –10 mV in control conditions (top), in 10 µM TRIP8b_nano dialyzed cell (middle) and in cells perfused with 10 µM TAT-TRIP8b_nano (bottom), before (black trace) and after (red trace) application of ISO 100 nM. The voltage-clamp protocol used for recordings is shown above current traces. (E) Mean normalized $I_{Ca,L}$ current density recorded at –10 mV in absence and in presence of 10 µM TRIP8b_nano in the patch pipette, before (open bars) and after (filled bars) 100 nM ISO perfusion. Data are presented as mean ± SEM. Number of cells (N) ≥ 8. Statistical analysis performed with two-way ANOVA test, followed by Sidak multiple comparisons test (\*\*\*p<0.001; \*\*\*\*p<0.0001). (F) Mean normalized $I_{Ca,L}$ current density recorded at –10 mV in control solution or in the solution containing 10 µM TAT-TRIP8b_nano, in absence (open bars) and in the presence (filled bars) of 100 nM ISO. Data are presented as mean ± SEM. Number of cells (N) ≥ 7. Statistical analysis performed with two-way ANOVA test, followed by Sidak multiple comparison test (\*\*p<0.01; \*\*\*p<0.001).

DOI: https://doi.org/10.7554/eLife.35753.017

The following figure supplement is available for figure 6:

**Figure supplement 1.** Effect of TAT- (SCRAMBLED) TRIP8b_nano on $I_f$ in mouse sinoatrial node (SAN) myocytes.

DOI: https://doi.org/10.7554/eLife.35753.018

functional data, which imply a competition between the two regulators (Han et al., 2011; DeBerg et al., 2015; Bankston et al., 2017). Another study, however, has indicated that a direct competition model cannot fully explain the mutually antagonistic effect of the two ligands (Hu et al., 2013). Specifically, the fact that the inhibitory effect of TRIP8b on channel activity persists even at saturating cAMP concentrations advocated an allosteric component in the regulation mechanism. Our data, showing that the N547D mutation in the N-bundle loop controls cAMP affinity in the binding pocket support the conclusion that the N-bundle loop allosterically controls cAMP binding. This is not surprising, given its crucial role of mechanically transducing to the pore the cAMP binding event within the CNBD (Saponaro et al., 2014).

The structural model also explains why a previously identified peptide selected by Lyman et al. (2017) failed to reproduce the binding affinity of the TRIP8b_core for the CNBD. This 37 aa long fragment is lacking one important contact residue, namely E240, which, in our model, forms a salt bridge with R650 of the CNBD. The loss of one crucial interaction is presumably the reason for the major decrease in affinity (about 20 times lower) reported for this peptide.

## TRIP8b_nano as a tool for modulating native $I_f$ currents

In functional assays, we showed that TRIP8b_nano binds the HCN channel CNBD with high affinity and fully abolishes the cAMP effect in all tested isoforms (HCN1, 2 and 4).

Given the small size of the peptide (<5 kDa), TRIP8b_nano is a good candidate for in vivo delivery into intact cells. As a proof of concept, we fused TRIP8b_nano to an internalization sequence, the TAT peptide (YGRKKRRQRRGG). This arginine-rich Cell Penetrating Peptide (CPP) from HIV has been used in several studies as a vehicle for the delivery of large molecules across the plasma membrane (Guidotti et al., 2017). In our case, the challenge was to construct a TAT-fusion protein, which would be efficiently delivered in the cell without compromising TRIP8b_nano function. Indeed, covalent conjugation of a CPP may negatively affect both the function of the cargo, and the cell-penetrating efficacy of the CPP-peptide chimera (Kristensen et al., 2016). The design of the construct was greatly supported by the detailed knowledge of the electrostatic interactions with the target protein CNBD, provided by the NMR model structure. This structure suggested that the polycationic TAT sequence would be best linked to the N-terminus of TRIP8b_nano to avoid interference with the cationic residues of CNBD, mainly located in the distal region of C-helix, which are crucial for the binding of the peptide. From test experiments with the TAT-TRIP8b_nano peptide in SAN myocytes, we can conclude that this strategy was successful in that: (i) the peptide is efficiently delivered inside the cells; (ii) it is kept in its active conformation; (iii) the TAT sequence did not damage cell membranes and did not interfere with the basic features of $I_f$ and $I_{Ca,L}$ currents; (iv) the modification did not affect the proteolytic stability of the TRIP8b_nano peptide at least in the time frame of our experiments (30 min to 1 hr).
In conclusion, we successfully used the miniaturized TRIP8b<sub>nano</sub> peptide to selectively control native I<sub>f</sub> currents and the rate of spontaneous firing in SAN myocytes. Unlike channel blockers, which inhibit ionic currents, the peptide only interferes with the cAMP-based regulation of HCN channels, while leaving basal HCN functions unaltered. In addition and in contrast to even the most selective blockers, it is selective for HCN and it does not interfere with other cAMP-modulated channels present in the SAN, such as L-type Ca<sup>2+</sup> channels. Collectively, this makes TRIP8b<sub>nano</sub> a promising tool in targeted therapeutic interventions.

## Materials and methods

### Key resources table

<table><thead><tr><td>Reagent type species</td><td>Designation</td><td>Source or reference</td><td>Identifiers</td><td>Additional information</td></tr></thead><tbody><tr><td>Gene (human)</td><td>HCN1</td><td>Xention Ltd. (Cambridge, UK)</td><td>NM_021072.3</td><td></td></tr><tr><td>Gene (mouse)</td><td>HCN2</td><td>PMID: 11331358</td><td>NM_008226.2</td><td>Laboratory of Steven A. Siegelbaum</td></tr><tr><td>Gene (rabbit)</td><td>HCN4</td><td>PMID: 10212270</td><td>NM_001082707.1</td><td></td></tr><tr><td>Gene (mouse)</td><td>TRIP8b</td><td>PMID: 19555649</td><td></td><td>Laboratory of Steven A. Siegelbaum</td></tr><tr><td>Strain, strain background (E. coli)</td><td>DH5α</td><td>Thermo Fisher Scientific</td><td></td><td></td></tr><tr><td>Strain, strain background (E. coli)</td><td>Stbl2</td><td>Thermo Fisher Scientific</td><td></td><td></td></tr><tr><td>Strain, strain background (Mus musculus)</td><td>Male or female C57BL/6J mice</td><td>The Jackson Laboratory</td><td>RRID:MGI:5653012</td><td></td></tr><tr><td>Strain, strain background (Oryctolagus cuniculus)</td><td>New Zealand white female rabbits</td><td>Envigo</td><td>ID strain:HsdOkd:NZW</td><td></td></tr><tr><td>Cell line (human)</td><td>HEK 293T</td><td>ATCC</td><td>RRID:CVCL_0063</td><td>Tested negative for mycoplasma</td></tr><tr><td>Biological sample (Mus musculus)</td><td>Isolated adult Sinoatrial node (SAN) cardiomyocytes</td><td>PMID: 11557233</td><td></td><td></td></tr><tr><td>Biological sample (Oryctolagus cuniculus)</td><td>Isolated adult Sinoatrial node (SAN) cardiomyocytes</td><td>PMID: 2432247</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>pET-52b (plasmid)</td><td>EMD Millipore</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>modified pET-24b (plasmid)</td><td>Laboratory of Daniel L. Minor, Jr.</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>pcDNA 3.1 (plasmid)</td><td>Clontech Laboratories</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>pCI (plasmid)</td><td>Promega</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (cDNA)</td><td>This paper</td><td></td><td>Made by PCR and cloning; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (E240R) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of TRIP8bnano wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (E241R) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of TRIP8bnano wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (E264K) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of TRIP8bnano wt; see Constructs</td></tr></tbody></table>

Continued on next page
Continued

<table><thead><tr><th>Reagent type species</th><th>Designation</th><th>Source or reference</th><th>Identifiers</th><th>Additional information</th></tr></thead><tbody><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (E265K) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of TRIP8bnano wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bnano (D252N) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of TRIP8bnano wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8bcore (cDNA)</td><td>PMID: 25197093</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN2 CNBD (cDNA)</td><td>PMID: 25197093</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN2 CNBD (N547D) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of human HCN2 CNBD wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN2 CNBD (K665E) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of human HCN2 CNBD wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN2 CNBD (K666E) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of human HCN2 CNBD wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN2 CNBD (R650E) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of human HCN2 CNBD wt; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>human HCN1 (cDNA)</td><td>Xention Ltd. (Cambridge, UK)</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>TRIP8b (1a4) (cDNA)</td><td>This paper</td><td></td><td>Made by PCR and cloning; see Constructs</td></tr><tr><td>Recombinant DNA reagent</td><td>mouse HCN2 (cDNA)</td><td>Laboratory of Steven A. Siegelbaum</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>rabbit HCN4 (cDNA)</td><td>PMID: 10212270</td><td></td><td></td></tr><tr><td>Recombinant DNA reagent</td><td>mouse HCN2 (N520D) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of mouse HCN2 wt; see Constructs</td></tr><tr><td>Sequence-based reagent</td><td>human HCN1 (R549E) (cDNA)</td><td>This paper</td><td></td><td>Made by site-directed mutagenesis of human HCN1 wt; see Constructs</td></tr><tr><td>Peptide, recombinant protein</td><td>TAT-TRIP8b<sub>nano</sub> (YGRKKRRQRRRGG-NHSLEEEFERAKAAAVESTEFWDKMQAE EWEEMARRNWISEN)</td><td>CASLO ApS</td><td></td><td></td></tr><tr><td>Peptide, recombinant protein</td><td>TAT-(SCRAMBLED)TRIP8b<sub>nano</sub> (YGRKKRRQRRGG-RNEAEAAE VAQKDMINERARTHEFEWESWE MWENLSESK)</td><td>CASLO ApS</td><td></td><td></td></tr><tr><td>Commercial assay or kit</td><td>QuickChange Lightning Site-Directed Mutagenesis Kit</td><td>Agilent</td><td></td><td></td></tr><tr><td>Commercial assay or kit</td><td>Thermo Scientific TurboFect Transfection Reagent</td><td>Thermo Fisher Scientific</td><td></td><td></td></tr><tr><td>Chemical compound, drug</td><td>Adenosine 3', 5'-cyclic monophosphate (cAMP)</td><td>SIGMA</td><td></td><td></td></tr><tr><td>Software, algorithm</td><td>Clampfit 10.5/10.7</td><td>Molecular Devices</td><td>RRID:SCR_011323</td><td></td></tr><tr><td>Software, algorithm</td><td>CYANA-2.1</td><td>L. A. Systems, Inc.</td><td></td><td></td></tr><tr><td>Software, algorithm</td><td>AMBER 12.0</td><td>http://pyenmr.cerm.unifi.it/access/index/amps-nmr</td><td></td><td></td></tr><tr><td>Software, algorithm</td><td>HADDOCK2.2</td><td>www.wenmr.eu</td><td></td><td></td></tr></tbody></table>
## Constructs

The cDNA fragment encoding residues 235–275 (TRIP8b_nano) of mouse TRIP8b (splice variant 1a4) was cloned into pET-52b (EMD Millipore) downstream of a Strep (II) tag sequence, while the cDNA fragment encoding residues 521–672 of human HCN2 (HCN2 CNBD) was cloned, in a previous study, into a modified pET-24b downstream of a double His<sub>6</sub>-maltose-binding protein (MBP) (*Saponaro et al., 2014*). The cDNA encoding full-length human HCN1 channel and mouse TRIP8b (1a4) were cloned into the mammalian expression vector pcDNA 3.1 (Clontech Laboratories), while mouse HCN2 channel and rabbit HCN4 channel were cloned into the mammalian expression vector pCI (Promega). Mutations were generated by site-directed mutagenesis (QuikChange site-directed mutagenesis kit; Agilent Technologies) and confirmed by sequencing.

## Preparation of proteins

The HCN2 CNBD WT and mutant proteins, as well as the TRIP8b<sub>core</sub> and TRIP8b<sub>nano</sub> proteins (WT and mutants) were produced and purified following the procedure previously described (*Saponaro et al., 2014*).

## Structure calculation of the cAMP-free human HCN2 CNBD in complex with TRIP8b_nano and vice versa

NMR experiments were acquired on Bruker Avance III 950, 700 and 500 MHz NMR spectrometers equipped with a TXI-cryoprobe at 298 K. The acquired triple resonance NMR experiments for the assignment of backbone resonances of cAMP-free HCN2 CNBD (CNBD hereafter) in complex with TRIP8b_nano and vice versa are summarized in Figure 3-source data 1. ¹⁵N, ¹³C', ¹³C<sub>α</sub>, ¹³C<sub>β</sub>, and H<sub>α</sub> chemical shifts were used to derive φ and ψ dihedral angles by TALOS + program (*Cornilescu et al., 1999*) for both CNBD and TRIP8b_nano. For TRIP8b_nano, CYANA-2.1 structure calculation (*Güntert and Buchner, 2015*) was performed using 68 φ and ψ dihedral angles and 40 backbone hydrogen bonds as input. For CNBD, CYANA-2.1 structure calculation was performed using 108 φ and ψ dihedral angles, combined with the NOEs obtained in our previous determination of the cAMP-free form of the CNBD (*Saponaro et al., 2014*) for those regions not affected by the interaction with TRIP8b_nano. The 10 conformers of TRIP8b_nano and CNBD with the lowest residual target function values were subjected to restrained energy minimization with AMBER 12.0 (*Case, 2012*) (http://pyenmr.cerm.unifi.it/access/index/amps-nmr) and used as input in docking calculations.

## Docking calculations

Docking calculations were performed with HADDOCK2.2 implemented in the WeNMR/West-Life GRID-enabled web portal (www.wenmr.eu). The docking calculations are driven by ambiguous interaction restraints (AIRs) between all residues involved in the intermolecular interactions (*Dominguez et al., 2003*). Active residues of the CNBD were defined as the surface exposed residues (at least 50% of solvent accessibility), which show chemical shift perturbation upon TRIP8b_nano binding.

The assignment of the CNBD bound to TRIP8b_nano allowed to highlight the residues of CNBD whose backbone featured appreciable Combined Chemical Shift Perturbation (CSP) (Figure 3—figure supplement 1). The combined CSP (Δ<sub>HN</sub>) is given by the equation Δ<sub>HN</sub>=((H<sub>Nfree</sub>-H<sub>Nbound</sub>)<sup>2</sup>+(-N<sub>free</sub>-N<sub>bound</sub>)/5)<sup>2</sup>/2)½ (*Garrett et al., 1997*).

Passive residues of CNBD were defined as the residues close in space to active residues and with at least 50% solvent accessibility.

In the case of TRIP8b_nano, the conserved stretch E<sub>239</sub>-E<sub>243</sub>, located in helix N, was defined as active region in a first docking calculation, while all the other solvent accessible residues of the peptide were defined as passive. This docking calculation generated several clusters. A post-docking filter step allowed us to select those clusters having an orientation of TRIP8b_nano bound to CNBD in agreement with a DEER study on the CNBD - TRIP8b_nano interaction (*DeBerg et al., 2015*). The selected clusters grouped in two classes on the basis of the orientation of helix N of TRIP8b_nano (N) relative to CNBD (Figure 3—figure supplement 2. A second docking calculation was subsequently performed introducing also residues E<sub>264</sub>-E<sub>265</sub>, located in helix C of TRIP8b_nano as active residues. The active residues for CNBD were the same used for the first calculation. For this second HADDOCK calculation, 14 clusters were obtained and ranked according to their HADDOCK score.
Among them only four clusters showed both an orientation of TRIP8b_nano bound to CNBD in agree-
ment with the DEER study (DeBerg et al., 2015) and the involvement of E<sub>239</sub>-E<sub>243</sub> stretch of TRIP8b-
nano in the binding to CNBD. These clusters were manually analyzed and subjected to a per-cluster
re-analysis following the protocol reported in http://www.bonvinlab.org/software/haddock2.2/analy-
sis/#reanal. From this analysis, it resulted that the top-ranking cluster, i.e. the one with the best ener-
getic and scoring functions, has a conformation in agreement with mutagenesis experiments
(Figure 3—figure supplement 3). Energy parameters (van der Waals energy, electrostatic energy,
desolvation energy, and the penalty energy due to violation of restraints) for this complex model are
reported in Figure 3—source data 2.

Both docking calculations were performed using 10 NMR conformers of both the CNBD and the TRIP8b_nano structures calculated as described above. In the TRIP8b_nano structures the unfolded N- and C-terminal regions were removed, while in the CNBD structures only the unfolded N-terminal region was removed. This is because the C-terminal region of the CNBD is known to comprise residues involved in TRIP8b_nano binding (Saponaro et al., 2014). Flexible regions of the proteins were defined based on the active and passive residues plus two preceding and following residues. The residue solvent accessibility was calculated with the program NACCESS (Hubbard and Hornton, 1993). In the initial rigid body docking calculation phase, 5000 structures of the complex were generated, and the best 400 in terms of total intermolecular energy were further submitted to the semi-flexible simulated annealing and a final refinement in water. Random removal of the restraints was turned off. The number of flexible refinement steps was increased from the default value of 500/500/1000/1000 to 2000/2000/2000/4000. The final 400 structures were then clustered using a cutoff of 5.0 Å of RMSD to take into consideration the smaller size of protein-peptide interface.

Electrophysiology of HEK 293 T cells

HEK 293 T cells were cultured in Dulbecco's modified Eagle's medium (Euroclone) supplemented with 10% fetal bovine serum (Euroclone), 1% Pen Strep (100 U/mL of penicillin and 100 µg/ml of streptomycin), and stored in a 37°C humidified incubator with 5% CO2. The plasmid containing cDNA of wild-type and mutant HCN1, HCN2 and HCN4 channels (1 µg) was co-transfected for transient expression into HEK 293 T cells with a plasmid containing cDNA of Green Fluorescent Protein (GFP) (1.3 µg). For co-expression with TRIP8b (1a-4), HEK 293 T cells were transiently transfected with wild-type (wt) and/or mutant human HCN1 cDNA (1 µg), wt TRIP8b (1a-4) cDNA (1 µg) and cDNA of Green Fluorescent Protein (GFP) (0.3 µg).

One day after transfection, GFP-expressing cells were selected for patch-clamp experiments in whole-cell configuration. The experiments were conducted at R.T. The pipette solution in whole cell experiments contained: 10 mM NaCl, 130 mM KCl, 1 mM egtazic acid (EGTA), 0.5 mM MgCl₂, 2 mM ATP (Mg salt) and 5 mM HEPES-KOH buffer (pH 7.4). The extracellular bath solution contained 110 mM NaCl, 30 mM KCl, 1.8 mM CaCl₂, 0.5 mM MgCl₂ and 5 mM HEPES-KOH buffer (pH 7.4).

TRIP8b_nano was added (10 µM) to the pipette solution. cAMP was added at different concentration to the pipette solution depending on the HCN isoform used: 0 µM for HCN1, 5 µM for HCN2 and 1 µM for HCN4.

Whole-cell measurements of HCN channels were performed using the following voltage clamp protocol depending on the HCN isoform measured: for HCN1, holding potential was –30 mV (1 s), with steps from –20 mV to –120 mV (10 mV interval, 3.5 s) and tail currents recorded at –40 mV (3 s); for HCN2, holding potential was –30 mV (1 s), with steps from –40 mV to –130 mV (10 mV interval, 5 s) and tail currents recorded at –40 mV (5 s); for HCN4, holding potential was –30 mV (1 s), steps from –30 mV to –165 mV (15 mV interval, 4.5 s) and tail currents were recorded at –40 mV (5 s). Current voltage relations and activation curves were obtained by the above activation and deactivation protocols and analyzed by the Boltzmann equation, see data analysis.

Isolation and electrophysiology of rabbit sinoatrial node cells

Animal protocols conformed to the guidelines of the care and use of laboratory animals established by Italian and European Directives (D. Lgs n° 2014/26, 2010/63/UE). New Zealand white female rabbits (0.8–1.2 kg) were anesthetized (xylazine 5 mg/Kg, i.m.), and euthanized with an overdose of sodium thiopental (i.v.); hearts were quickly removed, and the SAN region was isolated and cut in small pieces. Single SAN cardiomyocytes were isolated following an enzymatic and mechanical
procedure as previously described (DiFrancesco et al., 1986). Following isolation, cells were maintained at 4°C in Tyrode solution: 140 mM NaCl, 5.4 mM KCl, 1.8 mM CaCl₂, 1 mM MgCl₂, 5.5 mM D-glucose, 5 mM HEPES-NaOH (pH 7.4).

For patch clamp experiments cells were placed in a chamber on an inverted microscope and experiments were performed in the whole-cell configuration at 35 ± 0.5°C. The pipette solution contained: 10 mM NaCl, 130 mM KCl, 1 mM egtazic acid (EGTA), 0.5 mM MgCl₂, and 5 mM HEPES-KOH buffer (pH 7.2). The $I_f$ current was recorded from single cells superfused with Tyrode solution with 1 mM BaCl₂, and 2 mM MnCl₂.

$I_f$ activation curves were obtained using a two-step protocol in which test voltage steps (from -30 to -120 mV, 15 mV interval) were applied from a holding potential of -30 mV and were followed by a step to -125 mV. Test steps had variable durations so as to reach steady-state activation at all voltages. Analysis was performed with the Boltzmann equation (see data analysis).

In current-clamp studies, spontaneous action potentials were recorded from single cells superfused with Tyrode solution, and rate was measured from the interval between successive action potential. When indicated cAMP (1 µM) and/or nanoTRIP8b (10 µM) were added to the pipette solution.

## Isolation and electrophysiology of mouse sinoatrial node cells

Mice were killed by cervical dislocation under general anesthesia consisting of 0.01 mg/g xylazine (2% Rompun; Bayer AG), 0.1 mg/g ketamine (Imalgène; Merial) and 0.04 mg/g of Na-pentobarbital (Euthanasol VET, Laboratoire TVM, Lempdes, France), and beating hearts were quickly removed. The SAN region was excised in warmed (35°C) Tyrode's solution containing: 140 mM NaCl, 5.4 mM KCl, 1.8 mM CaCl₂, 1 mM MgCl₂, 1 mM Hepes-NaOH (pH = 7.4), and 5.5 mM D-glucose and cut in strips. Strips were then transferred into a 'low-Ca²⁺-low-Mg²⁺' solution containing: 140 mM NaCl; 5.4 mM KCl, 0.5 mM MgCl₂, 0.2 mM CaCl₂, 1.2 mM KH₂PO₄, 50 mM taurine, 5.5 mM D-glucose, 1 mg/ml bovine serum albumin (BSA), 5 mM Hepes-NaOH (pH = 6.9).

Tissue was digested by adding Liberase TH (0.15 mg/ml, Roche Diagnostics GmbH, Mannheim, Germany), elastase (1.9 U/ml, Worthington, Lakewood). Digestion was carried out for a variable time of 15-18 min at 35°C. Tissue strips were then washed and transferred into a modified 'Kraftbrühe' (KB) medium containing: 70 mM L-glutamic acid, 20 mM KCl, 80 mM KOH, (±) 10 mM D-b-OH-butyric acid; 10 mM KH₂PO₄, 10 mM taurine, 1 mg/ml BSA and 10 mM Hepes-KOH (pH = 7.4).

Single SAN cells were isolated by manual agitation in KB solution at 35°C for 30-50 s.

Cellular automaticity was recovered by re-adapting the cells to a physiological extracellular Ca²⁺ concentration by addition of a solution containing: 10 mM NaCl, 1.8 mM CaCl₂ and normal Tyrode solution containing BSA (1 mg/ml). The final storage solution contained: 100 mM NaCl, 35 mM KCl, 1.3 mM CaCl₂, 0.7 mM MgCl₂, 14 mM L-glutamic acid, (±) 2 mM D-b-OH-butyric acid, 2 mM KH₂PO₄, 2 mM taurine, 1 mg/ml BSA, (pH = 7.4). Cells were then stored at room temperature until use. All chemicals were from SIGMA (St Quentin Fallavier, France).

For electrophysiological recording, SAN cells in the storage solution were harvested in special custom-made recording plexiglas chambers with glass bottoms for proper cell attachment and mounted on the stage of an inverted microscope (Olympus IX71) and perfused with normal Tyrode solution. The recording temperature was 36°C. We used the whole-cell variation of the patch-clamp technique to record cellular ionic currents, by employing a Multiclamp 700B (Axon Instruments Inc., Foster USA) patch clamp amplifier. Recording electrodes were fabricated from borosilicate glass, by employing a WZ DMZ-Universal microelectrode puller (Zeitz-Instruments Vertriebs GmbH, Martinsried, Germany).

$I_f$ was recorded under standard whole-cell configuration during perfusion of standard Tyrode's containing 2 mM BaCl₂ to block $I_{K1}$. Patch-clamp pipettes were filled with an intracellular solution containing: 130 mM KCl, 10 mM NaCl, 1 mM EGTA, 0.5 mM MgCl₂ and 5 mM HEPES (pH 7.2).

For recording of L-type Ca²⁺ currents, pipette solution contained: 125 mM CsOH, 20 mM tetraethylammonium chloride (TEA-Cl), 1.2 mM CaCl₂, 5 mM Mg-ATP, 0.1 mM Li₂-GTP, 5 mM EGTA and 10 mM HEPES (pH 7.2 with aspartate). 30 µM TTX (Latoxan, Portes lès Valence, France) to block INa was added to external solution containing: 135 mM tetraethylammonium chloride (TEA-Cl), 4 mM CaCl₂, 10 mM 4-amino-pyridine, 1 mM MgCl₂, 10 mM HEPES and 1 mg/ml Glucose (pH 7.4 with TEA-OH).
Electrodes had a resistance of about 3 MΩ. Seal resistances were in the range of 2–5 GΩ. 10 μM TRIPb8<sub>nano</sub> was added to pipette solution. 10 μM TAT-TRIPb8<sub>nano</sub> was added in cell storage solution for at least 30 min before patch clamp recording.

## TAT-peptides

TAT-TRIP8b<sub>nano</sub> (YGRKKRRQRRRGG-NHSLEEFERAKAAVESTEFWDKMQAEWEEMARRNWISEN, TAT sequence is shown in bold type) and TAT-(SCRAMBLED)TRIP8b<sub>nano</sub> (YGRKKRRQRRRGG-RNEAEAAEVAQKDMINERARTHEFEWESWEMWENLSESK, TAT sequence is shown in bold type) were purchased from CASLO ApS. TAT-peptides were dissolved in Milliq water (1.5 mM) and added to the petri dish at the final concentration (10 μM) 30 min before current recordings. During the patch clamp experiments, cells were perfused with standard Tyrode with 2 mM BaCl<sub>2</sub> (see above) without the peptides. Recordings from the same petri dish were performed over a time window of 10 to 60 min in peptide-free solution

## Data analysis

Data were acquired at 1 kHz using an Axopatch 200B amplifier and pClamp10.5 or 10.7 software (Axon Instruments). Data were analyzed off-line using Clampfit 10.5 or 10.7 (Molecular Devices) and Origin 2015 or 16 (OriginLab Corp., Northampton MA). Activation curves were analyzed by the Boltzmann equation, y = 1/{1 + exp[(V-V<sub>1/2</sub>)/s]}, where y is fractional activation, V is voltage, V<sub>1/2</sub> half-activation voltage, and s the inverse slope factor (mV) (DiFrancesco, 1999). Mean activation curves were obtained by fitting individual curves from each cell to the Boltzmann equation and then averaging all curves obtained. Activation time constants (τ<sub>on</sub>) were obtained by fitting a single exponential function,

I=I<sub>0</sub> exp(-t/τ) to current traces recorded at the indicated voltages.

## Ethics statement

Experiments on rabbit SAN cells were performed using left-over cells obtained during experiments approved by the Animal Welfare Body of the University of Milan and by the Italian Ministry of Health (license n.1127/2015-PR). Animal procedures were conformed to the guidelines of the care and use of laboratory animals established by Italian and European Directives (D. Lgs n° 2014/26, 2010/63/UE).

Mouse primary pacemaker cells were isolated from adult C57BL/6J mice as previously described (Mangoni and Nargeot, Cardiovasc Res 2001), in accordance with the Guide for the Care and Use of Laboratory Animals (eighth edition, 2011), published by the US National Institute of Health and European directives (2010/63/EU). The protocol was approved by the ethical committee of the University of Montpellier and the French Ministry of Agriculture (protocol N°: 2017010310594939).

## Acknowledgements

This work has been supported by Fondazione CARIPLO grant 2014–0796 to AM, BS and LB, by 2016 Schaefer Research Scholars Program of Columbia University to AM, by European Research Council (ERC) 2015 Advanced Grant 495 (AdG) n. 695078 noMAGIC to AM and GT, by National Institutes for Health Grant R01 NS036658 to BS, by Instruct-ERIC and national member subscriptions to LB, and by Accademia Nazionale dei Lincei (Giuseppe Levi foundation) to AS. We specially thank the EU ESFRI Instruct Core Centre CERM-Italy.

## Additional information

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Accademia Nazionale dei Lincei</td><td>Postdoctoral Fellowship in Neurobiology research</td><td>Andrea Saponaro</td></tr><tr><td>H2020 European Research Council</td><td>2015 Advanced Grant 495 n. 695078 noMAGIC</td><td>Gerhard Thiel Anna Moroni</td></tr></tbody></table>
<table><tr><td>Fondazione Cariplo</td><td>2014-0796</td><td>Lucia Banci<br>Bina Santoro<br>Anna Moroni</td></tr><tr><td>Instruct-ERIC</td><td></td><td>Lucia Banci</td></tr><tr><td>National Institutes of Health</td><td>R01 NS036658</td><td>Bina Santoro</td></tr><tr><td>Columbia University</td><td>2016 Schaefer Research<br>Scholars Program</td><td>Anna Moroni</td></tr></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

## Author contributions

Andrea Saponaro, Conceptualization, Investigation, Methodology, Writing—original draft; Francesca Cantini, Formal analysis, Investigation, Methodology; Alessandro Porro, Annalisa Bucchi, Vincenzo Maione, Formal analysis, Investigation; Dario DiFrancesco, Methodology; Chiara Donadoni, Bianca Introini, Investigation; Pietro Mesirca, Data curation, Formal analysis, Investigation; Matteo E Mangoni, Conceptualization, Data curation, Writing—original draft; Gerhard Thiel, Conceptualization, Supervision, Writing—original draft; Lucia Banci, Conceptualization, Data curation, Supervision; Bina Santoro, Conceptualization, Supervision, Funding acquisition, Writing—original draft; Anna Moroni, Conceptualization, Resources, Supervision, Funding acquisition, Methodology, Writing—original draft

## Author ORCIDs

Andrea Saponaro [http://orcid.org/0000-0001-5035-5174](http://orcid.org/0000-0001-5035-5174)
Francesca Cantini [http://orcid.org/0000-0003-0526-6732](http://orcid.org/0000-0003-0526-6732)
Vincenzo Maione [http://orcid.org/0000-0002-8229-6612](http://orcid.org/0000-0002-8229-6612)
Matteo E Mangoni [https://orcid.org/0000-0002-8892-3373](https://orcid.org/0000-0002-8892-3373)
Bina Santoro [http://orcid.org/0000-0002-4277-1992](http://orcid.org/0000-0002-4277-1992)
Anna Moroni [http://orcid.org/0000-0002-1860-406X](http://orcid.org/0000-0002-1860-406X)

## Ethics

Animal experimentation: Experiments on rabbit SAN cells were performed using left-over cells obtained during experiments approved by the Animal Welfare Body of the University of Milan and by the Italian Ministry of Health (license n.1127/2015-PR). Animal procedures were conformed to the guidelines of the care and use of laboratory animals established by Italian and European Directives (D. Lgs no 2014/26, 2010/63/UE). Mouse primary pacemaker cells were isolated from adult C57BL/ 6J mice as previously described (Mangoni and Nargeot, Cardiovasc Res 2001), in accordance with the Guide for the Care and Use of Laboratory Animals (eighth edition, 2011), published by the US National Institute of Health and European directives (2010/63/EU). The protocol was approved by the ethical committee of the University of Montpellier and the French Ministry of Agriculture (protocol N°: 2017010310594939).

## Decision letter and Author response

Decision letter [https://doi.org/10.7554/eLife.35753.021](https://doi.org/10.7554/eLife.35753.021)
Author response [https://doi.org/10.7554/eLife.35753.022](https://doi.org/10.7554/eLife.35753.022)

## Additional files

### Supplementary files

* Transparent reporting form

DOI: [https://doi.org/10.7554/eLife.35753.019](https://doi.org/10.7554/eLife.35753.019)

### Data availability

All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figure 3.
## References

Bankston JR, DeBerg HA, Stoll S, Zagotta WN. 2017. Mechanism for the inhibition of the cAMP dependence of HCN ion channels by the auxiliary subunit TRIP8b. *Journal of Biological Chemistry* **292**:17794-17803.
DOI: https://doi.org/10.1074/jbc.M117.800722, PMID: 28864772

Case DA. 2012. AMBER 12. San Francisco: University of California.

Chen S, Wang J, Siegelbaum SA. 2001. Properties of hyperpolarization-activated pacemaker current defined by coassembly of HCN1 and HCN2 subunits and basal modulation by cyclic nucleotide. *The Journal of General Physiology* **117**:491-504. DOI: https://doi.org/10.1085/jgp.117.5.491, PMID: 11331358

Chen S, Wang J, Zhou L, George MS, Siegelbaum SA. 2007. Voltage sensor movement and cAMP binding allosterically regulate an inherently voltage-independent closed-open transition in HCN channels. *The Journal of General Physiology* **129**:175-188. DOI: https://doi.org/10.1085/jgp.200609585, PMID: 17261842

Cornilescu G, Delaglio F, Bax A. 1999. Protein backbone angle restraints from searching a database for chemical shift and sequence homology. *Journal of Biomolecular NMR* **13**:289-302. DOI: https://doi.org/10.1023/A:1008392405740, PMID: 10212987

Craven KB, Zagotta WN. 2004. Salt bridges and gating in the COOH-terminal region of HCN2 and CNGA1 channels. *The Journal of General Physiology* **124**:663-677. DOI: https://doi.org/10.1085/jgp.200409178, PMID: 15572346

DeBerg HA, Bankston JR, Rosenbaum JC, Brzovic PS, Zagotta WN, Stoll S. 2015. Structural mechanism for the regulation of HCN ion channels by the accessory protein TRIP8b. *Structure* **23**:734-744. DOI: https://doi.org/10.1016/j.str.2015.02.007, PMID: 25800552

DiFrancesco D, Ducouret P, Robinson RB. 1989. Muscarinic modulation of cardiac rate at low acetylcholine concentrations. *Science* **243**:669-671. DOI: https://doi.org/10.1126/science.2916119, PMID: 2916119

DiFrancesco D, Ferroni A, Mazzanti M, Tromba C. 1986. Properties of the hyperpolarizing-activated current (if) in cells isolated from the rabbit sino-atrial node. *The Journal of Physiology* **377**:61-88. DOI: https://doi.org/10.1113/jphysiol.1986.sp016177, PMID: 2432247

DiFrancesco D, Tortora P. 1991. Direct activation of cardiac pacemaker channels by intracellular cyclic AMP. *Nature* **351**:145-147. DOI: https://doi.org/10.1038/351145a0, PMID: 1709448

DiFrancesco D. 1993. Pacemaker mechanisms in cardiac tissue. *Annual Review of Physiology* **55**:455-472. DOI: https://doi.org/10.1146/annurev.ph.55.030193.002323, PMID: 7682045

DiFrancesco D. 1999. Dual allosteric modulation of pacemaker (f) channels by cAMP and voltage in rabbit SA node. *The Journal of Physiology* **515**:367-376. DOI: https://doi.org/10.1111/j.1469-7793.1999.367ac.x, PMID: 10050004

Dominguez C, Boelens R, Bonvin AM. 2003. HADDOCK: a protein-protein docking approach based on biochemical or biophysical information. *Journal of the American Chemical Society* **125**:1731-1737. DOI: https://doi.org/10.1021/ja026939x, PMID: 12580598

Dyson HJ, Wright PE. 2004. Unfolded proteins and protein folding studied by NMR. *Chemical Reviews* **104**: 3607-3622. DOI: https://doi.org/10.1021/cr030403s, PMID: 15303830

Emery EC, Young GT, McNaughton PA. 2012. HCN2 ion channels: an emerging role as the pacemakers of pain. *Trends in Pharmacological Sciences* **33**:456-463. DOI: https://doi.org/10.1016/j.tips.2012.04.004, PMID: 22613784

Fosgerau K, Hoffmann T. 2015. Peptide therapeutics: current status and future directions. *Drug Discovery Today* **20**:122-128. DOI: https://doi.org/10.1016/j.drudis.2014.10.003, PMID: 25450771

Garrett DS, Seok YJ, Peterkofsky A, Clore GM, Gronenborn AM. 1997. Identification by NMR of the binding surface for the histidine-containing phosphocarrier protein HPr on the N-terminal domain of enzyme I of the Escherichia coli phosphotransferase system. *Biochemistry* **36**:4393-4398. DOI: https://doi.org/10.1021/bi970221q, PMID: 9109646

Guidotti G, Brambilla L, Rossi D. 2017. Cell-Penetrating peptides: from basic research to clinics. *Trends in Pharmacological Sciences* **38**:406-424. DOI: https://doi.org/10.1016/j.tips.2017.01.003, PMID: 28209404

Güntert P, Buchner L. 2015. Combined automated NOE assignment and structure calculation with CYANA. *Journal of Biomolecular NMR* **62**:453-471. DOI: https://doi.org/10.1007/s10858-015-9924-9, PMID: 25801209

Han Y, Noam Y, Lewis AS, Gallagher JJ, Wadman WJ, Baram TZ, Chetkovich DM. 2011. Trafficking and gating of hyperpolarization-activated cyclic nucleotide-gated channels are regulated by interaction with tetratricopeptide repeat-containing Rab8b-interacting protein (TRIP8b) and cyclic AMP at distinct sites. *Journal of Biological Chemistry* **286**:20823-20834. DOI: https://doi.org/10.1074/jbc.M111.236125, PMID: 21504900

Henninot A, Collins JC, Nuss JM. 2018. The current state of peptide drug discovery: back to the future? *Journal of Medicinal Chemistry* **61**:1382-1414. DOI: https://doi.org/10.1021/acs.jmedchem.7b00318, PMID: 28737935

Herce HD, Garcia AE, Cardoso MC. 2014. Fundamental molecular mechanism for the cellular uptake of guanidinium-rich molecules. *Journal of the American Chemical Society* **136**:17459-17467. DOI: https://doi.org/10.1021/ja507790z, PMID: 25405895

Herrmann S, Rajab H, Christ I, Schirdewahn C, Höfler D, Fischer MJM, Bruno A, Fenske S, Gruner C, Kramer F, Wachsmann T, Wahl-Schott C, Stieber J, Biel M, Ludwig A. 2017. Protein kinase A regulates inflammatory pain sensitization by modulating HCN2 channel activity in nociceptive sensory neurons. *PAIN* **158**:2012-2024. DOI: https://doi.org/10.1097/j.pain.0000000000001005, PMID: 28767511
Hu L, Santoro B, Saponaro A, Liu H, Moroni A, Siegelbaum S. 2013. Binding of the auxiliary subunit TRIP8b to HCN channels shifts the mode of action of cAMP. *The Journal of General Physiology* **142**:599-612.
DOI: https://doi.org/10.1085/jgp.201311013, PMID: 24277603

Hubbard SJ, Hornton JM. 1993. NACCESS. University College London.

Kristensen M, Birch D, Mørck Nielsen H. 2016. Applications and challenges for use of Cell-Penetrating peptides as delivery vectors for peptide and protein cargos. *International Journal of Molecular Sciences* **17**:185.
DOI: https://doi.org/10.3390/ijms17020185

Lee CH, MacKinnon R. 2017. Structures of the Human HCN1 Hyperpolarization-Activated Channel. *Cell* **168**:111-120.
DOI: https://doi.org/10.1016/j.cell.2016.12.023, PMID: 28086084

Lolicato M, Nardini M, Gazzarrini S, Möller S, Bertinetti D, Herberg FW, Bolognesi M, Martin H, Fasolini M, Bertrand JA, Arrigoni C, Thiel G, Moroni A. 2011. Tetramerization dynamics of C-terminal domain underlies isoform-specific cAMP gating in hyperpolarization-activated cyclic nucleotide-gated channels. *Journal of Biological Chemistry* **286**:44811-44820.
DOI: https://doi.org/10.1074/jbc.M111.297606, PMID: 22006928

Lyman KA, Han Y, Heuermann RJ, Cheng X, Kurz JE, Lyman RE, Van Veldhoven PP, Chetkovich DM. 2017. Allostery between two binding sites in the ion channel subunit TRIP8b confers binding specificity to HCN channels. *Journal of Biological Chemistry* **292**:17718-17730.
DOI: https://doi.org/10.1074/jbc.M117.802256, PMID: 28887304

Puljung MC, Zagotta WN. 2013. A secondary structural transition in the C-helix promotes gating of cyclic nucleotide-regulated ion channels. *Journal of Biological Chemistry* **288**:12944-12956.
DOI: https://doi.org/10.1074/jbc.M113.464123, PMID: 23525108

Robinson RB, Siegelbaum SA. 2003. Hyperpolarization-activated cation currents: from molecules to physiological function. *Annual Review of Physiology* **65**:453-480.
DOI: https://doi.org/10.1146/annurev.physiol.65.092101.142734, PMID: 12471170

Santoro B, Hu L, Liu H, Saponaro A, Pian P, Piskorowski RA, Moroni A, Siegelbaum SA. 2011. TRIP8b regulates HCN1 channel trafficking and gating through two distinct C-terminal interaction sites. *Journal of Neuroscience* **31**:4074-4086.
DOI: https://doi.org/10.1523/JNEUROSCI.5707-10.2011, PMID: 21411649

Santoro B, Piskorowski RA, Pian P, Hu L, Liu H, Siegelbaum SA. 2009. TRIP8b splice variants form a family of auxiliary subunits that regulate gating and trafficking of HCN channels in the brain. *Neuron* **62**:802-813.
DOI: https://doi.org/10.1016/j.neuron.2009.05.009, PMID: 19555649

Saponaro A, Pauleta SR, Cantini F, Matzapetakis M, Hammann C, Donadoni C, Hu L, Thiel G, Banci L, Santoro B, Moroni A. 2014. Structural basis for the mutual antagonism of cAMP and TRIP8b in regulating HCN channel function. *PNAS* **111**:14577-14582.
DOI: https://doi.org/10.1073/pnas.1410389111, PMID: 25197093

Wainger BJ, DeGennaro M, Santoro B, Siegelbaum SA, Tibbs GR. 2001. Molecular mechanism of cAMP modulation of HCN pacemaker channels. *Nature* **411**:805-810.
DOI: https://doi.org/10.1038/35081088, PMID: 11459060

Weißgraeber S, Saponaro A, Thiel G, Hamacher K. 2017. A reduced mechanical model for cAMP-modulated gating in HCN channels. *Scientific Reports* **7**:40168.
DOI: https://doi.org/10.1038/srep40168, PMID: 28074902

Zagotta WN, Olivier NB, Black KD, Young EC, Olson R, Gouaux E. 2003. Structural basis for modulation and agonist specificity of HCN pacemaker channels. *Nature* **425**:200-205.
DOI: https://doi.org/10.1038/nature01922, PMID: 12968185

Zhou L, Siegelbaum SA. 2007. Gating of HCN channels by cyclic nucleotides: residue contacts that underlie ligand binding, selectivity, and efficacy. *Structure* **15**:655-670.
DOI: https://doi.org/10.1016/j.str.2007.04.012, PMID: 17562313

Zolles G, Wenzel D, Bildl W, Schulte U, Hofmann A, Müller CS, Thumfart JO, Vlachos A, Deller T, Pfeifer A, Fleischmann BK, Roeper J, Fakler B, Klöcker N. 2009. Association with the auxiliary subunit PEX5R/Trip8b controls responsiveness of HCN channels to cAMP and adrenergic stimulation. *Neuron* **62**:814-825.
DOI: https://doi.org/10.1016/j.neuron.2009.05.008, PMID: 19555650