Membrane anchoring and interaction between transmembrane domains are crucial for K+ channel function

Yeast complementation assays provide limited information on functional features of K+ channels

We investigate to what extent yeast complementation assays, which in principle can provide large amounts of training data for machine-learning models, yield quantitative correlations between growth rescue and single-channel recordings. If this were the case, yeast complementation results could be used as surrogate data for machine-learning-based channel design. Therefore, we mutated position L94 at the cavity entry of the model K+ channel KcvPBCV1 to all proteinogenic amino acids. The function of the wild-type channel and its mutants was investigated by reconstituting them in planar lipid bilayers and by their ability to rescue the growth of a yeast strain deficient in K+ uptake. The single-channel data show a distinct effect of mutations in this critical position on unitary conductance and open probability, with no apparent causal relationship between the two functional parameters. We also found that even conservative amino acid replacements can alter the unitary conductance and/or open probability and that most functional changes show no systematic relationship with the physicochemical nature of the amino acids. This emphasizes that the functional influence of an amino acid on channel function cannot be reduced to a single chemical property. Mutual comparison of single-channel data and yeast complementation results exhibit only a partial correlation between their electrical parameters and their potency of rescuing growth. Hence, complementation data alone are not sufficient for enabling functional channel design; they need to be complemented by additional parameters such as the number of channels in the plasma membrane.

Potassium viroporins as model systems for understanding eukaryotic ion channel behaviour

Ion channels are membrane proteins essential for a plethora of cellular functions including maintaining cell shape, ion homeostasis, cardiac rhythm and action potential in neurons. The complexity and often extensive structure of eukaryotic membrane proteins makes it difficult to understand their basic biological regulation. Therefore, this article suggests, viroporins - the miniature versions of eukaryotic protein homologs from viruses - might serve as model systems to provide insights into behaviour of eukaryotic ion channels in general. The structural requirements for correct assembly of the channel along with the basic functional properties of a K+ channel exist in the minimal design of the viral K+ channels from two viruses, Chlorella virus (Kcv) and Ectocarpus siliculosus virus (Kesv). These small viral proteins readily assemble into tetramers and they sort in cells to distinct target membranes. When these viruses-encoded channels are expressed into the mammalian cells, they utilise their protein machinery and hence can serve as excellent tools to study the cells protein sorting machinery. This combination of small size and robust function makes viral K+ channels a valuable model system for detection of basic structure-function correlations. It is believed that molecular and physiochemical analyses of these viroporins may serve as basis for the development of inhibitors or modulators to ion channel activity for targeting ion channel diseases - so called channelopathies. Therefore, it may provide a potential different scope for molecular pharmacology studies aiming at novel and innovative therapeutics associated with channel related diseases. This article reviews the structural and functional properties of Kcv and Kesv upon expression in mammalian cells and Xenopus oocytes. The mechanisms behind differential protein sorting in Kcv and Kesv are also thoroughly discussed.

Inferring functional units in ion channel pores via relative entropy

Coarse-grained protein models approximate the first-principle physical potentials. Among those modeling approaches, the relative entropy framework yields promising and physically sound results, in which a mapping from the target protein structure and dynamics to a model is defined and subsequently adjusted by an entropy minimization of the model parameters. Minimization of the relative entropy is equivalent to maximization of the likelihood of reproduction of (configurational ensemble) observations by the model. In this study, we extend the relative entropy minimization procedure beyond parameter fitting by a second optimization level, which identifies the optimal mapping to a (dimension-reduced) topology. We consider anisotropic network models of a diverse set of ion channels and assess our findings by comparison to experimental results.

Assessing the extent of the structural and dynamic modulation of membrane lipids due to pore forming toxins: insights from molecular dynamics simulations

