Upgraded molecular models of the human KCNQ1 potassium channel


Authors: Georg Kuenze aff001;  Amanda M. Duran aff001;  Hope Woods aff001;  Kathryn R. Brewer aff001;  Eli Fritz McDonald aff001;  Carlos G. Vanoye aff004;  Alfred L. George, Jr. aff004;  Charles R. Sanders aff001;  Jens Meiler aff001
Authors place of work: Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America aff001;  Department of Chemistry, Vanderbilt University, Nashville, Tennessee, United States of America aff002;  Department of Biochemistry, Vanderbilt University, Nashville, Tennessee, United States of America aff003;  Department of Pharmacology, Northwestern University Feinberg School of Medicine, Chicago, Illinois, United States of America aff004;  Department of Pharmacology, Vanderbilt University, Nashville, Tennessee, United States of America aff005
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0220415

Summary

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.

Keywords:

Research and analysis methods – Simulation and modeling – Microscopy – Electron microscopy – Electron cryo-microscopy – Database and informatics methods – Biological databases – Protein structure databases – Bioinformatics – Sequence analysis – Sequence alignment – Biology and life sciences – Biochemistry – Biochemical simulations – Proteins – Ion channels – Ion channel gating – Protein structure – Protein structure databases – Lipids – Computational biology – Biochemical simulations – Biophysics – Ion channels – Ion channel gating – Physiology – Electrophysiology – Ion channels – Ion channel gating – Neurophysiology – Ion channels – Ion channel gating – Neuroscience – Neurophysiology – Ion channels – Ion channel gating – Molecular biology – Macromolecular structure analysis – Protein structure – Protein structure databases – Physical sciences – Chemistry – Physical chemistry – Chemical bonding – Hydrogen bonding – Physics – Biophysics – Ion channels – Ion channel gating – Medicine and health sciences – Physiology – Electrophysiology – Ion channels – Ion channel gating – Neurophysiology – Ion channels – Ion channel gating

Introduction

Voltage-gated ion channels are ubiquitously expressed in human tissues and contribute to diverse physiological phenomena such as generation and modulation of the membrane potential in excitable cells, myocyte contraction, modulation of neurotransmitter and hormone release, and electrolyte transport in epithelia. KCNQ1 (KV7.1) is a voltage-gated potassium (KV) channel expressed in the heart and in epithelial cells of the inner ear, stomach, kidney and colon [1]. Like other KV channels KCNQ1 is comprised of four identical α-subunits, each containing six membrane-spanning segments (S1-S6) and a pore loop (P loop) that contributes to the ion selectivity filter in the homo-tetramers to create the KCNQ1 channel. The central pore domain (PD, S5-P-S6) is surrounded by four voltage-sensing domains (VSDs, S1-S4) that respond to membrane depolarization by undergoing a conformational change that triggers structural rearrangements in the PD (called electromechanical coupling), which opens the channel gate making the channel conductive [2, 3]. VSD activation occurs stepwise and proceeds from an initial resting VSD conformation in which the PD is closed (resting/closed, RC) to an activated VSD with an open PD (activated/open, AO) [4, 5]. In KCNQ1, these transitions also involve an experimentally resolvable intermediate state [6, 7]. KCNQ1 also contains a sizable C-terminus that interacts with both phosphatidyl-4,5-bisphosphate (PIP2) and calmodulin [8, 9]. While modeling its structure is beyond the scope of this study, it plays important roles in regulating KCNQ1 channel function and contains many disease mutation sites [10, 11]. Finally, we note that a hallmark of the KCNQ1 channel is its co-assembly with the KCNE1 auxiliary subunit in the heart to generate the channel complex that is responsible for the slow delayed rectifier current (IKS) necessary for myocardial repolarization [12, 13].

Heritable mutations in KCNQ1 are associated with several cardiac diseases including LQTS, atrial fibrillation, and short QT syndrome [14]. About 50% of the genetic cases of LQTS, which predisposes children and young adults to sudden cardiac death, are associated with dominant mutations in KCNQ1 (type 1 LQTS) [15]. While progress in the functional characterization of LQTS-associated mutations has been made [1619], the molecular mechanisms underlying channel dysfunction remain difficult to assess without the availability of high-accuracy structural data. A detailed molecular understanding is needed to improve decision making for new unclassified KCNQ1 mutations and to support the development of new anti-arrhythmic therapeutics.

The lack of an experimentally determined structure for human KCNQ1 has prompted molecular modeling efforts. The first structural model of the KCNQ1 channel containing all membrane-embedded regions S1-S6 was published in 2007 [20] and since then has stimulated numerous structural and functional studies that provided new insights into KCNQ1, e.g. its regulation by KCNE1 [2124] and other KCNE proteins [25, 26], binding of PIP2 [27, 28], the structural determinants and transitions in VSD activation [5, 29], and mechanisms of VSD-PD electromechanical coupling [30]. The recent determination of a cryo-electron microscopy (cryo-EM) structure of X. leavis KCNQ1 [31] which shares 78% sequence identity to human KCNQ1 provides a new template to upgrade the earlier Rosetta homology models of KCNQ1. The PD of this structure is believed to be in its closed state conformation, while its VSD occupies the activated state. In addition to this cryo-EM structure, recent experimentally determined structures of non-mammalian VSDs in resting conformations [32, 33] provide new templates for modeling the resting/closed channel, which previously was informed by purely computational models built by combining the bacterial KcsA template with de novo prediction methods [20, 34].

Here, using this new set of experimental structures leveraged in Rosetta multiple-template homology modeling, we report the development and structural validation of a second generation of Rosetta models of the human KCNQ1 channel in the closed and open states. These models have improved completeness and accuracy and thus represent useful tools for studying the structure-function relationships in KCNQ1 and the molecular mechanisms of mutation-induced changes in KCNQ1 phenotypes. We showcase the utility of the new models by (1) applying molecular dynamics (MD) simulation to our KCNQ1 models which generated useful insight into the VSD and PD conformations, channel-PIP2 interaction properties and VSD-PD coupling contacts, and (2) by applying Rosetta energy calculations to KCNQ1 VSD mutants to estimate the mutation-induced thermodynamic changes in protein stability, with results supporting the notion that protein destabilization is a common cause of mutation-induced KCNQ1 dysfunction. We anticipate that these Rosetta models will stimulate a new series of integrated structure-electrophysiology studies and personalized medicine approaches in KCNQ1 research.

Results

Model building and quality check of human KCNQ1 models

Structural models of human KCNQ1 (helix S0-S6) in the RC and AO state (Fig 1) were generated by Rosetta homology modeling guided by the multiple sequence alignments shown in S1 and S2 Figs. The template with the highest sequence identity (78%) to human KCNQ1 is the recently published cryo-EM structure of X. leavis KCNQ1 (PDB 5VMS) [31], which appears to represent a decoupled state with an activated VSD and a closed pore (referred to as the AC state). Thus, the two domains served as separate templates for the AO and RC states, respectively. The X. leavis KCNQ1 VSD provided a template for the human KCNQ1 AO model whereas its PD was used in RC state modeling. An additional template for the open pore conformation in the AO model was derived from the crystal structure of the chimeric KV1.2–2.1 channel (PDB 2R9R) [35]. Because no experimental resting state VSD structure of a voltage-gated potassium channel was available, templates of the VSD in the RC model were obtained from VSD structures of related non-mammalian proteins and another structural model: the resting VSD of the C. intestinalis voltage-sensing phosphatase (Ci-VSP) (PDB 4G7Y) [32], VSD2 in the A. thaliana two pore calcium channel protein 1 (TPC1) (PDB 5DQQ) [33], and the resting VSD conformation C3 in a model of the Shaker channel [36].

<h2>Overall architecture of the <i>X</i>. <i>leavis</i> KCNQ1 structure and human KCNQ1 homology models.</h2>
Fig. 1.

Overall architecture of the X. leavis KCNQ1 structure and human KCNQ1 homology models.


(A) Schematic cartoon depicting the functional states of the VSD (light green box) and PD (light yellow box). The S4 helix (with positive gating charges “+”) and S4-S5 linker helix are shown as green and blue cylinders, respectively. The X. leavis KCNQ1 structure is in a decoupled state with an activated VSD and closed PD. Homology models of human KCNQ1 were developed with the VSD / PD in the resting / closed (RC) and activated / open (AO) conformation. Black arrows illustrate conformational changes during voltage-dependent channel activation. (B) Side view of the KCNQ1 channel models. The approximate position of the membrane bilayer is indicated by vertical lines and the extracellular and intracellular sides are labeled EXT and INT, respectively. (C)+(D) View of the KCNQ1 channel models from the extracellular and intracellular side, respectively.

Rosetta homology modeling yielded an ensemble of low energy models that were assessed with Procheck [37] and Molprobity [38]. The RC and AO models with the best Molprobity statistics (compare with S1 Table) were selected as final models. They exhibited favorable Molprobity and clash scores that were equal or better than the 90th percentile of structures in the Molprobity database. No violations for backbone bonds or Cβ deviations were observed, and only four bond angles in both models fell outside of four times the standard deviation of dictionary values. In the AO model, all residues fell within favored (95.9%) or allowed (4.1%) regions of the Ramachandran plot. In the RC model, only 1% of the residues fell in disallowed regions with the remaining being in favored (92.8%) and allowed (6.2%) regions. Analysis of the side-chain conformations revealed 100% (AO model) and 99.5% (RC model) favored rotamers and no bad rotamers. Taken together, these values are well within the range of high-resolution structural models and a positive indicator of the quality of our KCNQ1 models. PDB files of these models are provided in the Supporting Information (S1 PDB and S2 PDB). Morphs between the RC and AO model are shown in supporting videos S1 Movie, S2 Movie and S3 Movie.

Structural validation of human KCNQ1 RC and AO models

Our KCNQ1 RC and AO models are displayed in Fig 1 and their topology is compared to that of X. leavis KCNQ1. KCNQ1 has a domain-swapped homotetrameric topology with the VSD of one subunit interacting with helix S5’ from a neighboring subunit. The VSD and PD in the AO model had a Cα-atom RMSD of 1.9 Å and 2.3 Å relative to the X. leavis KCNQ1 AC structure, respectively. For the RC model the Cα-RMSD for the VSD and PD were 4.1 Å and 1.0 Å, respectively. These values are considerably lower than those calculated for the first generation of human KCNQ1 models [20] which had RMSD values for the VSD and PD of 4.8 Å (3.5 Å when excluding the S2-S3 linker) / 3.0 Å (AO model) and 4.3 Å (3.6 Å) / 3.1 Å (RC model), respectively. Thus, the new KCNQ1 models are closer to the X. leavis KCNQ1 structure (i.e. the most closely related structural homolog available), which is a positive indicator for their improvement over the Rosetta models by Smith et al. [20] (see Discussion).

The KCNQ1 RC and AO models were validated rigorously against known experimental and structural data to assess model confidence. First, we calculated the channel pore radius and compared it to the closed PD conformation in X. leavis KCNQ1 (Fig 2). In the RC model, the positions with the narrowest constriction along the ion conduction pathway are at residues G345, S349 and L353 which correspond to G335, S339 and L343 in X. leavis KCNQ1. The pore radius at those residues (~0.5 Å) is clearly below the ionic radius of a K+ ion (1.38 Å) confirming that the channel gate in the RC model is closed. In contrast, the channel gate is open in the AO model. The radius at the narrowest point of the pore at A344 is not smaller than 2.0 Å, which will allow a K+ ion to pass through. Beyond that point, the pore radius increases and remains larger than the K+ ionic radius. This transition point in the KCNQ1 pore coincides with the functionally important PAG motif, which corresponds to PXP in other KV channels. The PAG motif serves as a hinge point allowing the S6 helix to bend and swing outward, which opens the channel gate (S4 Movie). Mutations of the PAG motif abolish current but not channel surface expression [39] consistent with the idea that this region is important for movement of the activation gate. These observations for our Rosetta models were reproduced by pore radius measurements made for these models when employed in MD simulations (S3 Fig, see explanations below), further confirming the validity and stability of the PD structure in our KCNQ1 models.

<h2>Pore conformation in the <i>X</i>. <i>leavis</i> KCNQ1 structure and human KCNQ1 RC and AO homology models.</h2>
Fig. 2.

Pore conformation in the X. leavis KCNQ1 structure and human KCNQ1 RC and AO homology models.


Left: surface representation of the channel pore. Subunits at the front and back are excluded for clarity. The selectivity filter is displayed as sticks. Amino acid residues facing the pore are labeled. Right: pore radius calculated from the X. leavis KCNQ1 structure and the human KCNQ1 homology models. The approximate radius of a K+ ion is indicated by a gray dashed line.

