Methods of Measuring Mitochondrial Potassium Channels: A Critical Assessment.

In this paper, the techniques used to study the function of mitochondrial potassium channels are critically reviewed. The majority of these techniques have been known for many years as a result of research on plasma membrane ion channels. Hence, in this review, we focus on the critical evaluation of techniques used in the studies of mitochondrial potassium channels, describing their advantages and limitations. Functional analysis of mitochondrial potassium channels in comparison to that of plasmalemmal channels presents additional experimental challenges. The reliability of functional studies of mitochondrial potassium channels is often affected by the need to isolate mitochondria and by functional properties of mitochondria such as respiration, metabolic activity, swelling capacity, or high electrical potential. Three types of techniques are critically evaluated: electrophysiological techniques, potassium flux measurements, and biochemical techniques related to potassium flux measurements. Finally, new possible approaches to the study of the function of mitochondrial potassium channels are presented. We hope that this review will assist researchers in selecting reliable methods for studying, e.g., the effects of drugs on mitochondrial potassium channel function. Additionally, this review should aid in the critical evaluation of the results reported in various articles on mitochondrial potassium channels.

Critical Assessment of Common Force Fields for Molecular Dynamics Simulations of Potassium Channels.

For the last two decades, the KcsA K+ channel has served as a case study to understand how potassium channels operate at the atomic scale, and molecular dynamics simulations have contributed significantly to the current knowledge of the atomic mechanisms of conduction, inactivation, and gating in this family of membrane proteins. Currently, microsecond trajectories are becoming the new standard in the field, and consequently, it is critical to assess and compare the performance of the classical force fields ordinarily used in simulations of biological systems as well as to quantitatively assess the agreement with experimental data for trajectories of this order of magnitude. To that extent, we performed classical molecular dynamics simulations with CHARMM36 and AMBER-ff14sb force fields using atomic models based on the experimental structure of the KcsA channel in the open/conductive state, at conditions of ionic concentrations and membrane potentials resembling the ones adopted in experiments. In simulations using the CHARMM force field, the experimental open/conductive structure of the channel exhibited high conformational plasticity and fast collapse toward an occluded state. In contrast, in an identical set of simulations using the AMBER force field, no major deviations from the experimental structure were recorded. Force field development is a complex process in which many approximations are typically required and adopted. The results presented here provide additional motivation to further improve the existing models and, crucially, alert practitioners about limitations.

Revealing kinetics and state-dependent binding properties of IKur-targeting drugs that maximize atrial fibrillation selectivity.

The KV1.5 potassium channel, which underlies the ultra-rapid delayed-rectifier current (IKur) and is predominantly expressed in atria vs. ventricles, has emerged as a promising target to treat atrial fibrillation (AF). However, while numerous KV1.5-selective compounds have been screened, characterized, and tested in various animal models of AF, evidence of antiarrhythmic efficacy in humans is still lacking. Moreover, current guidelines for pre-clinical assessment of candidate drugs heavily rely on steady-state concentration-response curves or IC50 values, which can overlook adverse cardiotoxic effects. We sought to investigate the effects of kinetics and state-dependent binding of IKur-targeting drugs on atrial electrophysiology in silico and reveal the ideal properties of IKur blockers that maximize anti-AF efficacy and minimize pro-arrhythmic risk. To this aim, we developed a new Markov model of IKur that describes KV1.5 gating based on experimental voltage-clamp data in atrial myocytes from patient right-atrial samples in normal sinus rhythm. We extended the IKur formulation to account for state-specificity and kinetics of KV1.5-drug interactions and incorporated it into our human atrial cell model. We simulated 1- and 3-Hz pacing protocols in drug-free conditions and with a [drug] equal to the IC50 value. The effects of binding and unbinding kinetics were determined by examining permutations of the forward (kon) and reverse (koff) binding rates to the closed, open, and inactivated states of the KV1.5 channel. We identified a subset of ideal drugs exhibiting anti-AF electrophysiological parameter changes at fast pacing rates (effective refractory period prolongation), while having little effect on normal sinus rhythm (limited action potential prolongation). Our results highlight that accurately accounting for channel interactions with drugs, including kinetics and state-dependent binding, is critical for developing safer and more effective pharmacological anti-AF options.

