Molecular Dynamics Simulation of the Cytosolic Mouth in Kcv-Type Potassium Channels

HCN1 mutation spectrum: from neonatal epileptic encephalopathy to benign generalized epilepsy and beyond

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels control neuronal excitability and their dysfunction has been linked to epileptogenesis but few individuals with neurological disorders related to variants altering HCN channels have been reported so far. In 2014, we described five individuals with epileptic encephalopathy due to de novo HCN1 variants. To delineate HCN1-related disorders and investigate genotype-phenotype correlations further, we assembled a cohort of 33 unpublished patients with novel pathogenic or likely pathogenic variants: 19 probands carrying 14 different de novo mutations and four families with dominantly inherited variants segregating with epilepsy in 14 individuals, but not penetrant in six additional individuals. Sporadic patients had epilepsy with median onset at age 7 months and in 36% the first seizure occurred during a febrile illness. Overall, considering familial and sporadic patients, the predominant phenotypes were mild, including genetic generalized epilepsies and genetic epilepsy with febrile seizures plus (GEFS+) spectrum. About 20% manifested neonatal/infantile onset otherwise unclassified epileptic encephalopathy. The study also included eight patients with variants of unknown significance: one adopted patient had two HCN1 variants, four probands had intellectual disability without seizures, and three individuals had missense variants inherited from an asymptomatic parent. Of the 18 novel pathogenic missense variants identified, 12 were associated with severe phenotypes and clustered within or close to transmembrane domains, while variants segregating with milder phenotypes were located outside transmembrane domains, in the intracellular N- and C-terminal parts of the channel. Five recurrent variants were associated with similar phenotypes. Using whole-cell patch-clamp, we showed that the impact of 12 selected variants ranged from complete loss-of-function to significant shifts in activation kinetics and/or voltage dependence. Functional analysis of three different substitutions altering Gly391 revealed that these variants had different consequences on channel biophysical properties. The Gly391Asp variant, associated with the most severe, neonatal phenotype, also had the most severe impact on channel function. Molecular dynamics simulation on channel structure showed that homotetramers were not conducting ions because the permeation path was blocked by cation(s) strongly complexed to the Asp residue, whereas heterotetramers showed an instantaneous current component possibly linked to deformation of the channel pore. In conclusion, our results considerably expand the clinical spectrum related to HCN1 variants to include common generalized epilepsy phenotypes and further illustrate how HCN1 has a pivotal function in brain development and control of neuronal excitability.

High bandwidth approaches in nanopore and ion channel recordings - A tutorial review

Transport processes through ion-channel proteins, protein pores, or solid-state nanopores are traditionally recorded with commercial patch-clamp amplifiers. The bandwidth of these systems is typically limited to 10 kHz by signal-to-noise-ratio (SNR) considerations associated with these measurement platforms. At high bandwidth, the input-referred current noise in these systems dominates, determined by the input-referred voltage noise of the transimpedance amplifier applied across the capacitance at the input of the amplifier. This capacitance arises from several sources: the parasitic capacitance of the amplifier itself; the capacitance of the lipid bilayer harboring the ion channel protein (or the membrane used to form the solid-state nanopore); and the capacitance from the interconnections between the electronics and the membrane. Here, we review state-of-the-art applications of high-bandwidth conductance recordings of both ion channels and solid-state nanopores. These approaches involve tightly integrating measurement electronics fabricated in complementary metal-oxide semiconductors (CMOS) technology with lipid bilayer or solid-state membranes. SNR improvements associated with this tight integration push the limits of measurement bandwidths, in some cases in excess of 10 MHz. Recent case studies demonstrate the utility of these approaches for DNA sequencing and ion-channel recordings. In the latter case, studies with extended bandwidth have shown the potential for providing new insights into structure-function relations of these ion-channel proteins as the temporal resolutions of functional recordings matches time scales achievable with state-of-the-art molecular dynamics simulations.

Characterizing hydration sites in protein-ligand complexes towards the design of novel ligands

Water is an essential part of protein binding sites and mediates interactions to ligands. Its displacement by ligand parts affects the free binding energy of resulting protein-ligand complexes. Therefore the characterization of solvation properties is important for design. Of particular interest is the propensity of localized water to be favorably displaced by a ligand. This review discusses two popular computational approaches addressing these questions, namely WaterMap based on statistical mechanics analysis of MD simulations and 3D RISM based on integral equation theory of liquids. The theoretical background and recent applications in structure-based design will be presented.

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.

