We explored the structural basis of voltage sensing in HCN1 hyperpolarization-activated channels by examining the relative orientation of the voltage-sensor and pore domains. The opening of channels engineered to contain single cysteine residues at the extracellular ends of the voltage-sensing S4 (V246C) and pore-forming S5 (C303) domains is inhibited by formation of disulfide or cysteine:Cd2+ bonds. As Cd2+ coordination is promoted by depolarization, the S4-S5 interaction occurs preferentially in the closed state. The failure of oxidation to catalyze dimer formation, as assayed by Western blotting, indicates the V246C:C303 interaction occurs within a subunit. Intriguingly, a similar interaction has been observed in depolarization-activated Shaker Kv channels at depolarized potentials but such an intrasubunit interaction is inconsistent with the X-ray crystal structure of Kv1.2, wherein S4 approaches S5 of an adjacent subunit. These findings suggest channels of opposite voltage-sensing polarity adopt a conserved S4-S5 orientation in the depolarized state that is distinct from that trapped upon crystallization.

The positively charged S4 transmembrane segment of voltage-gated channels is thought to function as the voltage sensor by moving charge through the membrane electric field in response to depolarization. Here we studied S4 movements in the mammalian HCN pacemaker channels. Unlike most voltage-gated channel family members that are activated by depolarization, HCN channels are activated by hyperpolarization. We determined the reactivity of the charged sulfhydryl-modifying reagent, MTSET, with substituted cysteine (Cys) residues along the HCN1 S4 segment. Using an HCN1 channel engineered to be MTS resistant except for the chosen S4 Cys substitution, we determined the reactivity of 12 S4 residues to external or internal MTSET application in either the closed or open state of the channel. Cys substitutions in the NH2-terminal half of S4 only reacted with external MTSET; the rates of reactivity were rapid, regardless of whether the channel was open or closed. In contrast, Cys substitutions in the COOH-terminal half of S4 selectively reacted with internal MTSET when the channel was open. In the open state, the boundary between externally and internally accessible residues was remarkably narrow (∼3 residues). This suggests that S4 lies in a water-filled gating canal with a very narrow barrier between the external and internal solutions, similar to depolarization-gated channels. However, the pattern of reactivity is incompatible with either classical gating models, which postulate a large translational or rotational movement of S4 within a gating canal, or with a recent model in which S4 forms a peripheral voltage-sensing paddle (with S3b) that moves within the lipid bilayer (the KvAP model). Rather, we suggest that voltage sensing is due to a rearrangement in transmembrane segments surrounding S4, leading to a collapse of an internal gating canal upon channel closure that alters the shape of the membrane field around a relatively static S4 segment.

Previous expression quantitative trait loci (eQTL) studies have performed genetic association studies for gene expression, but most of these studies examined lymphoblastoid cell lines from non-diseased individuals. We examined the genetics of gene expression in a relevant disease tissue from chronic obstructive pulmonary disease (COPD) patients to identify functional effects of known susceptibility genes and to find novel disease genes. By combining gene expression profiling on induced sputum samples from 131 COPD cases from the ECLIPSE Study with genomewide single nucleotide polymorphism (SNP) data, we found 4315 significant cis-eQTL SNP-probe set associations (3309 unique SNPs). The 3309 SNPs were tested for association with COPD in a genomewide association study (GWAS) dataset, which included 2940 COPD cases and 1380 controls. Adjusting for 3309 tests (p<1.5e-5), the two SNPs which were significantly associated with COPD were located in two separate genes in a known COPD locus on chromosome 15: CHRNA5 and IREB2. Detailed analysis of chromosome 15 demonstrated additional eQTLs for IREB2 mapping to that gene. eQTL SNPs for CHRNA5 mapped to multiple linkage disequilibrium (LD) bins. The eQTLs for IREB2 and CHRNA5 were not in LD. Seventy-four additional eQTL SNPs were associated with COPD at p<0.01. These were genotyped in two COPD populations, finding replicated associations with a SNP in PSORS1C1, in the HLA-C region on chromosome 6. Integrative analysis of GWAS and gene expression data from relevant tissue from diseased subjects has located potential functional variants in two known COPD genes and has identified a novel COPD susceptibility locus.