A model for stimulation of enzyme activity by a competitive inhibitor based on the interaction of terazosin and phosphoglycerate kinase 1.

The drug terazosin (TZ) binds to and can enhance the activity of the glycolytic enzyme phosphoglycerate kinase 1 (PGK1) and can increase ATP levels. That finding prompted studies of TZ in Parkinson's disease (PD) in which decreased neuronal energy metabolism is a hallmark feature. TZ was neuroprotective in cell-based and animal PD models and in large epidemiological studies of humans. However, how TZ might increase PGK1 activity has remained a perplexing question because structural data revealed that the site of TZ binding to PGK1 overlaps with the site of substrate binding, predicting that TZ would competitively inhibit activity. Functional data also indicate that TZ is a competitive inhibitor. To explore the paradoxical observation of a competitive inhibitor increasing enzyme activity under some conditions, we developed a mass action model of TZ and PGK1 interactions using published data on PGK1 kinetics and the effect of varying TZ concentrations. The model indicated that TZ-binding introduces a bypass pathway that accelerates product release. At low concentrations, TZ binding circumvents slow product release and increases the rate of enzymatic phosphotransfer. However, at high concentrations, TZ inhibits PGK1 activity. The model explains stimulation of enzyme activity by a competitive inhibitor and the biphasic dose-response relationship for TZ and PGK1 activity. By providing a plausible mechanism for interactions between TZ and PGK1, these findings may aid development of TZ or other agents as potential therapeutics for neurodegenerative diseases. The results may also have implications for agents that interact with the active site of other enzymes.

Atrial fibrillation and sinus node dysfunction in human ankyrin-B syndrome: a computational analysis.

Ankyrin-B is a multifunctional adapter protein responsible for localization and stabilization of select ion channels, transporters, and signaling molecules in excitable cells including cardiomyocytes. Ankyrin-B dysfunction has been linked with highly penetrant sinoatrial node (SAN) dysfunction and increased susceptibility to atrial fibrillation. While previous studies have identified a role for abnormal ion homeostasis in ventricular arrhythmias, the molecular mechanisms responsible for atrial arrhythmias and SAN dysfunction in human patients with ankyrin-B syndrome are unclear. Here, we develop a computational model of ankyrin-B dysfunction in atrial and SAN cells and tissue to determine the mechanism for increased susceptibility to atrial fibrillation and SAN dysfunction in human patients with ankyrin-B syndrome. Our simulations predict that defective membrane targeting of the voltage-gated L-type Ca(2+) channel Cav1.3 leads to action potential shortening that reduces the critical atrial tissue mass needed to sustain reentrant activation. In parallel, increased fibrosis results in conduction slowing that further increases the susceptibility to sustained reentry in the setting of ankyrin-B dysfunction. In SAN cells, loss of Cav1.3 slows spontaneous pacemaking activity, whereas defects in Na(+)/Ca(2+) exchanger and Na(+)/K(+) ATPase increase variability in SAN cell firing. Finally, simulations of the intact SAN reveal a shift in primary pacemaker site, SAN exit block, and even SAN failure in ankyrin-B-deficient tissue. These studies identify the mechanism for increased susceptibility to atrial fibrillation and SAN dysfunction in human disease. Importantly, ankyrin-B dysfunction involves changes at both the cell and tissue levels that favor the common manifestation of atrial arrhythmias and SAN dysfunction.

Ca2+/calmodulin-dependent protein kinase II-based regulation of voltage-gated Na+ channel in cardiac disease.

BACKGROUND: Human gene variants affecting ion channel biophysical activity and/or membrane localization are linked to potentially fatal cardiac arrhythmias. However, the mechanism for many human arrhythmia variants remains undefined despite more than a decade of investigation. Posttranslational modulation of membrane proteins is essential for normal cardiac function. Importantly, aberrant myocyte signaling has been linked to defects in cardiac ion channel posttranslational modifications and disease. We recently identified a novel pathway for posttranslational regulation of the primary cardiac voltage-gated Na(+) channel (Na(v)1.5) by Ca(2+)/calmodulin-dependent protein kinase II (CaMKII). However, a role for this pathway in cardiac disease has not been evaluated. METHODS AND RESULTS: We evaluated the role of CaMKII-dependent phosphorylation in human genetic and acquired disease. We report an unexpected link between a short motif in the Na(v)1.5 DI-DII loop, recently shown to be critical for CaMKII-dependent phosphorylation, and Na(v)1.5 function in monogenic arrhythmia and common heart disease. Experiments in heterologous cells and primary ventricular cardiomyocytes demonstrate that the human arrhythmia susceptibility variants (A572D and Q573E) alter CaMKII-dependent regulation of Na(v)1.5, resulting in abnormal channel activity and cell excitability. In silico analysis reveals that these variants functionally mimic the phosphorylated channel, resulting in increased susceptibility to arrhythmia-triggering afterdepolarizations. Finally, we report that this same motif is aberrantly regulated in a large-animal model of acquired heart disease and in failing human myocardium. CONCLUSIONS: We identify the mechanism for 2 human arrhythmia variants that affect Na(v)1.5 channel activity through direct effects on channel posttranslational modification. We propose that the CaMKII phosphorylation motif in the Na(v)1.5 DI-DII cytoplasmic loop is a critical nodal point for proarrhythmic changes to Na(v)1.5 in congenital and acquired cardiac disease.