As a second validation, we measured state-dependent distances in our models between amino acid residues suggested to be responsible for voltage-dependent activation of the VSD and compared them to expectations from known residue pairings inferred from electrophysiological studies on charge-reversal VSD mutants of KCNQ1 [5, 6]. Positively charged residues located along the S4 helix, commonly referred to as gating charges and labeled R1 through R6, undergo a series of stepwise transitions and successively pair with acidic residues on S2 and S3 that creates movement of S4 and confers voltage-sensitivity to the ion channel [40] (S5 Movie). Compared to the Shaker [36] and KV1.2 channel [41], it is currently assumed that KCNQ1 adopts a limited number of VSD states (resting, intermediate and activated) [42], possibly due to the lower number of positively charged S4 residues. Residues R3 and K5 in KV1.2 are replaced by Q3 and H5 in KCNQ1. Critical residues that coordinate S4 movement include E1 (E160/E150 in human/X. leavis KCNQ1), E2 (E170/E160) and D202 (D192 in X. leavis KCNQ1). In the resting VSD, E1 likely pairs with R1 (R228 in human KCNQ1) whereas the activated conformation is stabilized by an interaction between E1 and R4 (R237) [5, 6, 43]. In addition, residue Q3 (Q234/Q224 in human/X. leavis KCNQ1) is found close to E2 in the cryo-EM structure of KCNQ1, which likely represents an activated state of the VSD. The described state-dependent residue pairings are correctly captured by our KCNQ1 Rosetta models and persistently detected in MD simulations of the RC and AO model (Fig 3; see also S5 Movie for a morph between the resting and activated VSD structure). The structural models also suggest additional interactions occur between helix S4 and residues of the charge transfer center: in the RC model R2 and Q3 interact with E2 and D202 whereas in the AO model H5 and R6 are close to those two acidic residues. Moreover, another source of structural stabilization of the VSD in its resting and activated conformation may come from interactions of the gating charges with membrane phospholipids as observed for the VSD of other KV channels [4446]. This idea is supported by the observation of frequent hydrogen bond contacts of R6 (RC) and R1 (AO) with phospholipids in the lower and upper bilayer leaflet during MD simulation (Fig 3).

<h2>Voltage sensor conformation in the <i>X</i>. <i>leavis</i> KCNQ1 structure and human KCNQ1 RC and AO homology models.</h2>
Fig. 3.

Voltage sensor conformation in the X. leavis KCNQ1 structure and human KCNQ1 RC and AO homology models.


Left: VSD conformation; helices S0-S3 are colored pale green and helix S4 is colored blue. The side-chains of positively charged residues on S4 (yellow) and the gating charge transfer center residues E150/E160 (S2, yellow), F157/F167 (S2, magenta), E160/E170 (S2, yellow) and D192/D202 (S3, yellow) are depicted as sticks. Phosphate groups of phospholipids in the outer (top) and inner (bottom) bilayer leaflet are represented as spheres. Right: side-chain–side-chain distances between residues on S4 and gating charge transfer center residues. Values represent the average over four VSD copies, and were calculated from the cryo-EM structure of X. leavis KCNQ1 (after embedding in a POPC bilayer with CHARMM-GUI [47] and short minimization in Amber16 [48]) and from MD simulations of human KCNQ1 models. Distances were measured between the geometric centers of the side-chain atoms H2N = Cζ(NH2)-NεH-CδH2 (Arg), H3Nζ-CεH2 (Lys), Cγ-Nδ1-Cε1H-Nε2H-Cδ2H (His), HOOCγ-CβH2 (Asp), HOOCε-CγH2 (Glu) and the lipid phosphate group PO4-.

Overall, the measured channel pore dimensions and the good agreement of residue pairings in the VSD with experimental data lend confidence to the reliability of our KCNQ1 structural models. Further, they provide insight into the structural dynamics of VSD activation and channel opening and closing.

Refinement of KCNQ1 models by MD simulation and study of KCNQ1-lipid interactions

Computational simulations have become an indispensable tool in the study of voltage-gated ion channels and have provided insights into different aspects of KCNQ1 function and regulation, such as KCNE1 binding [21, 49, 50], drug [51] and lipid binding [52]. Here, we assessed the structural stability of our KCNQ1 AO and RC homology models in all-atom MD simulations with Amber [48]. We further validated that the simulations can reproduce known interactions of KCNQ1 with PIP2 previously revealed by experimental data. PIP2 is an important second messenger for cell signaling that binds to and regulates a wide variety of ion channels, including KCNQ1 [53]. Depletion of PIP2 from the membrane suppresses VSD-PD coupling in KCNQ1 [27] but leaves the VSD activation intact.

A total of four RC state and four AO state simulations starting from the selected KCNQ1 models as well as three additional conformations from the ensemble of Rosetta homology models were conducted, each running for 400 ns. Representative snapshots from the trajectory of AO and RC models are provided as supporting files S1 MDtraj and S2 MDtraj, respectively. Simulations were performed under constant temperature and pressure conditions in explicit POPC bilayers that contained 10 mol% (~28 molecules) of PIP2 in the cytosolic leaflet (Fig 4A). Overall, the KCNQ1 models remained stable in MD and deviated by not more than a Cα-RMSD of 4–5 Å from their respective starting conformations (Fig 4B and S4 Fig). The RC and AO models maintained their closed and open pore structure (Fig 2 and S3 Fig), respectively, as well as their characteristic gating charge pairings (Fig 3) as discussed above.

<h2>MD simulation of KCNQ1 homology models.</h2>
Fig. 4.

MD simulation of KCNQ1 homology models.


(A) RC and AO models embedded in a POPC bilayer that contained 10 mol% PIP2 in the inner leaflet. The fatty acid chains of POPC and PIP2 molecules are depicted by gray and cyan sticks, respectively, and the head group phosphates are represented as olive spheres. Potassium and chloride ions are displayed as gray and light blue spheres and water molecules are excluded for clarity. (B) Cα-atom RMSD relative to the starting structure over the course of an MD trajectory for the RC and AO model. RMSD plots for all four RC and AO state MD simulations, respectively, are shown in S4 Fig. (C) Snapshot of the interaction of two PIP2 molecules with KCNQ1 as observed in MD. PIP2 molecules are drawn as spheres with carbon atoms colored yellow or cyan, respectively, and oxygen and phosphorous atoms colored red and orange, respectively. Different KCNQ1 subunits are colored green, blue, purple and orange. Side-chains of nearby KCNQ1 residues are depicted as sticks with heteroatoms and polar hydrogen atoms colored by their chemical identity (hydrogen: white, oxygen: red, nitrogen: blue). Potential hydrogen bond contacts are indicated by dashed lines. (D) Average number of hydrogen bonds per MD frame between KCNQ1 residues and phosphate groups in PC and PIP2 head groups calculated over the course of all RC and AO state MD trajectories.

The number of hydrogen bonds with POPC and PIP2 was analyzed after letting the system equilibrate for 100 ns. Because the mean square displacement values for lateral diffusion of POPC and PIP2 in Amber lipid simulations are on the order of several ten Å2 per 10 ns [54], this time span was long enough to allow the initially randomly placed PIP2 molecules to diffuse towards KCNQ1 and engage in close protein contacts (Fig 4C). PIP2 was observed to bind preferentially in the cleft between two neighboring VSDs (Fig 4C). It made numerous hydrogen bond interactions with residues on the S2-S3 linker (S2-S3L) (R181, K183, R190, R192), the S4-S5 linker (S4-S5L) (R249, R259) and the cytoplasmic end of S6 (S6C) (K354, R360, K362) (Fig 4D) in agreement with the experimentally reported PIP2 binding sites [27, 28]. Previous studies suggested that interactions of PIP2 with these positively charged residues strengthen the coupling between the VSD and PD and are critical for pore opening [55]. Consistent with its high functional importance, a series of LQTS-associated mutations affect the PIP2 binding region: K183R [56], K183M [57], R190Q [58], R190L [59], R190W [60], R192H [61], R192P [60], R259C [62], R259H [61], R259L [63], R360M [59], R360T [60].

Comparison of the locations of POPC and PIP2 interaction sites on the cytosolic juxtamembrane side of KCNQ1 (Fig 4D) revealed some significant differences: POPC was unable to participate in hydrogen bonding with residues R181 and K183 in the short loop between S2 and the S2-S3 linker and with residues on S6c which are located within or below the water-interface region of the membrane. As a consequence of the small PC head group those residues are outside of the range of POPC but can still be spanned by PIP2, with its 1-phosphoinositol-4,5-diphospate headgroup.

In conclusion, all-atom MD simulations of our KCNQ1 models correctly identified experimentally established PIP2 binding sites in S2-S3L, S4, S4-S5L and S6C, and revealed the molecular details of KCNQ1-PIP2 interaction. Owing to its increased head group size, PIP2 can reach and bridge sites in the VSD (S2-S3L, S4) with sites in the PD (S4-S5L, S6C) to act as a molecular coupling element for KCNQ1.

High impact mutation sites map intra- and inter-subunit contact regions between the VSD and PD in KCNQ1 models

Using MD simulation, we further characterized the interactions between the VSD and PD in the RC and AO model and compared them to reported experimental data [30, 6466]. VSD-PD interactions translate the voltage-dependent movement of S4 to channel opening and closing. Inference from studies of other KV channels [6769] and mutational scanning experiments on KCNQ1 [30, 64] have pointed to the interface between the S4-S5 linker (S4-S5L) and the cytosolic end of S6 (S6C) as one important region for electromechanical coupling. In those studies, one group of mutations in S4-S5L [30] and S6C [64] reduced the channel opening rate and shifted channel activation to more depolarized voltages whereas another group of mutations, specifically at V254 in S4-S5L and at L353 in S6C, prevented the channel from closing and led to a constitutively open channel. These results were taken to suggest that S4-S5L binds to S6C to stabilize the closed state of the channel and that S4-S5L relocation in voltage-dependent activation releases tension on S6 that allows it to kink, promoting channel opening.

In our RC and AO models, we identified extended regions of intra- and inter-subunit contacts (displayed in Fig 5A) that tie the VSD and S4-S5L to the neighboring PD. Intra-subunit contacts were found between S4-S5L and S6C and between S5 and the N-terminal end of S6 (S6N). Inter-subunit contacts were observed for S4 and S4-S5L with S5’ from an adjacent KCNQ1 chain and between two S6 helices from neighboring chains that cross each other near the channel gate. Strikingly, mutation sites in S4-S5L, S5 and S6 with high impact on channel opening and closing [30, 64, 66] are clustered in those contact regions and have consistently high contact numbers (Fig 5B and 5D) indicating that those regions are important for channel function. For clarity, mutation sites in the selectivity filter and P helix were not marked in Fig 5B. Those elements share only a small contact area with the VSD and many mutations cause complete channel loss of function (LOF) because of the detrimental effect that these mutations have on the structural integrity of the selectivity filter and extracellular pore entrance gate.

<h2>High impact mutations and LQTS sites fall within intra- and inter-subunit contact regions in KCNQ1 channel models.</h2>
Fig. 5.

High impact mutations and LQTS sites fall within intra- and inter-subunit contact regions in KCNQ1 channel models.


(A) Extracellular view of the KCNQ1 RC and AO models. For clarity, only two neighboring subunits are displayed while the other two chains are transparent. Residues in S4 to S6 that make contacts within the same subunit or with residues in the neighboring subunit are shown as surface representation and are colored light blue or orange, respectively. A cutoff value of 1.0 and 0.6 for the normalized contact number (i.e. mean total number of heteroatom contacts within a 4 Å distance divided by the number of heteroatoms of a given amino acid residue) of intra- and inter-subunit contacts was used for making the surface representations. The normalized contact number was calculated as the average over the four RC or AO state MD simulations, respectively (see Methods). (B) Residues corresponding to high impact mutation sites in KCNQ1 [30, 64, 66] located in S4-S5L, S5 or S6, respectively, are plotted onto the KCNQ1 models and depicted as sticks. Circled areas indicate intra- and inter-subunit contact regions. (C) Location of LQTS sites mapped on the KCNQ1 models depicted as yellow spheres. (D) Normalized contact number of residues at high impact mutation sites in S4-S5L, S5 and the N- (S6N) and C-terminal (S6C) end of S6. (E) Receiver operating characteristic (ROC) curve for classifying a KCNQ1 variant as LQTS or non-LQTS based on a residue’s contact number. The area under the curve (AUC) when using the RC and AO model was 59% and 58%, respectively.

