Article Outline

Article Information

PubMed

Related Articles

  • …more

Copyright © 2007 The Biophysical Society. All rights reserved.
Biophysical Journal, Volume 92, Issue 5, 1559-1572, 1 March 2007

doi:10.1529/biophysj.106.092478

Channels, Receptors, and Electrical Signaling

Dual Stretch Responses of mHCN2 Pacemaker Channels: Accelerated Activation, Accelerated Deactivation

Wei LinUlrike LaitkoPeter F. Juranka and Catherine E. MorrisGo To Corresponding Author 

Neuroscience, Ottawa Health Research Institute, Ottawa Hospital, Ottawa, Ontario, Canada K1Y 4E9

Address reprint requests to Catherine E. Morris, Neuroscience, Ottawa Health Research Institute, Ottawa Hospital, Ottawa, Ontario, Canada K1Y 4E9. Tel.: 613-798-5555 ext. 18608; Fax: 613-761-5330.

Abstract

Mechanoelectric feedback in heart and smooth muscle is thought to depend on diverse channels that afford myocytes a mechanosensitive cation conductance. Voltage-gated channels (e.g., Kv1) are stretch sensitive, but the only voltage-gated channels that are cation permeant, the pacemaker or HCN (hyperpolarization-activated cyclic nucleotide-gated) channels, have not been tested. To assess if HCN channels could contribute to a mechanosensitive cation conductance, we recorded IHCN in cell-attached oocyte patches before, during, and after stretch for a range of voltage protocols. ImHCN2 has voltage-dependent and instantaneous components; only the former was stretch sensitive. Stretch reversibly accelerated hyperpolarization-induced ImHCN2 activation (likewise for IspHCN) and depolarization-induced deactivation. HCN channels (like Kv1 channels) undergo mode-switch transitions that render their activation midpoints voltage history dependent. The result, as seen from sawtooth clamp, is a pronounced hysteresis. During sawtooth clamp, stretch increased current magnitudes and altered the hysteresis pattern consistent with stretch-accelerated activation and deactivation. ImHCN2 responses to step protocols indicated that at least two transitions were mechanosensitive: an unspecified rate-limiting transition along the hyperpolarization-driven path, mode Iclosed→mode IIopen, and depolarization-induced deactivation (from mode Iopen and/or from mode IIopen). How might this affect cardiac rhythmicity? Since hysteresis patterns and “on” and “off” IHCN responses all changed with stretch, predictions are difficult. For an empirical overview, we therefore clamped patches to cyclic action potential waveforms. During the diastolic potential of sinoatrial node cell and Purkinje fiber waveforms, net stretch effects were frequency dependent. Stretch-inhibited (SI) ImHCN2 dominated at low frequencies and stretch-augmented (SA) ImHCN2 was progressively more important as frequency increased. HCN channels might therefore contribute to either SI or SA cation conductances that in turn contribute to stretch arrhythmias and other mechanoelectric feedback phenomena.

Introduction

Mechanical forces affect the electrophysiology of cardiac and smooth muscle cells. Prime candidates for force transduction in mechanoelectric feedback are channels whose open probability changes reversibly with membrane stretch, i.e., mechanosensitive (MS) channels. MS channels include (along with TRP, 2PDK, and others) voltage-gated channels 1. Here, we examine stretch responses of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, the only cation-selective voltage-gated channels 2. HCN tetramers form pacemaker channels and the subunits have a Kv-like arrangement (six transmembrane-spanning domains each with S1–S4 sensor modules and S5–S6 pore modules).

Highlighting the need to test pacemaker channel mechanosensitivity are the following observations: 1), cell swelling reversibly increases ImHCN23; 2), swelling of ventricular myocytes, cells now known to express HCN channels (e.g., 4,5), activates an unidentified inwardly rectifying Ication6; 3), hypothalamic osmosensory neurons have a swelling- and voltage-dependent Cs-blockable Ication7; 4), cholesterol is a bilayer mechanical reagent that modulates voltage-gated channels 8 and altered cholesterol levels in sinoatrial node (SAN) cells alter the kinetics of HCN4-based pacemaker channels 9; 5), inhibitors of another type of MS cation conductance that is a candidate mechanotransducer for mechanoelectric feedback, the MS TRP channels, inhibit neither mouse SAN cell MS Ication10 nor certain ventricular mechanoelectric feedback phenomena 11; 6), repeated attempts to record MS TRP-like events from adult cardiomyocytes have failed 12; and 7), many MS cells and MS cellular structures are abundantly endowed with HCN channels (Table 1).

Information on native cell MS voltage-gated channel responses is sparse, but in smooth muscle and SAN cells, voltage-gated calcium channel current reversibly increases with stretch 1 (in recombinant calcium channels this trait resides with the pore subunit 13). MS modulation of voltage-gated channels is best understood for Kv channels, where voltage sensing, slow inactivation, and a concerted step before pore opening have all been characterized as MS transitions 14,15. Reversible stretch responses in these Kv channels occur in both excised and cell-attached patches 16; and since patch excision destroys membrane skeleton organization, we infer that the bilayer transmits force. Why does Kv gating respond to bilayer stretch? We hypothesize that each conformation of a Kv channel’s extensive protein/bilayer interface 17 has a sufficiently different lateral pressure profile 18 that bilayer mechanical perturbations change the relative stabilities of the conformations. Lateral pressure profiles can be perturbed physically (i.e., bilayer expansion/compression) and/or chemically (e.g., depletion/enrichment of surface-active lipids, alcohols, anesthetics, and drugs). Since the “bilayer mechanics” of membrane stretch and the bilayer mechanics of cholesterol, lysophospholipids, short chain alcohols, phenothiazines, etc. are related forms of lipid stress 19, the fact that HCN channels can be modulated by stretch has wide potential physiological and pathophysiological implications.

Our principal subject here, mHCN2, is known to participate in SAN cell and neuronal rhythmicity through the phenotypes of a knockout mouse 20. We also briefly tested spHCN 21 (from sea urchin) since this is the only isoform for which gating current information has been obtained 22 and the one for which it was shown that HCN and Kv1 voltage sensors respond similarly to voltage 23. In HCN channels, however, it is hyperpolarization, not depolarization, that triggers channel opening. In Kv1 channels, depolarization and stretch facilitate the same step, namely, the independent depolarization-driven activation transition 15. On the face of it, if HCN channels respond like Kv1 channels, then stretch should slow IHCN activation and speed IHCN deactivation. However, that is not the full story. Kv1 channel mutants with identified rate-limiting transitions reveal additional MS transitions: stretch accelerates slow inactivation 15 and it decelerates a concerted voltage-dependent step just before pore opening 14. If comparable transitions occur in HCN channels (not yet established in the case of concerted transitions), they may be similarly MS.

