Upgraded molecular models of the human KCNQ1 potassium channel

The voltage-gated potassium channel KCNQ1 (KV7.1) assembles with the KCNE1 accessory protein to generate the slow delayed rectifier current, IKS, which is critical for membrane repolarization as part of the cardiac action potential. Loss-of-function (LOF) mutations in KCNQ1 are the most common cause of congenital long QT syndrome (LQTS), type 1 LQTS, an inherited genetic predisposition to cardiac arrhythmia and sudden cardiac death. A detailed structural understanding of KCNQ1 is needed to elucidate the molecular basis for KCNQ1 LOF in disease and to enable structure-guided design of new anti-arrhythmic drugs. In this work, advanced structural models of human KCNQ1 in the resting/closed and activated/open states were developed by Rosetta homology modeling guided by newly available experimentally-based templates: X. leavis KCNQ1 and various resting voltage sensor structures. Using molecular dynamics (MD) simulations, the capacity of the models to describe experimentally established channel properties including state-dependent voltage sensor gating charge interactions and pore conformations, PIP2 binding sites, and voltage sensor–pore domain interactions were validated. Rosetta energy calculations were applied to assess the utility of each model in interpreting mutation-evoked KCNQ1 dysfunction by predicting the change in protein thermodynamic stability for 50 experimentally characterized KCNQ1 variants with mutations located in the voltage-sensing domain. Energetic destabilization was successfully predicted for folding-defective KCNQ1 LOF mutants whereas wild type-like mutants exhibited no significant energetic frustrations, which supports growing evidence that mutation-induced protein destabilization is an especially common cause of KCNQ1 dysfunction. The new KCNQ1 Rosetta models provide helpful tools in the study of the structural basis for KCNQ1 function and can be used to generate hypotheses to explain KCNQ1 dysfunction.

Related papers:

title: Atomic structure of a voltage-dependent K+ channel in a lipid membrane-like environment
abstract: Voltage-dependent K+ (Kv) channels repolarize the action potential in neurons and muscle. This type of channel is gated directly by membrane voltage through protein domains known as voltage sensors, which are molecular voltmeters that read the membrane voltage and regulate the pore. Here we describe the structure of a chimaeric voltage-dependent K+ channel, which we call the ‘paddle-chimaera channel’, in which the voltage-sensor paddle has been transferred from Kv2.1 to Kv1.2. Crystallized in complex with lipids, the complete structure at 2.4 ångström resolution reveals the pore and voltage sensors embedded in a membrane-like arrangement of lipid molecules. The detailed structure, which can be compared directly to a large body of functional data, explains charge stabilization within the membrane and suggests a mechanism for voltage-sensor movements and pore gating.

title: Structural mechanism of voltage-dependent gating in an isolated voltage-sensing domain
abstract: The transduction of transmembrane electric fields into protein motion has an essential role in the generation and propagation of cellular signals. Voltage-sensing domains (VSDs) carry out these functions through reorientations of positive charges in the S4 helix. Here, we determined crystal structures of the Ciona intestinalis VSD (Ci-VSD) in putatively active and resting conformations. S4 undergoes an ~5-Å displacement along its main axis, accompanied by an ~60° rotation. This movement is stabilized by an exchange in countercharge partners in helices S1 and S3 that generates an estimated net charge transfer of ~1 eo. Gating charges move relative to a ''hydrophobic gasket' that electrically divides intra- and extracellular compartments. EPR spectroscopy confirms the limited nature of S4 movement in a membrane environment. These results provide an explicit mechanism for voltage sensing and set the basis for electromechanical coupling in voltage-dependent enzymes and ion channels.

title: Tracking a complete voltage-sensor cycle with metal-ion bridges
abstract: Voltage-gated ion channels open and close in response to changes in membrane potential, thereby enabling electrical signaling in excitable cells. The voltage sensitivity is conferred through four voltage-sensor domains (VSDs) where positively charged residues in the fourth transmembrane segment (S4) sense the potential. While an open state is known from the Kv1.2/2.1 X-ray structure, the conformational changes underlying voltage sensing have not been resolved. We present 20 additional interactions in one open and four different closed conformations based on metal-ion bridges between all four segments of the VSD in the voltage-gated Shaker K channel. A subset of the experimental constraints was used to generate Rosetta models of the conformations that were subjected to molecular simulation and tested against the remaining constraints. This achieves a detailed model of intermediate conformations during VSD gating. The results provide molecular insight into the transition, suggesting that S4 slides at least 12 Å along its axis to open the channel with a 310 helix region present that moves in sequence in S4 in order to occupy the same position in space opposite F290 from open through the three first closed states.

title: PROCHECK: a program to check the stereochemical quality of protein structures
abstract: The PROCHECK suite of programs provides a detailed check on the stereochemistry of a protein structure. Its outputs comprise a number of plots in PostScript format and a comprehensive residue-by-residue listing. These give an assessment of the overall quality of the structure as compared with well refined structures of the same resolution and also highlight regions that may need further investigation. The PROCHECK programs are useful for assessing the quality not only of protein structures in the process of being solved but also of existing structures and of those being modelled on known structures.

title: MolProbity: all-atom contacts and structure validation for proteins and nucleic acids
abstract: MolProbity is a general-purpose web server offering quality validation for 3D structures of proteins, nucleic acids and complexes. It provides detailed all-atom contact analysis of any steric problems within the molecules as well as updated dihedral-angle diagnostics, and it can calculate and display the H-bond and van der Waals contacts in the interfaces between components. An integral step in the process is the addition and full optimization of all hydrogen atoms, both polar and nonpolar. New analysis functions have been added for RNA, for interfaces, and for NMR ensembles. Additionally, both the web site and major component programs have been rewritten to improve speed, convenience, clarity and integration with other resources. MolProbity results are reported in multiple forms: as overall numeric scores, as lists or charts of local problems, as downloadable PDB and graphics files, and most notably as informative, manipulable 3D kinemage graphics shown online in the KiNG viewer. This service is available free to all users at http://molprobity.biochem.duke.edu.

title: Disease-associated mutations cause premature oligomerization of myelin proteolipid protein in the endoplasmic reticulum.
abstract: Pelizaeus-Merzbacher disease (PMD) is a dysmyelinating disease caused by mutations, deletions, or duplications of the proteolipid protein (PLP) gene. Mutant forms of PLP are retained in the endoplasmic reticulum (ER), and the resulting accumulation of mutant protein is thought to be a direct cause of oligodendrocyte cell death, which is the primary clinical feature of PMD. The molecular mechanisms underlying the toxicity of mutant PLP are however currently unknown. We report here that PMD-linked mutations of PLP are associated with the accelerated assembly of the protein into stable homooligomers that resemble mature, native PLP. Thus although WT PLP forms stable oligomers after an extended maturation period, most likely at the cell surface, mutant forms of PLP rapidly assemble into such oligomers at the ER. Using PLP mutants associated with diseases of varying severity, we show that the formation of stable oligomers correlates with the development of PMD. Based on these findings, we propose that the premature oligomerization of PLP in the ER of oligodendrocytes contributes to the pathology of PMD.

title: Structural models for the KCNQ1 voltage-gated potassium channel.
abstract: Mutations in the human voltage-gated potassium channel KCNQ1 are associated with predisposition to deafness and various cardiac arrhythmia syndromes including congenital long QT syndrome, familial atrial fibrillation, and sudden infant death syndrome. In this work 3-D structural models were developed for both the open and closed states of human KCNQ1 to facilitate structurally based hypotheses regarding mutation-phenotype relationships. The KCNQ1 open state was modeled using Rosetta in conjunction with Molecular Operating Environment software, and is based primarily on the recently determined open state structure of rat Kv1.2 (Long, S. B., et al. (2005) Science 309, 897-903). The closed state model for KCNQ1 was developed based on the crystal structures of bacterial potassium channels and the closed state model for Kv1.2 of Yarov-Yarovoy et al. ((2006) ). Using the new models for KCNQ1, we generated a database for the location and predicted residue-residue interactions for more than 85 disease-linked sites in both open and closed states. These data can be used to generate structure-based hypotheses for disease phenotypes associated with each mutation. The potential utility of these models and the database is exemplified by the surprising observation that four of the five known mutations in KCNQ1 that are associated with gain-of-function KCNQ1 defects are predicted to share a common interface in the open state structure between the S1 segment of the voltage sensor in one subunit and both the S5 segment and top of the pore helix from another subunit. This interface evidently plays an important role in channel gating.

title: Differential roles of S6 domain hinges in the gating of KCNQ potassium channels.
abstract: Voltage-gated K+ channel activation is proposed to result from simultaneous bending of all S6 segments away from the central axis, enlarging the aperture of the pore sufficiently to permit diffusion of K+ into the water-filled central cavity. The hinge position for the bending motion of each S6 segment is proposed to be a Gly residue and/or a Pro-Val-Pro motif in Kv1-Kv4 channels. The KCNQ1 (Kv7.1) channel has Ala-336 in the Gly-hinge position and Pro-Ala-Gly. Here we show that mutation of Ala-336 to Gly in KCNQ1 increased current amplitude and shifted the voltage dependence of activation to more negative potentials, consistent with facilitation of hinge activity that favors the open state. In contrast, mutation of Ala-336 to Cys or Thr shifted the voltage dependence of activation to more positive potentials and reduced current amplitude. Mutation of the putative Gly hinge to Ala in KCNQ2 (Kv7.2) abolished channel function. Mutation-dependent changes in current amplitude, but not kinetics, were found in heteromeric KCNQ1/KCNE1 channels. Mutation of the Pro or Gly of the Pro-Ala-Gly motif to Ala abolished KCNQ1 function and introduction of Gly in front of the Ala-mutations partially recovered channel function, suggesting that flexibility at the PAG is important for channel activation.

title: State-dependent electrostatic interactions of S4 arginines with E1 in S2 during Kv7.1 activation
abstract: The voltage-sensing domain of voltage-gated channels is comprised of four transmembrane helices (S1–S4), with conserved positively charged residues in S4 moving across the membrane in response to changes in transmembrane voltage. Although it has been shown that positive charges in S4 interact with negative countercharges in S2 and S3 to facilitate protein maturation, how these electrostatic interactions participate in channel gating remains unclear. We studied a mutation in Kv7.1 (also known as KCNQ1 or KvLQT1) channels associated with long QT syndrome (E1K in S2) and found that reversal of the charge at E1 eliminates macroscopic current without inhibiting protein trafficking to the membrane. Pairing E1R with individual charge reversal mutations of arginines in S4 (R1–R4) can restore current, demonstrating that R1–R4 interact with E1. After mutating E1 to cysteine, we probed E1C with charged methanethiosulfonate (MTS) reagents. MTS reagents could not modify E1C in the absence of KCNE1. With KCNE1, (2-sulfonatoethyl) MTS (MTSES)− could modify E1C, but [2-(trimethylammonium)ethyl] MTS (MTSET)+ could not, confirming the presence of a positively charged environment around E1C that allows approach by MTSES− but repels MTSET+. We could change the local electrostatic environment of E1C by making charge reversal and/or neutralization mutations of R1 and R4, such that MTSET+ modified these constructs depending on activation states of the voltage sensor. Our results confirm the interaction between E1 and the fourth arginine in S4 (R4) predicted from open-state crystal structures of Kv channels and reveal an E1–R1 interaction in the resting state. Thus, E1 engages in electrostatic interactions with arginines in S4 sequentially during the gating movement of S4. These electrostatic interactions contribute energetically to voltage-dependent gating and are important in setting the limits for S4 movement.

title: Domain–domain interactions determine the gating, permeation, pharmacology, and subunit modulation of the IKs ion channel
abstract: Voltage-gated ion channels generate electrical currents that control muscle contraction, encode neuronal information, and trigger hormonal release. Tissue-specific expression of accessory (β) subunits causes these channels to generate currents with distinct properties. In the heart, KCNQ1 voltage-gated potassium channels coassemble with KCNE1 β-subunits to generate the IKs current (Barhanin et al., 1996; Sanguinetti et al., 1996), an important current for maintenance of stable heart rhythms. KCNE1 significantly modulates the gating, permeation, and pharmacology of KCNQ1 (Wrobel et al., 2012; Sun et al., 2012; Abbott, 2014). These changes are essential for the physiological role of IKs (Silva and Rudy, 2005); however, after 18 years of study, no coherent mechanism explaining how KCNE1 affects KCNQ1 has emerged. Here we provide evidence of such a mechanism, whereby, KCNE1 alters the state-dependent interactions that functionally couple the voltage-sensing domains (VSDs) to the pore.

title: Spectroscopic mapping of voltage sensor movement in the Shaker potassium channel
abstract: Voltage-gated ion channels underlie the generation of action potentials and trigger neurosecretion and muscle contraction. These channels consist of an inner pore-forming domain, which contains the ion permeation pathway and elements of its gates, together with four voltage-sensing domains, which regulate the gates. To understand the mechanism of voltage sensing it is necessary to define the structure and motion of the S4 segment, the portion of each voltage-sensing domain that moves charged residues across the membrane in response to voltage change. We have addressed this problem by using fluorescence resonance energy transfer as a spectroscopic ruler to determine distances between S4s in the Shaker K+ channel in different gating states. Here we provide evidence consistent with S4 being a tilted helix that twists during activation. We propose that helical twist contributes to the movement of charged side chains across the membrane electric field and that it is involved in coupling voltage sensing to gating.