In the RC model, the S4-S5L –S6C contact surface is formed by residues L250, L251, V254, V255, H258 on S4-S5L and I346, L347, G350, F351, L353 and K354 on S6C. The side-chains of V254 and L353 were found to be proximal to each other in the model consistent with the results of combinations of mutations at those two sites [30] that were able to rescue channel function. Interestingly, L353 also engages in side-chain–backbone interactions with one of its symmetry mates in the neighboring subunit. As viewed from the cytosolic side (S4 Movie), the four symmetric copies of L353 are arranged circularly with the orientation of their side-chains directed clockwise. L353 is located one helical turn below the helix-helix crossing point at S349, which forms the narrowest constriction of the channel pore (compare with Fig 2). Hence, our structural model can provide a conclusive explanation for the constitutively open channel phenotype of KCNQ1 mutants with charged amino acid substitutions (K/E) at L353 [64]: steric repulsion between close local charges at the position of L353 pushes the S6 helices apart and prevents the channel from closing.

The contact region between S4-S5L and S6C is maintained in the AO model as well and comprises the same group of residues but is reoriented slightly when S4-S5L and S6C slide against each other during pore opening. The L353 sites, however, become separated from each other when S6 is kinked (S4 Movie) and are not in contact in the AO model.

The inter-subunit contacts span similar regions in the RC and AO model, respectively, but differ in their local structure reflecting the conformational changes that occur during channel activation as described below. Inter-subunit contacts are probably also the reason for the global motions in KCNQ1 that we observed, as characteristic regular patterns in the dynamic residue cross-correlation matrix (DCCM) (S5 Fig), and by performing principal component analysis (PCA) on the intra-membrane segments in KCNQ1 (S5 Fig). The global motions of KCNQ1 in MD simulation feature a rigid-body-like swing movement of the VSDs that includes both an ‘in-membrane-plane’ as well as an ‘out-of-membrane-plane’ component. Two neighboring VSDs were observed moving antiparallel to each other whereas the two domains on opposite sides of the tetrameric channel moved in the same direction (see S5 Fig for further details).

We further noticed that many mutation sites that have been associated with LQTS fall within intra- and inter-subunit contact regions in KCNQ1 (see Fig 5C). A simple classifier that distinguishes LQTS from non-LQTS sites based on residue contact number yields a receiver operating characteristic (ROC) with an area under the curve (AUC) of 59% and 58% for the RC and AO model, respectively (Fig 5E). This shows that LQTS sites are enriched in areas with high contact number and that mutations at those sites are critical to ensure channel function.

Taken together, we show that the VSD and PD in KCNQ1 share extended intra- and inter-subunit contact regions between S4, S4-S5L, S5 and S6. Experimentally tested high-impact mutation sites and sites associated with LQTS fall in those channel regions, which highlights their importance for channel function and renders them possible candidates for roles in electromechanical coupling.

Transition from the RC to the AO state triggers changes in inter-subunit contacts between the VSD and pore

Comparison of the RC and AO model revealed conformational changes at the VSD-PD interface and in the PD (Fig 6) (in addition to S4 in the VSD), which likely reflect the changes happening during depolarization-activated channel gating. Triggered by the gating charge movement in the VSD, S4 slides along S5’ and moves upwards by about two helical turns. Hydrophobic residues on the S5’-facing side of S4 swap their positions and come into contact with new residues on S5’ while S4 is moving (compare with left plot in Fig 6A). The S4-S5L is pulled upwards and sideward (Fig 6B and 6C) and engages in new interactions with residues located two to three turns higher on S5’. For instance, in the RC model, T247 makes contacts with the first residues in S5’ (Q260, E261) but is packed against T264, Y267 and I268 in the AO model. W248 makes few interactions with S5’ in the RC model but makes many putative stabilizing contacts with Y267, I268 and L271 in the AO model (Fig 6B). Moreover, movement of S4-S5L allows S6 to kink and move sideward cantilever-like. This movement breaks contacts between the ends of S6C and induces a crosswise sliding of the S6 helices (compare with right plot in Fig 6A) that widens the channel pore diameter. Specifically, residue F351 slides upwards by about two helical turns and becomes tightly packed in the AO model where it appears to bridge sites on S4-S5L with S6C of the same and S6’C of the neighboring subunit (Fig 6B). The observation that F351 is at a critical spot in KCNQ1 with many contacts to other residues is consistent with the fact that F351 is highly conserved in KV channels, and that it is a high impact mutation site [64]. Experimental studies [64] showed that the F351A mutant featured a rightward shifted current-voltage relationship and longer activation and deactivation times that mimicked the IKS-phenotype of KCNQ1, and was ascribed to a destabilization of the KCNQ1 intermediate/open state [6] similar to the effects of KCNE1.

<h2>Conformational changes in KCNQ1 during channel activation inferred from modeling.</h2>
Fig. 6.

Conformational changes in KCNQ1 during channel activation inferred from modeling.


(A) Change in average contact number between residues in two neighboring KCNQ1 chains A and B. The contact number was defined as number of heteroatom pairs within a 4 Å distance and averaged over the last 300 ns of MD and all four pairs of neighboring subunits in the KCNQ1 tetramer. The section of the contact matrix corresponding to the interface of S4 and S4-S5L with S5’ (left) and of two neighboring helices S6 and S6’ (right) is shown. A gray color means this contact is observed in the RC model whereas a red color denotes a contact formed in the AO model. Changes in specific residue contacts occurring while the channel transitions from the RC to the AO state are framed and labeled by their corresponding amino acid residue. The direction of the structural changes with channel activation is indicated by a blue arrow. (B) Cartoon representation of the inter-subunit interface in the RC (left) and AO (right) model, respectively. For clarity, only helices S4 to S6 are shown. Residues which are part of the inter-subunit interface and fall within regions of the contact matrix in (A) are depicted as sticks. Residues with drastic changes in their contact pattern as identified in (A) are labeled. (C) Movement of helix S4, S4-S5L, and S6C during transition from the RC to AO state. Only one KCNQ1 subunit is displayed for clarity with the superimposition optimized for the pore domain. Helix axes are marked by dashed lines and the direction of their movement is indicated by arrows. S4-S5L moves upward to the extracellular side and tilts sideward allowing S6C to kink.

Hou and coworkers [7] demonstrated that KCNQ1 when expressed alone can also undergo an opening transition once the VSD is in the intermediate or activated conformation. To characterize the structural changes associated with channel opening from an activated VSD, we further compared the AO model with an additional model for the KCNQ1 AC state. The AC model was generated by comparative modeling based on the X leavis cryo-EM structure [31], and is available as supporting file S3 PDB. Residue contacts were analyzed by conducting four independent 350 ns all-atom MD simulations in POPC bilayers. The AC model maintained its activated VSD and closed PD conformation during MD (S3 and S4 Figs). The S4 helix is in an upward position as in the AO model, but the intracellular channel gate is closed with the narrowest constrictions formed by residues S349 and L353 (S3 and S6 Figs). At its N-terminal site, S4 makes fewer contacts with S5’ in the AC model (S6A and S6B Fig) suggesting that the VSD-PD interface is less tightly packed in the AC state. As seen from within the membrane plane, S4-S5L adopts a similar position in the AC and AO model, but is oriented slightly different and tilted inwards and closer to S6 in the AC model (S6B and S6C Fig). Inward tilting of S4-S5L likely stabilizes the closed PD conformation in the AC state by holding neighboring S6 helices together and increasing the number of inter-subunit contacts (S6A and S6B Fig). This observation suggests that in the AC state the channel is primed to open but that an additional outward tilting of S4-S5L is needed before S6 can bend to facilitate widening of the channel pore. We speculate that the separation of S4 and S4-S5L movement as inferred from our modeling of the AC and AO states and supported by the X. leavis KCNQ1 structure [31] is the result of decreased electromechanical coupling in KCNQ1 which is strongly dependent on PIP2 binding [6, 27, 55].

In conclusion, we characterized VSD-PD interactions and their state-dependent changes in KCNQ1 models when the channel transitions from the closed (RC or AC) to the open (AO) state. Our observations provide a stepping stone for systematically testing protein-contact-mediated coupling mechanisms in the activated channel state.

Rosetta-predicted stability changes in KCNQ1 models correlate with experimental channel loss of function data

While progress in the functional characterization of LQTS-associated mutations has been made [1619, 7072] the mechanistic molecular basis of channel dysfunction for most KCNQ1 mutations is still unclear. The dearth of direct experimental data that can help to determine how mutations alter KCNQ1 structure and function has prompted computational modeling approaches [20, 73].

Here, we assessed the capability of our structural models to predict the change in thermodynamic stability (ΔΔG) in KCNQ1 variants and whether the predicted energy changes can account for alterations in channel function. We used a Rosetta ΔΔG prediction protocol [74] allowing for flexibility in the protein backbone and side-chains around the mutation site. Protein backbone flexibility was modeled by a series of backrub moves [75] that rotate short main-chain segments (three to twelve residues) as rigid bodies about an axis defined by the starting and ending Cα-atom of the segment. This aims at alleviating the prediction for mutations with a large change in amino acid size. Side-chain flexibility was modeled using discrete rotamer conformations and simulated annealing (i.e. side chain packing). The ΔΔG was calculated as mean difference between the three top-scoring models of wild type (WT) and mutant KCNQ1 and averaged over ten different starting conformations taken from the ensemble of Rosetta homology models as described in Methods.

ΔΔG predictions were made for a set of single-site KCNQ1 variants with the affected sites located in the VSD (shown in Fig 7A, and divided by functional classes in S7 Fig). These KCNQ1 VSD variants were previously biochemically, biophysically, and functionally characterized [19], which provided the mechanistic basis for their pathogenicity. Comparison with this variant dataset allowed us to validate our computational stability predictions and assess the capabilities and limitations of a structure-based approach to predict KCNQ1 dysfunction. Not included from this original set of 51 KCNQ1 variants was the deletion mutant ΔF167 because modeling of deletions and insertions is currently not supported by the Rosetta ΔΔG prediction protocol. Furthermore, ΔΔG predictions for mutations to proline were found to be off-scale due to incompatible backbone torsions yielding Ramachandran and ring closure penalties as noted earlier [76]. Therefore, their ΔΔG values were considered less reliable and also not included in the analysis leaving ΔΔG data for 44 single-site KCNQ1 variants.

<h2>Computationally predicted stability changes of KCNQ1 variants with mutations in the VSD.</h2>
Fig. 7.

Computationally predicted stability changes of KCNQ1 variants with mutations in the VSD.


(A) Location of mutation sites in the KCNQ1 VSD. Backbone sites are mapped by yellow spheres and the native amino acid residue is indicated by yellow sticks. (B) Distribution of Rosetta ΔΔG values for the six functionally distinct classes of KCNQ1 VSD variants [19] (i.e. WT-like and classes I to V for non-functional variants) calculated with the AO homology model. The median and average values are drawn as a black horizontal line and a green triangle, respectively. The median of class IV and V is compared to that one of WT-like variants using a Kruskal-Wallis H-test (* p < 0.05, ** p < 0.01, *** p < 0.001, nWT-like = 15, nIV = 5, nV = 11). ΔΔG values for mutations to proline are off-scale (ΔΔGL114P = 70.5 ± 4.1 REU, ΔΔGL131P = 76.0 ± 3.0 REU, ΔΔGL134P = 68.3 ± 4.1 REU, ΔΔGR195P = 56.8 ± 2.0 REU, ΔΔGQ234P = 71.0 ± 5.0 REU, ΔΔGL236P = 65.8 ± 6.0 REU) due to incompatible backbone torsions in the starting model yielding bad backbone and proline ring geometries and were not used in the analysis. (C) Correlation plots of total expression level, cell surface expression, trafficking efficiency and channel peak current density versus calculated Rosetta ΔΔG values (mean ± S.E.M.). KCNQ1 variant classes are indicated with different colors. Variants falling outside of or being close to the boundary of the 95% confidence interval for a linear regression model (gray shaded area) are labeled and their structural models are shown in S11 Fig in the supplement. The experimental data and mutant classification scheme are from reference [19] and are listed together with the computed ΔΔG values for RC and AO models in S2 Table. Class I variants traffic normally but exhibit channel LOF. Class II mutations induce mistrafficking, but do not alter the functionality of the small population of correctly-trafficked channel. Class III variants traffic normally and are conductive, but exhibit altered channel regulation. Class IV mutations cause mistrafficking and channel LOF. Class V variants are expression and trafficking-defective and have so low current that channel function is unmeasurable.

