Bayesian inference of kinetic schemes for ion channels by Kalman filtering

From Average Transient Transporter Currents to Microscopic Mechanism─A Bayesian Analysis

Electrophysiology studies of secondary active transporters have revealed quantitative mechanistic insights over many decades of research. However, the emergence of new experimental and analytical approaches calls for investigation of the capabilities and limitations of the newer methods. We examine the ability of solid-supported membrane electrophysiology (SSME) to characterize discrete-state kinetic models with >10 rate constants. We use a Bayesian framework applied to synthetic data for three tasks: to quantify and check (i) the precision of parameter estimates under different assumptions, (ii) the ability of computation to guide the selection of experimental conditions, and (iii) the ability of our approach to distinguish among mechanisms based on SSME data. When the general mechanism, i.e., event order, is known in advance, we show that a subset of kinetic parameters can be "practically identified" within ∼1 order of magnitude, based on SSME current traces that visually appear to exhibit simple exponential behavior. This remains true even when accounting for systematic measurement bias and realistic uncertainties in experimental inputs (concentrations) are incorporated into the analysis. When experimental conditions are optimized or different experiments are combined, the number of practically identifiable parameters can be increased substantially. Some parameters remain intrinsically difficult to estimate through SSME data alone, suggesting that additional experiments are required to fully characterize parameters. We also demonstrate the ability to perform model selection and determine the order of events when that is not known in advance, comparing Bayesian and maximum-likelihood approaches. Finally, our studies elucidate good practices for the increasingly popular but subtly challenging Bayesian calculations for structural and systems biology.

Fast functional mapping of ligand-gated ion channels

Ligand-gated ion channels are formed by three to five subunits that control the opening of the pore in a cooperative fashion. We developed a microfluidic chip-based technique for studying ion currents and fluorescence signals in either excised membrane patches or whole cells to measure activation and deactivation kinetics of the channels as well as ligand binding and unbinding when using confocal patch-clamp fluorometry. We show how this approach produces in a few seconds either unidirectional concentration-activation relationships at or near equilibrium and, moreover, respective time courses of activation and deactivation for a large number of freely designed steps of the ligand concentration. The short measuring period strongly minimizes the contribution of disturbing superimposing effects such as run-down phenomena and desensitization effects. To validate gating mechanisms, complex kinetic schemes are quantified without the requirement to have data at equilibrium. The new method has potential for functionally analyzing any ligand-gated ion channel and, beyond, also for other receptors.