title: Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations
abstract: The response of a membrane-bound Kv1.2 ion channel to an applied transmembrane potential has been studied using molecular dynamics simulations. Channel deactivation is shown to involve three intermediate states of the voltage sensor domain (VSD), and concomitant movement of helix S4 charges 10–15 Å along the bilayer normal; the latter being enabled by zipper-like sequential pairing of S4 basic residues with neighboring VSD acidic residues and membrane-lipid head groups. During the observed sequential transitions S4 basic residues pass through the recently discovered charge transfer center with its conserved phenylalanine residue, F233. Analysis indicates that the local electric field within the VSD is focused near the F233 residue and that it remains essentially unaltered during the entire process. Overall, the present computations provide an atomistic description of VSD response to hyperpolarization, add support to the sliding helix model, and capture essential features inferred from a variety of recent experiments.

title: Voltage-Dependent Gating: Novel Insights from KCNQ1 Channels.
abstract: Gating of voltage-dependent cation channels involves three general molecular processes: voltage sensor activation, sensor-pore coupling, and pore opening. KCNQ1 is a voltage-gated potassium (Kv) channel whose distinctive properties have provided novel insights on fundamental principles of voltage-dependent gating. 1) Similar to other Kv channels, KCNQ1 voltage sensor activation undergoes two resolvable steps; but, unique to KCNQ1, the pore opens at both the intermediate and activated state of voltage sensor activation. The voltage sensor-pore coupling differs in the intermediate-open and the activated-open states, resulting in changes of open pore properties during voltage sensor activation. 2) The voltage sensor-pore coupling and pore opening require the membrane lipid PIP2 and intracellular ATP, respectively, as cofactors, thus voltage-dependent gating is dependent on multiple stimuli, including the binding of intracellular signaling molecules. These mechanisms underlie the extraordinary KCNE1 subunit modification of the KCNQ1 channel and have significant physiological implications.

title: The Role of S4 Charges in Voltage-dependent and Voltage-independent KCNQ1 Potassium Channel Complexes
abstract: Voltage-gated potassium (Kv) channels extend their functional repertoire by coassembling with MinK-related peptides (MiRPs). MinK slows the activation of channels formed with KCNQ1 α subunits to generate the voltage-dependent IKs channel in human heart; MiRP1 and MiRP2 remove the voltage dependence of KCNQ1 to generate potassium “leak” currents in gastrointestinal epithelia. Other Kv α subunits interact with MiRP1 and MiRP2 but without loss of voltage dependence; the mechanism for this disparity is unknown. Here, sequence alignments revealed that the voltage-sensing S4 domain of KCNQ1 bears lower net charge (+3) than that of any other eukaryotic voltage-gated ion channel. We therefore examined the role of KCNQ1 S4 charges in channel activation using alanine-scanning mutagenesis and two-electrode voltage clamp. Alanine replacement of R231, at the N-terminal side of S4, produced constitutive activation in homomeric KCNQ1 channels, a phenomenon not observed with previous single amino acid substitutions in S4 of other channels. Homomeric KCNQ4 channels were also made constitutively active by mutagenesis to mimic the S4 charge balance of R231A-KCNQ1. Loss of single S4 charges at positions R231 or R237 produced constitutively active MinK-KCNQ1 channels and increased the constitutively active component of MiRP2-KCNQ1 currents. Charge addition to the CO2H-terminal half of S4 eliminated constitutive activation in MiRP2-KCNQ1 channels, whereas removal of homologous charges from KCNQ4 S4 produced constitutively active MiRP2-KCNQ4 channels. The results demonstrate that the unique S4 charge paucity of KCNQ1 facilitates its unique conversion to a leak channel by ancillary subunits such as MiRP2.

title: Phospholipids and the origin of cationic gating charges in voltage sensors
abstract: Cells communicate with their external environment through physical and chemical processes that take place in the cell-surrounding membrane. The membrane serves as a barrier as well as a special environment in which membrane proteins are able to carry out important processes. Certain membrane proteins have the ability to detect the membrane voltage and regulate ion conduction or enzyme activity. Such voltage-dependent processes rely on the action of protein domains known as voltage sensors, which are embedded inside the cell membrane and contain an excess of positively charged amino acids, which react to an electric field. How does the membrane create an environment suitable for voltage sensors? Here we show under a variety of conditions that the function of a voltage-dependent K+ channel is dependent on the negatively charged phosphodiester of phospholipid molecules. A non-voltage-dependent K+ channel does not exhibit the same dependence. The data lead us to propose that the phospholipid membrane, by providing stabilizing interactions between positively charged voltage-sensor arginine residues and negatively charged lipid phosphodiester groups, provides an appropriate environment for the energetic stability and operation of the voltage-sensing machinery. We suggest that the usage of arginine residues in voltage sensors is an adaptation to the phospholipid composition of cell membranes.

title: Solution structure and phospholipid interactions of the isolated voltage-sensor domain from KvAP.
abstract: Voltage-sensor domains (VSDs) are specialized transmembrane segments that confer voltage sensitivity to many proteins such as ion channels and enzymes. The activities of these domains are highly dependent on both the chemical properties and the physical properties of the surrounding membrane environment. To learn about VSD-lipid interactions, we used nuclear magnetic resonance spectroscopy to determine the structure and phospholipid interface of the VSD from the voltage-dependent K(+) channel KvAP (prokaryotic Kv from Aeropyrum pernix). The solution structure of the KvAP VSD solubilized within phospholipid micelles is similar to a previously determined crystal structure solubilized by a nonionic detergent and complexed with an antibody fragment. The differences observed include a previously unidentified short amphipathic α-helix that precedes the first transmembrane helix and a subtle rigid-body repositioning of the S3-S4 voltage-sensor paddle. Using (15)N relaxation experiments, we show that much of the VSD, including the pronounced kink in S3 and the S3-S4 paddle, is relatively rigid on the picosecond-to-nanosecond timescale. In contrast, the kink in S3 is mobile on the microsecond-to-millisecond timescale and may act as a hinge in the movement of the paddle during channel gating. We characterized the VSD-phospholipid micelle interactions using nuclear Overhauser effect spectroscopy and showed that the micelle uniformly coats the KvAP VSD and approximates the chemical environment of a phospholipid bilayer. Using paramagnetically labeled phospholipids, we show that bilayer-forming lipids interact with the S3 and S4 helices more strongly than with S1 and S2.

title: Phosphatidic acid modulation of Kv channel voltage sensor function
abstract: Membrane phospholipids can function as potent regulators of ion channel function. This study uncovers and investigates the effect of phosphatidic acid on Kv channel gating. Using the method of reconstitution into planar lipid bilayers, in which protein and lipid components are defined and controlled, we characterize two effects of phosphatidic acid. The first is a non-specific electrostatic influence on activation mediated by electric charge density on the extracellular and intracellular membrane surfaces. The second is specific to the presence of a primary phosphate group, acts only through the intracellular membrane leaflet and depends on the presence of a particular arginine residue in the voltage sensor. Intracellular phosphatidic acid accounts for a nearly 50 mV shift in the midpoint of the activation curve in a direction consistent with stabilization of the voltage sensor's closed conformation. These findings support a novel mechanism of voltage sensor regulation by the signaling lipid phosphatidic acid.  http://dx.doi.org/10.7554/eLife.04366.001

title: CHARMM‐GUI Membrane Builder toward realistic biological membrane simulations
abstract: CHARMM‐GUI Membrane Builder, http://www.charmm‐gui.org/input/membrane, is a web‐based user interface designed to interactively build all‐atom protein/membrane or membrane‐only systems for molecular dynamics simulations through an automated optimized process. In this work, we describe the new features and major improvements in Membrane Builder that allow users to robustly build realistic biological membrane systems, including (1) addition of new lipid types, such as phosphoinositides, cardiolipin (CL), sphingolipids, bacterial lipids, and ergosterol, yielding more than 180 lipid types, (2) enhanced building procedure for lipid packing around protein, (3) reliable algorithm to detect lipid tail penetration to ring structures and protein surface, (4) distance‐based algorithm for faster initial ion displacement, (5) CHARMM inputs for P21 image transformation, and (6) NAMD equilibration and production inputs. The robustness of these new features is illustrated by building and simulating a membrane model of the polar and septal regions of E. coli membrane, which contains five lipid types: CL lipids with two types of acyl chains and phosphatidylethanolamine lipids with three types of acyl chains. It is our hope that CHARMM‐GUI Membrane Builder becomes a useful tool for simulation studies to better understand the structure and dynamics of proteins and lipids in realistic biological membrane environments.

title: Structure of KCNE1 and implications for how it modulates the KCNQ1 potassium channel.
abstract: KCNE1 is a single-span membrane protein that modulates the voltage-gated potassium channel KCNQ1 (K V7.1) by slowing activation and enhancing channel conductance to generate the slow delayed rectifier current ( I Ks) that is critical for the repolarization phase of the cardiac action potential. Perturbation of channel function by inherited mutations in KCNE1 or KCNQ1 results in increased susceptibility to cardiac arrhythmias and sudden death with or without accompanying deafness. Here, we present the three-dimensional structure of KCNE1. The transmembrane domain (TMD) of KCNE1 is a curved alpha-helix and is flanked by intra- and extracellular domains comprised of alpha-helices joined by flexible linkers. Experimentally restrained docking of the KCNE1 TMD to a closed state model of KCNQ1 suggests that KCNE1 slows channel activation by sitting on and restricting the movement of the S4-S5 linker that connects the voltage sensor to the pore domain. We postulate that this is an adhesive interaction that must be disrupted before the channel can be opened in response to membrane depolarization.  Binding of KCNE1 to this "gain-of-function cleft" may explain how it increases conductance and stabilizes the open state. These working models for the KCNE1-KCNQ1 complexes may be used to formulate testable hypotheses for the molecular bases of disease phenotypes associated with the dozens of known inherited mutations in KCNE1 and KCNQ1.

title: Building KCNQ1/KCNE1 channel models and probing their interactions by molecular-dynamics simulations.
abstract: The slow delayed rectifier (I(KS)) channel is composed of KCNQ1 (pore-forming) and KCNE1 (auxiliary) subunits, and functions as a repolarization reserve in the human heart. Design of I(KS)-targeting anti-arrhythmic drugs requires detailed three-dimensional structures of the KCNQ1/KCNE1 complex, a task made possible by Kv channel crystal structures (templates for KCNQ1 homology-modeling) and KCNE1 NMR structures. Our goal was to build KCNQ1/KCNE1 models and extract mechanistic information about their interactions by molecular-dynamics simulations in an explicit lipid/solvent environment. We validated our models by confirming two sets of model-generated predictions that were independent from the spatial restraints used in model-building. Detailed analysis of the molecular-dynamics trajectories revealed previously unrecognized KCNQ1/KCNE1 interactions, whose relevance in I(KS) channel function was confirmed by voltage-clamp experiments. Our models and analyses suggest three mechanisms by which KCNE1 slows KCNQ1 activation: by promoting S6 bending at the Pro hinge that closes the activation gate; by promoting a downward movement of gating charge on S4; and by establishing a network of electrostatic interactions with KCNQ1 on the extracellular surface that stabilizes the channel in a pre-open activated state. Our data also suggest how KCNE1 may affect the KCNQ1 pore conductance.

title: How does KCNE1 regulate the Kv7.1 potassium channel? Model-structure, mutations, and dynamics of the Kv7.1-KCNE1 complex.
abstract: The voltage-gated potassium channel Kv7.1 and its auxiliary subunit KCNE1 are expressed in the heart and give rise to the major repolarization current. The interaction of Kv7.1 with the single transmembrane helix of KCNE1 considerably slows channel activation and deactivation, raises single-channel conductance, and prevents slow voltage-dependent inactivation. We built a Kv7.1-KCNE1 model-structure. The model-structure agrees with previous disulfide mapping studies and enables us to derive molecular interpretations of electrophysiological recordings that we obtained for two KCNE1 mutations. An elastic network analysis of Kv7.1 fluctuations in the presence and absence of KCNE1 suggests a mechanistic perspective on the known effects of KCNE1 on Kv7.1 function: slow deactivation is attributed to the low mobility of the voltage-sensor domains upon KCNE1 binding, abolishment of voltage-dependent inactivation could result from decreased fluctuations in the external vestibule, and amalgamation of the fluctuations in the pore region is associated with enhanced ion conductivity.

