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

Permeant cations modulate pore dynamics and gating of TRPV1 ion channels

The transient receptor vanilloid 1 (TRPV1) is a non-selective ion channel, which is activated by several chemical ligands and heat. We have previously shown that activation of TRPV1 by different ligands results in single-channel openings with different conductance, suggesting that the selectivity filter is highly dynamic. TRPV1 is weakly voltage dependent; here, we sought to explore whether the permeation of different monovalent ions could influence the voltage dependence of this ion channel. By using single-channel recordings, we show that TRPV1 channels undergo rapid transitions to closed states that are directly connected to the open state, which may result from structural fluctuations of their selectivity filter. Moreover, we demonstrate that the rates of these transitions are influenced by the permeant ion, suggesting that ion permeation regulates the voltage dependence of these channels. Our data could be the basis for more detailed MD simulations exploring the permeation mechanism and how the occupancy of different ions alters the three-dimensional structure of the pore of TRPV1 channels.

2D-dwell-time analysis with simulations of ion-channel gating using high-performance computing

Single-channel patch-clamp recordings allow observing the action of a single protein complex in real time and hence the deduction of the underlying conformational changes in the ion-channel protein. Commonly, recordings are modeled using hidden Markov chains, connecting open and closed states in the experimental data with protein conformations. The rates between states denote transition probabilities that could be modified by membrane voltage or ligand binding. Modeling algorithms have to deal with limited recording bandwidth and a very noisy background. It was previously shown that the fit of two-dimensional (2D)-dwell-time histograms with simulations is very robust in that regard. Errors introduced by the low-pass filter or noise cancel out to a certain degree when comparing experimental and simulated data. In addition, the topology of models (that is, the chain of open and closed states) could be inferred from 2D-histograms. However, the 2D-fit was never applied to its full potential. A major reason may be the extremely time-consuming and often unreliable fitting process, due to the stochastic variability in the simulations. We have now solved these issues by introducing a message-passing interface (MPI) allowing massive parallel computing on a high-performance computing (HPC) cluster and obtaining ensemble solutions. With ensembles, we have demonstrated how important ranked solutions are for difficult tasks related to a noisy background, fast gating events beyond the corner frequency of the low-pass filter, and topology estimation of the underlying Markov model. Finally, we have shown that, by combining the objective function of the 2D-fit with the deviation of the current amplitude distributions, automatic determination of the current level of the conducting state is possible, even with an apparent current reduction due to low-pass filtering. Making use of an HPC cluster, the power of 2D-dwell-time analysis can be used to its fullest with minor input of the experimenter.

Silicon Nitride-Based Micro-Apertures Coated with Parylene for the Investigation of Pore Proteins Fused in Free-Standing Lipid Bilayers

In this work, we present a microsystem setup for performing sensitive biological membrane translocation measurements. Thin free-standing synthetic bilayer lipid membranes (BLM) were constructed in microfabricated silicon nitride apertures (<100 µm in diameter), conformal coated with Parylene (Parylene-C or Parylene-AF4). Within these BLMs, electrophysiological measurements were conducted to monitor the behavior of different pore proteins. Two approaches to integrate pore-forming proteins into the membrane were applied: direct reconstitution and reconstitution via outer membrane vesicles (OMVs) released from Gram-negative bacteria. The advantage of utilizing OMVs is that the pore proteins remain in their native lipid and lipopolysaccharide (LPS) environment, representing a more natural state compared to the usage of fused purified pore proteins. Multiple aperture chips can be easily assembled in the 3d-printed holder to conduct parallel membrane transport investigations. Moreover, well defined microfabricated apertures are achievable with very high reproducibility. The presented microsystem allows the investigation of fast gating events (down to 1 ms), pore blocking by an antibiotic, and gating events of small pores (amplitude of approx. 3 pA).

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 Kcv_NTS. 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, Kcv_NTS 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.

Effects of ionic strength on gating and permeation of TREK-2 K2P channels

In addition to the classical voltage-dependent behavior mediated by the voltage-sensing-domains (VSD) of ion channels, a growing number of voltage-dependent gating behaviors are being described in channels that lack canonical VSDs. A common thread in their mechanism of action is the contribution of the permeating ion to this voltage sensing process. The polymodal K2P K⁺ channel, TREK2 responds to membrane voltage through a gating process mediated by the interaction of K⁺ with its selectivity filter. Recently, we found that this action can be modulated by small molecule agonists (e.g. BL1249) which appear to have an electrostatic influence on K⁺ binding within the inner cavity and produce an increase in the single-channel conductance of TREK-2 channels. Here, we directly probed this K⁺-dependent gating process by recording both macroscopic and single-channel currents of TREK-2 in the presence of high concentrations of internal K⁺. Surprisingly we found TREK-2 is inhibited by high internal K⁺ concentrations and that this is mediated by the concomitant increase in ionic-strength. However, we were still able to determine that the increase in single channel conductance in the presence of BL1249 was blunted in high ionic-strength, whilst its activatory effect (on channel open probability) persisted. These effects are consistent with an electrostatic mechanism of action of negatively charged activators such as BL1249 on permeation, but also suggest that their influence on channel gating is complex.

Analysis of patchclamp recordings: model-free multiscale methods and software