How to resolve microsecond current fluctuations in single ion channels: the power of beta distributions

A main ingredient for the understanding of structure/function correlates of ion channels is the quantitative description of single-channel gating and conductance. However, a wealth of information provided from fast current fluctuations beyond the temporal resolution of the recording system is often ignored, even though it is close to the time window accessible to molecular dynamics simulations. This kind of current fluctuations provide a special technical challenge, because individual opening/closing or blocking/unblocking events cannot be resolved, and the resulting averaging over undetected events decreases the single-channel current. Here, I briefly summarize the history of fast-current fluctuation analysis and focus on the so-called "beta distributions." This tool exploits characteristics of current fluctuation-induced excess noise on the current amplitude histograms to reconstruct the true single-channel current and kinetic parameters. A guideline for the analysis and recent applications demonstrate that a construction of theoretical beta distributions by Markov Model simulations offers maximum flexibility as compared to analytical solutions.

3D RISM theory with fast reciprocal-space electrostatics

The calculation of electrostatic solute-solvent interactions in 3D RISM ("three-dimensional reference interaction site model") integral equation theory is recast in a form that allows for a computational treatment analogous to the "particle-mesh Ewald" formalism as used for molecular simulations. In addition, relations that connect 3D RISM correlation functions and interaction potentials with thermodynamic quantities such as the chemical potential and average solute-solvent interaction energy are reformulated in a way that calculations of expensive real-space electrostatic terms on the 3D grid are completely avoided. These methodical enhancements allow for both, a significant speedup particularly for large solute systems and a smoother convergence of predicted thermodynamic quantities with respect to box size, as illustrated for several benchmark systems.

Engineering a Ca⁺⁺⁺-sensitive (bio)sensor from the pore-module of a potassium channel

Signals recorded at the cell membrane are meaningful indicators of the physiological vs. pathological state of a cell and will become useful diagnostic elements in nanomedicine. In this project we present a coherent strategy for the design and fabrication of a bio-nano-sensor that monitors changes in intracellular cell calcium concentration and allows an easy read out by converting the calcium signal into an electrical current in the range of microampere that can be easily measured by conventional cell electrophysiology apparatus.

A virus-encoded potassium ion channel is a structural protein in the chlorovirus Paramecium bursaria chlorella virus 1 virion

Most chloroviruses encode small K(+) channels, which are functional in electrophysiological assays. The experimental finding that initial steps in viral infection exhibit the same sensitivity to channel inhibitors as the viral K(+) channels has led to the hypothesis that the channels are structural proteins located in the internal membrane of the virus particles. This hypothesis was questioned recently because proteomic studies failed to detect the channel protein in virions of the prototype chlorovirus Paramecium bursaria chlorella virus 1 (PBCV-1). Here, we used a mAb raised against the functional K(+) channel from chlorovirus MA-1D to search for the viral K(+) channel in the virus particle. The results showed that the antibody was specific and bound to the tetrameric channel on the extracellular side. The antibody reacted in a virus-specific manner with protein extracts from chloroviruses that encoded channels similar to that from MA-1D. There was no cross-reactivity with chloroviruses that encoded more diverse channels or with a chlorovirus that lacked a K(+) channel gene. Together with electron microscopic imaging, which revealed labelling of individual virus particles with the channel antibody, these results establish that the viral particles contain an active K(+) channel, presumably located in the lipid membrane that surrounds the DNA in the mature virions.

Viral potassium channels as a robust model system for studies of membrane-protein interaction