Infections caused by many virulent bacterial strains are triggered by the release of pore forming toxins (PFTs), which form oligomeric transmembrane pore complexes on the target plasma membrane. The spatial extent of the perturbation to the surrounding lipids during pore formation is relatively unexplored. Using all-atom molecular dynamics simulations, we investigate the changes in the structure and dynamics of lipids in a 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipid bilayer in the presence of contrasting PFTs. Cytolysin A (ClyA), an α toxin with its inserted wedge shaped bundle of inserted α helices, induces significant asymmetry across the membrane leaflets in comparison with α hemolysin (AHL), a β toxin. Despite the differences in hydrophobic mismatch and uniquely different topologies of the two oligomers, perturbations to lipid order as reflected in the tilt angle and order parameters and membrane thinning are short ranged, lying within ∼2.5 nm from the periphery of either pore complex, and commensurate with distances typically associated with van der Waals forces. In contrast, the spatial extent of perturbations to the lipid dynamics extends outward to at least 4 nm for both proteins, and the continuous survival probabilities reveal the presence of a tightly bound shell of lipids in this region. Displacement probability distributions show long tails and the distinctly non-Gaussian features reflect the induced dynamic heterogeneity. A detailed profiling of the protein-lipid contacts with tyrosine, tryptophan, lysine and arginine residues shows increased non-polar contacts in the cytoplasmic leaflet for both PFTs, with a higher number of atomic contacts in the case of AHL in the extracellular leaflet due to the mushroom-like topology of the pore complex. The short ranged nature of the perturbations observed in this simple one component membrane suggests inherent plasticity of membrane lipids enabling the recovery of the structure and membrane fluidity even in the presence of these large oligomeric transmembrane protein assemblies. This observation has implications in membrane repair processes such as budding or vesicle fusion events used to mitigate PFT virulence, where the underlying lipid dynamics and fluidity in the vicinity of the pore complex are expected to play an important role.

Role of an Aromatic-Aromatic Interaction in the Assembly and Trafficking of the Zebrafish Panx1a Membrane Channel

Pannexin 1 (Panx1) is a ubiquitously expressed hexameric integral membrane protein known to function as an adenosine triphosphate (ATP) release channel. Panx1 proteins exist in unglycosylated core form (Gly0). They undergo critical post-translational modifications forming the high mannose glycosylation state (Gly1) in the endoplasmic reticulum (ER) and the complex glycosylation state (Gly2) in the Golgi apparatus. The regulation of transition from the ER to the cell membrane is not fully understood. Using site-specific mutagenesis, dye uptake assays, and interaction testing, we identified two conserved aromatic residues, Trp123 and Tyr205, in the transmembrane domains 2 and 3 of the zebrafish panx1a protein. Results suggest that both residues primarily govern the assembly of panx1a subunits into channels, with mutant proteins failing to interact. The results provide insight into a mechanism enabling regulation of Panx1 oligomerization, glycosylation, and trafficking.

Genes for Membrane Transport Proteins: Not So Rare in Viruses

Some viruses have genes encoding proteins with membrane transport functions. It is unknown if these types of proteins are rare or are common in viruses. In particular, the evolutionary origin of some of the viral genes is obscure, where other viral proteins have homologs in prokaryotic and eukaryotic organisms. We searched virus genomes in databases looking for transmembrane proteins with possible transport function. This effort led to the detection of 18 different types of putative membrane transport proteins indicating that they are not a rarity in viral genomes. The most abundant proteins are K⁺ channels. Their predicted structures vary between different viruses. With a few exceptions, the viral proteins differed significantly from homologs in their current hosts. In some cases the data provide evidence for a recent gene transfer between host and virus, but in other cases the evidence indicates a more complex evolutionary history.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

Exploring the Viral Channel KcvPBCV-1 Function via Computation

Viral potassium channels (Kcv) are homologous to the pore module of complex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {K}^+$$\end{document}K+-selective ion channels of cellular organisms. Due to their relative simplicity, they have attracted interest towards understanding the principles of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {K}^+$$\end{document}K+ conduction and channel gating. In this work, we construct a homology model of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {Kcv}}_{\text{PBCV-1}}$$\end{document}KcvPBCV-1 open state, which we validate by studying the binding of known blockers and by monitoring ion conduction through the channel. Molecular dynamics simulations of this model reveal that the re-orientation of selectivity filter carbonyl groups coincides with the transport of potassium ions, suggesting a possible mechanism for fast gating. In addition, we show that the voltage sensitivity of this mechanism can originate from the relocation of potassium ions inside the selectivity filter. We also explore the interaction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hbox {Kcv}}_{\text{PBCV-1}}$$\end{document}KcvPBCV-1 with the surrounding bilayer and observe the binding of lipids in the area between two adjacent subunits. The model is available to the scientific community to further explore the structure/function relationship of Kcv channels.