Analysis of patchclamp recordings is often a challenging issue. We give practical guidance how such recordings can be analyzed using the model-free multiscale idealization methodology JSMURF, JULES, and HILDE. We provide an operational manual how to use the accompanying software available as an R-package and as a graphical user interface. This includes selection of the right approach and tuning of parameters. We also discuss advantages and disadvantages of model-free approaches in comparison to hidden Markov model approaches and explain how they complement each other.

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.

Heterogeneous Idealization of Ion Channel Recordings - Open Channel Noise

We propose a new model-free segmentation method for idealizing ion channel recordings. This method is designed to deal with heterogeneity of measurement errors. This in particular applies to open channel noise which, in general, is particularly difficult to cope with for model-free approaches. Our methodology is able to deal with lowpass filtered data which provides a further computational challenge. To this end we propose a multiresolution testing approach, combined with local deconvolution to resolve the lowpass filter. Simulations and statistical theory confirm that the proposed idealization recovers the underlying signal very accurately at presence of heterogeneous noise, even when events are shorter than the filter length. The method is compared to existing approaches in computer experiments and on real data. We find that it is the only one which allows to identify openings of the PorB porine at two different temporal scales. An implementation is available as an R package.

Comparing Current Noise in Biological and Solid-State Nanopores

Nanopores bear great potential as single-molecule tools for bioanalytical sensing and sequencing, due to their exceptional sensing capabilities, high-throughput, and low cost. The detection principle relies on detecting small differences in the ionic current as biomolecules traverse the nanopore. A major bottleneck for the further progress of this technology is the noise that is present in the ionic current recordings, because it limits the signal-to-noise ratio (SNR) and thereby the effective time resolution of the experiment. Here, we review the main types of noise at low and high frequencies and discuss the underlying physics. Moreover, we compare biological and solid-state nanopores in terms of the SNR, the important figure of merit, by measuring translocations of a short ssDNA through a selected set of nanopores under typical experimental conditions. We find that SiNx solid-state nanopores provide the highest SNR, due to the large currents at which they can be operated and the relatively low noise at high frequencies. However, the real game-changer for many applications is a controlled slowdown of the translocation speed, which for MspA was shown to increase the SNR > 160-fold. Finally, we discuss practical approaches for lowering the noise for optimal experimental performance and further development of the nanopore technology.

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.

Site-specific ion occupation in the selectivity filter causes voltage-dependent gating in a viral K+ channel

Many potassium channels show voltage-dependent gating without a dedicated voltage sensor domain. This is not fully understood yet, but often explained by voltage-induced changes of ion occupation in the five distinct K+ binding sites in the selectivity filter. To better understand this mechanism of filter gating we measured the single-channel current and the rate constant of sub-millisecond channel closure of the viral K+ channel KcvNTS for a wide range of voltages and symmetric and asymmetric K+ concentrations in planar lipid membranes. A model-based analysis employed a global fit of all experimental data, i.e., using a common set of parameters for current and channel closure under all conditions. Three different established models of ion permeation and various relationships between ion occupation and gating were tested. Only one of the models described the data adequately. It revealed that the most extracellular binding site (S0) in the selectivity filter functions as the voltage sensor for the rate constant of channel closure. The ion occupation outside of S0 modulates its dependence on K+ concentration. The analysis uncovers an important role of changes in protein flexibility in mediating the effect from the sensor to the gate.

Fully Automatic Multiresolution Idealization for Filtered Ion Channel Recordings: Flickering Event Detection

We propose a new model-free segmentation method, JULES, which combines recent statistical multiresolution techniques with local deconvolution for idealization of ion channel recordings. The multiresolution criterion takes into account scales down to the sampling rate enabling the detection of flickering events, i.e., events on small temporal scales, even below the filter frequency. For such small scales the deconvolution step allows for a precise determination of dwell times and, in particular, of amplitude levels, a task which is not possible with common thresholding methods. This is confirmed theoretically and in a comprehensive simulation study. In addition, JULES can be applied as a preprocessing method for a refined hidden Markov analysis. Our new methodology allows us to show that gramicidin A flickering events have the same amplitude as the slow gating events. JULES is available as an R function jules in the package clampSeg.

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 Ca²⁺-activated intracellular Ca²⁺-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 Ca²⁺ 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.

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.

Sensing Single Molecule Penetration into Nanopores: Pushing the Time Resolution to the Diffusion Limit

To quantify small molecule penetration into and eventually permeation through nanopores, we applied an improved excess-noise analysis of the ion current fluctuation caused by entering molecules. The kinetic parameters of substrate entry and exit are derived from a two-state Markov model, analyzing the substrate concentration dependence of the average ion current and its variance. Including filter corrections allows one to detect the transition rates beyond the cutoff frequency, f_c, of the instrumental ion-current filter. As an application of the method, we performed an analysis of the single-channel ion current of Meropenem, an antibiotic of the carbapenem family, interacting with OmpF, the major general outer membrane channel of Escherichia coli bacteria. At 40 °C we detected the residence time of Meropenem inside OmpF of about 500 ns-more than 2 orders of magnitude smaller than f_c^-1 and close to the diffusion limit of few hundred nanoseconds. We also have established theoretical limit conditions under which the substrate-induced channel blockages can be detected and suggest that submicrosecond-scale gating kinetic parameters are accessible with existing experimental equipment.