Fig 7B shows the distribution of ΔΔG values for those variants calculated with the KCNQ1 AO model as a function of variant class. Variants were previously assigned one of six classes based on their channel phenotype: WT-like and classes I-V for variants that have perturbations in the expression level, surface trafficking, and/or channel current density (see also the caption for Fig 7 and the footnote of S2 Table for a description of the class phenotype). Roughly, the severity of the channel LOF increases from class I to V with variants in class IV and V having the most severe perturbations. The median ΔΔG value of these two classes was significantly higher (ΔΔG˜IV=10.7REU (pIV < 0.05), ΔΔG˜V=26.9REU (pV < 0.001)) than that one of WT-like variants (ΔΔG˜WT‐like=0.72REU) indicating the mutant protein is thermodynamically destabilized. It must be mentioned that Rosetta ΔΔG values (in Rosetta Energy Units, REU) do not necessarily have the same scale as experimentally measured protein stability values (in kcal/mol), but that their relative order (e.g. between different variant classes) is more instructive.

The predicted ΔΔG values were correlated with the expression, trafficking and current density data (Fig 7C) suggesting protein stability is an important factor that contributes to these channel properties, and protein destabilization is a main cause of channel LOF. This idea is supported by NMR measurements on the KCNQ1 voltage sensor [19], which for two of seven class IV variants and ten of 14 class V variants had revealed spectra that were characteristic for a folding-destabilized protein. Another four class V variants failed to express the voltage sensor in E. coli also indicating protein folding defects [19]. The channel state of the structural model (RC or AO) had only minor influence on the ΔΔG predictions. Similar results were also obtained with the RC model (S8 and S9 Figs), indicating that most of the investigated mutations affected both channel states to a similar degree. Examples of mutants with varying ΔΔG predictions in the RC and AO model are shown in supporting S10 Fig and further discussed in the caption to S10 Fig.

The correlation plots in Fig 7C seem to have a complex, non-linear dependency, specifically the plot for peak current density, suggesting additional factors other than protein stability affect the named properties and/or that some variants had incorrect ΔΔG predictions. Variants that fell outside the 95% confidence interval for a linear model included four class I variants (C122Y, A150T, K218E, I227L), three class V variants (G179S, W176R, H126L), A102S and A149V (both WT-like) as well as P197L (class II) and R109L (class III). Structural models for these variants are shown in S11 Fig in the supporting information and their experimental and ΔΔG data are listed in S2 Table. In particular, it was noticeable that class I variants had ΔΔG values (ΔΔG˜I=3.4REU or ΔΔG˜I=1.4REU with the AO or RC model, respectively) not significantly different from WT-like variants (ΔΔG˜WT‐like=0.7REU), which failed to explain why these channels have low current density. For example, C122Y was predicted a negative ΔΔG (i.e. more stabilized than WT) but this channel was nonconductive. Two other variants for which channel function was not correlated with thermodynamic stability were A102S and R109L which had normal (0.2 REU) or increased ΔΔG values (10.3 REU), respectively, but exhibited current densities more than two times larger than WT (i.e. these were gain-of-function mutations).

The other few cases where the stability prediction was inconsistent with the functional data were likely due to incorrect ΔΔG calculations, which was supported by visual inspection of the respective models (see S11 Fig). For instance, A149V was predicted to have a high positive ΔΔG but had only moderately reduced expression levels and a trafficking efficiency higher than WT. Visual inspection of the structural model of this variant (S11 Fig) indicated remaining steric clashes between the valine side-chain and the nearby protein backbone causing a high repulsive score. This suggests that the backrub moves that are performed in the Flex ddG protocol prior to amino acid mutation had optimized the packing of the native alanine residue with the neighboring S3/S4 loop and generated a backbone conformational ensemble that was not sufficiently diverse to accommodate the larger valine side-chain upon mutation. Mutants G179S, W176R and H126L had ΔΔG values that were lower than other class V variants, mostly due to score balancing interactions made by the substituted amino acid residue in the mutant model, inconsistent with their low expression levels, trafficking and conductance.

In summary, Rosetta energy calculations correctly predicted destabilizations for most of the folding-defective KCNQ1 LOF variants: 13/18 variants with low expression (<65%WT) had increased ΔΔG values above 10 REU. Our results support growing evidence [19] that mutation-induced protein destabilization is a common cause of KCNQ1 dysfunction. They further suggest other so far unknown mutational effects exist (e.g. in class I variants), which are not detected in protein stability calculations, but have an impact on KCNQ1 function.

Discussion

Improvements of upgraded KCNQ1 Rosetta models

Our ability to improve homology models of human KCNQ1 was largely informed by the new X. leavis KCNQ1 template, the closest structural homolog available. Consequently, the Rosetta models presented here bear high similarity to the X. leavis KCNQ1 cryo-EM structure and are closer to it than the Rosetta models by Smith et al. [20]. The backbone RMSD for the VSD / PD compared to X. leavis KCNQ1 is 1.9 Å / 2.3 Å (this study) and 4.8 Å / 2.9 Å (Smith et al. [20]), respectively. For the closed state, the RMSD values for the VSD / PD are 4.1 Å / 1.0 Å (this study) and 4.3 Å / 3.1 Å (Smith et al. [20]), respectively. This indicates considerable remodeling changes were applied to the Rosetta models: comparison of the models presented herein and in Smith el al. shows an RMSD for the VSD / PD of 5.3 Å / 3.2 Å (closed) and 5.2 Å / 2.4 Å (open). It is important to mention that the Rosetta models of this work also compare favorably with their other structural templates: the PD in the AO model with the PD in KV1.2–2.1 (1.9 Å) and the VSD in the RC model with the VSD in Ci-VSP (3.6 Å for sequence-aligned regions; 2.6 Å for S1+S2) and TPC1 (4.4 Å for sequence-aligned regions; 2.1 Å for S1+S2). These results, together with the favorable stereochemical properties of the new human KCNQ1 Rosetta models and their valid VSD and PD conformations as illustrated in Results, lends confidence that they represent a substantial improvement over the preceding models. It is also important to point out that the Rosetta models by Smith et al. did not include the S0 helix in the VSD, which is now appreciated to represent a crucial element for structural stability (see reference [19] and Results on KCNQ1 VSD mutations). The previous models further fail to satisfy the E1-R1 (E160-R228) gating charge interaction in the RC state which is an experimentally established contact stabilizing the resting state [5, 6], however, was not known for KCNQ1 at the time of the first model generation.

The improved Rosetta models should be regarded as working models that are subject to revision as more experimental structural data become available. Specifically, we need to point out a few limitations. Of lower confidence are those regions for which no coordinate information in the template structures was available, requiring de novo modeling. Furthermore, higher uncertainty applies to transition regions where templates from different sources were hybridized. Lower confidence regions include the loop between S4 and S4-S5L in the RC and AO models, the loop between S3 and S4 in the AO model and the loops between S0 and S1 and between S2-S3L and S3 in the RC model.

Implications for KCNQ1 structural biology questions

Given improved completeness and accuracy, the Rosetta models provide helpful tools to refine and advance our understanding of the structure-function relationships in KCNQ1. A hallmark of the KCNQ1 channel is its distinctive functional properties observed upon association with different members of the family of KCNE accessory proteins. While KCNE1 increases current density and delays channel opening and closing [12, 13], KCNE2 constitutively activates KCNQ1 and reduces current levels to approximately 50% of KCNQ1 alone [77]. KCNE3 also makes KCNQ1 constitutively open but increases current by ~10 fold [26, 78].

The mechanism by which KCNE1 modulates KCNQ1 is a longstanding yet still unsettled debate and different models have been put forward such as alteration of the VSD movement [72, 7983], perturbation of pore opening [84] and changes in VSD-PD coupling interactions [6, 81, 82]. A mapping onto our KCNQ1 Rosetta models of putative KCNE1 binding sites identified in cysteine crosslinking [22, 23, 8486], mutational scanning and double mutant cycle [25, 87] experiments (see S12 Fig) suggests very similar binding pockets for KCNE1 in the RC and AO state. Preliminary docking calculations with our Rosetta models also suggest that KCNE1 makes similar interactions with KCNQ1 in the RC and AO states and adopts similar orientations (S12 Fig), which is an unexpected result and was not observed in the original Rosetta models [21] indicating that a revision of the old KCNQ1/KCNE1 models is well-merited.

Of the many previous structure-function studies of KCNE1, we want to highlight here the model proposed by Zaydman et al. [6] who suggested that KCNE1 modulates KCNQ1 function by changing the VSD-PD coupling interactions. This idea stems from observations of KCNE1 inhibiting channel opening when the VSD is in the intermediate conformation between resting and activated state, as well as studies of F351A, a mutation that impairs coupling by weakening the contact between S4-S5L and S6C, mimicking the KCNE1-phenotype of KCNQ1. Combining electrophysiology and kinetic modeling, Zaydman et al. concluded that changes in state-dependent VSD-PD interactions may explain many of the effects of KCNE1 including its delay of channel opening, increased channel conductance and pharmacological properties. Consistent with this idea, we identified extended intra- and inter-subunit contact regions between the VSD and PD in our Rosetta models (Fig 5) that potentially contribute to electromechanical coupling in addition to the established S4-S5L –S6C contact. Several other studies [6, 7, 27, 55, 81] hinted at the role of protein-protein and protein-PIP2-dependent mechanisms in KCNQ1 electromechanical coupling. It appears to be intrinsically weak and requires PIP2 binding to strengthen the interaction between the VSD and PD [6, 27, 55]. This is evidenced by the cryo-EM structure of X. leavis KCNQ1 [31] which was determined in the absence of PIP2 and represents a decoupled state. These data led to the conclusion that PIP2 acts as important coupling element which is especially prominent in KCNQ1. Recently, Hou and coauthors [7] also proposed different protein-protein coupling mechanisms may exist in KCNQ1 for channel opening from the intermediate and activated state of the VSD, and that these differences underlie the inactivation phenotype of KCNQ1. We speculate that KCNE1 modifies these VSD-PD interactions, changing the properties of the open pore conformation and affecting the channel open probability. Previous docking studies [21, 49, 87] converged on a model in which KCNE1 is bound in a crevice between the VSD and PD of two neighboring subunits. The intimate interactions that KCNE1 makes in these models with several sites on S4, S4-S5L, S5 and S6 seems well suited to modulate electromechanical coupling interactions.

Explicit structure-function knowledge concerning KCNE2 binding and regulation of KCNQ1 is limited. Results from cysteine scanning [88], double mutant cycle analysis and computational docking [25] suggested a different interaction mode for KCNE2 than for KCNE1. However, a widely recognized molecular model for the KCNQ1-KCNE2 interaction has yet to be established. The KCNQ1 Rosetta models provide a means to devise new experimentally-testable hypotheses on KCNE2’s mode of action that should trigger further progress in this field.

For KCNE3, a series of functional studies [89, 90] have converged on the appreciation that it acts primarily by stabilizing the VSD, especially S4, in its fully activated state, thereby increasing the open probability at resting potential and making KCNQ1 a constitutively open channel. The stabilization of S4 arises from salt bridge interactions between R228 on S4 and residues D54 and D55 on KCNE3 [89] and possibly by another polar contact between Q244 on KCNQ1 and R83 on KCNE3 [26]. Our Rosetta model of the AO state supports this mechanistic explanation by showing R228 and Q244 at surface-accessible sites at the extracellular and cytosolic ends of S4 where those residues extend into the inter-subunit crevice formed between the VSD and PD. Binding of KCNE3 to R228 on the extracellular side and to Q244 on the cytosolic side is suggested to lock the S4 helix in place and stabilize the activated VSD conformation. Recent studies also suggest an additional effect of KCNE3 on the PD [81] which can be investigated with the new KCNQ1 models e.g. through docking and MD simulations.

Potential use of KCNQ1 models for personalized medicine

For a variety of genetic human disorders, substantial evidence suggests that a majority of disease-causing mutations act by destabilizing the folded conformation of the encoded protein [91, 92] thereby reducing the population of the functional form of the affected protein in cells. Predicting the phenotype of a mutation and interpreting variants of unknown significance discovered by genome sequencing data can greatly benefit from methods to quantitatively compare the stabilities of WT and mutant forms of proteins. For water-soluble proteins, computational algorithms are able to predict mutation-induced stability changes with reasonable accuracy [93]. However, these methods perform significantly worse when applied to membrane proteins [94], highlighting the pressing need to develop methods with improved performance for membrane proteins.