Assessment of Basilar Artery Reactivity in Stroke and Subarachnoid Hemorrhage Using Wire Myograph.

Blood flow regulation of normal cerebral arteries is a critical and important factor to supply the brain tissue with nutrients and oxygen. Stroke insult results in a disruption or reduction in cerebral arteries' blood flow with subsequent brain tissue damage. Hemorrhagic stroke is one type of stroke and accounts for about 13 % of all of stroke insults. In this type of stroke, the cerebral artery breaks open and causes bleeding in or surrounding the brain. Subsequently, this bleeding causes blood vessels to constrict in a process called vasospasm, in which the vessels narrow and impede the blood flow to brain tissue. Hemorrhagic stroke is the major cause of prolonged constriction of cerebral arteries. This leads to partial brain damage and sometimes death in patients with aneurysmal subarachnoid hemorrhage. Among the key delicate techniques to assess small blood vessel functionality is the wire myograph, which can be utilized in several cerebral injury models including stroke. The wire myograph is a device that provides information about the reactivity, stiffness, and elasticity of small blood vessels under isometric conditions. In this book chapter, we describe the techniques involved in wire myography assessment and the different measures and parameters recorded; we describe the utility of this technique in evaluating the effects of subarachnoid hemorrhage on basilar artery sensitivity to different agonists.

[Comparative assessment of different approaches for obtaining terminally differentiated muscle cells].

Relevant cell model is essential to study pathogenesis of muscle disorders. However, in the field of muscle research there is no ultimate cell line considered as a standard for studying muscular and neuromuscular diseases. Standard cell line claimed to be well differentiated in muscle lineage, be morphological and physiological similar to mature muscle cells and be easily genetically modified. Therefore, the goal of our study was to pick up available and fruitful cell model of muscle differentiation, that could be further applied for examination of muscular disorder pathogenesis in vitro. We characterized human mesenchymal stem cells (MSC), mature murine muscle fibers and primary murine satellite cells. It has been shown that MSC have very small capacity to myogenic differentiation; moreover, they were able to differentiate only in presence of C2C12 cells. Lentiviral transduction exhibited rather high toxic effect on primary myofibers, and positively transduced cells were not able to response to electrical stimulation, i. e. were functionally inactive. Satellite cells turned out to be the most fruitful cell model since they were easily transduced via lentiviruses and rapidly formed myotubes in differentiation media. Functional analysis of obtained myotubes has confirmed their ability to react to electrical and chemical stimulations; besides, potassium and calcium channels availability has been also demonstrated via patch-clump technique. Taken together, these results imply that satellite cells are the most promising cell line for further experiments aimed at exploring the molecular pathways of muscle pathologies.

In silico assessment of drug safety in human heart applied to late sodium current blockers.

