Modulating the voltage sensor of a cardiac potassium channel shows antiarrhythmic effects.

Cardiac arrhythmias are the most common cause of sudden cardiac death worldwide. Lengthening the ventricular action potential duration (APD), either congenitally or via pathologic or pharmacologic means, predisposes to a life-threatening ventricular arrhythmia, Torsade de Pointes. IKs (KCNQ1+KCNE1), a slowly activating K+ current, plays a role in action potential repolarization. In this study, we screened a chemical library in silico by docking compounds to the voltage-sensing domain (VSD) of the IKs channel. Here, we show that C28 specifically shifted IKs VSD activation in ventricle to more negative voltages and reversed the drug-induced lengthening of APD. At the same dosage, C28 did not cause significant changes of the normal APD in either ventricle or atrium. This study provides evidence in support of a computational prediction of IKs VSD activation as a potential therapeutic approach for all forms of APD prolongation. This outcome could expand the therapeutic efficacy of a myriad of currently approved drugs that may trigger arrhythmias.

A study of the outward background current conductance gK1, the pacemaker current conductance gf, and the gap junction conductance gj as determinants of biological pacing in single cells and in a two-cell syncytium using the dynamic clamp.

We previously demonstrated that a two-cell syncytium, composed of a ventricular myocyte and an mHCN2 expressing cell, recapitulated most properties of in vivo biological pacing induced by mHCN2-transfected hMSCs in the canine ventricle. Here, we use the two-cell syncytium, employing dynamic clamp, to study the roles of gf (pacemaker conductance), gK1 (background K+ conductance), and gj (intercellular coupling conductance) in biological pacing. We studied gf and gK1 in single HEK293 cells expressing cardiac sodium current channel Nav1.5 (SCN5A). At fixed gf, increasing gK1 hyperpolarized the cell and initiated pacing. As gK1 increased, rate increased, then decreased, finally ceasing at membrane potentials near EK. At fixed gK1, increasing gf depolarized the cell and initiated pacing. With increasing gf, rate increased reaching a plateau, then decreased, ceasing at a depolarized membrane potential. We studied gj via virtual coupling with two non-adjacent cells, a driver (HEK293 cell) in which gK1 and gf were injected without SCN5A and a follower (HEK293 cell), expressing SCN5A. At the chosen values of gK1 and gf oscillations initiated in the driver, when gj was increased synchronized pacing began, which then decreased by about 35% as gj approached 20 nS. Virtual uncoupling yielded similar insights into gj. We also studied subthreshold oscillations in physically and virtually coupled cells. When coupling was insufficient to induce pacing, passive spread of the oscillations occurred in the follower. These results show a non-monotonic relationship between gK1, gf, gj, and pacing. Further, oscillations can be generated by gK1 and gf in the absence of SCN5A.

MATLAB implementation of a dynamic clamp with bandwidth of >125 kHz capable of generating I Na at 37 °C.

We describe the construction of a dynamic clamp with a bandwidth of >125 kHz that utilizes a high-performance, yet low-cost, standard home/office PC interfaced with a high-speed (16 bit) data acquisition module. High bandwidth is achieved by exploiting recently available software advances (code-generation technology and optimized real-time kernel). Dynamic-clamp programs are constructed using Simulink, a visual programming language. Blocks for computation of membrane currents are written in the high-level MATLAB language; no programming in C is required. The instrument can be used in single- or dual-cell configurations, with the capability to modify programs while experiments are in progress. We describe an algorithm for computing the fast transient Na(+) current (I Na) in real time and test its accuracy and stability using rate constants appropriate for 37 °C. We then construct a program capable of supplying three currents to a cell preparation: I Na, the hyperpolarizing-activated inward pacemaker current (I f) and an inward-rectifier K(+) current (I K1). The program corrects for the IR drop due to electrode current flow and also records all voltages and currents. We tested this program on dual patch-clamped HEK293 cells where the dynamic clamp controls a current-clamp amplifier and a voltage-clamp amplifier controls membrane potential, and current-clamped HEK293 cells where the dynamic clamp produces spontaneous pacing behavior exhibiting Na(+) spikes in otherwise passive cells.

Transport efficiency and workload distribution in a mathematical model of the thick ascending limb.