In light of personalized treatment strategies for LQTS, the KCNQ1 Rosetta models will be broadly useful for interpreting the functional phenotype of KCNQ1 variants that we demonstrated here by predicting the energetic change in protein stability upon mutation compared to WT KCNQ1. Of the 44 single-site missense mutations studied here, 25 were LOF mutations (channel current <65%WT). Huang et al. [19] found 18 of these LOF mutations exhibited lower (<65%WT) cell surface expression than WT and 13 were severely folding-destabilized as indicated by their NMR spectrum or the inability of the VSD to express in E. coli. For these subsets of variants, we found the median ΔΔG was 17.7 REU and 24.2 REU, respectively, which translates to 13 of 18 and 11 of 13 variants having largely increased ΔΔG values above 10 REU (Fig 7C). In contrast, the median ΔΔG of the 19 non-LOF variants was significantly lower (0.7 REU, p<0.001), and for the 26 variants with WT-like cell surface expression the ΔΔG was 1.1 REU (p<0.001). Thus, the Rosetta ΔΔG predictions support the conclusion that the KCNQ1 LOF variants were energetically destabilized and that protein destabilization represents the likely cause for protein misfolding and low cell surface expression. We previously noticed [19] that many class IV and V variants, which are folding-destabilized, have mutation sites located in S0 and S0-contacting regions (see S7 Fig). The Rosetta models and MD simulations revealed a network of interactions of residues in S0 with residues in S2 and S4 that appear to glue the VSD helix bundle together. This observation illuminates the role of S0 as an important scaffolding element in KCNQ1 and potentially also in other KV channels.

Energetic destabilization and trafficking defects is a common disease mechanism associated with gene mutations that impact membrane proteins [95]. Examples include the HERG potassium channel and cardiac arrhythmia [96], the human peripheral membrane myelin protein 22 (PMP22) and Charcot-Marie Tooth disease [97], and a subset of mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) protein [98]. It is likely that mutation-induced destabilization is also a trigger for the disease-causative intracellular retention and/or degradation of other membrane proteins e.g. the myelin proteolipid protein that causes Pelizaeus-Merzbacher disease [99].

Noticeable exceptions of variants for which Rosetta-predicted changes in protein thermodynamic stability were not correlated with KCNQ1 loss- or gain-of-function were variants in class I (e.g. C122Y) as well as A102S (WT-like) and R109L (class III). C122Y is a LOF variant and is associated with LQTS [100]. For this case, the reason for LOF is unknown. In our KCNQ1 models, C122 is close to the putative KCNE1 binding site in KCNQ1 between S1 and S5’ of two neighboring subunits. One may speculate mutation-induced changes of the local structure around C122 could affect the presentation of side-chains on S1 towards KCNE1 and therefore affect KCNE1 binding. Alternatively, the C122Y mutation could affect S4 movement through an allosteric mechanism. Five other mutations in class I (with normal or higher expression and trafficking levels but low or no conductance) are located in regions that are critical for KCNE1 binding and channel gating: on S1 (L134P), S4 (I227L, Q234P) and on S4-contacting regions (A150, K218). L134 is lipid-exposed and its side-chain is pointed towards the putative KCNE1 binding site. I227 is one position before the first gating charge R1 while Q234 occupies the third gating charge position (R3 in other channels, Q3 in KCNQ1). A150 is located in the S1-S2 loop and K218 in the S3-S4 loop, both of which undergo varying interactions with S4 during VSD activation. Thus, the unusual behavior of these class I variants may be related to the location-specific role these residues play in voltage-dependent channel activation mechanism.

The reason for increased conductance of A102S and R109L is unknown. Both sites are located in the first half of S0 and are lipid exposed. A structural mechanistic explanation for them is currently hampered by the lack of structural data for the N-terminal domain before S0, which is not part of the Rosetta models.

Knowledge of KCNQ1 LOF mechanisms pertaining to a given patient may soon facilitate the design of personalized treatment strategies for LQTS. The observation that many LOF variants are energetically destabilized and prone to protein misfolding motivates the idea of developing targeted therapies that assist in folding of destabilized KCNQ1 variants and/or rescue trafficking-defective variants to reach the plasma membrane. The development of such tailored therapies for varying classes of channel defects has become reality for the cystic fibrosis transmembrane regulator (CFTR) protein. Many patient-associated CFTR mutations are trafficking-defective, some of which have been shown to be rescuable by a pharmacological approach [95, 101]. The structural models developed in this study provide a step toward realizing such tailored pharmacological approaches for KCNQ1.

Conclusion

The structural models of the human KCNQ1 channel in closed and open states were developed by comparative modeling based on the high-grade template structure from X. leavis KCNQ1 and other experimental template structures that together provided nearly complete coverage of all channel regions (S0-S6). The KCNQ1 models are experimentally validated and offer improved completeness and accuracy over earlier KCNQ1 Rosetta models. MD simulations demonstrated the capability of the models to describe experimentally established channel properties such as state-dependent VSD gating charge interactions and pore conformations, PIP2 binding, and VSD-PD coupling interactions. Rosetta energy calculations on KCNQ1 models successfully predicted energetic destabilizations in folding-defective KCNQ1 LOF mutants and distinguished them from mutants with normal phenotype. The results in this work highlight the utility of the KCNQ1 structural models in addressing open questions about KCNQ1 structure-function mechanisms and assisting in the diagnosis of disease-causative KCNQ1 mutations towards the development of new therapeutic approaches.

Materials and methods

Homology modeling of KCNQ1 RC state and AO state structures

Structural models of human KCNQ1 (helix S0-S6, residues 100–369) in the RC, AC and AO state were generated by homology modeling using the protein structure prediction software Rosetta (version 3.8) [102]. The X. leavis KCNQ1 cryo-EM structure (PDB 5VMS) [31] provided a template for the AC state as well as separate templates for the VSD in the AO model and for the pore domain in the RC model. Additional templates for the open pore conformation in AO state KCNQ1 were derived from the crystal structure of rat KV1.2–2.1 channel (PDB 2R9R) [35]. Templates for the VSD in RC state KCNQ1 were obtained from related VSD structures: the resting VSD in C. intestinalis voltage-sensing phosphatase (Ci-VSP) (PDB 4G7Y) [32], VSD2 in A. thaliana two pore calcium channel protein 1 (TPC1) (PDB 5DQQ) [33], and the resting VSD conformation C3 in a model of the Shaker K channel [36]. Hybridization and regularization of multiple templates was accomplished with the Rosetta comparative modeling (RosettaCM) protocol [103] guided by the RosettaMembrane energy function [104, 105]. Missing loop regions and connecting residues between the VSD and pore domain were modeled de novo through fragment insertion. Template threading in RosettaCM was guided by multiple sequence alignments created with ClustalW [106] (see S1 and S2 Figs). Fragment libraries were created with the Rosetta fragment picker [107] incorporating PSIPRED [108] secondary structure prediction. Transmembrane helix regions were predicted by OCTOPUS [109] and used to impose membrane-specific Rosetta energy terms on residues within the theoretical membrane bilayer. During homology modeling, C4 symmetry of the KCNQ1 tetramer was enforced by means of Rosetta symmetry definition files as described in DiMaio et al. [110]. Amino acid side chain positions were optimized by a simulated annealing protocol, referred to as rotamer packing in Rosetta, and models were refined in internal and Cartesian coordinate space by gradient-based minimization while applying a low harmonic restraint to the initial Cα-atom coordinates. Models were clustered based on Cα-atom RMSD, and the top-scoring models of the ten largest clusters as well as another 10–20 models with low Rosetta score and valid helix and loop conformations were selected as final ensemble. Model quality was additionally checked by MolProbity [38] analysis (see supporting S1 Table), and the model with the lowest Molprobity score was considered the final model.

MD simulations of KCNQ1 RC and AO models

MD simulations of KCNQ1 RC and AO models were performed in explicit phospholipid membranes at 310 K using Amber16 [48] and the Lipid17 force field (Gould, I.R., Skjevik A.A., Dickson, C.J., Madej, B.D., Walker, R.C.,"Lipid17: A Comprehensive AMBER Force Field for the Simulation of Zwitterionic and Anionic Lipids", 2018, manuscript in preparation). Representative starting conformations for MD were selected from the ensemble of KCNQ1 homology models and by hierarchical full-linkage clustering.

Models were aligned to the membrane normal using the PPM webserver [111] and embedded into bilayers of POPC (palmitoyloleoyl-phosphatidylcholine) and PIP2 (phosphatidyl-4,5-bisphosphate) (~280 lipids per leaflet) using the membrane builder tool of the CHARMM-GUI website [47]. A TIP3P water layer with 22.5 Å thickness containing 150 mM of KCl was added on either side of the membrane. In addition, four K+ ions were placed in the channel selectivity filter at positions inferred from the X-ray structure coordinates of PDB 2R9R. Bilayers contained 10 mol% of PIP2 in the inner leaflet which comprised equal numbers of C4-PO4- and C5-PO4-mono-protonated PIP2 molecules with stearoyl and arachidonoyl conjugations at the sn-1 and sn-2 position. The geometry of PIP2 was optimized with Gaussian 09 (Gaussian, Inc., Wallingford CT) on the B3LYP/6-31G* level of theory, and assignment of Amber atom types and calculation of RESP charges was done with Antechamber [112]. Bond and angle parameters of the protonated C4-PO4 or C5-PO4 group in PIP2 were adjusted to values previously reported for phosphorylated amino acids [113] to avoid simulation instabilities. No PIP2 was added to the MD system of the AC model because the decoupled state of KCNQ1 in the cryo-EM structure is assumed to be due to the absence of PIP2 and no conformational change should be induced by the membrane. SHAKE [114] bond length constraints were applied to all bonds involving hydrogen. Nonbonded interactions were evaluated with a 10 Å cutoff, and electrostatic interactions were calculated by the particle-mesh Ewald method [115].

Each MD system was first minimized for 15,000 steps using steepest descent followed by 15,000 steps of conjugate gradient minimization. With protein and ions restrained to their initial coordinates, the lipid and water were heated to 50 K over 1000 steps with a step size of 1 fs in the NVT ensemble using Langevin dynamics with a rapid collision frequency of 10,000 ps-1. The system was then heated to 100 K over 50,000 steps with a collision frequency of 1000 ps-1 and finally to 310 K over 200,000 steps and a collision frequency of 100 ps-1. After changing to the NPT ensemble, restraints on ions were gradually removed over 500 ps and the system was equilibrated for another 5 ns at 310 K with weak positional restraints (with a force constant of 1 kcal mol-1 Å-2) applied to protein Cα atoms. The protein restraints were then gradually removed over 10 ns, and production MD was conducted for 350 ns (AC model) or 400 ns (RC and AO model) using a step size of 2 fs, constant pressure periodic boundary conditions, anisotropic pressure scaling and Langevin dynamics. Four independent simulations were carried out for each channel model yielding 1.4–1.6 μs of MD data for each state.

Analysis of KCNQ1 MD simulations

Analysis of MD trajectories with CPPTRAJ (version 18.0) [116] included calculation of Cα-atom root-mean-square deviations (Cα-RMSD) and root-mean-square fluctuations (Cα-RMSF), enumeration of protein-protein and protein-lipid hydrogen bonds and KCNQ1 residue contacts as well as principal component analysis. Residue contact numbers were calculated by counting the number of heteroatoms within a 4 Å distance of any heteroatom in a given KCNQ1 residue averaged over the last 300 ns of production MD as well as over all four RC and AO state MD simulations, respectively. For calculating intra-subunit contact numbers, only residues within the same chain were used whereas inter-subunit contact numbers were calculated considering only residues located on two different, directly adjacent chains. Intra- and inter-subunit contacts numbers were averaged over all four chains (i.e. A–D) or pairs of adjacent chains (i.e. A+B, B+C, C+D, A+D) in the KCNQ1 tetramer, respectively. Pairs of diagonally opposite chains (i.e. A+C and B+D) were not included in the analysis because they make no contacts between the VSD and PD. Furthermore, only long-range contacts (|i-j| > 5 residues) were used in the analysis of intra-subunit interactions.

Measurement of the channel pore radius was carried out with the HOLE program [117] and using snapshots of KCNQ1 taken at 1 ns intervals during the last 300 ns of MD.

Stability calculations of KCNQ1 VSD variants

Free energy changes in 50 functionally characterized KCNQ1 variants [19] with mutations in the VSD were modeled using the Rosetta Flex ddG protocol [74] and the RosettaMembrane all-atom energy function [105]. In short, the Flex ddG protocol models mutation-induced conformational and energetic changes by a series of “backrub” moves (15,000 steps in this study) of the protein backbone together with side-chain repacking within an 8 Å shell around the mutation site, followed by global minimization of all protein backbone and side-chain torsion angles. Minimization is performed with harmonic Cα-atom pair distance restraints to avoid any large structural deviations from the input model. A mutation was modeled in each of the four channel subunits simultaneously by enforcing C4 symmetry by means of Rosetta symmetry definition files. 50 independent trajectories were carried out for each mutant model and the WT model, and the Rosetta energy change (ΔΔG) was calculated as average score difference between the three top-scoring mutant and WT models. Stability calculations were repeated on ten different input models of the RC as well as AO state taken from the selected ensemble of KCNQ1 homology models to calculate the average ΔΔG value (± S.E.M.) for each KCNQ1 variant.