Drug-induced action potential (AP) prolongation leading to Torsade de Pointes is a major concern for the development of anti-arrhythmic drugs. Nevertheless the development of improved anti-arrhythmic agents, some of which may block different channels, remains an important opportunity. Partial block of the late sodium current (I(NaL)) has emerged as a novel anti-arrhythmic mechanism. It can be effective in the settings of free radical challenge or hypoxia. In addition, this approach can attenuate pro-arrhythmic effects of blocking the rapid delayed rectifying K(+) current (I(Kr)). The main goal of our computational work was to develop an in-silico tool for preclinical anti-arrhythmic drug safety assessment, by illustrating the impact of I(Kr)/I(NaL) ratio of steady-state block of drug candidates on "torsadogenic" biomarkers. The O'Hara et al. AP model for human ventricular myocytes was used. Biomarkers for arrhythmic risk, i.e., AP duration, triangulation, reverse rate-dependence, transmural dispersion of repolarization and electrocardiogram QT intervals, were calculated using single myocyte and one-dimensional strand simulations. Predetermined amounts of block of I(NaL) and I(Kr) were evaluated. "Safety plots" were developed to illustrate the value of the specific biomarker for selected combinations of IC(50)s for I(Kr) and I(NaL) of potential drugs. The reference biomarkers at baseline changed depending on the "drug" specificity for these two ion channel targets. Ranolazine and GS967 (a novel potent inhibitor of I(NaL)) yielded a biomarker data set that is considered safe by standard regulatory criteria. This novel in-silico approach is useful for evaluating pro-arrhythmic potential of drugs and drug candidates in the human ventricle.

Biophysical model of ion transport across human respiratory epithelia allows quantification of ion permeabilities.

Lung health and normal mucus clearance depend on adequate hydration of airway surfaces. Because transepithelial osmotic gradients drive water flows, sufficient hydration of the airway surface liquid depends on a balance between ion secretion and absorption by respiratory epithelia. In vitro experiments using cultures of primary human nasal epithelia and human bronchial epithelia have established many of the biophysical processes involved in airway surface liquid homeostasis. Most experimental studies, however, have focused on the apical membrane, despite the fact that ion transport across respiratory epithelia involves both cellular and paracellular pathways. In fact, the ion permeabilities of the basolateral membrane and paracellular pathway remain largely unknown. Here we use a biophysical model for water and ion transport to quantify ion permeabilities of all pathways (apical, basolateral, paracellular) in human nasal epithelia cultures using experimental (Ussing Chamber and microelectrode) data reported in the literature. We derive analytical formulas for the steady-state short-circuit current and membrane potential, which are for polarized epithelia the equivalent of the Goldman-Hodgkin-Katz equation for single isolated cells. These relations allow parameter estimation to be performed efficiently. By providing a method to quantify all the ion permeabilities of respiratory epithelia, the model may aid us in understanding the physiology that regulates normal airway surface hydration.

A sodium-activated potassium channel supports high-frequency firing and reduces energetic costs during rapid modulations of action potential amplitude.

We investigated the ionic mechanisms that allow dynamic regulation of action potential (AP) amplitude as a means of regulating energetic costs of AP signaling. Weakly electric fish generate an electric organ discharge (EOD) by summing the APs of their electric organ cells (electrocytes). Some electric fish increase AP amplitude during active periods or social interactions and decrease AP amplitude when inactive, regulated by melanocortin peptide hormones. This modulates signal amplitude and conserves energy. The gymnotiform Eigenmannia virescens generates EODs at frequencies that can exceed 500 Hz, which is energetically challenging. We examined how E. virescens meets that challenge. E. virescens electrocytes exhibit a voltage-gated Na(+) current (I(Na)) with extremely rapid recovery from inactivation (τ(recov) = 0.3 ms) allowing complete recovery of Na(+) current between APs even in fish with the highest EOD frequencies. Electrocytes also possess an inwardly rectifying K(+) current and a Na(+)-activated K(+) current (I(KNa)), the latter not yet identified in any gymnotiform species. In vitro application of melanocortins increases electrocyte AP amplitude and the magnitudes of all three currents, but increased I(KNa) is a function of enhanced Na(+) influx. Numerical simulations suggest that changing I(Na) magnitude produces corresponding changes in AP amplitude and that K(Na) channels increase AP energy efficiency (10-30% less Na(+) influx/AP) over model cells with only voltage-gated K(+) channels. These findings suggest the possibility that E. virescens reduces the energetic demands of high-frequency APs through rapidly recovering Na(+) channels and the novel use of KNa channels to maximize AP amplitude at a given Na(+) conductance.

Markov models of use-dependence and reverse use-dependence during the mouse cardiac action potential.

