In silico modeling of the pore region of a KCNQ4 missense mutant from a patient with hearing loss
© Matsunaga et al; licensee BioMed Central Ltd. 2012
Received: 28 November 2011
Accepted: 15 March 2012
Published: 15 March 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.
KeywordsChannel Deafness Electrostatic Hearing loss Molecular Modeling Mutation Potassium Structure
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
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 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.
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.
- Smith RJ, Bale JF, White KR: Sensorineural hearing loss in children. Lancet. 2005, 365: 879-890. 10.1016/S0140-6736(05)71047-3.PubMedView ArticleGoogle Scholar
- Morton CC, Nance WE: Newborn hearing screening-a silent revolution. N Engl J Med. 2006, 354: 2151-2164. 10.1056/NEJMra050700.PubMedView ArticleGoogle Scholar
- Kochhar A, Hildebrand MS, Smith RJ: Clinical aspects of hereditary hearing loss. Genet Med. 2007, 9: 393-408. 10.1097/GIM.0b013e3180980bd0.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Robbins J: KCNQ potassium channels: physiology, pathophysiology, and pharmacology. Pharmacol Ther. 2001, 90: 1-19. 10.1016/S0163-7258(01)00116-4.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Peitsch MC: Protein modeling by E-mail. Bio/Technology. 1995, 13: 658-660. 10.1038/nbt0795-658.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Lüthy R, Bowie JU, Eisenberg D: Assessment of protein models with three-dimensional profiles. Nature. 1992, 356: 83-85. 10.1038/356083a0.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Hol WG: Effects of the α-helix dipole upon the functioning and structure of proteins and peptides. Adv Biophys. 1985, 19: 133-165.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Xu XH, Yeung ES: Long-range electrostatic trapping of single-protein molecules at a liquid-solid interface. Science. 1998, 281: 1650-1653.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Larsen AE, Grier DG: Like-charge attractions in metastable colloidal crystallites. Nature. 1995, 385: 230-233.View ArticleGoogle Scholar
- 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: e12722Google Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.