The viral channel KcvNTS belongs to the smallest K(+) channels known so far. A monomer of a functional homotetramer contains only 82 amino acids. As a consequence of the small size the protein is almost fully submerged into the membrane. This suggests that the channel is presumably sensitive to its lipid environment. Here we perform a comparative analysis for the function of the channel protein embedded in three different membrane environments. 1. Single-channel currents of KcvNTS were recorded with the patch clamp method on the plasma membrane of HEK293 cells. 2. They were also measured after reconstitution of recombinant channel protein into classical planar lipid bilayers and 3. into horizontal bilayers derived from giant unilamellar vesicles (GUVs). The recombinant channel protein was either expressed and purified from Pichia pastoris or from a cell-free expression system; for the latter a new approach with nanolipoprotein particles was used. The data show that single-channel activity can be recorded under all experimental conditions. The main functional features of the channel like a large single-channel conductance (80pS), high open-probability (>50%) and the approximate duration of open and closed dwell times are maintained in all experimental systems. An apparent difference between the approaches was only observed with respect to the unitary conductance, which was ca. 35% lower in HEK293 cells than in the other systems. The reason for this might be explained by the fact that the channel is tagged by GFP when expressed in HEK293 cells. Collectively the data demonstrate that the small viral channel exhibits a robust function in different experimental systems. This justifies an extrapolation of functional data from these systems to the potential performance of the channel in the virus/host interaction. This article is part of a Special Issue entitled: Viral Membrane Proteins-Channels for Cellular Networking.

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.

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.

Integral Equation Theory of Biomolecules and Electrolytes

The so-called three-dimensional version (3D-RISM) can be used to describe the interactions of solvent components (here we treat water and ions) with a chemical or biomolecular solute of arbitrary size and shape. Here we give an overview of the current status of such models, describing some aspects of “pure” electrolytes (water plus simple ions) and of ionophores, proteins and nucleic acids in the presence of water and salts. Here we focus primarily on interactions with water and dissolved salts; as a practical matter, the discussion is mostly limited to monovalent ions, since studies of divalent ions present many difficult problems that have not yet been addressed. This is not a comprehensive review, but covers a few recent examples that illustrate current issues.

Mechanism of function of viral channel proteins and implications for drug development

Viral channel-forming proteins comprise a class of viral proteins which, similar to their host companions, are made to alter electrochemical or substrate gradients across lipid membranes. These proteins are active during all stages of the cellular life cycle of viruses. An increasing number of proteins are identified as channel proteins, but the precise role in the viral life cycle is yet unknown for the majority of them. This review presents an overview about these proteins with an emphasis on those with available structural information. A concept is introduced which aligns the transmembrane domains of viral channel proteins with those of host channels and toxins to give insights into the mechanism of function of the viral proteins from potential sequence identities. A summary of to date investigations on drugs targeting these proteins is given and discussed in respect of their mode of action in vivo.

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

The small viral channel Kcv is a Kir-like K(+) channel of only 94 amino acids. With this simple structure, the tetramer of Kcv represents the pore module of all complex K(+) channels. To examine the structural contribution of the transmembrane domains (TMDs) to channel function, we performed Ala scanning mutagenesis of the two domains and tested the functionality of the mutants in a yeast complementation assay. The data reveal, in combination with computational models, that the upper halves of both TMDs, which face toward the external medium, are rather rigid, whereas the inner parts are more flexible. The rigidity of the outer TMD is conferred by a number of essential aromatic amino acids that face the membrane and probably anchor this domain in the bilayer. The inner TMD is intimately connected with the rigid part of the outer TMD via π···π interactions between a pair of aromatic amino acids. This structural principle is conserved within the viral K(+) channels and also present in Kir2.2, implying a general importance of this architecture for K(+) channel function.

Minimal art: or why small viral K(+) channels are good tools for understanding basic structure and function relations

Some algal viruses contain genes that encode proteins with the hallmarks of K(+) channels. One feature of these proteins is that they are less than 100 amino acids in size, which make them truly minimal for a K(+) channel protein. That is, they consist of only the pore module present in more complex K(+) channels. The combination of miniature size and the functional robustness of the viral K(+) channels make them ideal model systems for studying how K(+) channels work. Here we summarize recent structure/function correlates from these channels, which provide insight into functional properties such as gating, pharmacology and sorting in cells.

Fast and slow gating are inherent properties of the pore module of the K+ channel Kcv

Kcv from the chlorella virus PBCV-1 is a viral protein that forms a tetrameric, functional K⁺ channel in heterologous systems. Kcv can serve as a model system to study and manipulate basic properties of the K⁺ channel pore because its minimalistic structure (94 amino acids) produces basic features of ion channels, such as selectivity, gating, and sensitivity to blockers. We present a characterization of Kcv properties at the single-channel level. In symmetric 100 mM K⁺, single-channel conductance is 114 ± 11 pS. Two different voltage-dependent mechanisms are responsible for the gating of Kcv. “Fast” gating, analyzed by β distributions, is responsible for the negative slope conductance in the single-channel current–voltage curve at extreme potentials, like in MaxiK potassium channels, and can be explained by depletion-aggravated instability of the filter region. The presence of a “slow” gating is revealed by the very low (in the order of 1–4%) mean open probability that is voltage dependent and underlies the time-dependent component of the macroscopic current.