title: An Allosteric Mechanism for Drug Block of the Human Cardiac Potassium Channel KCNQ1
abstract: The intracellular aspect of the sixth transmembrane segment within the ion-permeating pore is a common binding site for many voltage-gated ion channel blockers. However, the exact site(s) at which drugs bind remain controversial. We used extensive site-directed mutagenesis coupled with molecular modeling to examine mechanisms in drug block of the human cardiac potassium channel KCNQ1. A total of 48 amino acid residues in the S6 segment, S4-S5 linker, and the proximal C-terminus of the KCNQ1 channel were mutated individually to alanine; alanines were mutated to cysteines. Residues modulating drug block were identified when mutant channels displayed <50% block on exposure to drug concentrations that inhibited wild-type current by ≥90%. Homology modeling of the KCNQ1 channel based on the Kv1.2 structure unexpectedly predicted that the key residue modulating drug block (F351) faces away from the permeating pore. In the open-state channel model, F351 lines a pocket that also includes residues L251 and V254 in S4-S5 linker. Docking calculations indicated that this pocket is large enough to accommodate quinidine. To test this hypothesis, L251A and V254A mutants were generated that display a reduced sensitivity to blockage with quinidine. Thus, our data support a model in which open state block of this channel occurs not via binding to a site directly in the pore but rather by a novel allosteric mechanism: drug access to a side pocket generated in the open-state channel configuration and lined by S6 and S4-S5 residues.

title: Regulation of ion channels by phosphatidylinositol 4,5-bisphosphate
abstract: Phosphatidylinositol 4,5-bisphosphate is a signaling phospholipid of the plasma membrane that has a dynamically changing concentration. In addition to being the precursor of inositol trisphosphate and diacylglycerol, it complexes with and regulates many cytoplasmic and membrane proteins. Recent work has characterized the regulation of a wide range of ion channels by phosphatidylinositol 4,5-bisphosphate, helping to redefine the role of this lipid in cells and in neurobiology. In most cases, phosphatidylinositol 4,5-bisphosphate increases channel activity, and its hydrolysis by phospholipase C reduces channel activity.

title: Kv7.1 ion channels require a lipid to couple voltage sensing to pore opening
abstract: Voltage-gated ion channels generate dynamic ionic currents that are vital to the physiological functions of many tissues. These proteins contain separate voltage-sensing domains, which detect changes in transmembrane voltage, and pore domains, which conduct ions. Coupling of voltage sensing and pore opening is critical to the channel function and has been modeled as a protein–protein interaction between the two domains. Here, we show that coupling in Kv7.1 channels requires the lipid phosphatidylinositol 4,5-bisphosphate (PIP2). We found that voltage-sensing domain activation failed to open the pore in the absence of PIP2. This result is due to loss of coupling because PIP2 was also required for pore opening to affect voltage-sensing domain activation. We identified a critical site for PIP2-dependent coupling at the interface between the voltage-sensing domain and the pore domain. This site is actually a conserved lipid-binding site among different K+ channels, suggesting that lipids play an important role in coupling in many ion channels.

title: Novel Kv7.1-Phosphatidylinositol 4,5-Bisphosphate Interaction Sites Uncovered by Charge Neutralization Scanning*
abstract: Phosphatidylinositol 4,5-bisphosphate (PIP2) regulates a variety of ion channels including Kv7 channels.  Several Kv7.1 mutations impair functional interaction with PIP2.  PIP2 interacts with several Kv7.1 amino acids, indicating the existence of at least two binding pockets.  Impaired molecular interaction between mutant Kv7.1 and PIP2 may underlie heart disease. Kv7.1 to Kv7.5 α-subunits belong to the family of voltage-gated potassium channels (Kv). Assembled with the β-subunit KCNE1, Kv7.1 conducts the slowly activating potassium current IKs, which is one of the major currents underlying repolarization of the cardiac action potential. A known regulator of Kv7 channels is the lipid phosphatidylinositol 4,5-bisphosphate (PIP2). PIP2 increases the macroscopic current amplitude by stabilizing the open conformation of 7.1/KCNE1 channels. However, knowledge about the exact nature of the interaction is incomplete. The aim of this study was the identification of the amino acids responsible for the interaction between Kv7.1 and PIP2. We generated 13 charge neutralizing point mutations at the intracellular membrane border and characterized them electrophysiologically in complex with KCNE1 under the influence of diC8-PIP2. Electrophysiological analysis of corresponding long QT syndrome mutants suggested impaired PIP2 regulation as the cause for channel dysfunction. To clarify the underlying structural mechanism of PIP2 binding, molecular dynamics simulations of Kv7.1/KCNE1 complexes containing two PIP2 molecules in each subunit at specific sites were performed. Here, we identified a subset of nine residues participating in the interaction of PIP2 and Kv7.1/KCNE1. These residues may form at least two binding pockets per subunit, leading to the stabilization of channel conformations upon PIP2 binding.

title: PIP2-dependent coupling is prominent in Kv7.1 due to weakened interactions between S4-S5 and S6
abstract: Among critical aspects of voltage-gated potassium (Kv) channels' functioning is the effective communication between their two composing domains, the voltage sensor (VSD) and the pore. This communication, called coupling, might be transmitted directly through interactions between these domains and, as recently proposed, indirectly through interactions with phosphatidylinositol-4,5-bisphosphate (PIP2), a minor lipid of the inner plasma membrane leaflet. Here, we show how the two components of coupling, mediated by protein-protein or protein-lipid interactions, both contribute in the Kv7.1 functioning. On the one hand, using molecular dynamics simulations, we identified a Kv7.1 PIP2 binding site that involves residues playing a key role in PIP2-dependent coupling. On the other hand, combined theoretical and experimental approaches have shown that the direct interaction between the segments of the VSD (S4–S5) and the pore (S6) is weakened by electrostatic repulsion. Finally, we conclude that due to weakened protein-protein interactions, the PIP2-dependent coupling is especially prominent in Kv7.1.

title: Long QT syndrome with compound mutations is associated with a more severe phenotype: a Japanese multicenter study.
abstract: Long QT syndrome (LQTS) can be caused by mutations in the cardiac ion channels. Compound mutations occur at a frequency of 4% to 11% among genotyped LQTS cases.



The purpose of this study was to determine the clinical characteristics and manner of onset of cardiac events in Japanese patients with LQTS and compound mutations.



Six hundred three genotyped LQTS patients (310 probands and 293 family members) were divided into two groups: those with a single mutation (n = 568) and those with two mutations (n = 35). Clinical phenotypes were compared between the two groups.



Of 310 genotyped probands, 26 (8.4%) had two mutations in the same or different LQTS-related genes (compound mutations). Among the 603 LQTS patients, compound mutation carriers had significantly longer QTc interval (510 ± 56 ms vs 478± 53 ms, P = .001) and younger age at onset of cardiac events (10 ± 8 years vs 18 ± 16 years, P = .043) than did single mutation carriers. The incidence rate of cardiac events before age 40 years and use of beta-blocker therapy among compound mutation carriers also were different than in single mutation carriers. Subgroup analysis showed more cardiac events in LQTS type 1 (LQT1) and type 2 (LQT2) compound mutations compared to single LQT1 and LQT2 mutations.



Compound mutation carriers are associated with a more severe phenotype than single mutation carriers.

title: Mutation analysis of cases of sudden unexplained death, 15 years after death: prompt genetic evaluation after resuscitation can save future lives.
abstract: The aim of this study is to use genetic mutation analysis to determine the cause of sudden unexpected death in young (SUDY) persons with normal autopsy findings, and to provide relatives with an identified cardiac mutation with suitable cardiovascular prevention.



We performed mutation analysis on blood samples from first-degree relatives of 25 cases with normal autopsy findings identified in the national Swedish study of sudden cardiac death in 15- to 35-year-olds from 1992 to 1999.



We found three families with long QT syndrome through mutation screening, and the mutations were verified in two of the deceased. Eight family members were found to be mutation carriers and have been provided with suitable cardiovascular prevention. Mutation screening also identified a number of common polymorphisms in the individuals screened. Clinical history revealed one family each with short QT syndrome and hypertrophic cardiomyopathy, respectively, but no mutations were found in the family members or in the deceased. Two SCDs each had occurred in two of the affected families.



Cardiac/genetic evaluation of relatives long after SUDY can reveal a diagnosis in 5/25 (20%) of cases. Since DNA extraction of formalin fixed paraffin embedded samples is unreliable, it is important that blood or tissue samples be stored at autopsy of such cases. This can facilitate establishing a diagnosis and thereby save lives in the future.

title: Positional cloning of a novel potassium channel gene: KVLQT1 mutations cause cardiac arrhythmias
abstract: Genetic factors contribute to the risk of sudden death from cardiac arrhythmias. Here, positional cloning methods establish KVLQT1 as the chromosome 11-linked LQT1 gene responsible for the most common inherited cardiac arrhythmia. KVLQT1 is strongly expressed in the heart and encodes a protein with structural features of a voltage-gated potassium channel. KVLQT1 mutations are present in affected members of 16 arrhythmia families, including one intragenic deletion and ten different missense mutations. These data define KVLQT1 as a novel cardiac potassium channel gene and show that mutations in this gene cause susceptibility to ventricular tachyarrhythmias and sudden death.

title: Spectrum and prevalence of mutations from the first 2,500 consecutive unrelated patients referred for the FAMILION long QT syndrome genetic test.
abstract: Long QT syndrome (LQTS) is a potentially lethal, highly treatable cardiac channelopathy for which genetic testing has matured from discovery to translation and now clinical implementation.



Here we examine the spectrum and prevalence of mutations found in the first 2,500 unrelated cases referred for the FAMILION LQTS clinical genetic test.



Retrospective analysis of the first 2,500 cases (1,515 female patients, average age at testing 23 +/- 17 years, range 0 to 90 years) scanned for mutations in 5 of the LQTS-susceptibility genes: KCNQ1 (LQT1), KCNH2 (LQT2), SCN5A (LQT3), KCNE1 (LQT5), and KCNE2 (LQT6).



Overall, 903 referral cases (36%) hosted a possible LQTS-causing mutation that was absent in >2,600 reference alleles; 821 (91%) of the mutation-positive cases had single genotypes, whereas the remaining 82 patients (9%) had >1 mutation in > or =1 gene, including 52 cases that were compound heterozygous with mutations in >1 gene. Of the 562 distinct mutations, 394 (70%) were missense, 428 (76%) were seen once, and 336 (60%) are novel, including 92 of 199 in KCNQ1, 159 of 226 in KCNH2, and 70 of 110 in SCN5A.



This cohort increases the publicly available compendium of putative LQTS-associated mutations by >50%, and approximately one-third of the most recently detected mutations continue to be novel. Although control population data suggest that the great majority of these mutations are pathogenic, expert interpretation of genetic test results will remain critical for effective clinical use of LQTS genetic test results.

title: Genetic testing in the long QT syndrome: development and validation of an efficient approach to genotyping in clinical practice.
abstract: In long QT syndrome (LQTS), disease severity and response to therapy vary according to the genetic loci. There exists a critical need to devise strategies to expedite genetic analysis.



To perform genetic screening in patients with LQTS to determine the yield of genetic testing, as well as the type and the prevalence of mutations.



We investigated whether the detection of a set of frequently mutated codons in the KCNQ1, KCNH2, and SCN5A genes may translate in a novel strategy for rapid efficient genetic testing of 430 consecutive patients referred to our center between June 1996 and June 2004. The entire coding regions of KCNQ1, KCNH2, SCN5A, KCNE1, and KCNE2 were screened by denaturing high-performance liquid chromatography and DNA sequencing. The frequency and the type of mutations were defined to identify a set of recurring mutations. A separate cohort of 75 consecutive probands was used as a validation group to quantify prospectively the prevalence of the recurring mutations identified in the primary LQTS population.



Development of a novel approach to LQTS genotyping.



We identified 235 different mutations, 138 of which were novel, in 310 (72%) of 430 probands (49% KCNQ1, 39% KCNH2, 10% SCN5A, 1.7% KCNE1, and 0.7% KCNE2). Fifty-eight percent of probands carried nonprivate mutations in 64 codons of KCNQ1, KCNH2, and SCN5A genes. A similar occurrence of mutations at these codons (52%) was confirmed in the prospective cohort of 75 probands and in previously published LQTS cohorts.



We have developed an approach to improve the efficiency of genetic screening for LQTS. This novel method may facilitate wider access to genotyping resulting in better risk stratification and treatment of LQTS patients.

title: Spectrum of pathogenic mutations and associated polymorphisms in a cohort of 44 unrelated patients with long QT syndrome
abstract: Long QT syndrome (LQTS) is a rare and clinically heterogeneous inherited disorder characterized by a long QT interval on the electrocardiogram, increased risk of syncope and sudden death caused by arrhythmias. This syndrome is mostly caused by mutations in genes encoding various cardiac ion channels. The clinical heterogeneity is usually attributed to variable penetrance. One of the reasons for this variability in expression could be the coexistence of common single nucleotide polymorphisms (SNPs) on LQTS‐causing genes and/or unknown genes. Some synonymous and nonsynonymous exonic SNPs identified in LQTS‐causing genes may have an effect on the cardiac repolarization process and modulate the clinical expression of a latent LQTS pathogenic mutation. We report the molecular pattern of 44 unrelated patients with LQTS using denaturing high‐performance liquid chromatography analysis of the KCNQ1, KCNH2, SCN5A, KCNE1 and KCNE2 genes. Forty‐five disease‐causing mutations (including 24 novel ones) were identified in this cohort. Most of our patients (84%) showed complex molecular pattern with one mutation (and even two for four patients) associated with several SNPs located in several LQTS genes.