Mode switching, a kinetic trait whose consequences become striking during prolonged voltage excursions, is a feature common to HCN and Kv channels 22. The Kv1 version is termed “C-type inactivation”. Depending on recent voltage history, Kv1 channels are in “active” or “slow-inactivated” mode; “active” has open and closed states, “slow-inactivated”, only closed states. In HCN channels, by contrast, both modes have open and closed states and hence discrete g (V) curves (Figure 1A). Mannikko et al. hypothesize 22 that once a HCN channel opens, a voltage-independent conformation change occurs at the pore/sensor domain interface that stabilizes the voltage sensor in its retracted position and thus gives rise to a depolarizing shift in the Q (V) and activation curves. Mammalian HCN1, 2, and 4 isoforms all exhibit mode-switch behavior which, for all of them, occurs on a 100-ms timescale even though activation kinetics vary by well over an order of magnitude (HCN1 is fastest, HCN4 is slowest) 24. The resulting hysteresis in the HCN channel current/voltage (I/V) relations is most prominent in HCN1 and least in HCN4 25. Prolonged depolarization favors (stabilizes) mode I and prolonged hyperpolarization favors mode II (Figure 1A) with the overall consequence that open HCN channels experience longer openings and closed ones, longer closings.

Display large version of this figure
Figure 1
Modes and a common kinetic framework for HCN and Kv channels. (A) After prolonged sojourns (∞) at hyperpolarized or depolarized voltages, HCN channels principally occupy mode II (gray) or mode I (black) 22, respectively. Activation and gating charge curves (g (V), Q (V)) for the modes are depicted. The vertical arrow depicts the consequence (at a given voltage) of undergoing a mode I to mode II transition. (B) Open (O) and closed (C) states connected in circular schemes 22 for HCN and Kv channels (d, depolarization; h, hyperpolarization). The x through the OI of Kv signifies that in Kv-WT, “OI”, unlike its HCN counterpart, is occluded (i.e., “inactivated-nonconducting”) (see Discussion). The dominant states for typical Vhold potentials (depolarized-HCN; hyperpolarized-Kv) are boxed. Asterisks and inverse asterisks signify candidate stretch-accelerated and stretch-decelerated transitions, respectively, that would represent a consistent pattern across subfamilies. Two additional candidates for stretch-accelerated transitions bear question marks (see Discussion). The stretch-decelerated (1/*) Kv transition 14 is the concerted depolarization-driven closed-closed transition just before the last (voltage-independent) pore opening step 67,68. Whether HCN channels have analogous concerted transitions (shown gray) into the open state(s) is unknown. The asterisks in mode II correspond to a known stretch-accelerated voltage sensor motion for Kv 15; evidence suggesting that the comparable mode I transition is also stretch-accelerated is presented here for mHCN2.

Cardiomyocytes express mixtures of isoforms HCN1, 2, and 4 26. Smooth muscle shows electrophysiological and reverse transcription-polymerase chain reaction evidence of HCN2, 3, and 4 27. In mammalian brain, all four isoforms are expressed, each with a unique pattern of distribution 28; immunogold electronmicroscopy shows that, as well as being expressed in neurons, HCN2 has a widespread expression in oligodendrocytes 29. Native HCN-based channels, with conductances of <10pS, pass Ca2+30 as well as Na+ and K+.

We expressed HCN channels in oocytes, monitoring macroscopic IHCN in cell-attached patches before, during, and after stretch, using a variety of voltage clamp protocols. In Kv1 channels depolarization elicits activation of closed channels, and the process accelerates with stretch 15. In HCN channels, depolarization elicits deactivation of open channels, and we predicted it too would accelerate with stretch since it involves comparable motions of the voltage sensors 23. Thinking of HCN channels as “inverted Kv1 channels” we therefore expected to observe “stretch inhibition” (SI) during hyperpolarization (due to stretch-enhanced deactivation). We were therefore surprised to observe a very unequivocal “stretch augmentation” (SA). Eventually we found that, indeed, deactivation was accelerated with stretch, but this strong effect was initially cryptic, since a), it did not result in a slower activation of IHCN during hyperpolarizations, and b), we did not immediately understand the workings of HCN channels protracted tail currents.

Having unexpectedly found that HCN current activation accelerated with stretch, a second surprise came when we tested currents during rhythmic activity. During slow action potential (AP) waveforms, the SI effects, not the SA effects, were what dominated. More particularly, during the diastolic part of APs, ImHCN2 showed SI at low and SA at high “beat” frequencies. It required the emerging understanding of HCN channel mode switch and the hysteresis that results from mode switch 22 to help clarify the time-dependent stretch responses of the mHCN2 pacemaker channels.


Methods

Oocyte preparation and cRNA injection

Xenopus laevis oocytes were injected with 5–55ng cRNA as described previously 14. Capped cRNA for mHCN2 (mHCN2-pGEM-HE (kindly provided by Steve Siegelbaum); linearized with SphI) 31 and spHCN 21 (SPIH-pGEM-HE (kindly provided by U. B. Kaupp 21); linearized with NheI) was produced by in vitro transcription using Ambion message machine.


Electrophysiology

Manually devitellinated oocytes were transferred to an inverted microscope rig. Thick-walled, sylgard-coated, fire-polished (using a soda glass-coated platinum filament) pipettes were prepared as previously 14 and had resistances of 2–4MΩ. Macroscopic currents measured from cell-attached patches (Axopatch 200B; Axon Instruments, Foster City, CA) at RT (20°C–22°C) were filtered at 5kHz (Axopatch 200B low-pass filter) and analyzed using WinASCD (Guy Droogmans) and Origin (Microcal, Northampton, MA).

Experimental voltage protocols (step, ramp, and waveform protocols) were controlled via pClamp 8, as indicated in Results. Vhold was 0mV unless otherwise noted.

Membranes were stretched reversibly via negative pipette pressure, controlled, and monitored by a DPM-1B transducer (Bio-Tek Instruments, Winooski, VT) or a high-speed pressure clamp (HSPC-1, ALA Scientific Instruments, Westbury, NY). MS responses were typically elicited using −30 or −40mmHg. Because a given pressure generates a different tension in different patches, several levels were sometimes tested. To obtain difference currents, protocols were structured as Before/During/After stretch (B/D/A) sets, with at least 15s between runs within a B/D/A set. Data were discarded when stretch-induced changes were irreversible.

Water-injected (in actuality, Tris buffer-injected) control oocytes were incubated for 2–3 days, patched and clamped using protocols like those used for Figure 2 and Figure 3 and Figure 4 and Figure 5 and Figure 6 and Figure 7 and Figure 8 (i.e., step, ramp, and AP clamp protocols) to make B/D/A sets (using −40mm Hg or more). The same lanthanum (1mM), high-K pipette solution was used for control oocytes as for HCN expressing oocytes. The control records (available online as Supplementary Material ) showed virtually no evidence of time- and voltage-dependent conductances.


Solutions

High-K solution (for bath and pipette) contained (in mM) 89 KCl, 0.4 CaCl2, 5 HEPES, 0.8 MgCl2 (pH 7.5 with KOH). Then 40μM GdCl3 or 1mM LaCl3 was added to pipettes to block the endogenous MS cation channels. Lanthanum (adopted partway through the project) was always fully effective. Gadolinium was sometimes fully effective, but sometimes MS cation channel activity made records unusable, despite adding degassed acidified stock solution (100mM GdCl3) to degassed pipette solution immediately before recording 32. Either gadolinium precipitated imperceptibly or oocytes sporadically express La-sensitive-Gd-insensitive IMScation channels; renal cells 33 reportedly have such channels.



Results

Monitoring recombinant HCN channels against an IMScation background

To illustrate why lanthanides were needed, Figure 2Ai shows a “worst-case” example of interference by endogenous IMScation. A hyperpolarizing ramp was applied before, during, and after stretch (B/D/A). Here and throughout, black/red/gray signify B/D/A. The B and A traces run along the x axis, and IMScation (i.e., the current elicited by stretch, the D trace) corresponds to ∼40 open MScation channels. In Figure 2Aii, as in 2 Ai, −100-mmHg suction was used, but there was no stretch current because 40μM gadolinium in the pipette inhibited IMScation. IHCN families in the presence of lanthanides are shown in Figure 2Biii (mHCN2; 40μM gadolinium and 1mM lanthanum chloride, respectively) and 2 Biii (spHCN; 40μM gadolinium chloride). A hyperpolarizing ramp (Figure 2C) reveals the minor ohmic component reversing at ∼0mV (comprising the instantaneous component 34IHCNinst, nonspecific Ileak, and any imperfectly inhibited IMScation) and IHCN (V,t) for mHCN2 channels, most of which would have been in mode I as they activated at ∼−100mV. The B/D/A ramp traces in Figure 2D show that inhibiting IMScation (as in Figure 2Bii) allowed for extraction of a stretch difference current, Idiff (i.e., D−[(B+A)/2]), shown in pink.

Display large version of this figure
Figure 2
Dealing with the endogenous MS cation channel. (A) Ramp (33mV/s) currents on patches from uninjected oocytes before, during, and after (B/D/A: black, red, gray) stretch without (i) or with (ii) 40μM Gd3+ in the pipette. Typically, stretch was applied using less intense suction than the −100mmHg used here. (B) spHCN1 (i) and ImHCN2 (ii and iii) current families, using 20-mV steps from a Vhold of 0mV (i) or +20mV (ii, iii) to −140mV (i and iii) or −160mV (ii), with 40μM Gd3+ (i, ii) or 1mm La3+ (iii) in the pipette; scales 200pA/0.5s; 100pA/0.25s; and 200pA/0.2s (iiii). (C) A typical inwardly rectifying ramp IHCN (33mV/s) for mHCN2; arrow indicates ramp direction. (D) Ramp currents (as in C) with a B/D/A protocol, using strong suction (−80mmHg); in the Idiff=D−[(B+A)/2], the small instantaneous component gets removed.

Capacitive transients interfered with early IHCN during voltage steps (e.g., Figure 3AA and Figure 4AA). Nevertheless, we used no online linear subtraction since that requires long excursions to mode-switch-inducing voltages and would obliterate any stretch changes in IHCNinst. Traces are either raw current, averaged replicates (2 or 3), or Idiff (i.e., D−[(B+A)/2] from a B/D/A protocol). “Stretch-augmented” (SA) and “stretch-inhibited” (SI) correspond to “inward” and “outward” Idiff, respectively.

Display large version of this figure
Figure 3
Stretch, oocyte patches, mHCN2, and the subcomponents of ImHCN2 40μM Gd3+ present. Vhold=−40mV. (A) Averaged ImHCN2 for two sequential B/D/A sets (D, −30mmHg). Notable are the large instantaneous current component (evident at −120mV above the dashed line) and the effect of stepping back to Vhold. Idiff for the averaged sets is plotted at right. The asterisk is mentioned in the text. (B) Five sequential B/D/A sets (D, as indicated in mmHg). The voltage protocol was to +40mV then to −120mV as in A. The unitary currents are discussed in the text.
Display large version of this figure
Figure 4
spHCN1 and stretch. (A). Three B/D/A sets recorded (same patch) in the sequence shown. The protocol was a simple step from 0mV to −80mV. The consistently smaller amplitude of the after traces (gray) within each B/D/A and the progressive decrease in the amplitudes of before traces (black) from one B/D/A indicates that slow rundown was occurring, but this did not prevent stretch acceleration of activation. (B) Expanded early part of the first B/D/A set (i.e., D=−30mmHg; black trace is (B+A)/2) with vertical red and black bars of equal amplitude. (C) Idiff for these B/D/A sets, with early section expanded. Scale bars (top to bottom); 100pA/150ms/, 30pA/400ms, 30pA/60ms.

A priori, stretch might alter unitary HCN channel conductance. Since smooth muscle 35 cells show large (>50pS) hyperpolarization- and stretch-activated Ication events, sporadic unitary events riding atop macroscopic ImHCN2 during stretch (e.g., Figure 3B) were of interest. However, brief depolarizations showed that such events never inwardly rectified, and since they never occurred with La3+ (see Methods), we take them to be endogenous currents.


Stretch increases IHCN during hyperpolarizing steps

Figure 3A shows that hyperpolarization-induced ImHCN2 was accelerated and augmented by stretch. As is common for recombinant IHCN, sporadic inward noise (e.g., as in red and gray traces) was sometimes seen at large hyperpolarizations but was unrelated to stretch. At Vhold=−40mV there was no IHCN (note I=0pA at the start of traces). After prolonged hyperpolarizations, however, stepping back to −40mV produced a sustained tail current. In other words, the 80-mV drop produced, as expected, an ohmic decrease in IHCN, but this was not followed by decay toward I=0. Channels stayed open. We return later to this manifestation that mode switching had occurred during the hyperpolarized sojourn.

As seen in Figure 3B, even with Gd3+ in the pipette, patches sometimes exhibited inward and outward unitary IMScation events. The events contaminating two of the stretch traces illustrate a stark difference between MS IHCN and endogenous IMScation: even at the few pA level, MS IHCN was “macroscopic”. Moreover, endogenous unitary IMScation was not time dependent (note the unitary events during the latent period for the voltage- and time-dependent component of ImHCN2) (see also Figure 5A (x, bottom)).

Display large version of this figure
Figure 5
Reversible, reproducible stretch changes in ImHCN2. (A) Four sequential B/D/A sets. Vhold=−40mV, followed by the indicated steps. Gadolinium was used. The boxed trace (bottom) expands “run 1, D” near the X, to show unitary endogenous events. For control run, D=0mmHg, for runs 1–3, D=−20mmHg. (B) Averaged traces and Idiff (expanded at right). Asterisk, nearly overlooked stretch-accelerated current decline. (C) For runs 1–3, expanded stretch (averaged D) and no-stretch (averaged (B+A)/2) segments from above show that ImHCN2inst was stretch insensitive and that for ∼500ms, >1/3 of the time-dependent ImHCN2 was due to stretch.

Stretch does not affect the instantaneous component of ImHCN2

Not only did stretch not induce inwardly rectifying high conductance states of the mHCN2 channels, it did not alter ImHCN2inst. This was directly evident from the insensitivity to stretch of ImHCN2inst in records where the voltage- and time-dependent component of ImHCN2 increased with stretch. This lack of effect is apparent by inspection of the latent periods for ImHCN2 in Figure 3B (see also Figure 5C and Figure 6AB) and by inspection of hyperpolarizing ramps without and with stretch in Figure 2D (see also Figure 7BC).

Display large version of this figure
Figure 6
Effect of stretch on ImHCN2 activation and deactivation. (A) B/D/A with early activation (scales, 4pA, 0.1s; red-black bar as in Figure 4A) and deactivation expanded; for ∼150ms after the hyperpolarizing step, ∼1/2 the time-dependent ImHCN2 was SA. The −50mmHg Idiff is shown, bottom right (with its B/D/A for reference) and above, the equivalent for a smaller stretch stimulus. The prolonged (2s) step to −40mV showed that after a latency, stretch accelerated the current decline. The Idiff shapes here clarify those from truncated protocols (Figure 3 and Figure 5). (B and C) B/D/A traces, as explained in the text, plus an Idiff for which an asterisk (Ication>0pA) indicates SI ImHCN2. (D) expanded Itail at −40mV for this series of B/D/A sets (labels identify the preceding hyperpolarization level). Vertical scales from top to bottom: 3pA, 2pA 5pA, 5pA; horizontal scale: 0.4s. Vhold=−40mV. Lanthanum was used here and in subsequent figures.
Display large version of this figure
Figure 7
Hysteresis and sawtooth ramp ImHCN2 (A) Top, a hyperpolarizing then depolarizing ramp set (with a 50-ms pause at −150mV) from +80mV to −150mV. A heavily filtered segment from this ramp set and from two others at higher speeds (same patch) is shown below (labels give the duration of each limb). (B) A B/D/A ramp set (D, −40mmHg, ramps as in A, 7.0s) and its Idiff. (C) i. The fast-then-slow ramp protocol for the B/D/A set (D=−40mmHg) seen in ii with fast (f) and slow (s) sets plotted separately alongside their respective Idiff, for which scales are 4pA (fast) and 10pA (slow). In iii, relative effects of stretch on rates in the two limbs are compared graphically for the slow B/D/A by normalizing (see arrows) the average no-stretch ((B+A)/2) and the stretch (D) I/V relations.

Activation of spHCN is stretch sensitive

In B/D/A sets, IspHCN (Figure 4A) exhibited stretch-accelerated activation (arrows). As a semiquantitative measure of the effect of stretch (using −30mmHg) note that, for more than 150ms, more than half of the time- and voltage-dependent component of IspHCN comprised SA IHCN (Figure 4B). Idiff peaked at ∼150ms for all three stretch intensities tested (Figure 4C). Qualitatively, the traces show unequivocal reversible acceleration of activation by stretch, but as the expanded Idiff sections show, the early difference currents were too nonuniform to allow for standardized curve fitting and assignment of rates (as was true, too, for mHCN2 activation currents).


Stretch-augmented and stretch-inhibited ImHCN2: current activation and tail currents

Figure 5A shows ImHCN2 for an unusually stable patch: first a control set, then three B/D/A sets (with D=−20mmHg; these are averaged in Figure 5B). (We did not attempt to study tail currents at positive depolarized voltages; the brief steps to +40mV were internal controls to ensure that endogenous MS channels were inhibited.) Stretch-accelerated activation resulted in an elevated steady-state ImHCN2. On stepping to −40mV, the magnitude of tail Idiff (pink arrow) reflects the stretch augmentation of current that had occurred at −110mV. The persistence of Idiff at −40mV (which, note was also Vhold) suggests that mode switch had occurred during the sojourn at −110mV (at hyperpolarized voltages, the channels prefer mode II with its right-shifted g (V)). Stretch evidently reversibly affected a rate-limiting transition along the pathway mode Iclosed→mode IIopen (i.e., the path that is suggested by the vertical arrow in Figure 1A). Longer tail currents than obtained here were needed to clarify the shape and polarities of Idiff (Figure 5B*).

At −40mV, tail currents required ∼2s for completion, as seen in Figure 6A. At right, the Idiff for the D=−50mmHg B/D/A set is shown below that for the D=−30mmHg set (acquired before the −50mmHg set). Tail currents exhibited an open-channel latency (∼100ms) after which stretch reversibly accelerated the decay of IHCN to zero. Note that tail Idiff (pink) was initially stretch augmented, but then crossed the zero current line and became stretch inhibited (likewise, Figure 6C*). The inflection (after ∼0.5s) in the Figure 6A tail Idiff simply indicates the approach to steady state (likewise for spHCN, Fig. 4, for inflections in the rising phase Idiff and for the rising phase Idiff here at −140mV).

Stretch acceleration of ImHCN2 deactivation was robust. In the four B/D/A sets of Figure 6BD, the hyperpolarizing prepulse was successively deepened (−80, −100, −120, −140mV). At the onset of the −40mV tail currents, the post −140mV case (Figure 6D, bottom) should have the greatest fraction of channels in mode IIopen. Whether or not the latent period, which is clearly evident as the flat portion in the Figure 6C tail Idiff, was stretch sensitive is uncertain. The shorter latency here at −50mm Hg than at −30mmHg in Figure 6A (see the Idiff traces) suggests that it might be. The latency presumably included “electrically silent” mode IIopen→mode Iopen transitions and (as per the expanded version of the Mannikko et al. model 22) independent S4 motions between the (multiple) open states of either mode. The deactivation phase, which would include mode Iopen→mode Iclosed and mode IIopen→mode IIclosed transitions, was, however, unequivocally accelerated by stretch. This is the response depicted by the asterisk in Figure 1B*, HCN.

Thus, stretch accelerated ImHCN2 activation during hyperpolarization and ImHCN2 deactivation during depolarization. In mHCN2, therefore, at least two types of transition (with opposing effects on current) were susceptible to the bilayer mechanical perturbation 36 produced by membrane stretch.

Our data showed that the “sign” (SA versus SI) of Idiff at a given voltage (and time) depended on voltage history. We emphasize this because it would be crucial in designing B/D/A stretch tests for any native Ipacemaker and in interpreting the outcomes. Consider Figure 6D, where a 0.7-s depolarization to −40mV had four different “histories”. During that 0.7s, Idiff was mostly SA post −80mV but mostly SI post −140mV. With mode switching plus stretch acceleration (on both the activation and deactivation pathways) taken into account, this makes sense. Where there is net SA during hyperpolarization, the Idiff must start out SA, but given time they all become SI current (post −80mV, 0.7s was clearly not sufficient time). Except for being prematurely truncated, the Figure 6D tail currents behave as in Figure 6A (likewise for Figure 3A* and Fig. 5*). Consider, in Figure 6A, the small SI Idiff at −80mV versus, in Figure 6B, the (eventual) SA Idiff at −80mV. This too is perfectly understandable in terms of mode switch plus dual stretch effects but would be perfectly confusing if voltage history was ignored. Specifically, after a sojourn (>1s) at −40mV, mode Iclosed would dominate, but (as the SI traces for Idiff at −30mmHg and −50mmHg Figure 6A reveal) stretch-accelerated deactivation caused further closure. Consequently, during stretch, a step to −80mV produced what was initially an SI Idiff. Given sufficient time at −80mV, as seen in Figure 6B, activation (presumably into mode IIopen far more than into mode Iopen) began and stretch-accelerated activation began to overwhelm the effect of stretch on deactivation, thus producing the net SA effect.


Voltage hysteresis: dual effects of stretch during sawtooth ramps

HCN voltage hysteresis 25 during sawtooth ramp clamp is typified by right-shifted, more S-shaped I/V relations in depolarizing limbs 22. For HCN channels, unlike Kv channels, moderate slowing of ramp speeds fails to make hyperpolarizing and depolarizing limbs converge 22. Though hysteresis is less pronounced in mHCN2 than in mHCN1 24,25, mHCN2 (with 1mM La in the pipette) nevertheless showed the characteristic pattern noted for mHCN1 (i.e., increased rather than decreased hysteresis as ramps are made moderately slower 22), as seen in Figure 7A. The hysteresis persisted during stretch (Figure 7B). For the first (hyperpolarizing) limb, applying stretch was essentially like running a slower ramp; ImHCN2 increased. For the (second) depolarizing limb, stretch and ramp speed (= voltage history) interacted strongly (Figure 7Ciii). At first glance, only SA current is apparent, but that reflects the voltage history of the ramp protocol: the sawtooth ramps, like typical HCN step protocols, started from depolarized voltages. Figure 7Ciii shows a normalized plot of stretch and no-stretch traces for the slower (i.e., closer to equilibrium) sawtooth pair. This plot demonstrates, in relative terms, clear evidence of stretch-accelerated deactivation of ImHCN2 (i.e., clear evidence of “SI behavior”) in the depolarizing limb. In some patches (not shown) the effect of stretch on deactivation was sufficient to yield outward difference current in the depolarizing limb. Thus, step and sawtooth protocols yielded a consistent pattern of SA plus SI effects.


Dual effects of stretch during AP waveforms

Rhythmic cardiac AP represent even more complex voltage histories than sawtooth ramps. Lacking appropriate kinetic data for cardiac modeling, we looked empirically at MS ImHCN2 during APs. ImHCN2 before, during, and after stretch was compared from patches clamped to cyclic AP waveforms. This was complicated by the fact that the amplitude of ImHCN2 changed systematically with AP frequency (it was smaller at lower frequencies) and by the fact that during low frequency runs, stretch often decreased ImHCN2 irreversibly. Nevertheless, when a given patch was successfully tested (i.e., when reversible stretch-induced changes were obtained) at more than one frequency, the outcome for SAN cell-like AP excursions was as is seen in Figure 8A: during a train of APs a small amount of SA occurred at high frequency (e.g., 533bpm (beats per minute)), whereas a larger SI occurred at low frequency (e.g., 15bpm). Figure 8BDF shows data for a single patch at three frequencies: 60bpm showed SI and 240bpm and 480bpm showed SA. Figure 8C shows higher resolution excerpts, and Figure 8D shows their Idiff (overlying the voltage waveform). Figure 8E (different patch) shows three 60bpm trials with the AP hyperpolarizing to three different extents. Reversible SI occurred in each case, with the effect accumulating over several APs then, poststretch, decaying over several APs.

Display large version of this figure
Figure 8
Rhythmic AP currents before during and after stretch. The data in A, E,and G, and B, C, D, and F are from four different patches clamped to a rhythmic SAN cell-like AP waveform (AF) or to a Purkinje fiber-like AP waveform (G). All the figures except F (which has I/V relations for AP currents before and during stretch at the indicated beat frequencies) and D (in which stretch difference currents overlie the voltage waveform) are raw current traces. In B, the vertical arrows locate the traces shown in C. Vertical scales (in B) included for the bottom segment of each of the three current traces are 20, 10, and 10pA for 60, 240, and 480bpm, respectively. In A and C, before and after traces overlap. Other details are given in the text.

The I/V relations (B/D/A; B+A averaged) for APs (e.g., Figure 8F) revealed that with the membrane potential driven by AP waveforms, mHCN2 channels furnished the membrane with a cation conductance whose degree of inward rectification was far weaker than expected from g (V) relations implied by voltage step families or by hyperpolarizing ramp clamp (such as in Figure 2BC). Well-behaved pacemaker current is generally seen as being “on” during the diastolic potential of an AP and “off” during most of the spike 37; indeed, HCN1 behaves in this way 25. However, the outcome we noted—namely that, with or without stretch, AP waveforms keep mHCN2 channels disequilibrated at relatively high open probability throughout the AP regime—accords with recent observations of Azene et al. (25; see their Figure 4B) for HCN2 and HCN4.

The SAN waveforms were nearly sinusoidal, so we checked for SI/SA effects in a more asymmetrical waveform. Figure 8G shows ImHCN2 for a patch driven to a Purkinje fiber type AP (−90 to +5mV) at 120 then at 480bpm (inset). In both cases the net effect of stretch was SI behavior; at 240bpm (not shown) stretch produced an intermediate effect in terms of the extent and rate of SI. For this waveform, reversible stretch effects were obtained for six patches tested at two or more frequencies (30–480bpm). SI always dominated at the lower frequencies, but the frequency at which the SI/SA balance tipped had a wide range. One patch, e.g., showed SI at 60bpm and a mixed response at 120bpm (initial SA subsiding to SI over the nine-beat stimulus).

These findings serve as a “proof of principle” demonstration that during rhythmic activity, a given population of HCN channels can mediate both an SI and an SA cation conductance. The kinetic particulars of homotetrameric channels operating with elevated extracellular K+ at room temperature in an oocyte bilayer without auxiliary subunits are not, of course, those of native channels. For native channels in situ, therefore, the SI/SA balance might play out very differently. In diverse situations native pacemaker channels may feel bilayer stretch during voltage excursions. These include the following: the SAN during excessive filling of the right atrium, hypothalamic osmosensory neurons during lowered blood osmolarity, smooth muscle cells in distended gut, blood vessels, and bladder. What our data indicate is that under these conditions, the SI/SA balance of the pacemaker channels can be expected to vary with the waveform of the cell in question and with the duration of the stretch stimulus.



Discussion

The MS responses of HCN channels

Using recombinant homotetrameric HCN channels expressed in oocytes at levels that yielded macroscopic current, we measured patch current before, during, and after membrane stretch, with pipette suction used to apply stretch. We observed no effect of stretch on the amplitude of the instantaneous component of mHCN2 current. Overall, it appeared that, as in Kv channels 15, stretch altered HCN channel gating kinetics without inducing any novel states. In Kv1 channels, the independent depolarization-driven motions of the voltage sensor accelerate with stretch 15, so for HCN channels we thought that stretch-accelerated HCN deactivation might impede hyperpolarization-induced activation of IHCN. But for neither of two distantly related HCN channels was that true. In both a sea urchin and a mouse isoform, stretch accelerated the activation of IHCN. Further characterization was done with the mouse isoform, mHCN2. During voltage steps, both ImHCN2 activation (hyperpolarization-induced) and deactivation (depolarization-induced) accelerated reversibly with stretch. Sawtooth ramp clamp (hyperpolarize-then-depolarize) revealed a voltage hysteresis that persisted with stretch.

To predict stretch responses of ImHCN2 during APs or to model them from our qualitative information about MS responses was not feasible. Instead, we clamped patches to rhythmic AP waveforms then applied stretch. Like others 25, we found that rhythmic AP waveforms kept HCN channels disequilibrated at a high open probability; operational I/V relations (ImHCN2 versus Vm during rhythmic AP clamp, e.g., +5mV to −100mV at 1Hz) reflected open channel permeation at least as much as voltage-dependent gating. These I/V curves contrast starkly with those from hyperpolarizing ramp clamp (i.e., passing through 0mV and heading to −150mV), which elicited no current before ∼−100mV. During APs, given patches exhibited net SI and net SA responses, depending on waveform; but overall, at lower frequencies, it was SI that dominated.

SI dominance in APs might seem surprising given the strong showing of SA effects with step and sawtooth ramp clamp. However, in those protocols, typical depolarized Vhold levels created a bias for the MS responses of hyperpolarization-induced transitions, which, it turned out, were SA effects. Cardiac excitability machinery (mimicked here by AP clamp) keeps Vm more hyperpolarized, so our quasiphysiological waveforms had the opposite bias from that of our steps and ramps. During slow AP depolarizations (diastole), deactivation of mHCN2 channels comes increasingly into play. If, as we suspect, stretch accelerates deactivation transitions of HCN (both modes), domination by SI during slow APs seems reasonable—likewise for slow sawtooths, where SI was apparent but only after depolarization was underway. Substantial delays precede deactivation of mHCN2, perhaps explaining why SI effects accumulated over the course of several APs (and decayed in the same fashion upon stimulus release). With higher frequency APs, the interplay of dV/dt and channel kinetics evidently made mHCN2 channels hover closer to their open states. The characteristic abruptness of the net SA effects seen at high frequencies (fast onset and offset with stretch) was striking. The next section discusses which channel motions might be candidates for MS transitions.


Possible MS transitions

Generically speaking, for a two-state channel, stretch-accelerated activation of current could arise from reduced stability of the closed state (accelerated opening transition) and stretch-accelerated deactivation could arise from reduced stability of an open state (accelerating closing transition). If stretch stabilized the transition state (like a catalyst, lowering the barrier between the two states), both rates would accelerate. Practically speaking, in ImHCN2, neither activation nor deactivation is a kinetically simple (e.g., first order) process, so we could not attribute MS responses to specific transitions. As a proxy, since HCN and Kv channels have similar voltage-induced gating motions, we consider identified MS transitions of Shaker Kv channels (Figure 1B). With both protein structure and gating motions being largely conserved, the forces that structure the protein-bilayer interface are likely to be broadly conserved. The unified model of Kv and HCN mode switch 22 provides an elegant framework for comparing MS responses of HCN and Kv channels.

In Figure 1B, in keeping with that framework, major HCN and Kv transitions are shown in circular four-state kinetic schemes. The two modes are labeled; mode I corresponds, in Kv, to the slow inactivated mode. The Mannikko et al. scheme 22 highlights 4 states (drawn from a 20-state model), but here, concerted transitions (two short arrows in mode II) are also suggested. They are depicted in mode II only, since it is for mode II in Kv that such transitions have been identified 14. For typical voltage clamp experiments, Vhold would be ∼−100mV and ∼0mV for Kv and HCN, respectively, keeping channels closed in mode II and mode I, respectively, as indicated. The existence of a Kv mutant with an “inactivated”-conducting state 38 (“inactivated” in that prolonged depolarization left shifts its Q (V) curve) helps validate and explain the common Kv/HCN framework.

Vertical arrows signify mode-switch transitions and horizontal ones (except the short open arrow), charge-displacing voltage sensor motions. Based on the left-shifted Q (V), Mannikko et al. 22 suggest HCN mode-switch motions are molecularly related to Kv slow inactivation motions. Kv activation and Kv slow inactivation accelerate to the same extent with stretch, as if the two processes have near identical S4/bilayer interactions and thus experience the same bilayer mechanical perturbation during stretch 15. Based on simultaneous gating (and/or ionic) currents and residue-specific optical measurements as well as more extensive optical scans 39,40,41, Kv slow inactivation (including the earliest P-inactivation step) is thought to involve lateral S4 motions, that is, a net motion in the plane of the bilayer. This makes the sensitivity of Kv slow inactivation to bilayer stretch 15 eminently reasonable. The voltage-independent stretch-sensitive motion presumably includes the initial (P type) slow-inactivation process (i.e., the mode IIopen→ mode Iinactivated switch). A question mark on the relevant mode switch transition in Figure 1B suggests that stretch acceleration here might explain stretch-accelerated slow inactivation in Kv channels. That picture would be even more appealing, however, if it could explain our HCN channel observations in a straightforward way, but it cannot. To explain the stretch-accelerated activation seen in HCN channels, stretch would have to accelerate opposite-going transitions, as indicated by the question mark at CI→CII for HCN. Given the proclivity of HCN channels for “inverted gating”, however, this is a possibility.

For HCN, stretch acceleration of a rate-limiting forward transition in the path modeIclosed→modeIIopen (up arrow, Figure 1A) would accelerate IHCN activation as would stretch deceleration of a fast back-reaction. These possibilities are not mutually exclusive. As explained in the Figure 1B legend, asterisks indicate the class of voltage-dependent transitions identified in Kv as stretch-and-depolarization accelerated 15. HCN tail current data showed that, after a latency, stretch speeded up the depolarization-induced deactivation of ImHCN2. The simplest explanation for this, namely that deactivation of HCN channels (both modes) is a stretch and depolarization accelerated process, would be consistent with the Kv activation picture. The inverse asterisk in Figure 1B marks a voltage-dependent transition that, in Kv, decelerates with stretch 14; this last depolarization-driven motion before pore opening is a concerted step. In a Kv mutant whose concerted opening steps happen to be rate limiting (Shaker ILT), this stretch-decelerated transition yields SI activation 14, but it is not known if HCN channels open via comparable (albeit hyperpolarization-driven and thermal) concerted motions (depicted in gray in Figure 1B).


Physiological and pharmacological implications of MS transitions in HCN channels

Cardiac mechanoelectric feedback phenomena are thought to require an IMScation10,12,42. TRP-based channels are prime candidates in neurons, as they are in osmosensory neurons 43, another rhythmically active cell, but based on our observations, it will also be important to ask if the HCN-based cation channels in these cells transduce physiologically relevant force signals. During stretch and swelling, guinea pig (low intrinsic beat rate) and mouse (high intrinsic beat rate) SAN preparations make opposite responses, the former showing increased and the latter decreased rhythmicity 10. Our core finding is that HCN channels can generate both SA and SI cation current. Thus native pacemaker channels are now candidate contributors to these and similar responses. Since the SI/SA balance in our experiments produced net SI behavior at slower beat rates and net SA behavior at faster ones, it will be important to test a variety of native HCN channels in different mammalian models subjected to swelling and stretching stimuli before speculating further. The recombinant channels, of course, lacked auxiliary proteins that are likely present in situ and may tune the channels differently 44. Commenting on work by Zicha et al. 45, Giles 46 points out that during atrial flutter, the slow rate of native HCN channel deactivation manifests as what could be considered a steady-state current (e.g., as in Figure 8F). If, in such circumstances, the slowness of deactivation is rate limiting, stretch-accelerated deactivation might have an appreciable impact.

In addition to osmosensory neurons (which exhibit an osmosensory SI cation conductance; see references in Liu et al. 7 and Naeini et al. 43) HCN channels are abundant in (see Table 1) baroreceptor nerve terminals, various primary afferent mechanosensory neurons, hair cells, and sperm (where, we note, spHCN expression is restricted to the high-curvature, mechanically dynamic flagellar membrane). In these locations, the mechanosensitivity of HCN channels might provide a tunable element for mechanotransduction and/or electromotility.

With two distantly related HCN channels now included among the MS voltage-gated channels, mechanosensitivity can be regarded as a family-wide trait for voltage-gated channels. The gating of all voltage-gated channels types 8,19, including HCN channels (Table 2), is susceptible to many bilayer mechanical 18 reagents. It is likely that the susceptibilities of voltage-gated channels to bilayer stretch and to such reagents are both responses to lipid stress perturbations, one being physical, the other, physicochemical. Classic models of S4-sensor motions kept the sensor sequestered in the channel protein, away from lipids, but now, some S4 residues are known to contact bilayer lipids 17,47. The stretch and voltage dependencies of specific Kv channel transitions 14,15 constitute proof that physically induced lateral pressure profile changes alter gating rates. We suggest that the MS kinetics of mHCN2 and Kv channels follow a common pattern (Figure 1B) because broad conservation of structure in voltage-gated channels dictates robust conservation of bilayer mechanics. If so, modulation of voltage-gated channels by broad-spectrum lipophilic reagents may involve the same voltage sensor and gate transitions that are perturbed by membrane stretch. If this holds, much pharmacological perplexity evaporates 14,19,48. There would be no need to invoke halothane, ethanol, propofol, cholesterol, and capsazepine “binding sites” in HCN channels to explain kinetic modulation by these bilayer mechanical reagents (see Table 2) any more than we need to invoke “stretch-activation and stretch-inhibition gating motifs” in HCN channels to explain their repertoire of stretch responses.



Acknowledgments

This work was supported by a grant from the Canadian Institutes of Health Research.

References

1. Morris, C., and Laitko, U. (2005). The mechanosensitivity of voltage-gated channels may contribute to cardiac mechano-electric feedback. (Philadelphia: W. B. Saunders). PubMed

2. DiFrancesco, D. (2006). Serious workings of the funny current. Prog. Biophys. Mol. Biol. 90, 13–25. CrossRef | PubMed

3. Calloe, K., Elmedyb, P., Olesen, S.P., Jorgensen, N.K., and Grunnet, M. (2005). Hypoosmotic cell swelling as a novel mechanism for modulation of cloned HCN2 channels. Biophys. J. 89, 2159–2169. Abstract | Full Text | PDF (294 kb) | CrossRef | PubMed

4. Arinsburg, S.S., Cohen, I.S., and Yu, H.G. (2006). Constitutively active Src tyrosine kinase changes gating of HCN4 channels through direct binding to the channel proteins. J. Cardiovasc. Pharmacol. 47, 578–586. PubMed

5. Qu, J., Altomare, C., Bucchi, A., DiFrancesco, D., and Robinson, R.B. (2002). Functional comparison of HCN isoforms expressed in ventricular and HEK 293 cells. Pflugers Arch. 444, 597–601. CrossRef | PubMed

6. Suchyna, T.M., Johnson, J.H., Hamer, K., Leykam, J.F., Gage, D.A., Clemo, H.F., Baumgarten, C.M., and Sachs, F. (2000). Identification of a peptide toxin from Grammostola spatulata spider venom that blocks cation-selective stretch-activated channels. J. Gen. Physiol. 115, 583–598. CrossRef | PubMed

7. Liu, X.H., Zhang, W., and Fisher, T.E. (2005). A novel osmosensitive voltage gated cation current in rat supraoptic neurones. J. Physiol. 568, 61–68. CrossRef | PubMed

8. Lundbaek, J.A., Birn, P., Tape, S.E., Toombes, G.E., Sogaard, R., Koeppe 2nd, R.E., Gruner, S.M., Hansen, A.J., and Andersen, O.S. (2005). Capsaicin regulates voltage-dependent sodium channels by altering lipid bilayer elasticity. Mol. Pharmacol. 68, 680–689. PubMed

9. Barbuti, A., Gravante, B., Riolfo, M., Milanesi, R., Terragni, B., and DiFrancesco, D. (2004). Localization of pacemaker channels in lipid rafts regulates channel kinetics. Circ. Res. 94, 1325–1331. CrossRef | PubMed

10. Cooper, P.J., and Kohl, P. (2005). Species- and preparation-dependence of stretch effects on sino-atrial node pacemaking. Ann. N. Y. Acad. Sci. 1047, 324–335. CrossRef | PubMed

11. Sung, D., Mills, R.W., Schettler, J., Narayan, S.M., Omens, J.H., and McCulloch, A.D. (2003). Ventricular filling slows epicardial conduction and increases action potential duration in an optical mapping study of the isolated rabbit heart. J. Cardiovasc. Electrophysiol. 14, 739–749. PubMed

12. Sachs, F. (2005). Stretch-activated channels in the heart. Cardiac Mechano-Electric Feedback and Arrhythmias: From Pipette to Patient. (Philadelphia: W. B. Saunders). PubMed

13. Calabrese, B., Tabarean, I.V., Juranka, P., and Morris, C.E. (2002). Mechanosensitivity of N-type calcium channel currents. Biophys. J. 83, 2560–2574. Abstract | Full Text | PDF (306 kb) | PubMed

14. Laitko, U., Juranka, P.F., and Morris, C.E. (2006). Membrane stretch slows the concerted step prior to opening in a Kv channel. J. Gen. Physiol. 127, 687–701. CrossRef | PubMed

15. Laitko, U., and Morris, C.E. (2004). Membrane tension accelerates rate-limiting voltage-dependent activation and slow inactivation steps in a Shaker channel. J. Gen. Physiol. 123, 135–154. CrossRef | PubMed

16. Gu, C.X., Juranka, P.F., and Morris, C.E. (2001). Stretch-activation and stretch-inactivation of Shaker-IR, a voltage-gated K+ channel. Biophys. J. 80, 2678–2693. Abstract | Full Text | PDF (366 kb) | PubMed

17. Long, S.B., Campbell, E.B., and Mackinnon, R. (2005). Voltage sensor of Kv1.2: structural basis of electromechanical coupling. Science 309, 903–908. CrossRef | PubMed

18. Cantor, R.S. (2002). Size distribution of barrel-stave aggregates of membrane peptides: influence of the bilayer lateral pressure profile. Biophys. J. 82, 2520–2525. Abstract | Full Text | PDF (123 kb) | PubMed

19. Morris, C.E., and Juranka, P.F. (2007). Lipid stress at play: mechanosensitivity of voltage-gated channels. Current Topics in Membranes 59, 297–338. PubMed

20. Ludwig, A., Budde, T., Stieber, J., Moosmang, S., Wahl, C., Holthoff, K., Langebartels, A., Wotjak, C., Munsch, T., Zong, X., et al. (2003). Absence epilepsy and sinus dysrhythmia in mice lacking the pacemaker channel HCN2. EMBO J. 22, 216–224. CrossRef | PubMed

21. Gauss, R., Seifert, R., and Kaupp, U.B. (1998). Molecular identification of a hyperpolarization-activated channel in sea urchin sperm. Nature 393, 583–587. CrossRef | PubMed

22. Mannikko, R., Pandey, S., Larsson, H.P., and Elinder, F. (2005). Hysteresis in the voltage dependence of HCN channels: conversion between two modes affects pacemaker properties. J. Gen. Physiol. 125, 305–326. CrossRef | PubMed

23. Mannikko, R., Elinder, F., and Larsson, H.P. (2002). Voltage-sensing mechanism is conserved among ion channels gated by opposite voltages. Nature 419, 837–841. CrossRef | PubMed

24. Elinder, F., Mannikko, R., Pandey, S., and Larsson, H.P. (2006). Mode shifts in the voltage gating of the mouse and human HCN2 and HCN4 channels. J. Physiol. 575, 417–431. CrossRef | PubMed

25. Azene, E.M., Xue, T., Marban, E., Tomaselli, G.F., and Li, R.A. (2005). Non-equilibrium behavior of HCN channels: insights into the role of HCN channels in native and engineered pacemakers. Cardiovasc. Res. 67, 263–273. CrossRef | PubMed

26. Baruscotti, M., Bucchi, A., and Difrancesco, D. (2005). Physiology and pharmacology of the cardiac pacemaker (“funny”) current. Pharmacol. Ther. 107, 59–79. CrossRef | PubMed

27. Greenwood, I.A., and Prestwich, S.A. (2002). Characteristics of hyperpolarization-activated cation currents in portal vein smooth muscle cells. Am. J. Physiol. Cell Physiol. 282, C744–C753. PubMed

28. Monteggia, L.M., Eisch, A.J., Tang, M.D., Kaczmarek, L.K., and Nestler, E.J. (2000). Cloning and localization of the hyperpolarization-activated cyclic nucleotide-gated channel family in rat brain. Brain Res. Mol. Brain Res. 81, 129–139. CrossRef | PubMed

29. Notomi, T., and Shigemoto, R. (2004). Immunohistochemical localization of Ih channel subunits, HCN1-4, in the rat brain. J. Comp. Neurol. 471, 241–276. CrossRef | PubMed

30. Yu, X., Duan, K.L., Shang, C.F., Yu, H.G., and Zhou, Z. (2004). Calcium influx through hyperpolarization-activated cation channels (I (h) channels) contributes to activity-evoked neuronal secretion. Proc. Natl. Acad. Sci. USA 101, 1051–1056. CrossRef | PubMed

31. Chen, S., Wang, J., and Siegelbaum, S.A. (2001). Properties of hyperpolarization-activated pacemaker current defined by coassembly of HCN1 and HCN2 subunits and basal modulation by cyclic nucleotide. J. Gen. Physiol. 117, 491–504. CrossRef | PubMed

32. Morris, C.E., Juranka, P.F., Lin, W., Morris, T.J., and Laitko, U. (2006). Studying the mechanosensitivity of voltage-gated channels using oocyte patches. Methods Mol. Biol. 322, 315–329. PubMed

33. Hurwitz, C.G., Hu, V.Y., and Segal, A.S. (2002). A mechanogated nonselective cation channel in proximal tubule that is ATP sensitive. Am. J. Physiol. Renal Physiol.