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.

Single-channel recordings of RyR1 at microsecond resolution in CMOS-suspended membranes.

Single-channel recordings are widely used to explore functional properties of ion channels. Typically, such recordings are performed at bandwidths of less than 10 kHz because of signal-to-noise considerations, limiting the temporal resolution available for studying fast gating dynamics to greater than 100 µs. Here we present experimental methods that directly integrate suspended lipid bilayers with high-bandwidth, low-noise transimpedance amplifiers based on complementary metal-oxide-semiconductor (CMOS) integrated circuits (IC) technology to achieve bandwidths in excess of 500 kHz and microsecond temporal resolution. We use this CMOS-integrated bilayer system to study the type 1 ryanodine receptor (RyR1), a Ca2+-activated intracellular Ca2+-release channel located on the sarcoplasmic reticulum. We are able to distinguish multiple closed states not evident with lower bandwidth recordings, suggesting the presence of an additional Ca2+ binding site, distinct from the site responsible for activation. An extended beta distribution analysis of our high-bandwidth data can be used to infer closed state flicker events as fast as 35 ns. These events are in the range of single-file ion translocations.

Identification of Intrahelical Bifurcated H-Bonds as a New Type of Gate in K+ Channels.

Gating of ion channels is based on structural transitions between open and closed states. To uncover the chemical basis of individual gates, we performed a comparative experimental and computational analysis between two K+ channels, KcvS and KcvNTS. These small viral encoded K+ channel proteins, with a monomer size of only 82 amino acids, resemble the pore module of all complex K+ channels in terms of structure and function. Even though both proteins share about 90% amino acid sequence identity, they exhibit different open probabilities with ca. 90% in KcvNTS and 40% in KcvS. Single channel analysis, mutational studies and molecular dynamics simulations show that the difference in open probability is caused by one long closed state in KcvS. This state is structurally created in the tetrameric channel by a transient, Ser mediated, intrahelical hydrogen bond. The resulting kink in the inner transmembrane domain swings the aromatic rings from downstream Phes in the cavity of the channel, which blocks ion flux. The frequent occurrence of Ser or Thr based helical kinks in membrane proteins suggests that a similar mechanism could also occur in the gating of other ion channels.

Cyclic dinucleotides bind the C-linker of HCN4 to control channel cAMP responsiveness.

cAMP mediates autonomic regulation of heart rate by means of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, which underlie the pacemaker current If. cAMP binding to the C-terminal cyclic nucleotide binding domain enhances HCN open probability through a conformational change that reaches the pore via the C-linker. Using structural and functional analysis, we identified a binding pocket in the C-linker of HCN4. Cyclic dinucleotides, an emerging class of second messengers in mammals, bind the C-linker pocket (CLP) and antagonize cAMP regulation of the channel. Accordingly, cyclic dinucleotides prevent cAMP regulation of If in sinoatrial node myocytes, reducing heart rate by 30%. Occupancy of the CLP hence constitutes an efficient mechanism to hinder ß-adrenergic stimulation on If. Our results highlight the regulative role of the C-linker and identify a potential drug target in HCN4. Furthermore, these data extend the signaling scope of cyclic dinucleotides in mammals beyond their first reported role in innate immune system.

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.

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.

Role of Ion Distribution and Energy Barriers for Concerted Motion of Subunits in Selectivity Filter Gating of a K+ Channel.

Most potassium channels have two main gate locations, hosting an inner gate at the cytosolic entrance and a filter gate in the selectivity filter; the function of these gates is in many channels coupled. To obtain exclusive insights into the molecular mechanisms that determine opening and closing of the filter gate, we use a combination of single-channel recordings and gating analysis in the minimal viral channel KcvNTS. This channel has no inner gate, and its fast closing at negative voltages can therefore be entirely assigned to the filter gate. We find that mutations of S42 in the pore helix severely slow down closing of this filter gate, an effect which is not correlated with hydrogen bond formation by the amino acid at this position. Hence, different from KcsA, which contains the critical E71 in the equivalent position forming a salt bridge, the coupling between selectivity filter and surrounding structures for filter gating must in KcvNTS rely on different modes of interaction. Quantitative analysis of concatemers carrying different numbers of S42T mutations reveals that each subunit contributes the same amount of âˆ¼ 0.4 kcal/mol to the energy barrier for filter closure indicating a concerted action of the subunits. Since the mutations have neither an influence on the unitary current nor on the voltage dependency of the gate, the data stress that the high subunit cooperativity is mediated through conformational changes rather than through changes in the ion occupation in the selectivity filter.

Asymmetric Interplay Between K+ and Blocker and Atomistic Parameters From Physiological Experiments Quantify K+ Channel Blocker Release.

Modulating the activity of ion channels by blockers yields information on both the mode of drug action and on the biophysics of ion transport. Here we investigate the interplay between ions in the selectivity filter (SF) of K+ channels and the release kinetics of the blocker tetrapropylammonium in the model channel KcvNTS. A quantitative expression calculates blocker release rate constants directly from voltage-dependent ion occupation probabilities in the SF. The latter are obtained by a kinetic model of single-channel currents recorded in the absence of the blocker. The resulting model contains only two adjustable parameters of ion-blocker interaction and holds for both symmetric and asymmetric ionic conditions. This data-derived model is corroborated by 3D reference interaction site model (3D RISM) calculations on several model systems, which show that the K+ occupation probability is unaffected by the blocker, a direct consequence of the strength of the ion-carbonyl attraction in the SF, independent of the specific protein background. Hence, KcvNTS channel blocker release kinetics can be reduced to a small number of system-specific parameters. The pore-independent asymmetric interplay between K+ and blocker ions potentially allows for generalizing these results to similar potassium channels.