title: Hypokalemia‐Induced Long QT Syndrome with an Underlying Novel Missense Mutation in S4‐S5 Linker of KCNQ1
abstract: Novel Mutation in Hypokalemia‐Induced LQTS. Congenital long QT syndrome (LQTS) is caused by mutations in at least five genes coding for cardiac potassium or sodium channels that regulate the duration of ventricular action potentials. Acquired LQTS often is associated with drugs or metabolic abnormalities. A 47‐year‐old woman who presented with marked QT prolongation (QTc = 620 msec1/2) and repeated episodes of torsades de pointes associated with hypokalemia (2.6 mEq/L) was screened for mutations in LQTS genes using polymerase chain reaction/single‐strand conformation polymorphism (PCR/SSCP). We identified a novel missense mutation in the intracellular linker of S4‐S5 domains of KCNQ1, resulting in an amino acid substitution of cysteine for arginine at position 259 (R259C). Whole cell, patch clamp experiments were conducted on COS7 cells transfected with wild‐type and/or R259C KCNQ1 with or without KCNE1. Functional analyses of the mutant KCNQ1 subunit on COS7 cells revealed its functional channels in the homozygous state, producing a significantly smaller current than the KCNQ1 channels and a less severe dominant‐negative effect on IKS‐ The novel KCNQl mutation R259C is the molecular basis for IKs dysfunction underlying an apparently sporadic case of hypokalemia‐induced LQTS, consistent with a mild mutation likely to disclose the clinical manifestation of LQTS in a context of severe hypokalemia. Our findings suggest that gene carriers with such mild mutations might not he so rare as commonly expected in patients with acquired LQTS, and stress the importance of mutational analysis for detecting either “silent” forms of congenital LQTS or de novo mutations.

title: Spectrum and Frequency of Cardiac Channel Defects in Swimming-Triggered Arrhythmia Syndromes
abstract: Swimming is a relatively genotype-specific arrhythmogenic trigger for type 1 long-QT syndrome (LQT1). We hypothesize that mimickers of concealed LQT1, namely catecholaminergic polymorphic ventricular tachycardia (CPVT), may also underlie swimming-triggered cardiac events. Between August 1997 and May 2003, 388 consecutive, unrelated patients were referred specifically for LQTS genetic testing. The presence of a personal and/or family history of a near-drowning or drowning was determined by review of the medical records and/or phone interviews and was blinded to genetic test results. Comprehensive mutational analysis of the 5 LQTS-causing channel genes, KCNQ1 (LQT1), KCNH2 (LQT2), SCN5A (LQT3), KCNE1 (LQT5), and KCNE2 (LQT6), along with KCNJ2 (Andersen-Tawil syndrome) and targeted analysis of 18 CPVT1-associated exons in RyR2, was performed with the use of denaturing high-performance liquid chromatography and direct DNA sequencing. Approximately 11% (43 of 388) of the index cases had a positive swimming phenotype. Thirty-three of these 43 index cases had a “Schwartz” score (≥4) suggesting high clinical probability of LQTS. Among this subset, 28 patients (85%) were LQT1, 2 patients (6%) were LQT2, and 3 were genotype negative. Among the 10 cases with low clinical probability for LQTS, 9 had novel, putative CPVT1-causing RyR2 mutations. In contrast to previous studies that suggested universal LQT1 specificity, genetic heterogeneity underlies channelopathies that are suspected chiefly because of a near-drowning or drowning. CPVT1 and strategic genotyping of RyR2 should be considered when LQT1 is excluded in the pathogenesis of a swimming-triggered arrhythmia syndrome.

title: Ion conduction pore is conserved among potassium channels
abstract: Potassium channels, a group of specialized membrane proteins, enable K+ ions to flow selectively across cell membranes. Transmembrane K+ currents underlie electrical signalling in neurons and other excitable cells. The atomic structure of a bacterial K+ channel pore has been solved by means of X-ray crystallography. To the extent that the prokaryotic pore is representative of other K+ channels, this landmark achievement has profound implications for our general understanding of K+ channels. But serious doubts have been raised concerning whether the prokaryotic K+ channel pore does actually represent those of eukaryotes. Here we have addressed this fundamental issue by substituting the prokaryotic pore into eukaryotic voltage-gated and inward-rectifier K+ channels. The resulting chimaeras retain the respective functional hallmarks of the eukaryotic channels, which indicates that the ion conduction pore is indeed conserved among K+ channels.

title: Voltage Sensor of Kv1.2: Structural Basis of Electromechanical Coupling
abstract: Voltage-dependent ion channels contain voltage sensors that allow them to switch between nonconductive and conductive states over the narrow range of a few hundredths of a volt. We investigated the mechanism by which these channels sense cell membrane voltage by determining the x-ray crystal structure of a mammalian Shaker family potassium ion (K+) channel. The voltage-dependent K+ channel Kv1.2 grew three-dimensional crystals, with an internal arrangement that left the voltage sensors in an apparently native conformation, allowing us to reach three important conclusions. First, the voltage sensors are essentially independent domains inside the membrane. Second, they perform mechanical work on the pore through the S4-S5 linker helices, which are positioned to constrict or dilate the S6 inner helices of the pore. Third, in the open conformation, two of the four conserved Arg residues on S4 are on a lipid-facing surface and two are buried in the voltage sensor. The structure offers a simple picture of how membrane voltage influences the open probability of the channel.

title: Coupling between Voltage Sensors and Activation Gate in Voltage-gated K+ Channels
abstract: Current through voltage-gated K+ channels underlies the action potential encoding the electrical signal in excitable cells. The four subunits of a voltage-gated K+ channel each have six transmembrane segments (S1–S6), whereas some other K+ channels, such as eukaryotic inward rectifier K+ channels and the prokaryotic KcsA channel, have only two transmembrane segments (M1 and M2). A voltage-gated K+ channel is formed by an ion-pore module (S5–S6, equivalent to M1–M2) and the surrounding voltage-sensing modules. The S4 segments are the primary voltage sensors while the intracellular activation gate is located near the COOH-terminal end of S6, although the coupling mechanism between them remains unknown. In the present study, we found that two short, complementary sequences in voltage-gated K+ channels are essential for coupling the voltage sensors to the intracellular activation gate. One sequence is the so called S4–S5 linker distal to the voltage-sensing S4, while the other is around the COOH-terminal end of S6, a region containing the actual gate-forming residues.

title: The S4-S5 Linker of KCNQ1 Channels Forms a Structural Scaffold with the S6 Segment Controlling Gate Closure*
abstract: In vivo, KCNQ1 α-subunits associate with the β-subunit KCNE1 to generate the slowly activating cardiac potassium current (IKs). Structurally, they share their topology with other Kv channels and consist out of six transmembrane helices (S1–S6) with the S1–S4 segments forming the voltage-sensing domain (VSD). The opening or closure of the intracellular channel gate, which localizes at the bottom of the S6 segment, is directly controlled by the movement of the VSD via an electromechanical coupling. In other Kv channels, this electromechanical coupling is realized by an interaction between the S4-S5 linker (S4S5L) and the C-terminal end of S6 (S6T). Previously we reported that substitutions for Leu353 in S6T resulted in channels that failed to close completely. Closure could be incomplete because Leu353 itself is the pore-occluding residue of the channel gate or because of a distorted electromechanical coupling. To resolve this and to address the role of S4S5L in KCNQ1 channel gating, we performed an alanine/tryptophan substitution scan of S4S5L. The residues with a “high impact” on channel gating (when mutated) clustered on one side of the S4S5L α-helix. Hence, this side of S4S5L most likely contributes to the electromechanical coupling and finds its residue counterparts in S6T. Accordingly, substitutions for Val254 resulted in channels that were partially constitutively open and the ability to close completely was rescued by combination with substitutions for Leu353 in S6T. Double mutant cycle analysis supported this cross-talk indicating that both residues come in close contact and stabilize the closed state of the channel.

title: Role of the S6 C‐terminus in KCNQ1 channel gating
abstract: Co‐assembly of KCNQ1 α‐subunits with KCNE1 β‐subunits results in the channel complex underlying the cardiac IKs current in vivo. Like other voltage‐gated K+ channels, KCNQ1 has a tetrameric configuration. The S6 segment of each subunit lines the ion channel pore with the lower part forming the activation gate. To determine residues involved in protein–protein interactions in the C‐terminal part of S6 (S6T), alanine and tryptophan perturbation scans were performed from residue 348–362 in the KCNQ1 channel. Several residues were identified to be relevant in channel gating, as substitutions affected the activation and/or deactivation process. Some mutations (F351A and V355W) drastically altered the gating characteristics of the resultant KCNQ1 channel, to the point of mimicking the IKs current. Furthermore, mutagenesis of residue L353 to an alanine or a charged residue impaired normal channel closure upon hyperpolarization, generating a constitutively open phenotype. This indicates that the L353 residue is essential for stabilizing the closed conformation of the channel gate. These findings together with the identification of several LQT1 mutations in the S6 C‐terminus of KCNQ1 underscore the relevance of this region in KCNQ1 and IKs channel gating.

title: Allosteric features of KCNQ1 gating revealed by alanine scanning mutagenesis.
abstract: Controlled opening and closing of an ion-selective pathway in response to changes of membrane potential is a fundamental feature of voltage-gated ion channels. In recent decades, various details of this process have been revealed with unprecedented precision based on studies of prototypic potassium channels. Though current scientific efforts are focused more on a thorough description of voltage-sensor movement, much less is known about the similarities and differences of the gating mechanisms among potassium channels. Here, we describe the peculiarities of the KCNQ1 gating process in parallel comparison to Shaker. We applied alanine scanning mutagenesis to the S4-S5 linker and pore region and followed the regularities of gating perturbations in KCNQ1. We found a fractional constitutive conductance for wild-type KCNQ1. This component increased significantly in mutants with considerably leftward-shifted steady-state activation curves. In contrast to Shaker, no correlation between V(1/2) and Z parameters was observed for the voltage-dependent fraction of KCNQ1. Our experimental findings are explained by a simple allosteric gating scheme with voltage-driven and voltage-independent transitions. Allosteric features are discussed in the context of extreme gating adaptability of KCNQ1 upon interaction with KCNE β-subunits.

title: Inactivation of KCNQ1 potassium channels reveals dynamic coupling between voltage sensing and pore opening
abstract: In voltage-activated ion channels, voltage sensor (VSD) activation induces pore opening via VSD-pore coupling. Previous studies show that the pore in KCNQ1 channels opens when the VSD activates to both intermediate and fully activated states, resulting in the intermediate open (IO) and activated open (AO) states, respectively. It is also well known that accompanying KCNQ1 channel opening, the ionic current is suppressed by a rapid process called inactivation. Here we show that inactivation of KCNQ1 channels derives from the different mechanisms of the VSD-pore coupling that lead to the IO and AO states, respectively. When the VSD activates from the intermediate state to the activated state, the VSD-pore coupling has less efficacy in opening the pore, producing inactivation. These results indicate that different mechanisms, other than the canonical VSD-pore coupling, are at work in voltage-dependent ion channel activation.KCNQ1 is a voltage-gated potassium channel that is important in cardiac and epithelial function. Here the authors present a mechanism for KCNQ1 activation and inactivation in which voltage sensor activation promotes pore opening more effectively in the intermediate open state than the fully open state, generating inactivation.

title: Microscopic mechanisms for long QT syndrome type 1 revealed by single-channel analysis of I(Ks) with S3 domain mutations in KCNQ1.
abstract: The slowly activating delayed rectifier current IKs participates in cardiac repolarization, particularly at high heart rates, and mutations in this K(+) channel complex underlie long QT syndrome (LQTS) types 1 and 5.



The purpose of this study was to determine biophysical mechanisms of LQT1 through single-channel kinetic analysis of IKs carrying LQT1 mutations in the S3 transmembrane region of the pore-forming subunit KCNQ1.



We analyzed cell-attached recordings from mammalian cells in which a single active KCNQ1 (wild type or mutant) and KCNE1 complex could be detected.



The S3 mutants of KCNQ1 studied (D202H, I204F, V205M, and S209F), with the exception of S209F, all led to a reduction in channel activity through distinct kinetic mechanisms. D202H, I204F, and V205M showed decreased open probability (Po) compared with wild type (0.07, 0.04, and 0.12 vs 0.2); increased first latency from 1.66 to >2 seconds at +60 mV (I204F, V205M); variable-to-severe reductions in open dwell times (≥50% in V205M); stabilization of closed states (D202H); and an inability of channels to reach full conductance levels (V205M, I204F). S209F is a kinetic gain-of-function mutation with a high Po (0.40) and long open-state dwell times.



S3 mutations in KCNQ1 cause diverse kinetic defects in I(Ks), affecting opening and closing properties, and can account for LQT1 phenotypes.