Conversion of an instantaneous activating K+ channel into a slow activating inward rectifier

The miniature channel, Kcv, is a structural equivalent of the pore of all K⁺ channels. Here, we follow up on a previous observation that a largely voltage-insensitive channel can be converted into a slow activating inward rectifier after extending the outer transmembrane domain by one Ala. This gain of rectification can be rationalized by dynamic salt bridges at the cytosolic entrance to the channel; opening is favored by voltage-sensitive formation of salt bridges and counteracted by their disruption. Such latent voltage sensitivity in the pore could be relevant for the understanding of voltage gating in complex Kv channels.

Viruses infecting marine picoplancton encode functional potassium ion channels

Phycodnaviruses are dsDNA viruses, which infect algae. Their large genomes encode many gene products, like small K(+) channels, with homologs in prokaryotes and eukaryotes. Screening for K(+) channels revealed their abundance in viruses from fresh-water habitats. Recent sequencing of viruses from marine algae or from salt water in Antarctica revealed sequences with the predicted characteristics of K(+) channels but with some unexpected features. Two genes encode either 78 or 79 amino acid proteins, which are the smallest known K(+) channels. Also of interest is an unusual sequence in the canonical α-helixes in K(+) channels. Structural prediction algorithms indicate that the new channels have the conserved α-helix folds but the algorithms failed to identify the expected transmembrane domains flanking the K(+) channel pores. In spite of these unexpected properties electophysiological studies confirmed that the new proteins are functional K(+) channels.

The sorting of a small potassium channel in mammalian cells can be shifted between mitochondria and plasma membrane

The two small and similar viral K(+) channels Kcv and Kesv are sorted in mammalian cells and yeast to different destinations. Analysis of the sorting pathways shows that Kcv is trafficking via the secretory pathway to the plasma membrane, while Kesv is inserted via the TIM/TOM complex to the inner membrane of mitochondria. Studies with Kesv mutants show that an N-terminal mitochondrial targeting sequence in this channel is neither necessary nor sufficient for sorting of Kesv the mitochondria. Instead the sorting of Kesv can be redirected from the mitochondria to the plasma membrane by an insertion of ≥2 amino acids in a position sensitive manner into the C-terminal transmembrane domain (TMD2) of this channel. The available data advocate the presence of a C-terminal sorting signal in TMD2 of Kesv channel, which is presumably not determined by the length of this domain.

Tethered protein display identifies a novel Kir3.2 (GIRK2) regulator from protein scaffold libraries

Use of randomized peptide libraries to evolve molecules with new functions provides a means for developing novel regulators of protein activity. Despite the demonstrated power of such approaches for soluble targets, application of this strategy to membrane systems, such as ion channels, remains challenging. Here, we have combined libraries of a tethered protein scaffold with functional selection in yeast to develop a novel activator of the G-protein-coupled mammalian inwardly rectifying potassium channel Kir3.2 (GIRK2). We show that the novel regulator, denoted N5, increases Kir3.2 (GIRK2) basal activity by inhibiting clearance of the channel from the cellular surface rather than affecting the core biophysical properties of the channel. These studies establish the tethered protein display strategy as a means to create new channel modulators and highlight the power of approaches that couple randomized libraries with direct selections for functional effects. Our results further underscore the possibility for the development of modulators that influence channel function by altering cell surface expression densities rather than by direct action on channel biophysical parameters. The use of tethered library selection strategies coupled with functional selection bypasses the need for a purified target and is likely to be applicable to a range of membrane protein systems.

Creation of a reactive oxygen species-insensitive Kcv channel