Distinct lipid bilayer compositions have general and protein-specific effects on K+ channel function.

It has become increasingly apparent that the lipid composition of cell membranes affects the function of transmembrane proteins such as ion channels. Here, we leverage the structural and functional diversity of small viral K+ channels to systematically examine the impact of bilayer composition on the pore module of single K+ channels. In vitro-synthesized channels were reconstituted into phosphatidylcholine bilayers ± cholesterol or anionic phospholipids (aPLs). Single-channel recordings revealed that a saturating concentration of 30% cholesterol had only minor and protein-specific effects on unitary conductance and gating. This indicates that channels have effective strategies for avoiding structural impacts of hydrophobic mismatches between proteins and the surrounding bilayer. In all seven channels tested, aPLs augmented the unitary conductance, suggesting that this is a general effect of negatively charged phospholipids on channel function. For one channel, we determined an effective half-maximal concentration of 15% phosphatidylserine, a value within the physiological range of aPL concentrations. The different sensitivity of two channel proteins to aPLs could be explained by the presence/absence of cationic amino acids at the interface between the lipid headgroups and the transmembrane domains. aPLs also affected gating in some channels, indicating that conductance and gating are uncoupled phenomena and that the impact of aPLs on gating is protein specific. In two channels, the latter can be explained by the altered orientation of the pore-lining transmembrane helix that prevents flipping of a phenylalanine side chain into the ion permeation pathway for long channel closings. Experiments with asymmetrical bilayers showed that this effect is leaflet specific and most effective in the inner leaflet, in which aPLs are normally present in plasma membranes. The data underscore a general positive effect of aPLs on the conductance of K+ channels and a potential interaction of their negative headgroup with cationic amino acids in their vicinity.

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.

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.

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.

A small viral potassium ion channel with an inherent inward rectification.

Some algal viruses have coding sequences for proteins with structural and functional characteristics of pore modules of complex K+ channels. Here we exploit the structural diversity among these channel orthologs to discover new basic principles of structure/function correlates in K+ channels. The analysis of three similar K+ channels with ≤ 86 amino acids (AA) shows that one channel (Kmpv1) generates an ohmic conductance in HEK293 cells while the other two (KmpvSP1, KmpvPL1) exhibit typical features of canonical Kir channels. Like Kir channels, the rectification of the viral channels is a function of the K+ driving force. Reconstitution of KmpvSP1 and KmpvPL1 in planar lipid bilayers showed rapid channel fluctuations only at voltages negative of the K+ reversal voltage. This rectification was maintained in KCl buffer with 1 mM EDTA, which excludes blocking cations as the source of rectification. This means that rectification of the viral channels must be an inherent property of the channel. The structural basis for rectification was investigated by a chimera between rectifying and non-rectifying channels as well as point mutations making the rectifier similar to the ohmic conducting channel. The results of these experiments exclude the pore with pore helix and selectivity filter as playing a role in rectification. The insensitivity of the rectifier to point mutations suggests that tertiary or quaternary structural interactions between the transmembrane domains are responsible for this type of gating.

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.

Reconstitution and functional characterization of ion channels from nanodiscs in lipid bilayers.

Recent studies have shown that membrane proteins can be efficiently synthesized in vitro before spontaneously inserting into soluble nanoscale lipid bilayers called nanodiscs (NDs). In this paper, we present experimental details that allow a combination of in vitro translation of ion channels into commercially available NDs followed by their direct reconstitution from these nanobilayers into standard bilayer setups for electrophysiological characterization. We present data showing that two model K+ channels, Kcv and KcsA, as well as a recently discovered dual-topology F- channel, Fluc, can be reliably reconstituted from different types of NDs into bilayers without contamination from the in vitro translation cocktail. The functional properties of Kcv and KcsA were characterized electrophysiologically and exhibited sensitivity to the lipid composition of the target DPhPC bilayer, suggesting that the channel proteins were fully exposed to the target membrane and were no longer surrounded by the lipid/protein scaffold. The single-channel properties of the three tested channels are compatible with studies from recordings of the same proteins in other expression systems. Altogether, the data show that synthesis of ion channels into NDs and their subsequent reconstitution into conventional bilayers provide a fast and reliable method for functional analysis of ion channels.

Extended beta distributions open the access to fast gating in bilayer experiments-assigning the voltage-dependent gating to the selectivity filter.

Lipid bilayers provide many benefits for ion channel recordings, such as control of membrane composition and transport molecules. However, they suffer from high membrane capacitance limiting the bandwidth and impeding analysis of fast gating. This can be overcome by fitting the deviations of the open-channel noise from the baseline noise by extended beta distributions. We demonstrate this analysis step-by-step by applying it to the example of viral K+  channels (Kcv), from the choice of the gating model through the fitting process, validation of the results, and what kinds of results can be obtained. These voltage sensor-less channels show profoundly voltage-dependent gating with dwell times in the closed state of about 50 µs. Mutations assign it to the selectivity filter.

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.

A simple recipe for setting up the flux equations of cyclic and linear reaction schemes of ion transport with a high number of states: The arrow scheme.

The calculation of flux equations or current-voltage relationships in reaction kinetic models with a high number of states can be very cumbersome. Here, a recipe based on an arrow scheme is presented, which yields a straightforward access to the minimum form of the flux equations and the occupation probability of the involved states in cyclic and linear reaction schemes. This is extremely simple for cyclic schemes without branches. If branches are involved, the effort of setting up the equations is a little bit higher. However, also here a straightforward recipe making use of so-called reserve factors is provided for implementing the branches into the cyclic scheme, thus enabling also a simple treatment of such cases.

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.