Supporting information

S1 Table [docx]
Molprobity statistic of KCNQ1 AO, RC and AC homology models.

S2 Table [docx]
Comparison of Rosetta predicted stability changes (ΔΔG) with expression levels, trafficking efficiencies and peak current densities of KCNQ1 VSD mutants.

S1 Fig [tpc1]
Multiple sequence alignment of human KCNQ1 for RC homology modeling.

S2 Fig [tiff]
Multiple sequence alignment of human KCNQ1 for AO homology modeling.

S3 Fig [tiff]
Pore radius of KCNQ1 in MD simulations of RC, AC and AO models.

S4 Fig [ac]
Backbone Cα-RMSD of KCNQ1 RC, AC and AO state MD trajectories.

S5 Fig [a]
Global domain motions in KCNQ1 channel models revealed by residue cross-correlation and principal component analysis.

S6 Fig [a]
Conformational changes in KCNQ1 during transition from the AC to AO state inferred from modeling.

S7 Fig [tiff]
Location of mutation sites in the KCNQ1 VSD by variant class.

S8 Fig [a]
Predicted stability changes of KCNQ1 VSD variants calculated with the RC model.

S9 Fig [tiff]
Comparison of Rosetta ΔΔG values (mean ± S.E.M.) for KCNQ1 variants calculated with the RC and AO model.

S10 Fig [a]
Comparison of Rosetta RC and AO models of selected KCNQ1 mutants.

S11 Fig [a]
Structural models of selected KCNQ1 variants.

S12 Fig [a]
Experimentally restrained docking of KCNE1 to KCNQ1 closed and open channel models.

S1 PDB [pdb]
Coordinates of the Rosetta KCNQ1 AO model in PDB format.

S2 PDB [pdb]
Coordinates of the Rosetta KCNQ1 RC model in PDB format.

S3 PDB [pdb]
Coordinates of the Rosetta KCNQ1 AC model in PDB format.

S1 Movie [mp4]
Morph between the Rosetta KCNQ1 RC and AO channel models as viewed parallel to the membrane plane.

S2 Movie [mp4]
Morph between the Rosetta KCNQ1 RC and AO channel models as viewed from the extracellular side.

S3 Movie [mp4]
Morph between the Rosetta KCNQ1 RC and AO channel models as viewed from the cytosolic side.

S4 Movie [mp4]
Morph between the closed and open pore domain of KCNQ1.

S5 Movie [mp4]
Morph between the resting and activated VSD of KCNQ1.

S1 Input Files and Scripts []
This directory contains input data and scripts that were used to create the results in this manuscript.

S1 MDtraj [nc]
Snapshots from 400ns MD trajectory of KCNQ1 AO model in Amber netcdf format.

S2 MDtraj [nc]
Snapshots from 400ns MD trajectory of KCNQ1 RC model in Amber netcdf format.

S3 MDtraj [nc]
Snapshots from 350ns MD trajectory of KCNQ1 AC model in Amber netcdf format.


Zdroje

1. Abbott GW. Biology of the KCNQ1 Potassium Channel. New Journal of Science. 2014;2014:26. doi: 10.1155/2014/237431

2. Bezanilla F, Stefani E. Voltage-dependent gating of ionic channels. Annu Rev Biophys Biomol Struct. 1994;23:819–46. doi: 10.1146/annurev.bb.23.060194.004131 7522668.

3. Jensen MO, Jogini V, Borhani DW, Leffler AE, Dror RO, Shaw DE. Mechanism of voltage gating in potassium channels. Science. 2012;336(6078):229–33. doi: 10.1126/science.1216533 22499946.

4. Bezanilla F. The voltage sensor in voltage-dependent ion channels. Physiol Rev. 2000;80(2):555–92. doi: 10.1152/physrev.2000.80.2.555 10747201.

5. Wu D, Delaloye K, Zaydman MA, Nekouzadeh A, Rudy Y, Cui J. State-dependent electrostatic interactions of S4 arginines with E1 in S2 during Kv7.1 activation. J Gen Physiol. 2010;135(6):595–606. doi: 10.1085/jgp.201010408 20479111; PubMed Central PMCID: PMC2888051.

6. Zaydman MA, Kasimova MA, McFarland K, Beller Z, Hou P, Kinser HE, et al. Domain-domain interactions determine the gating, permeation, pharmacology, and subunit modulation of the IKs ion channel. eLife. 2014;3:e03606–e. doi: 10.7554/eLife.03606 25535795

7. Hou P, Eldstrom J, Shi J, Zhong L, McFarland K, Gao Y, et al. Inactivation of KCNQ1 potassium channels reveals dynamic coupling between voltage sensing and pore opening. Nat Commun. 2017;8(1):1730. doi: 10.1038/s41467-017-01911-8 29167462; PubMed Central PMCID: PMC5700111.

8. Sachyani D, Dvir M, Strulovich R, Tria G, Tobelaim W, Peretz A, et al. Structural basis of a Kv7.1 potassium channel gating module: studies of the intracellular c-terminal domain in complex with calmodulin. Structure. 2014;22(11):1582–94. doi: 10.1016/j.str.2014.07.016 25441029.

9. Tobelaim WS, Dvir M, Lebel G, Cui M, Buki T, Peretz A, et al. Ca(2+)-Calmodulin and PIP2 interactions at the proximal C-terminus of Kv7 channels. Channels (Austin). 2017;11(6):686–95. doi: 10.1080/19336950.2017.1388478 28976808; PubMed Central PMCID: PMC5786183.

10. Zheng R, Thompson K, Obeng-Gyimah E, Alessi D, Chen J, Cheng H, et al. Analysis of the interactions between the C-terminal cytoplasmic domains of KCNQ1 and KCNE1 channel subunits. Biochem J. 2010;428(1):75–84. doi: 10.1042/BJ20090977 20196769; PubMed Central PMCID: PMC2888147.

11. Tobelaim WS, Dvir M, Lebel G, Cui M, Buki T, Peretz A, et al. Competition of calcified calmodulin N lobe and PIP2 to an LQT mutation site in Kv7.1 channel. Proc Natl Acad Sci U S A. 2017;114(5):E869–E78. doi: 10.1073/pnas.1612622114 28096388; PubMed Central PMCID: PMC5293103.

12. Barhanin J, Lesage F, Guillemare E, Fink M, Lazdunski M, Romey G. K(v)LQT1 and IsK (minK) proteins associate to form the I-Ks cardiac potassium current. Nature. 1996;384(6604):78–80. ISI:A1996VR21900058. doi: 10.1038/384078a0 8900282

13. Sanguinetti MC, Curran ME, Zou A, Shen J, Spector PS, Atkinson DL, et al. Coassembly of K(v)LQT1 and minK (IsK) proteins to form cardiac I-Ks potassium channel. Nature. 1996;384(6604):80–3. ISI:A1996VR21900059. doi: 10.1038/384080a0 8900283

14. Hedley PL, Jørgensen P, Schlamowitz S, Wangari R, Moolman-Smook J, Brink PA, et al. The genetic basis of long QT and short QT syndromes: a mutation update. Hum Mutat. 2009;30:1486–511. doi: 10.1002/humu.21106 19862833

15. Wu J, Ding WG, Horie M. Molecular pathogenesis of long QT syndrome type 1. J Arrhythm. 2016;32(5):381–8. doi: 10.1016/j.joa.2015.12.006 27761162; PubMed Central PMCID: PMC5063268.

16. Eldstrom J, Wang Z, Werry D, Wong N, Fedida D. Microscopic mechanisms for long QT syndrome type 1 revealed by single-channel analysis of I(Ks) with S3 domain mutations in KCNQ1. Heart Rhythm. 2015;12:386–94. doi: 10.1016/j.hrthm.2014.10.029 25444851

17. Vanoye CG, Desai RR, Fabre KL, Gallagher SL, Potet F, DeKeyser JM, et al. High-Throughput Functional Evaluation of KCNQ1 Decrypts Variants of Unknown Significance. (2574–8300 (Electronic)).

18. Li B, Mendenhall JL, Kroncke BM, Taylor KC, Huang H, Smith DK, et al. Predicting the Functional Impact of KCNQ1 Variants of Unknown Significance. Circ Cardiovasc Genet. 2017;10(5). doi: 10.1161/CIRCGENETICS.117.001754 29021305; PubMed Central PMCID: PMC5679743.

19. Huang H, Kuenze G, Smith JA, Taylor KC, Duran AM, Hadziselimovic A, et al. Mechanisms of KCNQ1 channel dysfunction in long QT syndrome involving voltage sensor domain mutations. Sci Adv. 2018;4(3):eaar2631. doi: 10.1126/sciadv.aar2631 29532034; PubMed Central PMCID: PMC5842040.

20. Smith JA, Vanoye CG, George AL, Meiler J, Sanders CR. Structural Models for the KCNQ1 Voltage-Gated Potassium Channel. Biochemistry. 2007;46(49):14141–52. Epub 2007/11/15. doi: 10.1021/bi701597s 17999538; PubMed Central PMCID: PMC2565492.

21. Kang C, Tian C, Sonnichsen FD, Smith JA, Meiler J, George AL Jr, et al. Structure of KCNE1 and implications for how it modulates the KCNQ1 potassium channel. Biochemistry. 2008;47(31):7999–8006. Epub 2008/07/10. doi: 10.1021/bi800875q 18611041; PubMed Central PMCID: PMC2580054.

22. Chung DY, Chan PJ, Bankston JR, Yang L, Liu G, Marx SO, et al. Location of KCNE1 relative to KCNQ1 in the I(KS) potassium channel by disulfide cross-linking of substituted cysteines. Proc Natl Acad Sci U S A. 2009;106:743–8. doi: 10.1073/pnas.0811897106 19131515

23. Chan PJ, Osteen JD, Xiong D, Bohnen MS, Doshi D, Sampson KJ, et al. Characterization of KCNQ1 atrial fibrillation mutations reveals distinct dependence on KCNE1. J Gen Physiol. 2012;139:135–44. doi: 10.1085/jgp.201110672 22250012

24. Nakajo K, Kubo Y. Steric hindrance between S4 and S5 of the KCNQ1/KCNE1 channel hampers pore opening. Nat Commun. 2014;5:4100–. doi: 10.1038/ncomms5100 24920132

25. Li P, Liu H, Lai C, Sun P, Zeng W, Wu F, et al. Differential modulations of KCNQ1 by auxiliary proteins KCNE1 and KCNE2. Sci Rep. 2014;4:4973. doi: 10.1038/srep04973 24827085; PubMed Central PMCID: PMC4021338.

26. Kroncke BM, Van Horn WD, Smith J, Kang C, Welch RC, Song Y, et al. Structural basis for KCNE3 modulation of potassium recycling in epithelia. Sci Adv. 2016;2:e1501228–e. doi: 10.1126/sciadv.1501228 27626070

27. Zaydman MA, Silva JR, Delaloye K, Li Y, Liang H, Larsson HP, et al. Kv7.1 ion channels require a lipid to couple voltage sensing to pore opening. Proc Natl Acad Sci U S A. 2013;110:13180–5. doi: 10.1073/pnas.1305167110 23861489

28. Eckey K, Wrobel E, Strutz-Seebohm N, Pott L, Schmitt N, Seebohm G. Novel Kv7.1-phosphatidylinositol 4,5-bisphosphate interaction sites uncovered by charge neutralization scanning. J Biol Chem. 2014;289(33):22749–58. doi: 10.1074/jbc.M114.589796 24947509; PubMed Central PMCID: PMC4132781.

29. Restier L, Cheng L, Sanguinetti MC. Mechanisms by which atrial fibrillation-associated mutations in the S1 domain of KCNQ1 slow deactivation of IKs channels. J Physiol. 2008;586:4179–91. doi: 10.1113/jphysiol.2008.157511 18599533

30. Labro AJ, Boulet IR, Choveau FS, Mayeur E, Bruyns T, Loussouarn G, et al. The S4-S5 linker of KCNQ1 channels forms a structural scaffold with the S6 segment controlling gate closure. J Biol Chem. 2011;286(1):717–25. doi: 10.1074/jbc.M110.146977 21059661; PubMed Central PMCID: PMC3013030.

31. Sun J, MacKinnon R. Cryo-EM Structure of a KCNQ1/CaM Complex Reveals Insights into Congenital Long QT Syndrome. Cell. 2017;169:1042–50.e9. doi: 10.1016/j.cell.2017.05.019 28575668

32. Li Q, Wanderling S, Paduch M, Medovoy D, Singharoy A, McGreevy R, et al. Structural mechanism of voltage-dependent gating in an isolated voltage-sensing domain. Nat Struct Mol Biol. 2014;21:244–52. doi: 10.1038/nsmb.2768 24487958