The fast component of the cardiac transient outward current, I(Ktof), is blocked by a number of drugs. The major molecular bases of I(Ktof) are Kv4.2/Kv4.3 voltage-gated potassium channels. Drugs with similar potencies but different blocking mechanisms have differing effects on action potential duration (APD). We used in silico analysis to determine the effect of I(Ktof)-blocking drugs with different blocking mechanisms on mouse ventricular myocytes. We used our existing mouse model of the action potential, and developed 4 new Markov formulations for I(Ktof), I(Ktos), I(Kur), I(Ks). We compared effects of theoretical I(Ktof)-specific channel blockers: (1) a closed state, and (2) an open channel blocker. At concentrations lower or close to IC(50), the drug which bound to the open state always had a much greater effect on APD than the drug which bound to the closed state. At concentrations much higher than IC(50), both mechanisms had similar effects at very low pacing rates. However, an open state binding drug had a greater effect on APD at faster pacing rates, particularly around 10 Hz. In summary, our data indicate that drug effects on APD are strongly dependent not only on IC(50), but also on the drug binding state.

Dynamic monitoring of beating periodicity of stem cell-derived cardiomyocytes as a predictive tool for preclinical safety assessment.

BACKGROUND AND PURPOSE: Cardiac toxicity is a major concern in drug development and it is imperative that clinical candidates are thoroughly tested for adverse effects earlier in the drug discovery process. In this report, we investigate the utility of an impedance-based microelectronic detection system in conjunction with mouse embryonic stem cell-derived cardiomyocytes for assessment of compound risk in the drug discovery process. EXPERIMENTAL APPROACH: Beating of cardiomyocytes was measured by a recently developed microelectronic-based system using impedance readouts. We used mouse stem cell-derived cardiomyocytes to obtain dose-response profiles for over 60 compounds, including ion channel modulators, chronotropic/ionotropic agents, hERG trafficking inhibitors and drugs known to induce Torsades de Pointes arrhythmias. KEY RESULTS: This system sensitively and quantitatively detected effects of modulators of cardiac function, including some compounds missed by electrophysiology. Pro-arrhythmic compounds produced characteristic profiles reflecting arrhythmia, which can be used for identification of other pro-arrhythmic compounds. The time series data can be used to identify compounds that induce arrhythmia by complex mechanisms such as inhibition of hERG channels trafficking. Furthermore, the time resolution allows for assessment of compounds that simultaneously affect both beating and viability of cardiomyocytes. CONCLUSIONS AND IMPLICATIONS: Microelectronic monitoring of stem cell-derived cardiomyocyte beating provides a high throughput, quantitative and predictive assay system that can be used for assessment of cardiac liability earlier in the drug discovery process. The convergence of stem cell technology with microelectronic monitoring should facilitate cardiac safety assessment.

Monte Carlo study of gating and selection in potassium channels.

The study of selection and gating in potassium channels is a very important issue in modern biology. Indeed, such structures are known in essentially all types of cells in all organisms where they play many important functional roles. The mechanism of gating and selection of ionic species is not clearly understood. In this paper we study a model in which gating is obtained via an affinity-switching selectivity filter. We discuss the dependence of selectivity and efficiency on the cytosolic ionic concentration and on the typical pore open state duration. We demonstrate that a simple modification in the way in which the selectivity filter is modeled yields larger channel efficiency.

An exact stochastic hybrid model of excitable membranes including spatio-temporal evolution.