The thick ascending limb (TAL) is a major NaCl reabsorbing site in the nephron. Efficient reabsorption along that segment is thought to be a consequence of the establishment of a strong transepithelial potential that drives paracellular Na(+) uptake. We used a multicell mathematical model of the TAL to estimate the efficiency of Na(+) transport along the TAL and to examine factors that determine transport efficiency, given the condition that TAL outflow must be adequately dilute. The TAL model consists of a series of epithelial cell models that represent all major solutes and transport pathways. Model equations describe luminal flows, based on mass conservation and electroneutrality constraints. Empirical descriptions of cell volume regulation (CVR) and pH control were implemented, together with the tubuloglomerular feedback (TGF) system. Transport efficiency was calculated as the ratio of total net Na(+) transport (i.e., paracellular and transcellular transport) to transcellular Na(+) transport. Model predictions suggest that 1) the transepithelial Na(+) concentration gradient is a major determinant of transport efficiency; 2) CVR in individual cells influences the distribution of net Na(+) transport along the TAL; 3) CVR responses in conjunction with TGF maintain luminal Na(+) concentration well above static head levels in the cortical TAL, thereby preventing large decreases in transport efficiency; and 4) under the condition that the distribution of Na(+) transport along the TAL is quasi-uniform, the tubular fluid axial Cl(-) concentration gradient near the macula densa is sufficiently steep to yield a TGF gain consistent with experimental data.

Fluid dilution and efficiency of Na(+) transport in a mathematical model of a thick ascending limb cell.

Thick ascending limb (TAL) cells are capable of reducing tubular fluid Na(+) concentration to as low as ~25 mM, and yet they are thought to transport Na(+) efficiently owing to passive paracellular Na(+) absorption. Transport efficiency in the TAL is of particular importance in the outer medulla where O(2) availability is limited by low blood flow. We used a mathematical model of a TAL cell to estimate the efficiency of Na(+) transport and to examine how tubular dilution and cell volume regulation influence transport efficiency. The TAL cell model represents 13 major solutes and the associated transporters and channels; model equations are based on mass conservation and electroneutrality constraints. We analyzed TAL transport in cells with conditions relevant to the inner stripe of the outer medulla, the cortico-medullary junction, and the distal cortical TAL. At each location Na(+) transport efficiency was computed as functions of changes in luminal NaCl concentration ([NaCl]), [K(+)], [NH(4)(+)], junctional Na(+) permeability, and apical K(+) permeability. Na(+) transport efficiency was calculated as the ratio of total net Na(+) transport to transcellular Na(+) transport. Transport efficiency is predicted to be highest at the cortico-medullary boundary where the transepithelial Na(+) gradient is the smallest. Transport efficiency is lowest in the cortex where luminal [NaCl] approaches static head.

Suppression of phosphoinositide 3-kinase signaling and alteration of multiple ion currents in drug-induced long QT syndrome.

Many drugs, including some commonly used medications, can cause abnormal heart rhythms and sudden death, as manifest by a prolonged QT interval in the electrocardiogram. Cardiac arrhythmias caused by drug-induced long QT syndrome are thought to result mainly from reductions in the delayed rectifier potassium ion (K(+)) current I(Kr). Here, we report a mechanism for drug-induced QT prolongation that involves changes in multiple ion currents caused by a decrease in phosphoinositide 3-kinase (PI3K) signaling. Treatment of canine cardiac myocytes with inhibitors of tyrosine kinases or PI3Ks caused an increase in action potential duration that was reversed by intracellular infusion of phosphatidylinositol 3,4,5-trisphosphate. The inhibitors decreased the delayed rectifier K(+) currents I(Kr) and I(Ks), the L-type calcium ion (Ca(2+)) current I(Ca,L), and the peak sodium ion (Na(+)) current I(Na) and increased the persistent Na(+) current I(NaP). Computer modeling of the canine ventricular action potential showed that the drug-induced change in any one current accounted for less than 50% of the increase in action potential duration. Mouse hearts lacking the PI3K p110α catalytic subunit exhibited a prolonged action potential and QT interval that were at least partly a result of an increase in I(NaP). These results indicate that down-regulation of PI3K signaling directly or indirectly via tyrosine kinase inhibition prolongs the QT interval by affecting multiple ion channels. This mechanism may explain why some tyrosine kinase inhibitors in clinical use are associated with increased risk of life-threatening arrhythmias.

Epicardial border zone overexpression of skeletal muscle sodium channel SkM1 normalizes activation, preserves conduction, and suppresses ventricular arrhythmia: an in silico, in vivo, in vitro study.

BACKGROUND: In depolarized myocardial infarct epicardial border zones, the cardiac sodium channel (SCN5A) is largely inactivated, contributing to low action potential upstroke velocity (V(max)), slow conduction, and reentry. We hypothesized that a fast inward current such as the skeletal muscle sodium channel (SkM1) operating more effectively at depolarized membrane potentials might restore fast conduction in epicardial border zones and be antiarrhythmic. METHODS AND RESULTS: Computer simulations were done with a modified Hund-Rudy model. Canine myocardial infarcts were created by coronary ligation. Adenovirus expressing SkM1 and green fluorescent protein or green fluorescent protein alone (sham) was injected into epicardial border zones. After 5 to 7 days, dogs were studied with epicardial mapping, programmed premature stimulation in vivo, and cellular electrophysiology in vitro. Infarct size was determined, and tissues were immunostained for SkM1 and green fluorescent protein. In the computational model, modest SkM1 expression preserved fast conduction at potentials as positive as -60 mV; overexpression of SCN5A did not. In vivo epicardial border zone electrograms were broad and fragmented in shams (31.5 +/- 2.3 ms) and narrower in SkM1 (22.6 +/- 2.8 ms; P=0.03). Premature stimulation induced ventricular tachyarrhythmia/fibrillation >60 seconds in 6 of 8 shams versus 2 of 12 SkM1 (P=0.02). Microelectrode studies of epicardial border zones from SkM1 showed membrane potentials equal to that of shams and V(max) greater than that of shams as membrane potential depolarized (P<0.01). Infarct sizes were similar (sham, 30 +/- 2.8%; SkM1, 30 +/- 2.6%; P=0.86). SkM1 expression in injected epicardium was confirmed immunohistochemically. CONCLUSIONS: SkM1 increases V(max) of depolarized myocardium and reduces the incidence of inducible sustained ventricular tachyarrhythmia/fibrillation in canine infarcts. Gene therapy to normalize activation by increasing V(max) at depolarized potentials may be a promising antiarrhythmic strategy.

