One-shot design elevates functional expression levels of a voltage-gated potassium channel.

Membrane proteins play critical physiological roles as receptors, channels, pumps, and transporters. Despite their importance, however, low expression levels often hamper the experimental characterization of membrane proteins. We present an automated and web-accessible design algorithm called mPROSS (https://mPROSS.weizmann.ac.il), which uses phylogenetic analysis and an atomistic potential, including an empirical lipophilicity scale, to improve native-state energy. As a stringent test, we apply mPROSS to the Kv1.2-Kv2.1 paddle chimera voltage-gated potassium channel. Four designs, encoding 9-26 mutations relative to the parental channel, were functional and maintained potassium-selective permeation and voltage dependence in Xenopus oocytes with up to 14-fold increase in whole-cell current densities. Additionally, single-channel recordings reveal no significant change in the channel-opening probability nor in unitary conductance, indicating that functional expression levels increase without impacting the activity profile of individual channels. Our results suggest that the expression levels of other dynamic channels and receptors may be enhanced through one-shot design calculations.

Effect of Terahertz Electromagnetic Field on the Permeability of Potassium Channel Kv1.2.

In this paper, the influence of external terahertz electromagnetic fields with different frequencies of 4 THz, 10 THz, 15 THz, and 20 THz on the permeability of the Kv1.2 voltage-gated potassium ion channel on the nerve cell membrane was studied using the combined model of the "Constant Electric Field-Ion Imbalance" method by molecular dynamics. We found that although the applied terahertz electric field does not produce strong resonance with the -C=O groups of the conservative sequence T-V-G-Y-G amino acid residue of the selective filter (SF) of the channel, it would affect the stability of the electrostatic bond between potassium ions and the carbonyl group of T-V-G-Y-G of SF, and it would affect the stability of the hydrogen bond between water molecules and oxygen atoms of the hydroxyl group of the 374THR side chain at the SF entrance, changing the potential and occupied states of ions in the SF and the occurrence probability of the permeation mode of ions and resulting in the change in the permeability of the channel. Compared with no external electric field, when the external electric field with 15 THz frequency is applied, the lifetime of the hydrogen bond is reduced by 29%, the probability of the "soft knock on" mode is decreased by 46.9%, and the ion flux of the channel is activated by 67.7%. Our research results support the view that compared to "direct knock-on", "soft knock-on" is a slower permeation mode.

A Novel KCNA2 Variant in a Patient with Non-Progressive Congenital Ataxia and Epilepsy: Functional Characterization and Sensitivity to 4-Aminopyridine.

Kv1.2 channels, encoded by the KCNA2 gene, are localized in the central and peripheral nervous system, where they regulate neuronal excitability. Recently, heterozygous mutations in KCNA2 have been associated with a spectrum of symptoms extending from epileptic encephalopathy, intellectual disability, and cerebellar ataxia. Patients are treated with a combination of antiepileptic drugs and 4-aminopyridine (4-AP) has been recently trialed in specific cases. We identified a novel variant in KCNA2, E236K, in a Serbian proband with non-progressive congenital ataxia and early onset epilepsy, treated with sodium valproate. To ascertain the pathogenicity of E236K mutation and to verify its sensitivity to 4-AP, we transfected HEK 293 cells with Kv1.2 WT or E236K cDNAs and recorded potassium currents through the whole-cell patch-clamp. In silico analysis supported the electrophysiological data. E236K channels showed voltage-dependent activation shifted towards negative potentials and slower kinetics of deactivation and activation compared with Kv1.2 WT. Heteromeric Kv1.2 WT+E236K channels, resembling the condition of the heterozygous patient, confirmed a mixed gain- and loss-of-function (GoF/LoF) biophysical phenotype. 4-AP inhibited both Kv1.2 and E236K channels with similar potency. Homology modeling studies of mutant channels suggested a reduced interaction between the residue K236 in the S2 segment and the gating charges at S4. Overall, the biophysical phenotype of E236K channels correlates with the mild end of the clinical spectrum reported in patients with GoF/LoF defects. The response to 4-AP corroborates existing evidence that KCNA2-disorders could benefit from variant-tailored therapeutic approaches, based on functional studies.

Folding and misfolding of potassium channel monomers during assembly and tetramerization.

The dynamics and folding of potassium channel pore domain monomers are connected to the kinetics of tetramer assembly. In all-atom molecular dynamics simulations of Kv1.2 and KcsA channels, monomers adopt multiple nonnative conformations while the three helices remain folded. Consistent with this picture, NMR studies also find the monomers to be dynamic and structurally heterogeneous. However, a KcsA construct with a disulfide bridge engineered between the two transmembrane helices has an NMR spectrum with well-dispersed peaks, suggesting that the monomer can be locked into a native-like conformation that is similar to that observed in the folded tetramer. During tetramerization, fluoresence resonance energy transfer (FRET) data indicate that monomers rapidly oligomerize upon insertion into liposomes, likely forming a protein-dense region. Folding within this region occurs along separate fast and slow routes, with τfold ∼40 and 1,500 s, respectively. In contrast, constructs bearing the disulfide bond mainly fold via the faster pathway, suggesting that maintaining the transmembrane helices in their native orientation reduces misfolding. Interestingly, folding is concentration independent despite the tetrameric nature of the channel, indicating that the rate-limiting step is unimolecular and occurs after monomer association in the protein-dense region. We propose that the rapid formation of protein-dense regions may help with the assembly of multimeric membrane proteins by bringing together the nascent components prior to assembly. Finally, despite its name, the addition of KcsA's C-terminal "tetramerization" domain does not hasten the kinetics of tetramerization.

Multifaceted Regulation of Potassium-Ion Channels by Graphene Quantum Dots.

Graphene quantum dots (GQDs) are emerging as a versatile nanomaterial with numerous proposed biomedical applications. Despite the explosion in potential applications, the molecular interactions between GQDs and complex biomolecular systems, including potassium-ion (K+) channels, remain largely unknown. Here, we use molecular dynamics (MD) simulations and electrophysiology to study the interactions between GQDs and three representative K+ channels, which participate in a variety of physiological processes and are closely related to many disease states. Using MD simulations, we observed that GQDs adopt distinct contact poses with each of the three structurally distinct K+ channels. Our electrophysiological characterization of the effects of GQDs on channel currents revealed that GQDs interact with the extracellular voltage-sensing domain (VSD) of a Kv1.2 channel, augmenting current by left-shifting the voltage dependence of channel activation. In contrast, GQDs form a "lid" cluster over the extracellular mouth of inward rectifier Kir3.2, blocking the channel pore and decreasing the current in a concentration-dependent manner. Meanwhile, GQDs accumulate on the extracellular "cap domain" of K2P2 channels and have no apparent impact on the K+ flux through the channel. These results reveal a surprising multifaceted regulation of K+ channels by GQDs, which might help de novo design of nanomaterial-based channel probe openers/inhibitors that can be used to further discern channel function.

Structural and molecular insight into the pH-induced low-permeability of the voltage-gated potassium channel Kv1.2 through dewetting of the water cavity.

Understanding the gating mechanism of ion channel proteins is key to understanding the regulation of cell signaling through these channels. Channel opening and closing are regulated by diverse environmental factors that include temperature, electrical voltage across the channel, and proton concentration. Low permeability in voltage-gated potassium ion channels (Kv) is intimately correlated with the prolonged action potential duration observed in many acidosis diseases. The Kv channels consist of voltage-sensing domains (S1-S4 helices) and central pore domains (S5-S6 helices) that include a selectivity filter and water-filled cavity. The voltage-sensing domain is responsible for the voltage-gating of Kv channels. While the low permeability of Kv channels to potassium ion is highly correlated with the cellular proton concentration, it is unclear how an intracellular acidic condition drives their closure, which may indicate an additional pH-dependent gating mechanism of the Kv family. Here, we show that two residues E327 and H418 in the proximity of the water cavity of Kv1.2 play crucial roles as a pH switch. In addition, we present a structural and molecular concept of the pH-dependent gating of Kv1.2 in atomic detail, showing that the protonation of E327 and H418 disrupts the electrostatic balance around the S6 helices, which leads to a straightening transition in the shape of their axes and causes dewetting of the water-filled cavity and closure of the channel. Our work offers a conceptual advancement to the regulation of the pH-dependent gating of various voltage-gated ion channels and their related biological functions.

ImKTx96, a peptide blocker of the Kv1.2 ion channel from the venom of the scorpion Isometrus maculates.

Scorpion venom contains diverse bioactive peptides that can recognize and interact with membrane proteins such as ion channels. These natural toxins are believed to be useful tools for exploring the structure and function of ion channels. In this study, we characterized a K+-channel toxin gene, ImKTx96, from the venom gland cDNA library of the scorpion Isometrus maculates. The peptide deduced from the ImKTx96 precursor nucleotide sequence contains a signal peptide of 27 amino acid residues and a mature peptide of 29 residues with three disulfide bridges. Multiple sequence alignment indicated that ImKTx96 is similar with the scorpion toxins that typically target K+-channels. The recombined ImKTx96 peptide (rImKTx96) was expressed in the Escherichia coli system, and purified by GST-affinity chromatography and RP-HPLC. Results from whole-cell patch-clamp experiments revealed that rImKTx96 can inhibit the current of the Kv1.2 ion channel expressed in HEK293 cells. The 3D structure of ImKTx96 was constructed by molecular modeling, and the complex formed by ImKTx96 interacting with the Kv1.2 ion channel was obtained by molecular docking. Based on its structural features and pharmacological functions, ImKTx96 was identified as one member of K+-channel scorpion toxin α-KTx10 group and may be useful as a molecular probe for investigating the structure and function of the Kv1.2 ion channel.

Quantum Calculation of Proton and Other Charge Transfer Steps in Voltage Sensing in the Kv1.2 Channel.

Ion channels in cell membranes control entry and exit of ions; their gating (opening and closing) is key to their functioning. It is known that protons can pass through the voltage-sensing domain (VSD) of channels such as Kv1.2. Quantum calculations for a section of the VSD show the steps protons take in responding to voltage and show no major displacement of the protein backbone with voltage change; 70 amino acids are included, 42 with side chains (9 directly in the proton path), 28 as backbone only, and 24 water molecules. Protons provide much of the gating current, the capacitative current immediately preceding channel opening with significant additional contributions from charge transfer to other groups. Most gating models, in contrast, require major protein displacement during gating. Energy terms without classical analogues (exchange plus correlation energy, which are greater than thermal energy) show that quantum calculations are required. Energy as a function of voltage for a key proton transfer leads to, approximately, the correct voltage for channel opening. Calculated total charge transfer (not only protons) for gating is reasonable compared to experimental values. We are also able to account, at least qualitatively, for two mutations, one with the gating current curve left-shifted and one right-shifted, and show the alternate proton paths that are required to account for these.

Computational Protein-Protein Docking Reveals the Therapeutic Potential of Kunitz-type Venom against hKv1.2 Binding Sites.

BACKGROUND & OBJECTIVE: Kunitz-type venoms are bioactive proteins isolated from a wide variety of venomous animals. These venoms are involved in protease inhibitory activity or potassium channel blocking activity. Therefore, they are reported as an important source for lead drug candidates towards protease or channel associated diseases like neurological, metabolic and cardiovascular disorders. METHODS: This study aimed to check the inhibitory action of Kunitz-type venoms against potassium channels using computational tools. RESULTS: Among potassium channels, Human Voltage-Gated Potassium Channel 1.2 (hKv1.2) was used as a receptor whereas Kunitz-type peptides from the venoms of various species were selected as ligand dataset. CONCLUSION: This study helped in finding the binding interface between the receptor and ligand dataset for their potential therapeutic use in treating potassium channelopathies.

Molecular basis of Tityus stigmurus alpha toxin and potassium channel kV1.2 interactions.

The Tityus stigmurus scorpion is widely distributed in the Northeast of Brazil and is the main causal agent of human envenoming. The venom produced by this scorpion includes neurotoxins, which are peptides belonging to Family 2 toxins and are able to interact with ion channels. The KTx subfamily displays selectivity and affinity for Kv channel subtypes and the result of this interaction is the blockade of potassium channels, impairing vital functions. We report the optimized structural model of a transcript encoding a potassium channel blocker toxin from T. stigmurus. LC-MS analysis confirmed the presence of the toxin in the venom and the three-dimensional structure was obtained by computational homology modeling and refined by molecular dynamic simulations. Furthermore, docking simulations were performed using a Shaker kV-1.2 potassium channel from rats as receptor model and the contacts were identified revealing which amino acid residues and interactions could be involved in its blockade. These residues were mapped and their contact and electrostatic interactions were evaluated revealing the influence of positive lysine residues and the additional contribution of an asparagine to the stabilization of the complex, bringing new insights into the mechanism of action of this toxin.

Molecular Mechanism of Depolarization-Dependent Inactivation in W366F Mutant of Kv1.2.

Voltage-gated potassium channels play crucial roles in regulating membrane potential. They are activated by membrane depolarization, allowing the selective permeation of K+ ions across the plasma membrane, and enter a nonconducting state after lasting depolarization, a process known as inactivation. Inactivation in voltage-activated potassium channels occurs through two distinct mechanisms, N-type and C-type inactivation. C-type inactivation is caused by conformational changes in the extracellular mouth of the channel, whereas N-type inactivation is elicited by changes in the cytoplasmic mouth of the protein. The W434F-mutated Shaker channel is known as a nonconducting mutant and is in a C-type inactivation state at a depolarizing membrane potential. To clarify the structural properties of C-type inactivated protein, we performed molecular dynamics simulations of the wild-type and W366F (corresponding to W434F in Shaker) mutant of the Kv1.2-2.1 chimera channel. The W366F mutant was in a nearly nonconducting state with a depolarizing voltage and recovered from inactivation with a reverse voltage. Our simulations and three-dimensional reference interaction site model analysis suggested that structural changes in the selectivity filter upon membrane depolarization trap K+ ions around the inner mouth of the selectivity filter and prevent ion permeation. This pore restriction is involved in the molecular mechanism of C-type inactivation.

Single-particle cryo-EM structure of a voltage-activated potassium channel in lipid nanodiscs.

Voltage-activated potassium (Kv) channels open to conduct K+ ions in response to membrane depolarization, and subsequently enter non-conducting states through distinct mechanisms of inactivation. X-ray structures of detergent-solubilized Kv channels appear to have captured an open state even though a non-conducting C-type inactivated state would predominate in membranes in the absence of a transmembrane voltage. However, structures for a voltage-activated ion channel in a lipid bilayer environment have not yet been reported. Here we report the structure of the Kv1.2-2.1 paddle chimera channel reconstituted into lipid nanodiscs using single-particle cryo-electron microscopy. At a resolution of ~3 Å for the cytosolic domain and ~4 Å for the transmembrane domain, the structure determined in nanodiscs is similar to the previously determined X-ray structure. Our findings show that large differences in structure between detergent and lipid bilayer environments are unlikely, and enable us to propose possible structural mechanisms for C-type inactivation.

Gating interaction maps reveal a noncanonical electromechanical coupling mode in the Shaker K+ channel.

Membrane potential regulates the activity of voltage-dependent ion channels via specialized voltage-sensing modules, but the mechanisms involved in coupling voltage-sensor movement to pore opening remain unclear owing to a lack of resting state structures and robust methods to identify allosteric pathways. Here, using a newly developed interaction-energy analysis, we probe the interfaces of the voltage-sensing and pore modules in the Drosophila Shaker K+ channel. Our measurements reveal unexpectedly strong equilibrium gating interactions between contacts at the S4 and S5 helices in addition to those between S6 and the S4-S5 linker. Network analysis of MD trajectories shows that the voltage-sensor and pore motions are linked by two distinct pathways: a canonical pathway through the S4-S5 linker and a hitherto unknown pathway akin to rack-and-pinion coupling involving the S4 and S5 helices. Our findings highlight the central role of the S5 helix in electromechanical transduction in the voltage-gated ion channel (VGIC) superfamily.

Effect of Intercalated Water on Potassium Ion Transport through Kv1.2 Channels Studied via On-the-Fly Free-Energy Parametrization.

We introduce a two-dimensional version of the method called on-the-fly free energy parametrization (OTFP) to reconstruct free-energy surfaces using Molecular Dynamics simulations, which we name OTFP-2D. We first test the new method by reconstructing the well-known dihedral angles free energy surface of solvated alanine dipeptide. Then, we use it to investigate the process of K+ ions translocation inside the Kv1.2 channel. By comparing a series of two-dimensional free energy surfaces for ion movement calculated with different conditions on the intercalated water molecules, we first recapitulate the widely accepted knock-on mechanism for ion translocation and then confirm that permeation occurs with water molecules alternated among the ions, in accordance with the latest experimental findings. From a methodological standpoint, our new OTFP-2D algorithm demonstrates the excellent sampling acceleration of temperature-accelerated molecular dynamics and the ability to efficiently compute 2D free-energy surfaces. It will therefore be useful in large variety complex biomacromolecular simulations.

The role of entropic potential in voltage activation and K+ transport through Kv 1.2 channels.

We analyze the entropic effects of inner pore geometry changes of Kv 1.2 channel during membrane depolarization and their implications for the rate of transmembrane transport of potassium ions. We base this on the idea that spatial confinements within the channel pore give rise to entropic barriers which can both effectively affect the stability of open macroconformation and influence channel's ability to conduct the potassium ions through the membrane. First, we calculate the differences in entropy between voltage-activated and resting states of the channel. As a template, we take a set of structures of channel pore in an open state at different membrane potentials generated in our previous research. The obtained results indicate that tendency to occupy open states at membrane depolarization is entropy facilitated. Second, we describe the differences in rates of K+ transport through the channel pore at different voltages based on the results of appropriate random walk simulations in entropic and electric potentials. The simulated single channel currents (I) suggest that the geometry changes during membrane depolarization are an important factor contributing to the observed flow of potassium ions through the channel. Nevertheless, the charge distribution within the channel pore (especially at the extracellular entrance) seems most prominent for the observed I/Imax relation at a qualitative level at analyzed voltages.

Molecular Dynamics Simulation Reveals Specific Interaction Sites between Scorpion Toxins and Kv1.2 Channel: Implications for Design of Highly Selective Drugs.

The Kv1.2 channel plays an important role in the maintenance of resting membrane potential and the regulation of the cellular excitability of neurons, whose silencing or mutations can elicit neuropathic pain or neurological diseases (e.g., epilepsy and ataxia). Scorpion venom contains a variety of peptide toxins targeting the pore region of this channel. Despite a large amount of structural and functional data currently available, their detailed interaction modes are poorly understood. In this work, we choose four Kv1.2-targeted scorpion toxins (Margatoxin, Agitoxin-2, OsK-1, and Mesomartoxin) to construct their complexes with Kv1.2 based on the experimental structure of ChTx-Kv1.2. Molecular dynamics simulation of these complexes lead to the identification of hydrophobic patches, hydrogen-bonds, and salt bridges as three essential forces mediating the interactions between this channel and the toxins, in which four Kv1.2-specific interacting amino acids (D353, Q358, V381, and T383) are identified for the first time. This discovery might help design highly selective Kv1.2-channel inhibitors by altering amino acids of these toxins binding to the four channel residues. Finally, our results provide new evidence in favor of an induced fit model between scorpion toxins and K⁺ channel interactions.

Extracellular redox sensitivity of Kv1.2 potassium channels.

Kv1.2 is a prominent potassium channel subtype in the nervous system and serves as an important structural template for investigation of ion channel function. However, Kv1.2 voltage-dependence exhibits dramatic cell-to-cell variability due to a gating mode shift that is regulated by an unknown mechanism. We report that this variable behavior is regulated by the extracellular redox environment. Exposure to reducing agents promotes a shift in gating properties towards an 'inhibited' gating mode that resists opening, and causes channels to exhibit pronounced use-dependent activation during trains of repetitive depolarizations. This sensitivity to extracellular redox potential is absent in other Kv1 channels, but is apparent in heteromeric channels containing Kv1.2 subunits, and overlaps with the reported physiological range of extracellular redox couples. Mutagenesis of candidate cysteine residues fails to abolish redox sensitivity. Therefore, we suggest that an extrinsic, redox-sensitive binding partner imparts these properties.

Concentration-Dependent Binding of Small Ligands to Multiple Saturable Sites in Membrane Proteins.

Membrane proteins are primary targets for most therapeutic indications in cancer and neurological diseases, binding over 50% of all known small molecule drugs. Understanding how such ligands impact membrane proteins requires knowledge on the molecular structure of ligand binding, a reasoning that has driven relentless efforts in drug discovery and translational research. Binding of small ligands appears however highly complex involving interaction to multiple transmembrane protein sites featuring single or multiple occupancy states. Within this scenario, looking for new developments in the field, we investigate the concentration-dependent binding of ligands to multiple saturable sites in membrane proteins. The study relying on docking and free-energy perturbation provides us with an extensive description of the probability density of protein-ligand states that allows for computation of thermodynamic properties of interest. It also provides one- and three-dimensional spatial descriptions for the ligand density across the protein-membrane system which can be of interest for structural purposes. Illustration and discussion of the results are shown for binding of the general anesthetic sevoflurane against Kv1.2, a mammalian ion channel for which experimental data are available.

Gating Charge Calculations by Computational Electrophysiology Simulations.

Electrical cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the molecular details of voltage sensing remain insufficiently understood. Here, we present a molecular dynamics simulation-based method to determine the underlying charge movement across the membrane-the gating charge-by measuring electrical capacitor properties of membrane-embedded proteins. We illustrate the approach by calculating the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against experimental data. We then use the gating charge analysis to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with experiments, and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.