In this paper, we present a mathematical description for excitable biological membranes, in particular neuronal membranes. We aim to model the (spatio-) temporal dynamics, e.g., the travelling of an action potential along the axon, subject to noise, such as ion channel noise. Using the framework of Piecewise Deterministic Processes (PDPs) we provide an exact mathematical description-in contrast to pseudo-exact algorithms considered in the literature-of the stochastic process one obtains coupling a continuous time Markov chain model with a deterministic dynamic model of a macroscopic variable, that is coupling Markovian channel dynamics to the time-evolution of the transmembrane potential. We extend the existing framework of PDPs in finite dimensional state space to include infinite-dimensional evolution equations and thus obtain a stochastic hybrid model suitable for modelling spatio-temporal dynamics. We derive analytic results for the infinite-dimensional process, such as existence, the strong Markov property and its extended generator. Further, we exemplify modelling of spatially extended excitable membranes with PDPs by a stochastic hybrid version of the Hodgkin-Huxley model of the squid giant axon. Finally, we discuss the advantages of the PDP formulation in view of analytical and numerical investigations as well as the application of PDPs to structurally more complex models of excitable membranes.

Fluorescence detection of the movement of single KcsA subunits reveals cooperativity.

The prokaryotic KcsA channel is gated at the helical bundle crossing by intracellular protons and inactivates at the extracellular selectivity filter. The C-terminal transmembrane helix has to undergo a conformational change for potassium ions to access the central cavity. Whereas a partial opening of the tetrameric channel is suggested to be responsible for subconductance levels of ion channels, including KcsA, a cooperative opening of the 4 subunits is postulated as the final opening step. In this study, we used single-channel fluorescence spectroscopy of KcsA to directly observe the movement of each subunit and the temporal correlation between subunits. Purified KcsA channels labeled at the C terminus near the bundle crossing have been inserted into supported lipid bilayer, and the fluorescence traces analyzed by means of a cooperative or independent Markov model. The analysis revealed that the 4 subunits do not move fully independently but instead showed a certain degree of cooperativity. However, the 4 subunits do not simply open in 1 concerted step.

Modeling subunit cooperativity in opening of tetrameric ion channels.

Most potassium channels are tetramers of four homologous polypeptides (subunits). During channel gating, each subunit undergoes several conformational changes independent of the state of other subunits before reaching a permissive state, from which the channel can open. However, transition from the permissive states to the open state involves a concerted movement of all subunits. This cooperative transition must be included in Markov models of channel gating. Previously, it was implemented by considering all possible combinations of four subunit states in a much larger expanded model of channel states (e.g., 27,405 channel states versus 64 subunit states), which complicates modeling and is computationally intense, especially when accurate modeling requires a large number of subunit states. To overcome these complexities and retain the tetrameric molecular structure, a modeling approach was developed to incorporate the cooperative transition directly from the subunit models. In this approach, the open state is separated from the subunit models and represented by the net flux between the open state and the permissive states. Dynamic variations of the probability of state residencies computed using this direct approach and the expanded model were identical. Implementation of the direct approach is simple and its computational time is orders-of-magnitude shorter than the equivalent expanded model.

Optimization of chromone-2-carboxamide melanin concentrating hormone receptor 1 antagonists: assessment of potency, efficacy, and cardiovascular safety.

Evaluation of multiple structurally distinct series of melanin concentrating hormone receptor 1 antagonists in an anesthetized rat cardiovascualar assay led to the identification of a chromone-2-carboxamide series as having excellent safety against the chosen cardiovascular endpoints at high drug concentrations in the plasma and brain. Optimization of this series led to considerable improvements in affinity, functional potency, and pharmacokinetic profile. This led to the identification of a 7-fluorochromone-2-carboxamide (22) that was orally efficacious in a diet-induced obese mouse model, retained a favorable cardiovascular profile in rat, and demonstrated dramatic improvement in effects on mean arterial pressure in our dog cardiovascular model compared to other series reported by our group. However, this analogue also led to prolongation of the QT interval in the dog that was linked to affinity for hERG channel and unexpectedly potent functional blockade of this ion channel.

An overview of QT interval assessment in safety pharmacology.