title: High Throughput Functional Evaluation of KCNQ1 Decrypts Variants of Unknown Significance
abstract: The explosive growth in known human gene variation presents enormous challenges to current approaches for variant classification that impact diagnosis and treatment of many genetic diseases. For disorders caused by mutations in cardiac ion channels, such as congenital long-QT syndrome (LQTS), in vitro electrophysiological evidence has high value in discriminating pathogenic from benign variants, but these data are often lacking because assays are cost-, time- and labor-intensive.  We implemented a strategy for performing high throughput, functional evaluations of ion channel variants that repurposed an automated electrophysiology platform developed previously for drug discovery. We demonstrated success of this approach by evaluating 78 variants in KCNQ1, a major LQTS gene. We benchmarked our results with traditional electrophysiological approaches and observed a high level of concordance. Our results provided functional data useful for classifying ~70% of previously unstudied KCNQ1 variants annotated with uninformative descriptions in the public database ClinVar. Further, we show that rare and ultra-rare KCNQ1 variants in the general population exhibit functional properties ranging from normal to severe loss-of-function indicating that allele frequency is not a reliable predictor of channel function.  Our results illustrate an efficient and high throughput paradigm linking genotype to function for a human cardiac channelopathy that will enable data-driven classification of large numbers of variants and create new opportunities for precision medicine.

title: Predicting the Functional Impact of KCNQ1 Variants of Unknown Significance
abstract: An emerging standard-of-care for long-QT syndrome uses clinical genetic testing to identify genetic variants of the KCNQ1 potassium channel. However, interpreting results from genetic testing is confounded by the presence of variants of unknown significance for which there is inadequate evidence of pathogenicity.  In this study, we curated from the literature a high-quality set of 107 functionally characterized KCNQ1 variants. Based on this data set, we completed a detailed quantitative analysis on the sequence conservation patterns of subdomains of KCNQ1 and the distribution of pathogenic variants therein. We found that conserved subdomains generally are critical for channel function and are enriched with dysfunctional variants. Using this experimentally validated data set, we trained a neural network, designated Q1VarPred, specifically for predicting the functional impact of KCNQ1 variants of unknown significance. The estimated predictive performance of Q1VarPred in terms of Matthew’s correlation coefficient and area under the receiver operating characteristic curve were 0.581 and 0.884, respectively, superior to the performance of 8 previous methods tested in parallel. Q1VarPred is publicly available as a web server at http://meilerlab.org/q1varpred.  Although a plethora of tools are available for making pathogenicity predictions over a genome-wide scale, previous tools fail to perform in a robust manner when applied to KCNQ1. The contrasting and favorable results for Q1VarPred suggest a promising approach, where a machine-learning algorithm is tailored to a specific protein target and trained with a functionally validated data set to calibrate informatics tools.

title: Mechanisms of KCNQ1 Channel Dysfunction in Long QT Syndrome Involving Voltage Sensor Domain Mutations
abstract: Loss-of-function (LOF) mutations in human KCNQ1 are responsible for susceptibility to a life-threatening heart rhythm disorder, the congenital long-QT syndrome (LQTS). Hundreds of KCNQ1 mutations have been identified, but the molecular mechanisms responsible for impaired function are poorly understood. Here, we investigated the impact of 51 KCNQ1 variants located within the voltage sensor domain (VSD), with an emphasis on elucidating effects on cell surface expression, protein folding and structure. For each variant, the efficiency of trafficking to the plasma membrane, the impact of proteasome inhibition, and protein stability were assayed. The results of these experiments, combined with channel functional data, provided the basis for classifying each mutation into one of 6 mechanistic categories. More than half of the KCNQ1 LOF mutations destabilize the structure of the VSD, resulting in mistrafficking and degradation by the proteasome, an observation that underscores the growing appreciation that mutation-induced destabilization of membrane proteins may be a common human disease mechanism. Finally, we observed that 5 of the folding-defective LQTS mutants are located in the VSD S0 helix, where they interact with a number of other LOF mutation sites in other segments of the VSD. These observations reveal a critical role for the S0 helix as a central scaffold to help organize and stabilize the KCNQ1 VSD and, most likely, the corresponding domain of many other ion channels. One Sentence  Long QT syndrome-associated mutations in KCNQ1 most often destabilize the protein, leading to mistrafficking and degradation.

title: Single channel kinetic analysis of the cAMP effect on IKs mutants, S209F and S27D/S92D
abstract: The IKs current is important in the heart’s response to sympathetic stimulation. β-adrenergic stimulation increases the amount of IKs and creates a repolarization reserve that shortens the cardiac action potential duration. We have recently shown that 8-CPT-cAMP, a membrane-permeable cAMP analog, changes the channel kinetics and causes it to open more quickly and more often, as well as to higher subconductance levels, which produces an increase in the IKs current. The mechanism proposed to underlie these kinetic changes is increased activation of the voltage sensors. The present study extends our previous work and shows detailed subconductance analysis of the effects of 8-CPT-cAMP on an enhanced gating mutant (S209F) and on a double pseudo-phosphorylated IKs channel (S27D/S92D). 8-CPT-cAMP still produced kinetic changes in S209F + KCNE1, further enhancing gating, while S27D/S92D + KCNE1 showed no significant response to 8-CPT-cAMP, suggesting that these last two mutations fully recapitulate the effect of channel phosphorylation by cAMP.

title: Single-channel basis for the slow activation of the repolarizing cardiac potassium current, IKs
abstract: Coassembly of potassium voltage-gated channel, KQT-like subfamily, member 1 (KCNQ1) with potassium voltage-gated channel, Isk-related family, member 1 (KCNE1) the delayed rectifier potassium channel IKs. Its slow activation is critically important for membrane repolarization and for abbreviating the cardiac action potential, especially during sympathetic activation and at high heart rates. Mutations in either gene can cause long QT syndrome, which can lead to fatal arrhythmias. To understand better the elementary behavior of this slowly activating channel complex, we quantitatively analyzed direct measurements of single-channel IKs. Single-channel recordings from transiently transfected mouse ltk− cells confirm a channel that has long latency periods to opening (1.67 ± 0.073 s at +60 mV) but that flickers rapidly between multiple open and closed states in non-deactivating bursts at positive membrane potentials. Channel activity is cyclic with periods of high activity followed by quiescence, leading to an overall open probability of only ∼0.15 after 4 s under our recording conditions. The mean single-channel conductance was determined to be 3.2 pS, but unlike any other known wild-type human potassium channel, long-lived subconductance levels coupled to activation are a key feature of both the activation and deactivation time courses of the conducting channel complex. Up to five conducting levels ranging from 0.13 to 0.66 pA could be identified in single-channel recordings at 60 mV. Fast closings and overt subconductance behavior of the wild-type IKs channel required modification of existing Markov models to include these features of channel behavior.

title: A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential
abstract: Ion-channel function is determined by its gating movement. Yet, molecular dynamics and electrophysiological simulations were never combined to link molecular structure to function. We performed multiscale molecular dynamics and continuum electrostatics calculations to simulate a cardiac K+ channel (IKs) gating and its alteration by mutations that cause arrhythmias and sudden death. An all-atom model of the IKs α-subunit KCNQ1, based on the recent Kv1.2 structure, is used to calculate electrostatic energies during gating. Simulations are compared with experiments where varying degrees of positive charge—added via point mutation—progressively reduce current. Whole-cell simulations show that mutations cause action potential and ECG QT interval prolongation, consistent with clinical phenotypes. This framework allows integration of multiscale observations to study the molecular basis of excitation and its alteration by disease.

title: Flex ddG: Rosetta ensemble-based estimation of changes in protein-protein binding affinity upon mutation
abstract: Computationally modeling changes in binding free energies upon mutation (interface ΔΔG) allows large-scale prediction and perturbation of protein-protein interactions. Additionally, methods that consider and sample relevant conformational plasticity should be able to achieve higher prediction accuracy over methods that do not. To test this hypothesis, we developed a method within the Rosetta macromolecular modeling suite (flex ddG) that samples conformational diversity using “backrub” to generate an ensemble of models, then applying torsion minimization, side chain repacking and averaging across this ensemble to estimate interface ΔΔG values. We tested our method on a curated benchmark set of 1240 mutants, and found the method outperformed existing methods that sampled conformational space to a lesser degree. We observed considerable improvements with flex ddG over existing methods on the subset of small side chain to large side chain mutations, as well as for multiple simultaneous non-alanine mutations, stabilizing mutations, and mutations in antibody-antigen interfaces. Finally, we applied a generalized additive model (GAM) approach to the Rosetta energy function; the resulting non-linear reweighting model improved agreement with experimentally determined interface DDG values, but also highlights the necessity of future energy function improvements.

title: The backrub motion: how protein backbone shrugs when a sidechain dances.
abstract: Surprisingly, the frozen structures from ultra-high-resolution protein crystallography reveal a prevalent, but subtle, mode of local backbone motion coupled to much larger, two-state changes of sidechain conformation. This "backrub" motion provides an influential and common type of local plasticity in protein backbone. Concerted reorientation of two adjacent peptides swings the central sidechain perpendicular to the chain direction, changing accessible sidechain conformations while leaving flanking structure undisturbed. Alternate conformations in sub-1 angstroms crystal structures show backrub motions for two-thirds of the significant Cbeta shifts and 3% of the total residues in these proteins (126/3882), accompanied by two-state changes in sidechain rotamer. The Backrub modeling tool is effective in crystallographic rebuilding. For homology modeling or protein redesign, backrubs can provide realistic, small perturbations to rigid backbones. For large sidechain changes in protein dynamics or for single mutations, backrubs allow backbone accommodation while maintaining H bonds and ideal geometry.

title: An Integrated Framework Advancing Membrane Protein Modeling and Design
abstract: Membrane proteins are critical functional molecules in the human body, constituting more than 30% of open reading frames in the human genome. Unfortunately, a myriad of difficulties in overexpression and reconstitution into membrane mimetics severely limit our ability to determine their structures. Computational tools are therefore instrumental to membrane protein structure prediction, consequently increasing our understanding of membrane protein function and their role in disease. Here, we describe a general framework facilitating membrane protein modeling and design that combines the scientific principles for membrane protein modeling with the flexible software architecture of Rosetta3. This new framework, called RosettaMP, provides a general membrane representation that interfaces with scoring, conformational sampling, and mutation routines that can be easily combined to create new protocols. To demonstrate the capabilities of this implementation, we developed four proof-of-concept applications for (1) prediction of free energy changes upon mutation; (2) high-resolution structural refinement; (3) protein-protein docking; and (4) assembly of symmetric protein complexes, all in the membrane environment. Preliminary data show that these algorithms can produce meaningful scores and structures. The data also suggest needed improvements to both sampling routines and score functions. Importantly, the applications collectively demonstrate the potential of combining the flexible nature of RosettaMP with the power of Rosetta algorithms to facilitate membrane protein modeling and design.

title: Coassembly of K(V)LQT1 and minK (IsK) proteins to form cardiac I(Ks) potassium channel.
abstract: The slowly activating delayed-rectifier K+ current, I(Ks), modulates the repolarization of cardiac action potentials. The molecular structure of the I(Ks) channel is not known, but physiological data indicate that one component of the I(Ks), channel is minK, a 130-amino-acid protein with a single putative transmembrane domain. The size and structure of this protein is such that it is unlikely that minK alone forms functional channels. We have previously used positional cloning techniques to define a new putative K+-channel gene, KVLQT1. Mutations in this gene cause long-QT syndrome, an inherited disorder that increases the risk of sudden death from cardiac arrhythmias. Here we show that KVLQT1 encodes a K+ channel with biophysical properties unlike other known cardiac currents. We considered that K(V)LQT1 might coassemble with another subunit to form functional channels in cardiac myocytes. Coexpression of K(V)LQT1 with minK induced a current that was almost identical to cardiac I(Ks). Therefore, K(V)LQT1 is the subunit that coassembles with minK to form I(Ks) channels and I(Ks) dysfunction is a cause of cardiac arrhythmia.

title: KCNE2 confers background current characteristics to the cardiac KCNQ1 potassium channel
abstract: Mutations in HERG and KCNQ1 (or KVLQT1) genes cause the life‐threatening Long QT syndrome. These genes encode K+ channel pore‐forming subunits that associate with ancillary subunits from the KCNE family to underlie the two components, IKr and IKs, of the human cardiac delayed rectifier current IK. The KCNE family comprises at least three members. KCNE1 (IsK or MinK) recapitulates IKs when associated with KCNQ1, whereas it augments the amplitude of an IKr‐like current when co‐expressed with HERG. KCNE3 markedly changes KCNQ1 as well as HERG current properties. So far, KCNE2 (MirP1) has only been shown to modulate HERG current. Here we demonstrate the interaction of KCNE2 with the KCNQ1 subunit, which results in a drastic change of KCNQ1 current amplitude and gating properties. Furthermore, KCNE2 mutations also reveal their specific functional consequences on KCNQ1 currents. KCNQ1 and HERG appear to share unique interactions with KCNE1, 2 and 3 subunits. With the exception of KCNE3, mutations in all these partner subunits have been found to lead to an increased propensity for cardiac arrhythmias.