Parameter estimation for mathematical models of NKCC2 cotransporter isoforms.

An optimization problem, formulated using a nonlinear least-squares approach, was used to estimate parameters for kinetic models of the three isoforms of the kidney-specific Na-K-2Cl (NKCC2) cotransporter. Specifically, the optimization problem estimates the magnitude of model parameters (i.e., off-binding and translocation rate constants) by minimizing the distance between model unidirectional fluxes and published unidirectional (86)Rb(+) uptake curves for the A, B, and F isoforms of the NKCC2 cotransporter obtained in transfected Xenopus oocytes. By using different symmetry assumptions, NKCC2 models with five, six, seven, or eight parameters were evaluated. The optimization method identified parameter sets that yielded computed unidirectional fluxes consistent with the uptake data. However, the parameter values were not unique, in that systematic exploration of the parameter space revealed alternative parameter sets that fit the data with similar accuracy. Finally, we demonstrate that the optimization method can identify parameter sets for the three transporter isoforms that differ only in ion binding affinities, a result that is consistent with a published mutagenesis analysis of the molecular and structural bases for the differences in (86)Rb(+) uptake among the A, B, and F isoforms. These NKCC2 cotransporter models will facilitate the development of larger scale models of ion transport by thick ascending limb cells.

A mathematical model of pH(i) regulation in central CO2- chemoreception.

In CO2 chemosensitive neurons, an increase in CO2 (hypercapnia) leads to a maintained reduction in intracellular pH (pH(i)) while in non-chemosensitive neurons pH(i) recovery is observed. The precise mechanisms for the differential regulation of pH(i) recovery between these cell populations remain to be identified; however, studies have begun to explore the role of Na+/H+ exchange (NHE). Here, we compare the results of two different formulations of a mathematical model to begin to explore pH(i) regulation in central CO2 chemoreception.

Differences in pHi recovery in CO2-chemosensitive and non-chemosensitive cells: predictions from a mathematical model.

In this paper, we present a mathematic model designed to identify potential mechanisms responsible for the observed differences in pHi recovery in CO(2)-chemosensitive versus non-chemosensitive cells. The model suggests that differences in pHi regulation may be dependent upon differences in the activation set-point of the internal modifier site of the Na(+)/H(+) exchanger (NHE).

Field demonstration of DNAPL dehalogenation using emulsified zero-valent iron.

This paper describes the results of the first field-scale demonstration conducted to evaluate the performance of nanoscale emulsified zero-valent iron (EZVI) injected into the saturated zone to enhance in situ dehalogenation of dense, nonaqueous phase liquids (DNAPLs) containing trichloroethene (TCE). EZVI is an innovative and emerging remediation technology. EZVI is a surfactant-stabilized, biodegradable emulsion that forms emulsion droplets consisting of an oil-liquid membrane surrounding zero-valent iron (ZVI) particles in water. EZVI was injected over a five day period into eight wells in a demonstration test area within a larger DNAPL source area at NASA's Launch Complex 34 (LC34) using a pressure pulse injection method. Soil and groundwater samples were collected before and after treatment and analyzed for volatile organic compounds (VOCs) to evaluate the changes in VOC mass, concentration and mass flux. Significant reductions in TCE soil concentrations (>80%) were observed at four of the six soil sampling locations within 90 days of EZVI injection. Somewhat lower reductions were observed at the other two soil sampling locations where visual observations suggest that most of the EZVI migrated up above the target treatment depth. Significant reductions in TCE groundwater concentrations (57 to 100%) were observed at all depths targeted with EZVI. Groundwater samples from the treatment area also showed significant increases in the concentrations of cis-1,2-dichloroethene (cDCE), vinyl chloride (VC) and ethene. The decrease in concentrations of TCE in soil and groundwater samples following treatment with EZVI is believed to be due to abiotic degradation associated with the ZVI as well as biodegradation enhanced by the presence of the oil and surfactant in the EZVI emulsion.