The current of the minimal viral K(+) channel Kcv(PCBV-1) heterologously expressed in Xenopus oocytes is strongly inhibited by reactive oxygen species (ROS) like H(2)O(2). Possible targets for the ROS effect are two cysteines (C53 and C79) and four methionines (M1, M15, M23, and M26). The C53A/C79A and M23L/M26L double mutations maintained the same ROS sensitivity as the wild type. However, M15L as a single mutant or in combination with C53A/C79A and/or M23L/M26L caused a complete loss of sensitivity to H(2)O(2). These results indicate a prominent role of M15 at the cytosolic end of the outer transmembrane helix for gating of Kcv(PCBV-1). Furthermore, even though the channel lacks a canonical voltage sensor, it exhibits a weak voltage sensitivity, resulting in a slight activation in the millisecond range after a voltage step to negative potentials. The M15L mutation inverts this kinetics into an inactivation, underlining the critical role of this residue for gating. The negative slope of the I-V curves of M15L is the same as in the wild type, indicating that the selectivity filter is not involved.

Tuning the tetraethylammonium sensitivity of potassium channel Kcv by subunit combination

Tetraethylammonium (TEA) is a potassium (K⁺) channel inhibitor that has been extensively used as a molecular probe to explore the structure of channels’ ion pathway. In this study, we identified that Leu70 of the virus-encoded potassium channel Kcv is a key amino acid that plays an important role in regulating the channel’s TEA sensitivity. Site-directed mutagenesis of Leu70 can change the TEA sensitivity by 1,000-fold from ∼100 µM to ∼100 mM. Because no compelling trends exist to explain this amino acid’s specific interaction with TEA, the role of Leu70 at the binding site is likely to ensure an optimal conformation of the extracellular mouth that confers high TEA affinity. We further assembled the subunits of mutant and wt-Kcv into a series of heterotetramers. The differences in these heterochannels suggest that all of the four subunits in a Kcv channel additively participate in the TEA binding, and each of the four residues at the binding site independently contributes an equal binding energy. We therefore can present a series of mutant/wild-type tetramer combinations that can probe TEA over three orders of magnitude in concentration. This study may give insight into the mechanism for the interaction between the potassium channel and its inhibitor.

Relevance of lysine snorkeling in the outer transmembrane domain of small viral potassium ion channels

Transmembrane domains (TMDs) are often flanked by Lys or Arg because they keep their aliphatic parts in the bilayer and their charged groups in the polar interface. Here we examine the relevance of this so-called "snorkeling" of a cationic amino acid, which is conserved in the outer TMD of small viral K(+) channels. Experimentally, snorkeling activity is not mandatory for Kcv(PBCV-1) because K29 can be replaced by most of the natural amino acids without any corruption of function. Two similar channels, Kcv(ATCV-1) and Kcv(MT325), lack a cytosolic N-terminus, and neutralization of their equivalent cationic amino acids inhibits their function. To understand the variable importance of the cationic amino acids, we reanalyzed molecular dynamics simulations of Kcv(PBCV-1) and N-terminally truncated mutants; the truncated mutants mimic Kcv(ATCV-1) and Kcv(MT325). Structures were analyzed with respect to membrane positioning in relation to the orientation of K29. The results indicate that the architecture of the protein (including the selectivity filter) is only weakly dependent on TMD length and protonation of K29. The penetration depth of Lys in a given protonation state is independent of the TMD architecture, which leads to a distortion of shorter proteins. The data imply that snorkeling can be important for K(+) channels; however, its significance depends on the architecture of the entire TMD. The observation that the most severe N-terminal truncation causes the outer TMD to move toward the cytosolic side suggests that snorkeling becomes more relevant if TMDs are not stabilized in the membrane by other domains.

Functional HAK/KUP/KT-like potassium transporter encoded by chlorella viruses

Chlorella viruses are a source of interesting membrane transport proteins. Here we examine a putative K(+) transporter encoded by virus FR483 and related chlorella viruses. The protein shares sequence and structural features with HAK/KUP/KT-like K(+) transporters from plants, bacteria and fungi. Yeast complementation assays and Rb(+) uptake experiments show that the viral protein, termed HAKCV (high-affinity K(+) transporter of chlorella virus), is functional, with transport characteristics that are similar to those of known K(+) transporters. Expression studies revealed that the protein is expressed as an early gene during viral replication, and proteomics data indicate that it is not packaged in the virion. The function of HAKCV is unclear, but the data refute the hypothesis that the transporter acts as a substitute for viral-encoded K(+) channels during virus infection.