Medicinal products that prolong cardiac repolarization unintentionally, as assessed in terms of prolongation of the QT interval of the electrocardiogram, may trigger a potentially fatal arrhythmia called torsade de pointe (TDP). This lethal risk necessitates a detailed preclinical evaluation before initiating clinical trials. There are two different and complementary approaches to assess the potential of drugs to cause QT interval prolongation. The in vivo approach provides information on the potential of the compound to prolong the QT interval under near-physiological conditions. It is mostly descriptive and not explanatory in terms of mechanisms of action. The in vitro approach provides much more mechanistic information, but is far removed from the clinical situation. While both approaches appear to possess reasonable predictive value, the results may depend largely on the experimental conditions employed. This unit reviews these issues and discusses a strategy aimed at understanding the problems associated with this cardiovascular risk.

In vitro and in vivo evaluation of novel glibenclamide derivatives as imaging agents for the non-invasive assessment of the pancreatic islet cell mass in animals and humans.

Pancreatic islet cell mass (PICM) is a major determinant of the insulin secretory capacity in humans. Currently, the only method for accurate assessment of the PICM is an autopsy study. Thus, development of a technique allowing the non-invasive quantification of PICM is of great interest. The aim of this study was to develop such a non-invasive technique featuring novel fluorine- and (99m)Tc-labelled glibenclamide derivatives. Despite the structural modifications necessary to introduce fluorine into the glibenclamide molecule, all derivatives retained insulin stimulating capacity as well as high affinity binding to human SUR1 when compared to the original glibenclamide. Contrastingly, the lipophilicity of the fluorine-labelled derivatives was altered depending on the particular modification. In the human PET-study a constant but weak radioactive signal could be detected in the pancreas using a fluorine-labelled glibenclamide derivative. However, a reliable assessment and visualisation of the PICM could not be obtained. It can be assumed that the high uptake of the fluorine-labelled tracer e.g. into the the liver and the high plasma protein binding leads to a relatively low signal-to-noise ratio. In case of the presented fluorine-labelled glibenclamide based compounds this could be the result of their invariably high lipophilicity. The development of a (99 m)Tc-labelled glibenclamide derivative with a lower lipophilicity and differing in vivo behaviour, glibenclamide based compounds for non-invasive imaging of the pancreatic islet cell mass may be possible.

Correlation character of ionic current fluctuations: analysis of ion current through a voltage-dependent potassium single channel.

The gating of ion channels has widely been modeled by assuming the transition between open and closed states is a memoryless process. Nevertheless, the statistical analysis of an ionic current signal recorded from voltage dependence K(+) single channel is presented. Calculating the sample auto-correlation function of the ionic current based on the digitized signals, rather than the sequence of open and closed states duration time. The results provide evidence for the existence of memory. For different voltages, the ion channel current fluctuation has different correlation attributions. The correlations in data generated by simulation of two Markov models, on one hand, auto-correlation function of the ionic current shows a weaker memory, after a delayed period of time, the attribute of memory does not exist; on the other hand, the correlation depends on the number of states in the Markov model. For V(p)=-60 mV pipette potential, spectral analysis of ion channel current was conducted, the result indicates that the spectrum is not a flat spectrum, the data set from ionic current fluctuations shows considerable variability with a broad 1/f -like spectrum, alpha=1.261+/-0.24. Thus the ion current fluctuations give information about the kinetics of the channel protein, the results suggest the correlation character of ion channel protein nonlinear kinetics regardless of whether the channel is in open or closed state.

Model selection and parameter estimation for ion channel recordings with an application to the K+ outward-rectifier in barley leaf.

We present a statistical method, and its accompanying algorithms, for the selection of a mathematical model of the gating mechanism of an ion channel and for the estimation of the parameters of this model. The method assumes a hidden Markov model that incorporates filtering, colored noise and state-dependent white excess noise for the recorded data. The model selection and parameter estimation are performed via a Bayesian approach using Markov chain Monte Carlo. The method is illustrated by its application to single-channel recordings of the K(+) outward-rectifier in barley leaf.