Estimating sufficient statistics in co-evolutionary analysis by mutual information

Mutual information (MI) is a standard measure in information theory to observe and quantify correlated signals and events in both, empirical data sets and theoretical models. In the field of computational biology the MI turned out to be particularly useful in studies on co-evolutionary signals of sites within biomolecules. A key issue in the applicability of the MI is, however, a correct reference system or null model to understand finite-size effects in the underlying, finite data set. Although some bioinformatics studies exist with rigorous results for theoretical, well-designed random distributions, data from real-world proteins was never used to quantify the effect of finite-size samples. The impact of real-world statistics is, however, most relevant for researchers in all fields concerned with detecting evolutionary signals within biological sequences. We present results on such effects in finite-sized biological data sets and point to future research directions. We are most of all concerned with bacterial, ribosomal proteins as a prototypical example in molecular evolution. We compare to previous published suggestions, give an empirical formula, and propose a protocol to guide future research projects.

Model generation of viral channel forming 2B protein bundles from polio and coxsackie viruses

2B is a 99 amino acid membrane protein encoded by enteroviruses such as polio and coxsackie viruses with two transmembrane domains. The protein is found to make membranes of infected cells permeable. Using a computational approach which positions the models and assesses stability by molecular dynamics (MD) simulations a putative tetrameric bundle model of 2B is generated. The bundles show a pore lining motif of three lysines followed by a serine. The bundle is discussed in terms of different possible orientations of the helices in the membrane and the consequences this has on the in vivo activity of 2B.

Using a five-state model for fitting amplitude histograms from MaxiK channels: beta-distributions reveal more than expected

Fast gating of ion channels with rate constants higher than the corner frequency of the recording set-up can be evaluated by fitting so-called beta distributions to measured amplitude histograms. Up to now, this was preferentially done for O-C Markov sub-models with one open and one closed state. Here, a fit of the amplitude histograms from MaxiK (BK) single-channel records was achieved with a five-state model with two open and three closed states including three open-close transitions with rate constants higher than the corner frequency (20 kHz) of the inevitable low-pass filter of the recording system. The numerical values of the rate constants of these transitions enabled a nearly one-to-one relationship between typical regions of the histograms and the reactions in the Markov model. These characteristic features are the width of the peak at the apparent single-channel current, the side slopes at the open and at the closed peak, and the depth of the valley between the two peaks. However, the simplex routine alone was incapable of finding the solution but could do so if guided by hand along a suggested strategy.

Interference of shot noise of open-channel current with analysis of fast gating: patchers do not (Yet) have to care

Microsecond gating of ion channels can be evaluated by fitting beta distributions to amplitude histograms of measured time series. The shape of these histograms is determined not only by the rate constants of the gating process (in relation to the filter frequency) but also by baseline noise and shot noise, resulting from the stochastic nature of ion flow. Under normal temporal resolution, the small shot noise can be ignored. This simplification may no longer be legitimate when rate constants reach the range above 1 mus(-1). Here, the influence of shot noise is studied by means of simulated time series for several values of single-channel current of the fully open state and baseline noise. Under realistic optimal conditions (16 pA current, 1 pA noise, 50 kHz bandwidth), ignoring the shot noise leads to an underestimation of the rate constants above 1 mus(-1) by a factor of about 2.5. However, in that range, the scatter of the evaluated rate constants is at least of the same magnitude, obscuring the systematic error. The incorporation of shot noise into the analysis will become more important when amplifiers with significantly reduced noise become available.

Mimicking biopolymers on a molecular scale: nano(bio)technology based on engineered proteins

Proteins are ubiquitous biopolymers that adopt distinct three-dimensional structures and fulfil a multitude of elementary functions in organisms. Recent systematic studies in molecular biology and biotechnology have improved the understanding of basic functional and architectural principles of proteins, making them attractive candidates as concept generators for technological development in material science, particularly in biomedicine and nano(bio)technology. This paper highlights the potential of molecular biomimetics in mimicking high-performance proteins and provides concepts for applications in four case studies, i.e. spider silk, antifreeze proteins, blue mussel adhesive proteins and viral ion channels.