title: A constitutively open potassium channel formed by KCNQ1 and KCNE3
abstract: Mutations in all four known KCNQ potassium channel α-subunit genes lead to human diseases. KCNQ1 (KvLQT1) interacts with the β-subunit KCNE1 (IsK, minK) to form the slow, depolarization-activated potassium current IKs that is affected in some forms of cardiac arrhythmia. Here we show that the novel β-subunit KCNE3 markedly changes KCNQ1 properties to yield currents that are nearly instantaneous and depend linearly on voltage. It also suppresses the currents of KCNQ4 and HERG potassium channels. In the intestine, KCNQ1 and KCNE3 messenger RNAs colocalized in crypt cells. This localization and the pharmacology, voltage-dependence and stimulation by cyclic AMP of KCNQ1/KCNE3 currents indicate that these proteins may assemble to form the potassium channel that is important for cyclic AMP-stimulated intestinal chloride secretion and that is involved in secretory diarrhoea and cystic fibrosis.

title: KCNQ1 channel modulation by KCNE proteins via the voltage‐sensing domain
abstract: The gating of the KCNQ1 potassium channel is drastically regulated by auxiliary subunit KCNE proteins. KCNE1, for example, slows the activation kinetics of KCNQ1 by two orders of magnitude. Like other voltage‐gated ion channels, the opening of KCNQ1 is regulated by the voltage‐sensing domain (VSD; S1–S4 segments). Although it has been known that KCNE proteins interact with KCNQ1 via the pore domain, some recent reports suggest that the VSD movement may be altered by KCNE. The altered VSD movement of KCNQ1 by KCNE proteins has been examined by site‐directed mutagenesis, the scanning cysteine accessibility method (SCAM), voltage clamp fluorometry (VCF) and gating charge measurements. These accumulated data support the idea that KCNE proteins interact with the VSDs of KCNQ1 and modulate the gating of the KCNQ1 channel. In this review, we will summarize recent findings and current views of the KCNQ1 modulation by KCNE via the VSD. In this context, we discuss our recent findings that KCNE1 may alter physical interactions between the S4 segment (VSD) and the S5 segment (pore domain) of KCNQ1. Based on these findings from ourselves and others, we propose a hypothetical mechanism for how KCNE1 binding alters the VSD movement and the gating of the channel.

title: IKs channels open slowly because KCNE1 accessory subunits slow the movement of S4 voltage sensors in KCNQ1 pore-forming subunits
abstract: Significance E1 and Q1 protein subunits assemble to form IKslow channels in the heart and ear. Inherited mutations in either subunit that decrease protein level or alter function can cause life-threatening cardiac arrhythmias and deafness. The mechanism by which E1 slows channel opening has been the subject of active debate. Here, we use gating current measurements and simultaneous recordings of ionic currents and changes in fluorescence of a probe on the Q1 voltage sensors to demonstrate that E1 slows the movement of sensors in a manner that is both necessary and sufficient to determine the slow activation time course of IKs channels. Human IKs channels activate slowly with the onset of cardiac action potentials to repolarize the myocardium. IKs channels are composed of KCNQ1 (Q1) pore-forming subunits that carry S4 voltage-sensor segments and KCNE1 (E1) accessory subunits. Together, Q1 and E1 subunits recapitulate the conductive and kinetic properties of IKs. How E1 modulates Q1 has been unclear. Investigators have variously posited that E1 slows the movement of S4 segments, slows opening and closing of the conduction pore, or modifies both aspects of electromechanical coupling. Here, we show that Q1 gating current can be resolved in the absence of E1, but not in its presence, consistent with slowed movement of the voltage sensor. E1 was directly demonstrated to slow S4 movement with a fluorescent probe on the Q1 voltage sensor. Direct correlation of the kinetics of S4 motion and ionic current indicated that slowing of sensor movement by E1 was both necessary and sufficient to determine the slow-activation time course of IKs.

title: KCNE1 and KCNE3 modulate KCNQ1 channels by affecting different gating transitions
abstract: Regulatory β-subunits associate with voltage-gated K+ channels to modulate their biophysical properties and physiological roles. KCNE1 and KCNE3 β-subunits turn voltage-dependent KCNQ1 channels into delayed activating KCNQ1/KCNE1 channels and apparent voltage-independent KCNQ1/KCNE3 channels, respectively, which are important for cardiac action potentials and transport of water and salts across epithelial cells. Mutations in KCNQ1/KCNE1 and KCNQ1/KCNE3 channels are associated with diseases, such as cardiac arrhythmias, congenital deafness, secretory diarrhea, and tinnitus. Therefore, KCNQ1, KCNE1, and KCNE3 are potential drug targets. We here propose a model for how KCNE1 and KCNE3 differentially modulate KCNQ1 that will allow for a better understanding of how mutations in KCNQ1, KCNE1, and KCNE3 cause diseases and how to design drugs to treat these diseases. KCNE β-subunits assemble with and modulate the properties of voltage-gated K+ channels. In the heart, KCNE1 associates with the α-subunit KCNQ1 to generate the slowly activating, voltage-dependent potassium current (IKs) in the heart that controls the repolarization phase of cardiac action potentials. By contrast, in epithelial cells from the colon, stomach, and kidney, KCNE3 coassembles with KCNQ1 to form K+ channels that are voltage-independent K+ channels in the physiological voltage range and important for controlling water and salt secretion and absorption.

title: IKs ion-channel pore conductance can result from individual voltage sensor movements
abstract: A potassium current in the heart (IKs) is important in regulating the heartbeat. These channels, composed of KCNQ1 and KCNE1 subunits, have four voltage-sensitive components and a pore component, which opens to allow current to pass through. The number of voltage-sensitive components required to activate to allow the pore to open has been debated extensively. While some believe all four are needed, others believe the pore can open with movement of a single voltage-sensitive component. Here, we show that the latter idea is correct, with movement of each component resulting in increasingly more current. This makes for a uniquely flexible channel, tunable for a broad range of responses. The IKs current has an established role in cardiac action potential repolarization, and provides a repolarization reserve at times of stress. The underlying channels are formed from tetramers of KCNQ1 along with one to four KCNE1 accessory subunits, but how these components together gate the IKs complex to open the pore is controversial. Currently, either a concerted movement involving all four subunits of the tetramer or allosteric regulation of open probability through voltage-dependent subunit activation is thought to precede opening. Here, by using the E160R mutation in KCNQ1 or the F57W mutation in KCNE1 to prevent or impede, respectively, voltage sensors from moving into activated conformations, we demonstrate that a concerted transition of all four subunits after voltage sensor activation is not required for the opening of IKs channels. Tracking voltage sensor movement, via [2-(trimethylammonium)ethyl]methanethiosulfonate bromide (MTSET) modification and fluorescence recordings, shows that E160R-containing voltage sensors do not translocate upon depolarization. E160R, when expressed in all four KCNQ1 subunits, is nonconducting, but if one, two, or three voltage sensors contain the E160R mutation, whole-cell and single-channel currents are still observed in both the presence and absence of KCNE1, and average conductance is reduced proportional to the number of E160R voltage sensors. The data suggest that KCNQ1 + KCNE1 channels gate like KCNQ1 alone. A model of independent voltage sensors directly coupled to open states can simulate experimental changes in IKs current kinetics, including the nonlinear depolarization of the conductance–voltage (G–V) relationship, and tail current acceleration as the number of nonactivatable E160R subunits is increased.

title: Location and orientation of minK within the I(Ks) potassium channel complex.
abstract: The slowly activating cardiac potassium current (I(Ks)) is generated by a heteromultimeric potassium channel complex consisting of pore-forming (KvLQT1) and accessory (minK) subunits belonging to the KCNQ and KCNE gene families, respectively. Evidence indicating that minK residues line the I(Ks) pore originates from the observation that two minK cysteine mutants (G55C and F54C) render I(Ks) Cd2+-sensitive. We have identified a single cysteine residue in the KvLQT1 S6 segment (Cys-331) that contributes to Cd2+ coordination in conjunction with cysteine residues engineered into the minK transmembrane domain. This observation indicates that minK resides in close proximity to S6 in the I(Ks) channel complex. On the basis of homology modeling that compares the KvLQT1 S6 segment with the structure of the bacterial potassium channel KcsA, we predict that the sulfhydryl side chain of Cys-331 projects away from the central axis of the KvLQT1 pore and suggest that minK resides outside of the permeation pathway. A preliminary model illustrating the orientation of minK with S6 was validated by successful prediction of a novel Cd2+ binding site created within the I(Ks) channel complex by engineering additional cysteine residues into both subunits. Our results indicate the location and orientation of minK within the I(Ks) channel complex and further suggest that Cd2+ exerts its effect on I(Ks) through an allosteric mechanism rather than direct pore blockade.

title: Location of KCNE1 relative to KCNQ1 in the IKS potassium channel by disulfide cross-linking of substituted cysteines
abstract: The cardiac-delayed rectifier K+ current (IKS) is carried by a complex of KCNQ1 (Q1) subunits, containing the voltage-sensor domains and the pore, and auxiliary KCNE1 (E1) subunits, required for the characteristic IKS voltage dependence and kinetics. To locate the transmembrane helix of E1 (E1-TM) relative to the Q1 TM helices (S1–S6), we mutated, one at a time, the first four residues flanking the extracellular ends of S1–S6 and E1-TM to Cys, coexpressed all combinations of Q1 and E1 Cys-substituted mutants in CHO cells, and determined the extents of spontaneous disulfide-bond formation. Cys-flanking E1-TM readily formed disulfides with Cys-flanking S1 and S6, much less so with the S3-S4 linker, and not at all with S2 or S5. These results imply that the extracellular flank of the E1-TM is located between S1 and S6 on different subunits of Q1. The salient functional effects of selected cross-links were as follows. A disulfide from E1 K41C to S1 I145C strongly slowed deactivation, and one from E1 L42C to S6 V324C eliminated deactivation. Given that E1-TM is between S1 and S6 and that K41C and L42C are likely to point approximately oppositely, these two cross-links are likely to favor similar axial rotations of E1-TM. In the opposite orientation, a disulfide from E1 K41C to S6 V324C slightly slowed activation, and one from E1 L42C to S1 I145C slightly speeded deactivation. Thus, the first E1 orientation strongly favors the open state, while the approximately opposite orientation favors the closed state.

title: Characterization of KCNQ1 atrial fibrillation mutations reveals distinct dependence on KCNE1
abstract: The IKs potassium channel, critical to control of heart electrical activity, requires assembly of α (KCNQ1) and β (KCNE1) subunits. Inherited mutations in either IKs channel subunit are associated with cardiac arrhythmia syndromes. Two mutations (S140G and V141M) that cause familial atrial fibrillation (AF) are located on adjacent residues in the first membrane-spanning domain of KCNQ1, S1. These mutations impair the deactivation process, causing channels to appear constitutively open. Previous studies suggest that both mutant phenotypes require the presence of KCNE1. Here we found that despite the proximity of these two mutations in the primary protein structure, they display different functional dependence in the presence of KCNE1. In the absence of KCNE1, the S140G mutation, but not V141M, confers a pronounced slowing of channel deactivation and a hyperpolarizing shift in voltage-dependent activation. When coexpressed with KCNE1, both mutants deactivate significantly slower than wild-type KCNQ1/KCNE1 channels. The differential dependence on KCNE1 can be correlated with the physical proximity between these positions and KCNE1 as shown by disulfide cross-linking studies: V141C forms disulfide bonds with cysteine-substituted KCNE1 residues, whereas S140C does not. These results further our understanding of the structural relationship between KCNE1 and KCNQ1 subunits in the IKs channel, and provide mechanisms for understanding the effects on channel deactivation underlying these two atrial fibrillation mutations.

title: KCNQ1 and KCNE1 in the IKs Channel Complex Make State-dependent Contacts in their Extracellular Domains
abstract: KCNQ1 and KCNE1 (Q1 and E1) associate to form the slow delayed rectifier IKs channels in the heart. A short stretch of eight amino acids at the extracellular end of S1 in Q1 (positions 140–147) harbors six arrhythmia-associated mutations. Some of these mutations affect the Q1 channel function only when coexpressed with E1, suggesting that this Q1 region may engage in the interaction with E1 critical for the IKs channel function. Identifying the Q1/E1 contact points here may provide new insights into how the IKs channel operates. We focus on Q1 position 145 and E1 positions 40–43. Replacing all native cysteine (Cys) in Q1 and introducing Cys into the above Q1 and E1 positions do not significantly perturb the Q1 channel function or Q1/E1 interactions. Immunoblot experiments on COS-7 cells reveal that Q1 145C can form disulfide bonds with E1 40C and 41C, but not E1 42C or 43C. Correspondingly, voltage clamp experiments in oocytes reveal that Q1 145C coexpressed with E1 40C or E1 41C manifests unique gating behavior and DTT sensitivity. Our data suggest that E1 40C and 41C come close to Q1 145C in the activated and resting states, respectively, to allow disulfide bond formation. These data and those in the literature lead us to propose a structural model for the Q1/E1 channel complex, in which E1 is located between S1, S4, and S6 of three separate Q1 subunits. We propose that E1 is not a passive partner of the Q1 channel, but instead can engage in molecular motions during IKs gating.