33. Kintzer AF, Stroud RM. Structure, inhibition and regulation of two-pore channel TPC1 from Arabidopsis thaliana. Nature. 2016;531:258–62. doi: 10.1038/nature17194 26961658

34. Yarov-Yarovoy V, Baker D, Catterall WA. Voltage sensor conformations in the open and closed states in ROSETTA structural models of K(+) channels. Proc Natl Acad Sci U S A. 2006;103(19):7292–7. doi: 10.1073/pnas.0602350103 16648251; PubMed Central PMCID: PMC1464335.

35. Long SB, Tao X, Campbell EB, MacKinnon R. Atomic structure of a voltage-dependent K+ channel in a lipid membrane-like environment. Nature. 2007;450(7168):376–82. doi: 10.1038/nature06265 18004376.

36. Henrion U, Renhorn J, Borjesson SI, Nelson EM, Schwaiger CS, Bjelkmar P, et al. Tracking a complete voltage-sensor cycle with metal-ion bridges. Proc Natl Acad Sci U S A. 2012;109(22):8552–7. doi: 10.1073/pnas.1116938109 22538811; PubMed Central PMCID: PMC3365220.

37. Laskowski RA, MacArthur MW, Moss DS, Thornton JM. PROCHECK: a program to check the stereochemical quality of protein structures. Journal of Applied Crystallography. 1993;26(2):283–91. doi: 10.1107/S0021889892009944

38. Davis IW, Leaver-Fay A, Chen VB, Block JN, Kapral GJ, Wang X, et al. MolProbity: all-atom contacts and structure validation for proteins and nucleic acids. Nucleic Acids Res. 2007;35(Web Server issue):W375–83. doi: 10.1093/nar/gkm216 17452350; PubMed Central PMCID: PMC1933162.

39. Seebohm G, Strutz-Seebohm N, Ureche ON, Baltaev R, Lampert A, Kornichuk G, et al. Differential roles of S6 domain hinges in the gating of KCNQ potassium channels. Biophys J. 2006;90(6):2235–44. doi: 10.1529/biophysj.105.067165 16326905; PubMed Central PMCID: PMC1386802.

40. Glauner KS, Mannuzzu LM, Gandhi CS, Isacoff EY. Spectroscopic mapping of voltage sensor movement in the Shaker potassium channel. Nature. 1999;402(6763):813–7. doi: 10.1038/45561 10617202.

41. Delemotte L, Tarek M, Klein ML, Amaral C, Treptow W. Intermediate states of the Kv1.2 voltage sensor from atomistic molecular dynamics simulations. Proc Natl Acad Sci U S A. 2011;108:6109–14. doi: 10.1073/pnas.1102724108 21444776

42. Cui J. Voltage-Dependent Gating: Novel Insights from KCNQ1 Channels. Biophys J. 2016;110(1):14–25. doi: 10.1016/j.bpj.2015.11.023 26745405; PubMed Central PMCID: PMC4805877.

43. Panaghie G, Abbott GW. The role of S4 charges in voltage-dependent and voltage-independent KCNQ1 potassium channel complexes. J Gen Physiol. 2007;129(2):121–33. doi: 10.1085/jgp.200609612 17227916; PubMed Central PMCID: PMC2154355.

44. Schmidt D, Jiang QX, MacKinnon R. Phospholipids and the origin of cationic gating charges in voltage sensors. Nature. 2006;444(7120):775–9. doi: 10.1038/nature05416 17136096.

45. Butterwick JA, MacKinnon R. Solution structure and phospholipid interactions of the isolated voltage-sensor domain from KvAP. J Mol Biol. 2010;403(4):591–606. doi: 10.1016/j.jmb.2010.09.012 20851706; PubMed Central PMCID: PMC2971526.

46. Hite RK, Butterwick JA, MacKinnon R. Phosphatidic acid modulation of Kv channel voltage sensor function. Elife. 2014;3. doi: 10.7554/eLife.04366 25285449; PubMed Central PMCID: PMC4212207.

47. Wu EL, Cheng X, Jo S, Rui H, Song KC, Davila-Contreras EM, et al. CHARMM-GUI Membrane Builder toward realistic biological membrane simulations. J Comput Chem. 2014;35(27):1997–2004. doi: 10.1002/jcc.23702 25130509; PubMed Central PMCID: PMC4165794.

48. Case DA, Betz RM, Cerutti DS, Cheatham Iii TE, Darden TA, Duke RE, et al. AMBER 2016. San Francisco: University of California; 2016.

49. Xu Y, Wang Y, Meng X-Y, Zhang M, Jiang M, Cui M, et al. Building KCNQ1/KCNE1 channel models and probing their interactions by molecular-dynamics simulations. Biophys J. 2013;105:2461–73. doi: 10.1016/j.bpj.2013.09.058 24314077

50. Gofman Y, Shats S, Attali B, Haliloglu T, Ben-Tal N. How does KCNE1 regulate the Kv7.1 potassium channel? Model-structure, mutations, and dynamics of the Kv7.1-KCNE1 complex. Structure. 2012;20:1343–52. doi: 10.1016/j.str.2012.05.016 22771213

51. Yang T, Smith JA, Leake BF, Sanders CR, Meiler J, Roden DM. An allosteric mechanism for drug block of the human cardiac potassium channel KCNQ1. Mol Pharmacol. 2013;83(2):481–9. Epub 2012/11/30. doi: 10.1124/mol.112.081513 23193163; PubMed Central PMCID: PMC3558809.

52. Chen L, Zhang Q, Qiu Y, Li Z, Chen Z, Jiang H, et al. Migration of PIP2 lipids on voltage-gated potassium channel surface influences channel deactivation. Sci Rep. 2015;5:15079. doi: 10.1038/srep15079 26469389; PubMed Central PMCID: PMC4606798.

53. Suh BC, Hille B. Regulation of ion channels by phosphatidylinositol 4,5-bisphosphate. Curr Opin Neurobiol. 2005;15(3):370–8. doi: 10.1016/j.conb.2005.05.005 15922587.

54. Dickson CJ, Madej BD, Skjevik AA, Betz RM, Teigen K, Gould IR, et al. Lipid14: The Amber Lipid Force Field. J Chem Theory Comput. 2014;10(2):865–79. doi: 10.1021/ct4010307 24803855; PubMed Central PMCID: PMC3985482.

55. Kasimova MA, Zaydman MA, Cui J, Tarek M. PIP2-dependent coupling is prominent in Kv7.1 due to weakened interactions between S4-S5 and S6. Sci Rep. 2015;5:7474–. doi: 10.1038/srep07474 25559286

56. Itoh H, Shimizu W, Hayashi K, Yamagata K, Sakaguchi T, Ohno S, et al. Long QT syndrome with compound mutations is associated with a more severe phenotype: a Japanese multicenter study. Heart Rhythm. 2010;7(10):1411–8. doi: 10.1016/j.hrthm.2010.06.013 20541041.

57. Wisten A, Bostrom IM, Morner S, Stattin EL. Mutation analysis of cases of sudden unexplained death, 15 years after death: prompt genetic evaluation after resuscitation can save future lives. Resuscitation. 2012;83(10):1229–34. doi: 10.1016/j.resuscitation.2012.05.015 22659597.

58. Wang Q, Curran ME, Splawski I, Burn TC, Millholland JM, VanRaay TJ, et al. Positional cloning of a novel potassium channel gene: KVLQT1 mutations cause cardiac arrhythmias. Nat Genet. 1996;12:17–23. doi: 10.1038/ng0196-17 8528244

59. Kapplinger JD, Tester DJ, Salisbury BA, Carr JL, Harris-Kerr C, Pollevick GD, et al. Spectrum and prevalence of mutations from the first 2,500 consecutive unrelated patients referred for the FAMILION long QT syndrome genetic test. Heart Rhythm. 2009;6(9):1297–303. doi: 10.1016/j.hrthm.2009.05.021 19716085; PubMed Central PMCID: PMC3049907.

60. Napolitano C, Priori SG, Schwartz PJ, Bloise R, Ronchetti E, Nastoli J, et al. Genetic testing in the long QT syndrome: development and validation of an efficient approach to genotyping in clinical practice. JAMA. 2005;294(23):2975–80. doi: 10.1001/jama.294.23.2975 16414944.

61. Millat G, Chevalier P, Restier-Miron L, Da Costa A, Bouvagnet P, Kugener B, et al. Spectrum of pathogenic mutations and associated polymorphisms in a cohort of 44 unrelated patients with long QT syndrome. Clin Genet. 2006;70(3):214–27. doi: 10.1111/j.1399-0004.2006.00671.x 16922724.

62. Kubota T, Shimizu W, Kamakura S, Horie M. Hypokalemia-induced long QT syndrome with an underlying novel missense mutation in S4-S5 linker of KCNQ1. J Cardiovasc Electrophysiol. 2000;11(9):1048–54. 11021476.

63. Choi G, Kopplin LJ, Tester DJ, Will ML, Haglund CM, Ackerman MJ. Spectrum and frequency of cardiac channel defects in swimming-triggered arrhythmia syndromes. Circulation. 2004;110(15):2119–24. doi: 10.1161/01.CIR.0000144471.98080.CA 15466642.

64. Boulet IR, Labro AJ, Raes AL, Snyders DJ. Role of the S6 C-terminus in KCNQ1 channel gating. J Physiol. 2007;585(Pt 2):325–37. doi: 10.1113/jphysiol.2007.145813 17932138; PubMed Central PMCID: PMC2375480.

65. Choveau FS, Rodriguez N, Abderemane Ali F, Labro AJ, Rose T, Dahimene S, et al. KCNQ1 channels voltage dependence through a voltage-dependent binding of the S4-S5 linker to the pore domain. J Biol Chem. 2011;286(1):707–16. doi: 10.1074/jbc.M110.146324 20940310; PubMed Central PMCID: PMC3013029.

66. Ma LJ, Ohmert I, Vardanyan V. Allosteric features of KCNQ1 gating revealed by alanine scanning mutagenesis. Biophys J. 2011;100(4):885–94. doi: 10.1016/j.bpj.2010.12.3726 21320432; PubMed Central PMCID: PMC3037572.

67. Lu Z, Klem AM, Ramu Y. Ion conduction pore is conserved among potassium channels. Nature. 2001;413(6858):809–13. doi: 10.1038/35101535 11677598.

68. Long SB, Campbell EB, Mackinnon R. Voltage sensor of Kv1.2: structural basis of electromechanical coupling. Science. 2005;309(5736):903–8. doi: 10.1126/science.1116270 16002579.

69. Lu Z, Klem AM, Ramu Y. Coupling between voltage sensors and activation gate in voltage-gated K+ channels. J Gen Physiol. 2002;120(5):663–76. doi: 10.1085/jgp.20028696 12407078.

70. Thompson E, Eldstrom J, Fedida D. Single channel kinetic analysis of the cAMP effect on IKs mutants, S209F and S27D/S92D. Channels (Austin). 2018;12(1):276–83. doi: 10.1080/19336950.2018.1499369 30027808.

71. Zhou X, Bueno-Orovio A, Schilling RJ, Kirkby C, Denning C, Rajamohan D, et al. Investigating the Complex Arrhythmic Phenotype Caused by the Gain-of-Function Mutation KCNQ1-G229D. Front Physiol. 2019;10:259. doi: 10.3389/fphys.2019.00259 30967788; PubMed Central PMCID: PMC6430739.

72. Werry D, Eldstrom J, Wang Z, Fedida D. Single-channel basis for the slow activation of the repolarizing cardiac potassium current, I(Ks). Proc Natl Acad Sci U S A. 2013;110:E996–1005.

73. Silva JR, Pan H, Wu D, Nekouzadeh A, Decker KF, Cui J, et al. A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential. Proc Natl Acad Sci U S A. 2009;106(27):11102–6. doi: 10.1073/pnas.0904505106 19549851; PubMed Central PMCID: PMC2700153.

74. Barlow KA, S OC, Thompson S, Suresh P, Lucas JE, Heinonen M, et al. Flex ddG: Rosetta Ensemble-Based Estimation of Changes in Protein-Protein Binding Affinity upon Mutation. J Phys Chem B. 2018;122(21):5389–99. doi: 10.1021/acs.jpcb.7b11367 29401388; PubMed Central PMCID: PMC5980710.

75. Davis IW, Arendall WB 3rd, Richardson DC, Richardson JS. The backrub motion: how protein backbone shrugs when a sidechain dances. Structure. 2006;14(2):265–74. Epub 2006/02/14. S0969-2126(06)00040-2 [pii] doi: 10.1016/j.str.2005.10.007 16472746.