A computational investigation of cardiac caveolae as a source of persistent sodium current.

Recent studies of cholesterol-rich membrane microdomains, called caveolae, reveal that caveolae are reservoirs of "recruitable" sodium ion channels. Caveolar channels constitute a substantial and previously unrecognized source of sodium current in cardiac cells. In this paper we model for the first time caveolar sodium currents and their contributions to cardiac action potential morphology. We show that the ß-agonist-induced opening of caveolae may have substantial impacts on peak overshoot, maximum upstroke velocity, and ultimately conduction velocity. Additionally, we show that prolonged action potentials and the formation of potentially arrhythmogenic afterdepolarizations, can arise if caveolae open intermittently throughout the action potential. Our simulations suggest that caveolar sodium current may constitute a route, which is independent of channelopathies, to delayed repolarization and the arrhythmias associated with such delays.

Defining new insight into atypical arrhythmia: a computational model of ankyrin-B syndrome.

Normal cardiac excitability depends on the coordinated activity of specific ion channels and transporters within specialized domains at the plasma membrane and sarcoplasmic reticulum. Ion channel dysfunction due to congenital or acquired defects has been linked to human cardiac arrhythmia. More recently, defects in ion channel-associated proteins have been associated with arrhythmia. Ankyrin-B is a multifunctional adapter protein responsible for targeting select ion channels, transporters, cytoskeletal proteins, and signaling molecules in excitable cells, including neurons, pancreatic ß-cells, and cardiomyocytes. Ankyrin-B dysfunction has been linked to cardiac arrhythmia in human patients and ankyrin-B heterozygous (ankyrin-B(+/-)) mice with a phenotype characterized by sinus node dysfunction, susceptibility to ventricular arrhythmias, and sudden death ("ankyrin-B syndrome"). At the cellular level, ankyrin-B(+/-) cells have defects in the expression and membrane localization of the Na(+)/Ca(2+) exchanger and Na(+)-K(+)-ATPase, Ca(2+) overload, and frequent afterdepolarizations, which likely serve as triggers for lethal cardiac arrhythmias. Despite knowledge gathered from mouse models and human patients, the molecular mechanism responsible for cardiac arrhythmias in the setting of ankyrin-B dysfunction remains unclear. Here, we use mathematical modeling to provide new insights into the cellular pathways responsible for Ca(2+) overload and afterdepolarizations in ankyrin-B(+/-) cells. We show that the Na(+)/Ca(2+) exchanger and Na(+)-K(+)-ATPase play related, yet distinct, roles in intracellular Ca(2+) accumulation, sarcoplasmic reticulum Ca(2+) overload, and afterdepolarization generation in ankyrin-B(+/-) cells. These findings provide important insights into the molecular mechanisms underlying a human disease and are relevant for acquired human arrhythmia, where ankyrin-B dysfunction has recently been identified.

A two-current model for the dynamics of cardiac membrane.

In this paper we introduce and study a model for electrical activity of cardiac membrane which incorporates only an inward and an outward current. This model is useful for three reasons: (1) Its simplicity, comparable to the FitzHugh-Nagumo model, makes it useful in numerical simulations, especially in two or three spatial dimensions where numerical efficiency is so important. (2) It can be understood analytically without recourse to numerical simulations. This allows us to determine rather completely how the parameters in the model affect its behavior which in turn provides insight into the effects of the many parameters in more realistic models. (3) It naturally gives rise to a one-dimensional map which specifies the action potential duration as a function of the previous diastolic interval. For certain parameter values, this map exhibits a new phenomenon--subcritical alternans--that does not occur for the commonly used exponential map.

Precision of neural timing: effects of convergence and time-windowing.

We study the improvement in timing accuracy in a neural system having n identical input neurons projecting to one target neuron. The n input neurons receive the same stimulus but fire at stochastic times selected from one of four specified probability densities, f, each with standard deviation 1.0 msec. The target cell fires if and when it receives m inputs within a time window of epsilon msec. Let sigma(n,m,epsilon) denote the standard deviation of the time of firing of the target neuron (i.e. the standard deviation of the target neuron's latency relative to the arrival time of the stimulus). Mathematical analysis shows that sigma(n,m,epsilon) is a very complicated function of n, m, and epsilon. Typically, sigma(n,m,epsilon) is a non-monotone function of m and epsilon and the improvement of timing accuracy is highly dependent of the shape of the probability density for the time of firing of the input neurons. For appropriate choices of m, epsilon, and f, the standard deviation sigma(n,m,epsilon) may be as low as 1/n. Thus, depending on these variables, remarkable improvements in timing accuracy of such a stochastic system may occur.