title: Gating-Related Molecular Motions in the Extracellular Domain of the IKs Channel: Implications for IKs Channelopathy
abstract: Cardiac slow delayed rectifier (IKs) channel complex consists of KCNQ1 channel and KCNE1 auxiliary subunits. The extracellular juxtamembranous region of KCNE1 is an unstructured loop that contacts multiple KCNQ1 positions in a gating-state-dependent manner. Congenital arrhythmia-related mutations have been identified in the extracellular S1–S2 linker of KCNQ1. These mutations manifest abnormal phenotypes only when coexpressed with KCNE1, pointing to the importance of proper KCNQ1/KCNE1 interactions here in IKs channel function. We investigate the interactions between the KCNE1 loop (positions 36–47) and KCNQ1 S1–S2 linker (positions 140–148) by means of disulfide trapping and voltage clamp techniques. During transitions among the resting-state conformations, KCNE1 positions 36–43 make contacts with KCNQ1 positions 144, 145, and 147 in a parallel fashion. During conformational changes in the activated state, KCNE1 position 40 can make contacts with all three KCNQ1 positions, while the neighboring KCNE1 positions (36, 38, 39, and 41) can make contact with KCNQ1 position 147. Furthermore, KCNQ1 positions 143 and 146 are high-impact positions that cannot tolerate cysteine substitution. To maintain the proper IKs channel function, position 143 requires a small side chain with a hydroxyl group, and position 146 requires a negatively charged side chain. These data and the proposed molecular motions provide insights into the mechanisms by which mutations in the extracellular juxtamembranous region of the IKs channel impair its function.

title: Differential Modulations of KCNQ1 by Auxiliary Proteins KCNE1 and KCNE2
abstract: KCNQ1 channels play vital roles in cardiovascular, gastric and other systems. The conductance and dynamics of KCNQ1 could be modulated by different single transmembrane helical auxiliary proteins (such as KCNE1, KCNE2 and others). In this study, detail KCNQ1 function modulations by different regions of KCNE1 or KCNE2 were examined using combinational methods of electrophysiology, immunofluorescence, solution NMR and related backbone flexibility analysis. In the presence of KCNE2 N-terminus, decreased surface expression and consequent low activities of KCNQ1 were observed. The transmembrane domains (TMDs) of KCNE1 and KCNE2 were illustrated to associate with the KCNQ1 channel in different modes: Ile64 in KCNE2-TMD interacting with Phe340 and Phe275 in KCNQ1, while two pairs of interacting residues (Phe340-Thr58 and Ala244-Tyr65) in the KCNQ1/KCNE1 complex. The KCNE1 C-terminus could modulate gating property of KCNQ1, whereas KCNE2 C-terminus had only minimal influences on KCNQ1. All of the results demonstrated different KCNQ1 function modulations by different regions of the two auxiliary proteins.

title: Structural Basis of Slow Activation Gating in the Cardiac IKs Channel Complex
abstract: Accessory β-subunits of the KCNE gene family modulate the function of various cation channel α-subunits by the formation of heteromultimers. Among the most dramatic changes of biophysical properties of a voltage-gated channel by KCNEs are the effects of KCNE1 on KCNQ1 channels. KCNQ1 and KCNE1 are believed to form nativeIKs channels. Here, we characterize molecular determinants of KCNE1 interaction with KCNQ1 channels by scanning mutagenesis, double mutant cycle analysis, and molecular dynamics simulations. Our findings suggest that KCNE1 binds to the outer face of the KCNQ1 channel pore domain, modifies interactions between voltage sensor, S4-S5 linker and the pore domain, leading to structural modifications of the selectivity filter and voltage sensor domain. Molecular dynamics simulations suggest a stable interaction of the KCNE1 transmembrane α-helix with the pore domain S5/S6 and part of the voltage sensor domain S4 of KCNQ1 in a putative pre-open channel state. Formation of this state may induce slow activation gating, the pivotal characteristic of native cardiac IKs channels. This new KCNQ1-KCNE1 model may become useful for dynamic modeling of disease-associated mutant IKs channels.

title: Probing the structural basis for differential KCNQ1 modulation by KCNE1 and KCNE2
abstract: KCNE1 associates with KCNQ1 to increase its current amplitude and slow the activation gating process, creating the slow delayed rectifier channel that functions as a “repolarization reserve” in human heart. The transmembrane domain (TMD) of KCNE1 plays a key role in modulating KCNQ1 pore conductance and gating kinetics, and the extracellular juxtamembrane (EJM) region plays a modulatory role by interacting with the extracellular surface of KCNQ1. KCNE2 is also expressed in human heart and can associate with KCNQ1 to suppress its current amplitude and slow the deactivation gating process. KCNE1 and KCNE2 share the transmembrane topology and a high degree of sequence homology in TMD and surrounding regions. The structural basis for their distinctly different effects on KCNQ1 is not clear. To address this question, we apply cysteine (Cys) scanning mutagenesis to TMDs and EJMs of KCNE1 and KCNE2. We analyze the patterns of functional perturbation to identify high impact positions, and probe disulfide formation between engineered Cys side chains on KCNE subunits and native Cys on KCNQ1. We also use methanethiosulfonate reagents to probe the relationship between EJMs of KCNE subunits and KCNQ1. Our data suggest that the TMDs of both KCNE subunits are at about the same location but interact differently with KCNQ1. In particular, the much closer contact of KCNE2 TMD with KCNQ1, relative to that of KCNE1, is expected to impact the allosteric modulation of KCNQ1 pore conductance and may explain their differential effects on the KCNQ1 current amplitude. KCNE1 and KCNE2 also differ in the relationship between their EJMs and KCNQ1. Although the EJM of KCNE1 makes intimate contacts with KCNQ1, there appears to be a crevice between KCNQ1 and KCNE2. This putative crevice may perturb the electrical field around the voltage-sensing domain of KCNQ1, contributing to the differential effects of KCNE2 versus KCNE1 on KCNQ1 gating kinetics.

title: KCNE3 acts by promoting voltage sensor activation in KCNQ1
abstract: The association of KCNE3 beta subunits to KCNQ1 channels turns voltage-dependent KCNQ1 channels into apparent voltage-independent KCNQ1/KCNE3 channels that are important for the transport of water and salts across epithelial cell layers. Because KCNQ1/KCNE3 channels are necessary for water and salt secretion in the colon, KCNQ1/KCNE3 channels are a potential drug target in the treatment of secretory diarrhea. Mutations in KCNE3 have also been associated with diseases such as cardiac arrhythmia and tinnitus. We here propose a model for how KCNE3 turns KCNQ1 into a voltage-independent channel. Our model will allow for a better understanding of how mutations in KCNQ1 and KCNE3 cause diseases and how to design drugs to treat these diseases. KCNE β-subunits assemble with and modulate the properties of voltage-gated K+ channels. In the colon, stomach, and kidney, KCNE3 coassembles with the α-subunit KCNQ1 to form K+ channels important for K+ and Cl− secretion that appear to be voltage-independent. How KCNE3 subunits turn voltage-gated KCNQ1 channels into apparent voltage-independent KCNQ1/KCNE3 channels is not completely understood. Different mechanisms have been proposed to explain the effect of KCNE3 on KCNQ1 channels. Here, we use voltage clamp fluorometry to determine how KCNE3 affects the voltage sensor S4 and the gate of KCNQ1. We find that S4 moves in KCNQ1/KCNE3 channels, and that inward S4 movement closes the channel gate. However, KCNE3 shifts the voltage dependence of S4 movement to extreme hyperpolarized potentials, such that in the physiological voltage range, the channel is constitutively conducting. By separating S4 movement and gate opening, either by a mutation or PIP2 depletion, we show that KCNE3 directly affects the S4 movement in KCNQ1. Two negatively charged residues of KCNE3 (D54 and D55) are found essential for the effect of KCNE3 on KCNQ1 channels, mainly exerting their effects by an electrostatic interaction with R228 in S4. Our results suggest that KCNE3 primarily affects the voltage-sensing domain and only indirectly affects the gate.

title: KCNE1 and KCNE3 Stabilize and/or Slow Voltage Sensing S4 Segment of KCNQ1 Channel
abstract: KCNQ1 is a voltage-dependent K+ channel whose gating properties are dramatically altered by association with auxiliary KCNE proteins. For example, KCNE1, which is mainly expressed in heart and inner ear, markedly slows the activation kinetics of KCNQ1. Whether the voltage-sensing S4 segment moves differently in the presence of KCNE1 is not yet known, however. To address that question, we systematically introduced cysteine mutations, one at a time, into the first half of the S4 segment of human KCNQ1. A226C was found out as the most suited mutant for a methanethiosulfonate (MTS) accessibility analysis because it is located at the N-terminal end of S4 segment and its current was stable with repetitive stimuli in the absence of MTS reagent. MTS accessibility analysis revealed that the apparent second order rate constant for modification of the A226C mutant was state dependent, with faster modification during depolarization, and was 13 times slower in the presence of KCNE1 than in its absence. In the presence of KCNE3, on the other hand, the second order rate constant for modification was not state dependent, indicating that the C226 residue was always exposed to the extracellular milieu, even at the resting membrane potential. Taken together, these results suggest that KCNE1 stabilizes the S4 segment in the resting state and slows the rate of transition to the active state, while KCNE3 stabilizes the S4 segment in the active state. These results offer new insight into the mechanism of KCNQ1 channel modulation by KCNE1 and KCNE3.

title: Molecular mechanisms of disease-causing missense mutations.
abstract: Genetic variations resulting in a change of amino acid sequence can have a dramatic effect on stability, hydrogen bond network, conformational dynamics, activity and many other physiologically important properties of proteins. The substitutions of only one residue in a protein sequence, so-called missense mutations, can be related to many pathological conditions and may influence susceptibility to disease and drug treatment. The plausible effects of missense mutations range from affecting the macromolecular stability to perturbing macromolecular interactions and cellular localization. Here we review the individual cases and genome-wide studies that illustrate the association between missense mutations and diseases. In addition, we emphasize that the molecular mechanisms of effects of mutations should be revealed in order to understand the disease origin. Finally, we report the current state-of-the-art methodologies that predict the effects of mutations on protein stability, the hydrogen bond network, pH dependence, conformational dynamics and protein function.

title: Combining Structural Modeling with Ensemble Machine Learning to Accurately Predict Protein Fold Stability and Binding Affinity Effects upon Mutation
abstract: Advances in sequencing have led to a rapid accumulation of mutations, some of which are associated with diseases. However, to draw mechanistic conclusions, a biochemical understanding of these mutations is necessary. For coding mutations, accurate prediction of significant changes in either the stability of proteins or their affinity to their binding partners is required. Traditional methods have used semi-empirical force fields, while newer methods employ machine learning of sequence and structural features. Here, we show how combining both of these approaches leads to a marked boost in accuracy. We introduce ELASPIC, a novel ensemble machine learning approach that is able to predict stability effects upon mutation in both, domain cores and domain-domain interfaces. We combine semi-empirical energy terms, sequence conservation, and a wide variety of molecular details with a Stochastic Gradient Boosting of Decision Trees (SGB-DT) algorithm. The accuracy of our predictions surpasses existing methods by a considerable margin, achieving correlation coefficients of 0.77 for stability, and 0.75 for affinity predictions. Notably, we integrated homology modeling to enable proteome-wide prediction and show that accurate prediction on modeled structures is possible. Lastly, ELASPIC showed significant differences between various types of disease-associated mutations, as well as between disease and common neutral mutations. Unlike pure sequence-based prediction methods that try to predict phenotypic effects of mutations, our predictions unravel the molecular details governing the protein instability, and help us better understand the molecular causes of diseases.

title: Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability
abstract: There is a compelling and growing need to accurately predict the impact of amino acid mutations on protein stability for problems in personalized medicine and other applications. Here the ability of 10 computational tools to accurately predict mutation-induced perturbation of folding stability (ΔΔG) for membrane proteins of known structure was assessed. All methods for predicting ΔΔG values performed significantly worse when applied to membrane proteins than when applied to soluble proteins, yielding estimated concordance, Pearson, and Spearman correlation coefficients of <0.4 for membrane proteins. Rosetta and PROVEAN showed a modest ability to classify mutations as destabilizing (ΔΔG < −0.5 kcal/mol), with a 7 in 10 chance of correctly discriminating a randomly chosen destabilizing variant from a randomly chosen stabilizing variant. However, even this performance is significantly worse than for soluble proteins. This study highlights the need for further development of reliable and reproducible methods for predicting thermodynamic folding stability in membrane proteins.

title: Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis
abstract: Advances over the past 25 years have revealed much about how the structural properties of membranes and associated proteins are linked to the thermodynamics and kinetics of membrane protein (MP) folding. At the same time biochemical progress has outlined how cellular proteostasis networks mediate MP folding and manage misfolding in the cell. When combined with results from genomic sequencing, these studies have established paradigms for how MP folding and misfolding are linked to the molecular etiologies of a variety of diseases. This emerging framework has paved the way for the development of a new class of small molecule “pharmacological chaperones” that bind to and stabilize misfolded MP variants, some of which are now in clinical use. In this review, we comprehensively outline current perspectives on the folding and misfolding of integral MPs as well as the mechanisms of cellular MP quality control. Based on these perspectives, we highlight new opportunities for innovations that bridge our molecular understanding of the energetics of MP folding with the nuanced complexity of biological systems. Given the many linkages between MP misfolding and human disease, we also examine some of the exciting opportunities to leverage these advances to address emerging challenges in the development of therapeutics and precision medicine.