76. Alford RF, Koehler Leman J, Weitzner BD, Duran AM, Tilley DC, Elazar A, et al. An Integrated Framework Advancing Membrane Protein Modeling and Design. PLoS Comput Biol. 2015;11(9):e1004398. doi: 10.1371/journal.pcbi.1004398 26325167; PubMed Central PMCID: PMC4556676.

77. Tinel N, Diochot S, Borsotto M, Lazdunski M, Barhanin J. KCNE2 confers background current characteristics to the cardiac KCNQ1 potassium channel. EMBO J. 2000;19(23):6326–30. doi: 10.1093/emboj/19.23.6326 11101505; PubMed Central PMCID: PMC305874.

78. Schroeder BC, Waldegger S, Fehr S, Bleich M, Warth R, Greger R, et al. A constitutively open potassium channel formed by KCNQ1 and KCNE3. Nature. 2000;403(6766):196–9. doi: 10.1038/35003200 10646604.

79. Nakajo K, Kubo Y. KCNQ1 channel modulation by KCNE proteins via the voltage-sensing domain. J Physiol. 2015;593:2617–25. doi: 10.1113/jphysiol.2014.287672 25603957

80. Ruscic KJ, Miceli F, Villalba-Galea CA, Dai H, Mishina Y, Bezanilla F, et al. IKs channels open slowly because KCNE1 accessory subunits slow the movement of S4 voltage sensors in KCNQ1 pore-forming subunits. Proc Natl Acad Sci U S A. 2013;110(7):E559–66. doi: 10.1073/pnas.1222616110 23359697; PubMed Central PMCID: PMC3574954.

81. Barro-Soria R, Ramentol R, Liin SI, Perez ME, Kass RS, Larsson HP. KCNE1 and KCNE3 modulate KCNQ1 channels by affecting different gating transitions. Proc Natl Acad Sci U S A. 2017;114(35):E7367–E76. doi: 10.1073/pnas.1710335114 28808020; PubMed Central PMCID: PMC5584458.

82. Westhoff M, Eldstrom J, Murray CI, Thompson E, Fedida D. I Ks ion-channel pore conductance can result from individual voltage sensor movements. Proc Natl Acad Sci U S A. 2019;116(16):7879–88. doi: 10.1073/pnas.1811623116 30918124; PubMed Central PMCID: PMC6475427.

83. Barro-Soria R, Rebolledo S, Liin SI, Perez ME, Sampson KJ, Kass RS, et al. KCNE1 divides the voltage sensor movement in KCNQ1/KCNE1 channels into two steps. Nat Commun. 2014;5:3750. doi: 10.1038/ncomms4750 24769622; PubMed Central PMCID: PMC4019390.

84. Tapper AR, George AL. Location and orientation of minK within the I-Ks potassium channel complex. Journal of Biological Chemistry. 2001;276(41):38249–54. ISI:000171526500072. doi: 10.1074/jbc.M103956200 11479291

85. Xu X, Jiang M, Hsu K-L, Zhang M, Tseng G-N. KCNQ1 and KCNE1 in the IKs channel complex make state-dependent contacts in their extracellular domains. J Gen Physiol. 2008;131:589–603. doi: 10.1085/jgp.200809976 18504315

86. Wang YH, Jiang M, Xu XL, Hsu K-L, Zhang M, Tseng G-N. Gating-related molecular motions in the extracellular domain of the IKs channel: implications for IKs channelopathy. J Membr Biol. 2011;239:137–56. doi: 10.1007/s00232-010-9333-7 21152909

87. Strutz-Seebohm N, Pusch M, Wolf S, Stoll R, Tapken D, Gerwert K, et al. Structural basis of slow activation gating in the cardiac I Ks channel complex. Cell Physiol Biochem. 2011;27:443–52. doi: 10.1159/000329965 21691061

88. Wang Y, Zhang M, Xu Y, Jiang M, Zankov DP, Cui M, et al. Probing the structural basis for differential KCNQ1 modulation by KCNE1 and KCNE2. J Gen Physiol. 2012;140(6):653–69. doi: 10.1085/jgp.201210847 23183700; PubMed Central PMCID: PMC3514736.

89. Barro-Soria R, Perez ME, Larsson HP. KCNE3 acts by promoting voltage sensor activation in KCNQ1. Proc Natl Acad Sci U S A. 2015;112(52):E7286–92. doi: 10.1073/pnas.1516238112 26668384; PubMed Central PMCID: PMC4703023.

90. Nakajo K, Kubo Y. KCNE1 and KCNE3 stabilize and/or slow voltage sensing S4 segment of KCNQ1 channel. J Gen Physiol. 2007;130(3):269–81. doi: 10.1085/jgp.200709805 17698596; PubMed Central PMCID: PMC2151641.

91. Stefl S, Nishi H, Petukh M, Panchenko AR, Alexov E. Molecular mechanisms of disease-causing missense mutations. J Mol Biol. 2013;425(21):3919–36. doi: 10.1016/j.jmb.2013.07.014 23871686; PubMed Central PMCID: PMC3796015.

92. Kroncke BM, Vanoye CG, Meiler J, George AL Jr., Sanders CR. Personalized biochemistry and biophysics. Biochemistry. 2015;54(16):2551–9. doi: 10.1021/acs.biochem.5b00189 25856502; PubMed Central PMCID: PMC4415889.

93. Berliner N, Teyra J, Colak R, Garcia Lopez S, Kim PM. Combining structural modeling with ensemble machine learning to accurately predict protein fold stability and binding affinity effects upon mutation. PLoS One. 2014;9(9):e107353. doi: 10.1371/journal.pone.0107353 25243403; PubMed Central PMCID: PMC4170975.

94. Kroncke BM, Duran AM, Mendenhall JL, Meiler J, Blume JD, Sanders CR. Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability. Biochemistry. 2016;55(36):5002–9. doi: 10.1021/acs.biochem.6b00537 27564391; PubMed Central PMCID: PMC5024705.

95. Marinko JT, Huang H, Penn WD, Capra JA, Schlebach JP, Sanders CR. Folding and Misfolding of Human Membrane Proteins in Health and Disease: From Single Molecules to Cellular Proteostasis. Chem Rev. 2019. doi: 10.1021/acs.chemrev.8b00532 30608666.

96. Smith JL, Anderson CL, Burgess DE, Elayi CS, January CT, Delisle BP. Molecular pathogenesis of long QT syndrome type 2. J Arrhythm. 2016;32(5):373–80. doi: 10.1016/j.joa.2015.11.009 27761161; PubMed Central PMCID: PMC5063260.

97. Schlebach JP, Narayan M, Alford C, Mittendorf KF, Carter BD, Li J, et al. Conformational Stability and Pathogenic Misfolding of the Integral Membrane Protein PMP22. J Am Chem Soc. 2015;137(27):8758–68. doi: 10.1021/jacs.5b03743 26102530; PubMed Central PMCID: PMC4507940.

98. Farinha CM, Matos P, Amaral MD. Control of cystic fibrosis transmembrane conductance regulator membrane trafficking: not just from the endoplasmic reticulum to the Golgi. FEBS J. 2013;280(18):4396–406. doi: 10.1111/febs.12392 23773658.

99. Swanton E, Holland A, High S, Woodman P. Disease-associated mutations cause premature oligomerization of myelin proteolipid protein in the endoplasmic reticulum. Proc Natl Acad Sci U S A. 2005;102(12):4342–7. doi: 10.1073/pnas.0407287102 15753308; PubMed Central PMCID: PMC555485.

100. Tester DJ, Will ML, Haglund CM, Ackerman MJ. Compendium of cardiac channel mutations in 541 consecutive unrelated patients referred for long QT syndrome genetic testing. Heart Rhythm. 2005;2(5):507–17. doi: 10.1016/j.hrthm.2005.01.020 15840476.

101. Fajac I, De Boeck K. New horizons for cystic fibrosis treatment. Pharmacol Ther. 2017;170:205–11. doi: 10.1016/j.pharmthera.2016.11.009 27916649.

102. Bender BJ, Cisneros A 3rd, Duran AM, Finn JA, Fu D, Lokits AD, et al. Protocols for Molecular Modeling with Rosetta3 and RosettaScripts. Biochemistry. 2016;55(34):4748–63. doi: 10.1021/acs.biochem.6b00444 27490953; PubMed Central PMCID: PMC5007558.

103. Song Y, DiMaio F, Wang RY, Kim D, Miles C, Brunette T, et al. High-resolution comparative modeling with RosettaCM. Structure. 2013;21(10):1735–42. doi: 10.1016/j.str.2013.08.005 24035711; PubMed Central PMCID: PMC3811137.

104. Yarov-Yarovoy V, Schonbrun J, Baker D. Multipass membrane protein structure prediction using Rosetta. Proteins. 2006;62(4):1010–25. doi: 10.1002/prot.20817 16372357; PubMed Central PMCID: PMC1479309.

105. Barth P, Schonbrun J, Baker D. Toward high-resolution prediction and design of transmembrane helical protein structures. Proc Natl Acad Sci U S A. 2007;104(40):15682–7. doi: 10.1073/pnas.0702515104 17905872; PubMed Central PMCID: PMC2000396.

106. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8. Epub 2007/09/12. doi: 10.1093/bioinformatics/btm404 [pii] 17846036.

107. Gront D, Kulp DW, Vernon RM, Strauss CE, Baker D. Generalized fragment picking in Rosetta: design, protocols and applications. PLoS One. 2011;6(8):e23294. doi: 10.1371/journal.pone.0023294 21887241; PubMed Central PMCID: PMC3160850.

108. Jones DT. Protein Secondary Structure Prediction Based on Position-specific Scoring Matrices. J Mol Biol. 1999;292(2):195–202. doi: 10.1006/jmbi.1999.3091 10493868

109. Viklund H, Elofsson A. OCTOPUS: improving topology prediction by two-track ANN-based preference scores and an extended topological grammar. Bioinformatics. 2008;24(15):1662–8. Epub 2008/05/14. doi: 10.1093/bioinformatics/btn221 18474507.

110. Dimaio F, Leaver-Fay A, Bradley P, Baker D, Andre I. Modeling symmetric macromolecular structures in rosetta3. PLoS One. 2011;6(6):e20450. Epub 2011/07/07. doi: 10.1371/journal.pone.0020450 PONE-D-11-04312 [pii]. 21731614; PubMed Central PMCID: PMC3120754.

111. Lomize MA, Pogozheva ID, Joo H, Mosberg HI, Lomize AL. OPM database and PPM web server: resources for positioning of proteins in membranes. Nucleic Acids Res. 2012;40(Database issue):D370–6. doi: 10.1093/nar/gkr703 21890895; PubMed Central PMCID: PMC3245162.

112. Wang J, Wang W, Kollman PA, Case DA. Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model. 2006;25(2):247–60. doi: 10.1016/j.jmgm.2005.12.005 16458552.

113. Homeyer N, Horn AH, Lanig H, Sticht H. AMBER force-field parameters for phosphorylated amino acids in different protonation states: phosphoserine, phosphothreonine, phosphotyrosine, and phosphohistidine. J Mol Model. 2006;12(3):281–9. doi: 10.1007/s00894-005-0028-4 16240095.

114. Ryckaert JP, Ciccotti G, Berendsen HJC. Numerical-Integration of Cartesian Equations of Motion of a System with Constraints—Molecular-Dynamics of N-Alkanes. J Comput Phys. 1977;23(3):327–41. doi: 10.1016/0021-9991(77)90098-5 WOS:A1977CZ25300007.

115. Darden T, York D, Pedersen L. Particle Mesh Ewald—an N.Log(N) Method for Ewald Sums in Large Systems. Journal of Chemical Physics. 1993;98(12):10089–92. doi: 10.1063/1.464397 WOS:A1993LG10100091.

116. Roe DR, Cheatham TE 3rd. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J Chem Theory Comput. 2013;9(7):3084–95. doi: 10.1021/ct400341p 26583988.

117. Smart OS, Neduvelil JG, Wang X, Wallace BA, Sansom MS. HOLE: a program for the analysis of the pore dimensions of ion channel structural models. J Mol Graph. 1996;14(6):354–60, 76. 9195488.

118. Ortiz AR, Strauss CEM, Olmea O. MAMMOTH (Matching molecular models obtained from theory): An automated method for model comparison. Protein Sci. 2002;11:2606–11. doi: 10.1110/ps.0215902 12381844

119. Duran AM, Meiler J. Computational design of membrane proteins using RosettaMembrane. Protein Sci. 2018;27(1):341–55. doi: 10.1002/pro.3335 29090504; PubMed Central PMCID: PMC5734395.


Článek vyšel v časopise

PLOS One


2019 Číslo 9