- Research article
- Open Access
In silico modeling of the pore region of a KCNQ4 missense mutant from a patient with hearing loss
BMC Research Notesvolume 5, Article number: 145 (2012)
Mutation of the voltage-gated potassium channel KCNQ4 causes DFNA2-type nonsyndromic autosomal dominant sensorineural hearing loss. KCNQ4 is expressed predominantly in the auditory sensory outer hair cells, which are critical for sound amplification.
We sequenced KCNQ4 from Japanese patients with sensorineural hearing loss, and identified a novel missense mutation encoding a Tyr270His located at the N-terminus of the highly conserved pore helix sequence. As this patient was not accessible to us and information about them was limited, we used molecular modeling to investigate whether this novel mutation is hypothetically pathogenic. A careful examination of an in silico structural model of the KCNQ4 pore region revealed that the Tyr270His mutation caused an alteration in the electrostatic surface potential of the pore helix.
We propose two possible means by which the Tyr270His mutation causes hearing loss: a positively charged His270 side chain might enhance the helix dipole moment of the pore helix, thereby destabilizing the helix and/or the pore region, or it might disturb transport of K+ through the channel by electrostatic repulsion.
Deafness is a frequently inherited sensory disorder. For every 1000 newborns, more than 1% have bilateral sensorineural hearing loss (SNHL), and 50-70% of the cases are monogenic disorders [1, 2]. Hereditary hearing loss is classified as syndromic and nonsyndromic . The genetic causes of nonsyndromic hearing loss are autosomal dominant, autosomal recessive, X-chromosome linked, and mitochondrial in nature. To date, 25, 40, 3, and 6 genes have been identified that are responsible for autosomal dominant-, autosomal recessive-, X-chromosome linked-, and mitochondrial hearing loss, respectively http://hereditaryhearingloss.org/. In particular, autosomal dominant nonsyndromic sensorineural deafness type 2 (DFNA2) hearing loss affects the ability of children to hear high frequencies, and results in hearing loss at all frequencies later in life . KCNQ4 is the causative gene  and encodes the membrane protein KQT-like, subfamily member 4, which contains 695 amino acid residues in its longest isoform.
KCNQ4 is a member of the five muscarinic receptor-regulated and voltage-gated potassium (M-type K+ channel) channels present in ear and epithelia . In the inner ear, KCNQ4 subunits assemble to form homomeric channels that help with the recycling of potassium . Because patients with KCNQ4 mutations show progressive hearing loss, development of drugs that improve KCNQ4 function should attenuate hearing loss.
While conducting an ongoing genetic study of hereditary hearing loss, we identified a mutated KCNQ4 from a Japanese patient in which the encoded Tyr270 residue was replaced with His. Structural analysis of mutations using molecular modeling is a practical technology to investigate their pathogenicity in patients who are not readily available for analysis. To determine how the mutated protein (denoted Tyr270His) induced hearing loss, we built models of the KCNQ4 and Tyr270His pore regions. By examining models that included electrostatic potentials, we were able to propose two molecular mechanisms that could have caused the patient's hearing loss.
Results and discussion
Prior to this report, 14 missense mutations, a splice-site mutation, and 2 small deletion mutations in KCNQ4 were reported to be associated with hearing loss [7–9]. Because 11 of the missense mutations were found in the pore region, it was considered a pathogenic "hot spot". We therefore sequenced KCNQ4 exons 5, 6, and 7 (which encode the pore region) of Japanese patients with bilateral nonsyndromic hearing loss. A new heterozygotically presented KCNQ4 missense mutation was found at position 808, in which the wild-type thymine (T) was replaced by cytosine (C). This results in a Tyr270His mutation in the protein (Figure 1).
The patient who carried the Tyr270His mutation was a 28-year-old female with congenital, progressive hearing loss. An audiometric examination found mild hearing loss in her left ear and moderate hearing loss in her right ear, with steeply sloping audiometric configurations in both ears (Figure 2). The patient had no inner ear malformations, according to computational tomography. None of her family members had hearing loss, and no environmental factors that might have caused the hearing loss were reported (e.g., exposure to noise or ototoxic drugs.) DNA samples were not obtained from her family.
Similar to other voltage-gated K+ (Kv) channels, KCNQ4 probably contains six transmembrane α-helices, (S1-S6), a pore helix (PH), a pore-loop (P-loop), a short N-terminal region, and a long C-terminal region (Figure 3A). The pore regions of Kv channels are formed by the contiguous regions S5, PH, P-loop, and S6. Functional and structural characterizations of Kv channels have demonstrated that the PH and P-loop are responsible for selective K+ transport, whereas S4, acting as a voltage sensor, regulates the open/closed state of the channel [10, 11]. A modeling study of wild-type and mutant KCNQ3 channels suggested a critical stabilizing interaction between the PH and the S4-selectivity filter that is responsible for rearrangements of the P-loop architecture induced by the presence of a hydroxyl-containing residue at the 315 position, "unlocking" the channels into a conductive conformation . Tyr270 is predicted to be the N-terminal residue of the PH. Available mammalian KCNQ4 pore-region sequences are identical, and those of non-mammals (e.g., birds, fruit flies, tunicates, and nematodes; Figure 3B) are very similar to mammalian sequences. In addition, the PH and P-loop sequences of non-KCNQ4 Kv channels (e.g., Kv2.1 from Homo sapiens is identical in 14/24 amino acid residues, Kv1.2 from Rattus norvegicus in 11/24, and KvAP, a prokaryotic Kv channel from Aeropyrum pernix in 13/24) are somewhat conserved, suggesting that the pore residues, including Tyr270, are necessary for the proper functioning of the channel (Figure 3C).
Given the patient's clinical data and that the Tyr270His mutation is within the highly conserved pore region, we built structural models of wild-type KCNQ4 and the Tyr270His pore regions in an attempt to understand how the mutation affected K+ transport and, consequently, the patient's hearing. To identify a template crystal structure, we searched the Protein Data Bank using Gapped BLAST  and PDBsum , to obtain a three-dimensional structure of a KCNQ4 homolog. Among the channels with available three-dimensional structures, the Kv1.2 sequence was most similar to that of KCNQ4 (27.5% identity for the transmembrane sequence Ser32 to Thr417). Consequently, structural models for the transmembrane regions (Tyr80 to Gln329) of wild-type KCNQ4 and Tyr270His were built based on the template given by the crystal structure of Kv1.2 using SWISS-MODEL http://swissmodel.expasy.org/[15–17]. The reliabilities of the models were evaluated using Verify_3D http://nihserver.mbi.ucla.edu/Verify_3D/[18, 19], and were found to be trustworthy (see Additional file 1: Figure S1).
The α-carbon tracing of the KCNQ4 model was visualized using Chimera http://plato.cgl.ucsf.edu/chimera/. Tyr270, at the N-terminus of PH, was found at the same position as the corresponding Kv1.2 residue (Ile361) (Figure 4A). We therefore assessed the electrostatic characteristics of the two models to determine if a positively charged His270 side chain might be responsible for the patient's hearing loss.
The ribbon model of wild-type KCNQ4 with the electrostatic surface potential superimposed, suggests that the side chain of Tyr270 (colored white in Figure 4B) should be electrostatically neutral. In contrast, for Tyr270His, the side chain of His270 (colored blue in Figure 4C) is predicted to retain at least a partial positive charge, which reflects a standard p Ka value (6.5) for Histidine. Moreover, His270 is surrounded by the negatively charged residues Asp272, Asp266, Asp262, Glu260, and the polar residues Ser269, Ser268, and Ser265, which are capable of hydrogen bonding (all the aforementioned residues are within 10 Å of His270, Figure 4D). In hemoglobin (PDB ID: 2hhb), the side-chain p Ka of His97, which is located in the loop domain between helix3 and helix4 and surrounded by negatively charged residues of Asp94 (located contiguous to the C-terminus of helix 3) and Asp99 (N-terminus of helix 4) within 5 Å distance, is shifted to an abnormally high value (pKa = 7.8) . For these reasons, the side chain of His270 may also have a p K value that is greater than the standard value, and, therefore, it was predicted that more than half of the histidine side chains at position 270 carry a positive charge near physiological pH (Figure 4C).
The alteration of the electrostatic surface potential of a single helical residue in the pore region might affect the structural stability of the channel as a consequence of an alteration in the helix dipole moment. Substantial dipole moments with positive and negative unit charges at the N- and C-terminals, respectively, of α-helices are often found . Comparison of three Kv channel structures suggests that the dipoles of the transmembrane helices orient the helices and are responsible, at least in part, for the structure of the pore [22–24]. The electrically neutral Tyr270, which, as noted above, is at the N-terminus of the PH, is located within 9 Å of Ala263, the S5 C-terminus (Figure 4B, C). The presence of a positively charged His adjacent to, or at the N-terminal of, an α-helix can destabilize the helix dipole . Replacement of Tyr270 with His may therefore increase the dipole moment of the PH.
Alternatively, or in conjunction with the destabilization effect, His270 may impede K+ transport. The electrostatic repulsion between a positively charged His270 and K+ would be stronger than that between the electrically neutral Tyr270 and K+. In the Tyr270His model, the distance between His270 and the center pore of the channel was ~20 Å (Figure 4E, F). This distance is small enough to affect the electrostatic interaction between two charged molecules in a non-polar environment (e.g., the interior of a membrane protein) because the electrostatic force is strong and affects charges separated by 500-1800 Å [26–28]. Additionally, the dielectric constantsin the interior of a protein and in a lipid membrane are assumed to be 5 and 2, respectively, whereas that for bulk water is ~80 at room temperature . The electrostatic potential energy for two charges in protein or lipid membrane interiors is therefore between 16 and 40 times greater than that in water. Consequently, the long-range electrostatic repulsive force between the positively charged His270 and a K+ might possibly impede K+ passage through the channel , whereas the force on K+ would be smaller in the extracellular region.
Several mutations (e.g., Leu274His and Trp276Ser) in the KCNQ4 PH have been correlated with hearing loss . The clinical features, congenital, progressive high-frequency hearing loss without substantial loss of speech recognition during the first decade of life, were the same for the patients with the aforementioned mutated KCNQ4s and the patient with Tyr270His [8, 31]. Interestingly, ectopic expression of the Trp276Ser mutant and wild-type KCNQ4 in a cultured cell line caused a reduction in the channel current, whereas expression of only the mutant caused impaired trafficking of the protein to the cellular membrane . The similar clinical symptoms and locations of the mutations within KCNQ4 support the proposal that Tyr270His is a pathogenic mutant associated with progressive SNHL. Further physiological and three-dimensional structural characterization of Tyr270His may identify which of our working hypotheses, structural distortion of the channel caused by a change in the PH dipole moment or electrostatic impediment of K+ transport, or both, caused hearing loss in the patient, and provide insight into how to reverse hearing loss caused by KCNQ4 mutants.
In summary, we identified a new KCNQ4 mutant, in which Tyr270 was replaced with His in a patient with progressive SNHL. Structural models of wild-type and Tyr270His pore regions revealed that His270 caused an alteration in the electrostatic surface potential of the PH, which, in turn, suggests that the molecular pathology of Tyr270His might cause the hearing loss in two ways. First, the increased dipole moment of the PH caused by the positively charged His270 might destabilize the PH structure, as well as that of the pore region. In parallel, the electrostatic repulsive force between His270 and K+ might impede K+ transport through the channel.
The Research Ethics Committee of the National Tokyo Medical Center approved this study. Subjects with bilateral sensorineural hearing loss were recruited by the National Tokyo Medical Center and collaborating hospitals. Informed, written consent was obtained from the patients and their parents before removing genetic material. Medical histories and physical examinations were acquired for 127 patients whose hearing loss was not a syndromic symptom or a consequence of environmental factors (e.g., perinatal complication, meningitis, prenatal or postnatal drug ototoxicity, or acoustic trauma). Before sequencing the pore regions of the patients' KCNQ4, their mitochondrial DNA sequences (m.1555A > G or m.3243A > G) and GJB2 were determined to be normal (mutations in these genomic regions are known to be the most frequent causes of hereditary hearing loss [33, 34]). Each patient's hearing level was examined by pure-tone audiometry. Hearing loss was classified according to pure-tone average values at 0.5, 1, 2, and 4 kHz http://hereditaryhearingloss.org/.
Genomic DNA was extracted from blood samples using the reagents of the Gentra Puregene Blood kit (QIAGEN, Hamburg, Germany). For polymerase-chain-reaction (PCR) amplification of KCNQ4 exons 5, 6, and 7, the primer sets 5'-GAGATGGGGGACCTTTATCC-3' and 5'-AGCCCTACAAAGACCCTCAC-3'; 5'-GACCAGTCCTGCCTGTAACC-3' and 5'-AACTGAGCAGGAGGCAACTC-3'; and 5'-ACCCTTGCAGCCTCTTACTG-3' and 5'-CTGCTCCTAGGGCTTCTTCC-3; respectively, were used in conjunction with AmpliTaq Gold DNA polymerase (Applied Biosystems, Foster City, CA) and a PC-818 Program Temp Control System (ASTEC, Shizuoka, Japan). The PCR program was: 95°C for 4 min; 30 three-step cycles of 95°C for 30 s, 59.4°C for 30 s, and 72°C for 45 s; and then 72°C for 5 min. Primer sets for all KCNQ4 exons (Additional file 2: Table S1) were used to characterize the DNA encoding Tyr270His. These exons were amplified using PrimeSTAR HS DNA polymerase (Takara Bio, Shiga, Japan). All other PCR-amplification conditions were the same as mentioned above. PCR products were sequenced using an ABI 3730 DNA sequence analyzer and the reagents of the ABI Prism Big Dye Terminator Cycle Sequencing kit (Applied Biosystems). Characterization of the sequences was undertaken using SeqScape software v2.6 (Applied Biosystems) and DNASIS Pro (Hitachisoft, Tokyo, Japan), with the KCNQ4 sequence (NG_008139, NCBI Build37.1) as the reference. Control DNA was obtained from 96 Japanese subjects with normal hearing.
Molecular modeling of KCNQ4
To find a voltage-gated K+ channel with a sequence homologous to that of KCNQ4 that could be used as the structural template for the modeling exercise, we searched the Protein Data Bank using Gapped BLAST  and PDBsum (European Bioinformatics Institute: http://www.ebi.ac.uk/pdbsum/. The transmembrane helices and the voltage sensor and pore regions of KCNQ4 and Tyr270His were then modeled using SWISS-MODEL [15–17], with the crystal structure of Kv1.2 (PDB ID: 3LUT, chain B) as the template . Both KCNQ4 models were built using the automatic modeling mode and default parameters. The qualities of the models were evaluated using the Verify_3D Structure Evaluation Server (see Additional file 2: Table S1) [18, 19]. Finally, the models were each superimposed onto Kv1.2 using Chimera  to visualize the modeled α-carbon frames and ribbon models, with the electrostatic surface potentials of KCNQ4 or Tyr270His overlaid.
Nonsyndromic autosomal dominant sensorineural deafness type 2
Potassium voltage-gated channel
- KQT-like subfamily:
Voltage-gated potassium channel
Polymerase chain reaction
Sensorineural hearing loss.
Smith RJ, Bale JF, White KR: Sensorineural hearing loss in children. Lancet. 2005, 365: 879-890. 10.1016/S0140-6736(05)71047-3.
Morton CC, Nance WE: Newborn hearing screening-a silent revolution. N Engl J Med. 2006, 354: 2151-2164. 10.1056/NEJMra050700.
Kochhar A, Hildebrand MS, Smith RJ: Clinical aspects of hereditary hearing loss. Genet Med. 2007, 9: 393-408. 10.1097/GIM.0b013e3180980bd0.
Coucke P, Van Camp G, Djoyodiharjo B, Smith SD, Frants RR, Padberg GW, Darby JK, Huizing EH, Cremers CW, Kimberling WJ: Linkage of autosomal dominant hearing loss to the short arm of chromosome 1 in two families. N Engl J Med. 1994, 331: 425-431. 10.1056/NEJM199408183310702.
Kubisch C, Schroeder BC, Friedrich T, Lütjohann B, El-Amraoui A, Marlin S, Petit C, Jentsch TJ: KCNQ4, a novel potassium channel expressed in sensory outer hair cells, is mutated in dominant deafness. Cell. 1999, 96: 437-446. 10.1016/S0092-8674(00)80556-5.
Robbins J: KCNQ potassium channels: physiology, pathophysiology, and pharmacology. Pharmacol Ther. 2001, 90: 1-19. 10.1016/S0163-7258(01)00116-4.
Topsakal V, Pennings RJ, te Brinke H, Hamel B, Huygen PL, Kremer H, Cremers CW: Phenotype determination guides swift genotyping of a DFNA2/KCNQ4 family with a hot spot mutation (W276S). Otol Neurotol. 2005, 26: 52-58. 10.1097/00129492-200501000-00009.
Van Hauwe P, Coucke PJ, Ensink RJ, Huygen P, Cremers CW, Van Camp G: Mutations in the KCNQ4 K+ channel gene, responsible for autosomal dominant hearing loss, cluster in the channel pore region. Am J Med Genet. 2000, 93: 184-187. 10.1002/1096-8628(20000731)93:3<184::AID-AJMG4>3.0.CO;2-5.
Coucke PJ, Van Hauwe P, Kelley PM, Kunst H, Schatteman I, Van Velzen D, Meyers J, Ensink RJ, Verstreken M, Declau F, Marres H, Kastury K, Bhasin S, McGuirt WT, Smith RJ, Cremers CW, Van de Heyning P, Willems PJ, Smith SD, Van Camp G: Mutations in the KCNQ4 gene are responsible for autosomal dominant deafness in four DFNA2 families. Hum Mol Genet. 1999, 8: 1321-1328. 10.1093/hmg/8.7.1321.
Chen X, Wang Q, Ni F, Ma J: Structure of the full-length Shaker potassium channel Kv1.2 by normal-mode-based X-ray crystallographic refinement. Proc Natl Acad Sci USA. 2010, 107: 11352-11357. 10.1073/pnas.1000142107.
Long SB, Tao X, Campbell EB, MacKinnon R: Atomic structure of a voltage-dependent K + channel in a lipid membrane-like environment. Nature. 2007, 450: 376-382. 10.1038/nature06265.
Zaika O, Hernandez CC, Bal M, Tolstykh GP, Shapiro MS: Determinants within the turret and pore-loop domains of KCNQ3 K + channels governing functional activity. Biophys J. 2008, 95: 5121-5137. 10.1529/biophysj.108.137604.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Laskowski RA: Enhancing the functional annotation of PDB structures in PDBsum using key figures extracted from the literature. Bioinformatics. 2007, 23: 1824-1827. 10.1093/bioinformatics/btm085.
Arnold K, Bordoli L, Kopp J, Schwede T: The SWISS-MODEL Workspace: A web-based environment for protein structure homology modeling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.
Kiefer F, Arnold K, Künzli M, Bordoli L, Schwede T: The SWISS-MODEL Repository and associated resources. Nucleic Acids Res. 2009, 37: 387-392. 10.1093/nar/gkn750.
Peitsch MC: Protein modeling by E-mail. Bio/Technology. 1995, 13: 658-660. 10.1038/nbt0795-658.
Bowie JU, Lüthy R, Eisenberg D: A method to identify protein sequences that fold into a known three-dimensional structure. Science. 1991, 253: 164-170. 10.1126/science.1853201.
Lüthy R, Bowie JU, Eisenberg D: Assessment of protein models with three-dimensional profiles. Nature. 1992, 356: 83-85. 10.1038/356083a0.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, Ferrin TE: UCSF Chimera-A Visualization System for Exploratory Research and Analysis. J Comput Chem. 2004, 25: 1605-1612. 10.1002/jcc.20084.
Hol WG: Effects of the α-helix dipole upon the functioning and structure of proteins and peptides. Adv Biophys. 1985, 19: 133-165.
Doyle DA, Morais Cabral J, Pfuetzner RA, Kuo A, Gulbis JM, Cohen SL, Chait BT, MacKinnon R: The structure of the potassium channel: molecular basis of K+ conduction and selectivity. Science. 1998, 280: 69-77. 10.1126/science.280.5360.69.
Roux B, MacKinnon R: The cavity and pore helices in the KcsA K+ channel: electrostatic stabilization of monovalent cations. Science. 1999, 285: 100-102. 10.1126/science.285.5424.100.
Dutzler R, Campbell EB, Cadene M, Chait BT, MacKinnon R: X-ray structure of a ClC chloride channel at 3.0 A reveals the molecular basis of anion selectivity. Nature. 2002, 415: 287-294. 10.1038/415287a.
Armstrong KM, Baldwin RL: Charged histidine affects alpha-helix stability at all positions in the helix by interacting with the backbone charges. 1993, 90: 11337-11340.
Xu XH, Yeung ES: Long-range electrostatic trapping of single-protein molecules at a liquid-solid interface. Science. 1998, 281: 1650-1653.
Kang SH, Yeung ES: Dynamics of single-protein molecules at a liquid/solid interface: implications in capillary electrophoresis and chromatography. Anal Chem. 2002, 74: 6334-6339. 10.1021/ac0261202.
Larsen AE, Grier DG: Like-charge attractions in metastable colloidal crystallites. Nature. 1995, 385: 230-233.
Callenberg KM, Choudhary OP, de Forest GL, Gohara DW, Baker NA, Grabe M: APBSmem: a graphical interface for electrostatic calculations at the membrane. PLoS ONE. 2010, 5: pii: e12722
Li N, Tang H, Gai H, Dong X, Wang Q, Yeung ES: Determination of protein surface excess on a liquid/solid interface by single-molecule counting. Anal Bioanalytical Chem. 2009, 394: 1879-1885. 10.1007/s00216-009-2888-4.
De Leenheer EM, Huygen PL, Coucke PJ, Admiraal RJ, van Camp G, Cremers CW: Longitudinal and cross-sectional phenotype analysis in a new, large Dutch DFNA2/KCNQ4 family. Ann Otol Rhinol Laryngol. 2002, 111: 267-274.
Kim HJ, Lv P, Sihn CR, Yamoah EN: Cellular and molecular mechanisms of autosomal dominant form of progressive hearing loss, DFNA2. J Biol Chem. 2011, 286: 1517-1527. 10.1074/jbc.M110.179010.
Matsunaga T, Hirota E, Bito S, Niimi S, Usami S: Clinical course of hearing and language development in GJB2 and non-GJB2 deafness following habilitation with hearing aids. Audiol Neuro. 2006, 11: 59-68. 10.1159/000089607.
Matsunaga T: Value of genetic testing in determination of therapy for auditory disorders. Keio J Med. 2009, 58: 216-222. 10.2302/kjm.58.216.
This study was supported by a Health Science Research Grant (H16-kankakuki-006 to TM) from the Ministry of Health, Labor, and Welfare of Japan, a Grant-in-Aid for Clinical Research to TM from the National Hospital Organization, and by a Grant-in-Aid for Scientific Research (19592001 to TM) from the Ministry of Education, Culture, Sports, Science and Technology of Japan.
The authors declare that they have no competing interests.
KN was involved in the design of the study, performed the experiments and data analysis, and drafted the manuscript. HM was involved in the conception of the study and polished the manuscript. HK was involved in the conception of the study and gave considerations about structural information. SH provided information about the patient with sensorineural hearing loss. TM conceived the study, and participated in its coordination, design, data analysis and production of the manuscript. All authors read and approved the final manuscript.
Kazunori Namba, Hideki Mutai contributed equally to this work.