title: Conformational Stability and Pathogenic Misfolding of the Integral Membrane Protein PMP22
abstract: Despite broad biochemical relevance, our understanding of the physiochemical reactions that limit the assembly and cellular trafficking of integral membrane proteins remains superficial. In this work, we report the first experimental assessment of the relationship between the conformational stability of a eukaryotic membrane protein and the degree to which it is retained by cellular quality control in the secretory pathway. We quantitatively assessed both the conformational equilibrium and cellular trafficking of 12 variants of the α-helical membrane protein peripheral myelin protein 22 (PMP22), the intracellular misfolding of which is known to cause peripheral neuropathies associated with Charcot–Marie–Tooth disease (CMT). We show that the extent to which these mutations influence the energetics of Zn(II)-mediated PMP22 folding is proportional to the observed reduction in cellular trafficking efficiency. Strikingly, quantitative analyses also reveal that the reduction of motor nerve conduction velocities in affected patients is proportional to the extent of the mutagenic destabilization. This finding provides compelling evidence that the effects of these mutations on the energetics of PMP22 folding lie at the heart of the molecular basis of CMT. These findings highlight conformational stability as a key factor governing membrane protein biogenesis and suggest novel therapeutic strategies for CMT.

title: Control of cystic fibrosis transmembrane conductance regulator membrane trafficking: not just from the endoplasmic reticulum to the Golgi
abstract: Biogenesis of cystic fibrosis transmembrane conductance regulator (CFTR) starts with its cotranslational insertion into the membrane of the endoplasmic reticulum (ER) and core glycosylation. These initial events are followed by a complex succession of steps with the main goal of checking the overall quality of CFTR conformation in order to promote its exit from the ER through the secretory pathway. Failure to pass the various checkpoints of the ER quality control targets the most frequent disease‐causing mutant protein (F508del‐CFTR) for premature degradation. For wild‐type CFTR that exits the ER, trafficking through the Golgi is the major site for glycan processing, although nonconventional trafficking pathways have also been described for CFTR. Once CFTR is at the cell surface, its stability is also controlled by multiple protein interactors, including Rab proteins, Rho small GTPases, and PDZ proteins. These regulate not only anterograde trafficking to the cell surface, but also endocytosis and recycling, thus achieving fine and tight modulation of CFTR plasma membrane levels. Exciting recent data have related autophagy and epithelial differentiation to the regulation of CFTR trafficking. Herein, we review the various checkpoints of the complex quality control along the secretory trafficking pathway and the associated pathways that are starting to be explored for the benefit of cystic fibrosis patients.

title: Compendium of cardiac channel mutations in 541 consecutive unrelated patients referred for long QT syndrome genetic testing.
abstract: The purpose of this study was to determine the spectrum and prevalence of cardiac channel mutations among a large cohort of consecutive, unrelated patients referred for long QT syndrome (LQTS) genetic testing.



Congenital LQTS is a primary cardiac channelopathy. More than 300 mutations have been identified in five genes encoding key ion channel subunits. Until the recent release of the commercial clinical genetic test, LQTS genetic testing had been performed in research laboratories during the past decade.



A cardiac channel gene screen for LQTS-causing mutations in KCNQ1 (LQT1), KCNH2 (LQT2), SCN5A (LQT3), KCNE1 (LQT5), and KCNE2 (LQT6) was performed for 541 consecutive, unrelated patients (358 females, average age at diagnosis 24 +/- 16 years, average QTc 482 +/- 57 ms) referred to Mayo Clinic's Sudden Death Genomics Laboratory for LQTS genetic testing between August 1997 and July 2004. A comprehensive open reading frame and splice site analysis of the 60 protein-encoding exons was conducted using polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing.



Overall, 211 putative pathogenic mutations in KCNQ1 (88), KCNH2 (89), SCN5A (32), KCNE1 (1), and KCNE2 (1) were found in 272 unrelated patients (50%). Among the genotype positive patients (N = 272), 243 had single pathogenic mutations (LQT1: n = 120 patients; LQT2: n = 93; LQT3: n = 26; LQT5: n = 3; LQT6: n = 1), and 29 patients (10% of genotype-positive patients and 5% overall) had two LQTS-causing mutations. The majority of mutations were missense mutations (154/210 [73%]), singletons (identified in only a single unrelated patient: 165/210 [79%]), and novel (125/211 [59%]). None of the mutations identified were seen in more than 1,500 reference alleles. Those patients harboring multiple mutations were younger at diagnosis (15 +/- 11 years vs 24 +/- 16 years, P = .003).



In this comprehensive cardiac channel gene screen of the largest cohort of consecutive, unrelated patients referred for LQTS genetic testing, half of the patients had an identifiable mutation. The majority of mutations continue to represent novel singletons that expand the published compendium of LQTS-causing mutations by 35%. These observations should facilitate diagnostic interpretation of the clinical genetic test for LQTS.

title: Generalized Fragment Picking in Rosetta: Design, Protocols and Applications
abstract: The Rosetta de novo structure prediction and loop modeling protocols begin with coarse grained Monte Carlo searches in which the moves are based on short fragments extracted from a database of known structures. Here we describe a new object oriented program for picking fragments that greatly extends the functionality of the previous program (nnmake) and opens the door for new approaches to structure modeling. We provide a detailed description of the code design and architecture, highlighting its modularity, and new features such as extensibility, total control over the fragment picking workflow and scoring system customization. We demonstrate that the program provides at least as good building blocks for ab-initio structure prediction as the previous program, and provide examples of the wide range of applications that are now accessible.

title: OCTOPUS: improving topology prediction by two-track ANN-based preference scores and an extended topological grammar
abstract: As alpha-helical transmembrane proteins constitute roughly 25% of a typical genome and are vital parts of many essential biological processes, structural knowledge of these proteins is necessary for increasing our understanding of such processes. Because structural knowledge of transmembrane proteins is difficult to attain experimentally, improved methods for prediction of structural features of these proteins are important.



OCTOPUS, a new method for predicting transmembrane protein topology is presented and benchmarked using a dataset of 124 sequences with known structures. Using a novel combination of hidden Markov models and artificial neural networks, OCTOPUS predicts the correct topology for 94% of the sequences. In particular, OCTOPUS is the first topology predictor to fully integrate modeling of reentrant/membrane-dipping regions and transmembrane hairpins in the topological grammar.


AVAILABILITY
OCTOPUS is available as a web server at http://octopus.cbr.su.se.

title: Modeling Symmetric Macromolecular Structures in Rosetta3
abstract: Symmetric protein assemblies play important roles in many biochemical processes. However, the large size of such systems is challenging for traditional structure modeling methods. This paper describes the implementation of a general framework for modeling arbitrary symmetric systems in Rosetta3. We describe the various types of symmetries relevant to the study of protein structure that may be modeled using Rosetta's symmetric framework. We then describe how this symmetric framework is efficiently implemented within Rosetta, which restricts the conformational search space by sampling only symmetric degrees of freedom, and explicitly simulates only a subset of the interacting monomers. Finally, we describe structure prediction and design applications that utilize the Rosetta3 symmetric modeling capabilities, and provide a guide to running simulations on symmetric systems.

title: OPM database and PPM web server: resources for positioning of proteins in membranes
abstract: The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topology and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approximately 350 organisms that represent approximately 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielectric constant and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calculating spatial positions in membranes of newly determined proteins structures or theoretical models.

title: Automatic atom type and bond type perception in molecular mechanical calculations.
abstract: In molecular mechanics (MM) studies, atom types and/or bond types of molecules are needed to determine prior to energy calculations. We present here an automatic algorithm of perceiving atom types that are defined in a description table, and an automatic algorithm of assigning bond types just based on atomic connectivity. The algorithms have been implemented in a new module of the AMBER packages. This auxiliary module, antechamber (roughly meaning "before AMBER"), can be applied to generate necessary inputs of leap-the AMBER program to generate topologies for minimization, molecular dynamics, etc., for most organic molecules. The algorithms behind the manipulations may be useful for other molecular mechanical packages as well as applications that need to designate atom types and bond types.

title: AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine
abstract: AbstractWe report a consistent set of AMBER force-field parameters for the most common phosphorylated amino acids, phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine in different protonation states. The calculation of atomic charges followed the original restrained electrostatic potential fitting procedure used to determine the charges for the parm94/99 parameter set, taking α-helical and β-strand conformations of the corresponding ACE-/NME-capped model peptide backbone into account. Missing force-field parameters were taken directly from the general AMBER force field (gaff) and the parm99 data set with minor modifications, or were newly generated based on ab initio calculations for model systems. Final parameters were validated by geometry optimizations and molecular-dynamics simulations. Template libraries for the phosphorylated amino acids in Leap format and corresponding frcmod parameter files are made available.
 FigureSchematic illustration of the systems used for parameter generation. Acid hydrogens are shown in red

title: Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes
abstract: A numerical algorithm integrating the 3N Cartesian equations of motion of a system of N points subject to holonomic constraints is formulated. The relations of constraint remain perfectly fulfilled at each step of the trajectory despite the approximate character of numerical integration. The method is applied to a molecular dynamics simulation of a liquid of 64 n-butane molecules and compared to a simulation using generalized coordinates. The method should be useful for molecular dynamics calculations on large molecules with internal degrees of freedom.

title: Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems
abstract: An N⋅log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolutions using fast Fourier transforms. Timings and accuracies are presented for three large crystalline ionic systems.

title: PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data.
abstract: We describe PTRAJ and its successor CPPTRAJ, two complementary, portable, and freely available computer programs for the analysis and processing of time series of three-dimensional atomic positions (i.e., coordinate trajectories) and the data therein derived. Common tools include the ability to manipulate the data to convert among trajectory formats, process groups of trajectories generated with ensemble methods (e.g., replica exchange molecular dynamics), image with periodic boundary conditions, create average structures, strip subsets of the system, and perform calculations such as RMS fitting, measuring distances, B-factors, radii of gyration, radial distribution functions, and time correlations, among other actions and analyses. Both the PTRAJ and CPPTRAJ programs and source code are freely available under the GNU General Public License version 3 and are currently distributed within the AmberTools 12 suite of support programs that make up part of the Amber package of computer programs (see http://ambermd.org ). This overview describes the general design, features, and history of these two programs, as well as algorithmic improvements and new features available in CPPTRAJ.

title: HOLE: a program for the analysis of the pore dimensions of ion channel structural models.
abstract: A method (HOLE) that allows the analysis of the dimensions of the pore running through a structural model of an ion channel is presented. The algorithm uses a Monte Carlo simulated annealing procedure to find the best route for a sphere with variable radius to squeeze through the channel. Results can be displayed in a graphical fashion or visualized with most common molecular graphical packages. Advances include a method to analyze the anisotropy within a pore. The method can also be used to predict the conductance of channels using a simple empirically corrected ohmic model. As an example the program is applied to the cholera toxin B-subunit pentamer. The compatibility of the crystal structure and conductance data is established.

title: MAMMOTH (Matching molecular models obtained from theory): An automated method for model comparison
abstract: Advances in structural genomics and protein structure prediction require the design of automatic, fast, objective, and well benchmarked methods capable of comparing and assessing the similarity of low‐resolution three‐dimensional structures, via experimental or theoretical approaches. Here, a new method for sequence‐independent structural alignment is presented that allows comparison of an experimental protein structure with an arbitrary low‐resolution protein tertiary model. The heuristic algorithm is given and then used to show that it can describe random structural alignments of proteins with different folds with good accuracy by an extreme value distribution. From this observation, a structural similarity score between two proteins or two different conformations of the same protein is derived from the likelihood of obtaining a given structural alignment by chance. The performance of the derived score is then compared with well established, consensus manual‐based scores and data sets. We found that the new approach correlates better than other tools with the gold standard provided by a human evaluator. Timings indicate that the algorithm is fast enough for routine use with large databases of protein models. Overall, our results indicate that the new program (MAMMOTH) will be a good tool for protein structure comparisons in structural genomics applications. MAMMOTH is available from our web site at http://physbio.mssm.edu/∼ortizg/.

title: Computational design of membrane proteins using RosettaMembrane
abstract: Computational membrane protein design is challenging due to the small number of high‐resolution structures available to elucidate the physical basis of membrane protein structure, multiple functionally important conformational states, and a limited number of high‐throughput biophysical assays to monitor function. However, structural determination of membrane proteins has made tremendous progress in the past years. Concurrently the field of soluble computational design has made impressive inroads. These developments allow us to tackle the formidable challenge of designing functional membrane proteins. Herein, Rosetta is benchmarked for membrane protein design. We evaluate strategies to cope with the often reduced quality of experimental membrane protein structures. Further, we test the usage of symmetry in design protocols, which is particularly important as many membrane proteins exist as homo‐oligomers. We compare a soluble scoring function with a scoring function optimized for membrane proteins, RosettaMembrane. Both scoring functions recovered around half of the native sequence when completely redesigning membrane proteins. However, RosettaMembrane recovered the most native‐like amino acid property composition. While leucine was overrepresented in the inner and outer‐hydrophobic regions of RosettaMembrane designs, it resulted in a native‐like surface hydrophobicity indicating that it is currently the best option for designing membrane proteins with Rosetta.