Salt bridges in the miniature viral channel Kcv are important for function

The viral potassium channel Kcv comprises only 94 amino acids, which represent the pore module of more complex K(+) channels. As for Kir-type channels, Kcv also has a short N-terminal helix exposed to the cytoplasm, upstream of the first transmembrane domain. Here we show that this helix is relevant for Kcv function. The presence of charged amino acids, which form dynamic inter- and intra-subunit salt bridges is crucial. Electrophysiological measurements, yeast rescue experiments and molecular dynamics simulations show that mutants in which the critical salt bridge formation is impaired have no or reduced channel activity. We conclude that these salt bridges destabilise the complexation of K(+) ions by negative charges on the inner transmembrane domain at the entrance into the cavity. This feature facilitates a continuous and coordinated transfer of ions between the cavity and the cytoplasm for channels without the canonical bundle crossing.

Model development for the viral Kcv potassium channel

A computational model for the open state of the short viral Kcv potassium channel was created and tested based on homology modeling and extensive molecular-dynamics simulation in a membrane environment. Particular attention was paid to the structure of the highly flexible N-terminal region and to the protonation state of membrane-exposed lysine residues. Data from various experimental sources, NMR spectroscopy, and electrophysiology, as well as results from three-dimensional reference interaction site model integral equation theory were taken into account to select the most reasonable model among possible variants. The final model exhibits spontaneous ion transitions across the complete pore, with and without application of an external field. The nonequilibrium transport events could be induced reproducibly without abnormally large driving potential and without the need to place ions artificially at certain key positions along the transition path. The transport mechanism through the filter region corresponds to the classic view of single-file motion, which in our case is coupled to frequent exchange of ions between the innermost filter position and the cavity.

The Phycodnaviridae: the story of how tiny giants rule the world

The family Phycodnaviridae encompasses a diverse and rapidly expanding collection of large icosahedral, dsDNA viruses that infect algae. These lytic and lysogenic viruses have genomes ranging from 160 to 560 kb. The family consists of six genera based initially on host range and supported by sequence comparisons. The family is monophyletic with branches for each genus, but the phycodnaviruses have evolutionary roots that connect them with several other families of large DNA viruses, referred to as the nucleocytoplasmic large DNA viruses (NCLDV). The phycodnaviruses have diverse genome structures, some with large regions of noncoding sequence and others with regions of ssDNA. The genomes of members in three genera in the Phycodnaviridae have been sequenced. The genome analyses have revealed more than 1000 unique genes, with only 14 homologous genes in common among the three genera of phycodnaviruses sequenced to date. Thus, their gene diversity far exceeds the number of so-called core genes. Not much is known about the replication of these viruses, but the consequences of these infections on phytoplankton have global affects, including influencing geochemical cycling and weather patterns.

Saturation and microsecond gating of current indicate depletion-induced instability of the MaxiK selectivity filter

Patch clamp experiments on single MaxiK channels expressed in HEK293 cells were performed with a high temporal resolution (50-kHz filter) in symmetrical solutions with 50, 150, or 400 mM KCl and 2.5 mM CaCl(2) and 2.5 mM MgCl(2). At membrane potentials >+100 mV, the single-channel current showed a negative slope resistance, concomitantly with a flickery block, which was not influenced by Ca(2+) or Mg(2+). The analysis of the amplitude histograms by beta distributions revealed that current in this voltage range was reduced by two effects: rate limitation at the cytosolic side of the pore and gating with rate constants 10-20-fold higher than the cutoff frequency of the filter (i.e., dwell times in the microsecond range). The data were analyzed in terms of a model that postulates a coupling between both effects; if the voltage over the selectivity filter withdraws ions from the cavity at a higher rate than that of refilling from the cytosol, the selectivity filter becomes instable because of ion depletion, and current is interrupted by the resulting flickering. The fit of the IV curves revealed a characteristic voltage of 35 mV. In contrast, the voltage dependence of the gating factor R, i.e., the ratio between true and apparent single-channel current, could be fitted by exponentials with a characteristic voltage of 60 mV, suggesting that only part of the transmembrane potential is felt by the flux through the selectivity filter.