Co-expression of calcium and hERG potassium channels reduces the incidence of proarrhythmic events

Relationship between ion currents and membrane capacitance in canine ventricular myocytes

Current density, the membrane current value divided by membrane capacitance (Cm), is widely used in cellular electrophysiology. Comparing current densities obtained in different cell populations assume that Cm and ion current magnitudes are linearly related, however data is scarce about this in cardiomyocytes. Therefore, we statistically analyzed the distributions, and the relationship between parameters of canine cardiac ion currents and Cm, and tested if dividing original parameters with Cm had any effect. Under conventional voltage clamp conditions, correlations were high for IK1, moderate for IKr and ICa,L, while negligible for IKs. Correlation between Ito1 peak amplitude and Cm was negligible when analyzing all cells together, however, the analysis showed high correlations when cells of subepicardial, subendocardial or midmyocardial origin were analyzed separately. In action potential voltage clamp experiments IK1, IKr and ICa,L parameters showed high correlations with Cm. For INCX, INa,late and IKs there were low-to-moderate correlations between Cm and these current parameters. Dividing the original current parameters with Cm reduced both the coefficient of variation, and the deviation from normal distribution. The level of correlation between ion currents and Cm varies depending on the ion current studied. This must be considered when evaluating ion current densities in cardiac cells.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

Assessing drug safety by identifying the axis of arrhythmia in cardiomyocyte electrophysiology

Many classes of drugs can induce fatal cardiac arrhythmias by disrupting the electrophysiology of cardiomyocytes. Safety guidelines thus require all new drugs to be assessed for pro-arrhythmic risk prior to conducting human trials. The standard safety protocols primarily focus on drug blockade of the delayed-rectifier potassium current (IKr). Yet the risk is better assessed using four key ion currents (IKr, ICaL, INaL, IKs). We simulated 100,000 phenotypically diverse cardiomyocytes to identify the underlying relationship between the blockade of those currents and the emergence of ectopic beats in the action potential. We call that relationship the axis of arrhythmia. It serves as a yardstick for quantifying the arrhythmogenic risk of any drug from its profile of multi-channel block alone. We tested it on 109 drugs and found that it predicted the clinical risk labels with an accuracy of 88.1-90.8%. Pharmacologists can use our method to assess the safety of novel drugs without resorting to animal testing or unwieldy computer simulations.

Targeted activation of human ether-à-go-go-related gene channels rescues electrical instability induced by the R56Q+/- long QT syndrome variant

AIMS

Long QT syndrome type 2 (LQTS2) is associated with inherited variants in the cardiac human ether-à-go-go-related gene (hERG) K+ channel. However, the pathogenicity of hERG channel gene variants is often uncertain. Using CRISPR-Cas9 gene-edited hiPSC-derived cardiomyocytes (hiPSC-CMs), we investigated the pathogenic mechanism underlying the LQTS-associated hERG R56Q variant and its phenotypic rescue by using the Type 1 hERG activator, RPR260243.

METHODS AND RESULTS

The above approaches enable characterization of the unclear causative mechanism of arrhythmia in the R56Q variant (an N-terminal PAS domain mutation that primarily accelerates channel deactivation) and translational investigation of the potential for targeted pharmacologic manipulation of hERG deactivation. Using perforated patch clamp electrophysiology of single hiPSC-CMs, programmed electrical stimulation showed that the hERG R56Q variant does not significantly alter the mean action potential duration (APD90). However, the R56Q variant increases the beat-to-beat variability in APD90 during pacing at constant cycle lengths, enhances the variance of APD90 during rate transitions, and increases the incidence of 2:1 block. During paired S1-S2 stimulations measuring electrical restitution properties, the R56Q variant was also found to increase the variability in rise time and duration of the response to premature stimulations. Application of the hERG channel activator, RPR260243, reduces the APD variance in hERG R56Q hiPSC-CMs, reduces the variability in responses to premature stimulations, and increases the post-repolarization refractoriness.

CONCLUSION

Based on our findings, we propose that the hERG R56Q variant leads to heterogeneous APD dynamics, which could result in spatial dispersion of repolarization and increased risk for re-entry without significantly affecting the average APD90. Furthermore, our data highlight the antiarrhythmic potential of targeted slowing of hERG deactivation gating, which we demonstrate increases protection against premature action potentials and reduces electrical heterogeneity in hiPSC-CMs.

Pairwise biosynthesis of ion channels stabilizes excitability and mitigates arrhythmias

Coordinated expression of ion channels is crucial for cardiac rhythms, neural signaling, and cell cycle progression. Perturbation of this balance results in many disorders including cardiac arrhythmias. Prior work revealed association of mRNAs encoding cardiac NaV1.5 (SCN5A) and hERG1 (KCNH2), but the functional significance of this association was not established. Here, we provide a more comprehensive picture of KCNH2, SCN5A, CACNA1C, and KCNQ1 transcripts collectively copurifying with nascent hERG1, NaV1.5, CaV1.2, or KCNQ1 channel proteins. Single-molecule fluorescence in situ hybridization (smFISH) combined with immunofluorescence reveals that the channel proteins are synthesized predominantly as heterotypic pairs from discrete molecules of mRNA, not as larger cotranslational complexes. Puromycin disrupted colocalization of mRNA with its encoded protein, as expected, but remarkably also pairwise mRNA association, suggesting that transcript association relies on intact translational machinery or the presence of the nascent protein. Targeted depletion of KCHN2 by specific shRNA resulted in concomitant reduction of all associated mRNAs, with a corresponding reduction in the encoded channel currents. This co-knockdown effect, originally described for KCNH2 and SCN5A, thus appears to be a general phenomenon among transcripts encoding functionally related proteins. In multielectrode array recordings, proarrhythmic behavior arose when IKr was reduced by the selective blocker dofetilide at IC50 concentrations, but not when equivalent reductions were mediated by shRNA, suggesting that co-knockdown mitigates proarrhythmic behavior expected from the selective reduction of a single channel species. We propose that coordinated, cotranslational association of functionally related ion channel mRNAs confers electrical stability by co-regulating complementary ion channels in macromolecular complexes.

Qi-Po-Sheng-Mai granule ameliorates Ach-CaCl2 -induced atrial fibrillation by regulating calcium homeostasis in cardiomyocytes

BACKGROUND

Atrial fibrillation (AF) is one of the most common arrhythmias encountered in clinical settings. Currently, the pathophysiology of AF remains unclear, which severely limits the effectiveness and safety of medical therapies. The Chinese herbal formula Qi-Po-Sheng-Mai Granule (QPSM) has been widely used in China to treat AF. However, its pharmacological and molecular mechanisms remain unknown.

PURPOSE

The purpose of this study was to investigate the molecular mechanisms and potential targets of QPSM for AF.

STUDY DESIGN AND METHODS

The AF model was induced by Ach (66 μg/ml) and CaCl2 (10 mg/kg), and the dose of 0.1 ml/100 g was injected into the tail vein for 5 weeks. QPSM was administered daily at doses of 4.42 and 8.84 g/kg, and amiodarone (0.18 g/kg) was used as the positive control. The effect of QPSM on AF was assessed by electrocardiogram, echocardiography, and histopathological analysis. Then, we employed network pharmacology with single nucleus RNA sequencing (snRNA-Seq) to investigate the molecular mechanisms and potential targets of QPSM for AF. Furthermore, high performance liquid chromatography (HPLC) method was used for component analysis of QPSM, and molecular docking was used to verify the potential targets. Using the IonOptix single cell contraction and ion synchronization test equipment, single myocyte length and calcium ion variations were observed in real time. The expression levels of calcium Transporter-related proteins were detected by western blot and immunohistochemistry.

RESULTS

Based on an Ach-CaCl2-induced AF model, we found that QPSM treatment significantly reduced atrial electrical remodeling-related markers, such as AF inducibility and duration, and attenuated atrial dilation and fibrosis. Network pharmacology identified 52 active ingredients and 119 potential targets for QPSM in the treatment of AF, and 45 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were enriched, among which calcium pathway had the greatest impact. Using single nucleus sequencing (snRNA-seq), we identified cardiomyocytes as the most differentially expressed in response to drug treatment, with nine differentially expressed genes enriched in calcium signaling pathways. High performance liquid chromatography and molecular docking confirmed that the core components of QPSM strongly bind to the key factors in the calcium signaling pathway. Additional experiments have shown that QPSM increases calcium transients (CaT) and contractility in the individual cardiomyocyte. This was accomplished by increasing the expression of CACNA1C and SERCA2a and decreasing the expression of CAMK2B and NCX1.

CONCLUSION

The present study has systematically elucidated the role of QPSM in maintaining calcium homeostasis in cardiomyocytes through the regulation of calcium transporters, which could lead to new drug development ideas for AF.

Diversity of cells and signals in the cardiovascular system

This white paper is the outcome of the seventh UC Davis Cardiovascular Research Symposium on Systems Approach to Understanding Cardiovascular Disease and Arrhythmia. This biannual meeting aims to bring together leading experts in subfields of cardiovascular biomedicine to focus on topics of importance to the field. The theme of the 2022 Symposium was 'Cell Diversity in the Cardiovascular System, cell-autonomous and cell-cell signalling'. Experts in the field contributed their experimental and mathematical modelling perspectives and discussed emerging questions, controversies, and challenges in examining cell and signal diversity, co-ordination and interrelationships involved in cardiovascular function. This paper originates from the topics of formal presentations and informal discussions from the Symposium, which aimed to develop a holistic view of how the multiple cell types in the cardiovascular system integrate to influence cardiovascular function, disease progression and therapeutic strategies. The first section describes the major cell types (e.g. cardiomyocytes, vascular smooth muscle and endothelial cells, fibroblasts, neurons, immune cells, etc.) and the signals involved in cardiovascular function. The second section emphasizes the complexity at the subcellular, cellular and system levels in the context of cardiovascular development, ageing and disease. Finally, the third section surveys the technological innovations that allow the interrogation of this diversity and advancing our understanding of the integrated cardiovascular function and dysfunction.

Initiation of ventricular arrhythmia in the acquired long QT syndrome

AIMS

Long QT syndrome (LQTS) carries a risk of life-threatening polymorphic ventricular tachycardia (Torsades de Pointes, TdP) and is a major cause of premature sudden cardiac death. TdP is induced by R-on-T premature ventricular complexes (PVCs), thought to be generated by cellular early-afterdepolarisations (EADs). However, EADs in tissue require cellular synchronisation, and their role in TdP induction remains unclear. We aimed to determine the mechanism of TdP induction in rabbit hearts with acquired LQTS (aLQTS).

METHODS AND RESULTS

Optical mapping of action potentials (APs) and intracellular Ca2+ was performed in Langendorff-perfused rabbit hearts (n = 17). TdP induced by R-on-T PVCs was observed during aLQTS (50% K+/Mg++ & E4031) conditions in all hearts (P < 0.0001 vs. control). Islands of AP prolongation bounded by steep voltage gradients (VGs) were consistently observed before arrhythmia and peak VGs were more closely related to the PVC upstroke than EADs, both temporally (7 ± 5 ms vs. 44 ± 27 ms, P < 0.0001) and spatially (1.0 ± 0.7 vs. 3.6 ± 0.9 mm, P < 0.0001). PVCs were initiated at estimated voltages of ∼ -40 mV and had upstroke dF/dtmax and Vm-Ca2+ dynamics compatible with ICaL activation. Computational simulations demonstrated that PVCs could arise directly from VGs, through electrotonic triggering of ICaL. In experiments and the model, sub-maximal L-type Ca2+ channel (LTCC) block (200 nM nifedipine and 90% gCaL, respectively) abolished both PVCs and TdP in the continued presence of aLQTS.

CONCLUSION

These data demonstrate that ICaL activation at sites displaying steep VGs generates the PVCs which induce TdP, providing a mechanism and rationale for LTCC blockers as a novel therapeutic approach in LQTS.

Modeling incomplete penetrance in long QT syndrome type 3 through ion channel heterogeneity: an in silico population study

Many cardiac diseases are characterized by an increased late sodium current, including heart failure, hypertrophic cardiomyopathy, and inherited long QT syndrome type 3 (LQT3). The late sodium current in LQT3 is caused by a gain-of-function mutation in the voltage-gated sodium channel Nav1.5. Despite a well-defined genetic cause of LQT3, treatment remains inconsistent because of incomplete penetrance of the mutation and variability of antiarrhythmic efficacy. Here, we investigate the relationship between LQT3-associated mutation incomplete penetrance and variability in ion channel expression, simulating a population of 1,000 individuals with the O'Hara-Rudy model of the human ventricular myocyte. We first simulate healthy electrical activity (i.e., in the absence of a mutation) and then incorporate heterozygous expression for three LQT3-associated mutations (Y1795C, I1768V, and ΔKPQ), to directly compare the effects of each mutation on individuals across a diverse population. For all mutations, we find that susceptibility, defined by either the presence of an early afterdepolarization (EAD) or prolonged action potential duration (APD), primarily depends on the balance between the conductance of IKr and INa, for which individuals with a higher IKr-to-INa ratio are less susceptible. Furthermore, we find distinct differences across the population, observing individuals susceptible to zero, one, two, or all three mutations. Individuals tend to be less susceptible with an appropriate balance of repolarizing currents, typically via increased IKs or IK1. Interestingly, the more critical repolarizing current is mutation specific. We conclude that the balance between key currents plays a significant role in mutant-specific presentation of the disease phenotype in LQT3.NEW & NOTEWORTHY An in silico population approach investigates the relationship between variability in ion channel expression and gain-of-function mutations in the voltage-gated sodium channel associated with the congenital disorder long QT syndrome type 3 (LQT3). We find that ion channel variability can contribute to incomplete penetrance of the mutation, with mutant-specific differences in ion channel conductances leading to susceptibility to proarrhythmic action potential duration prolongation or early afterdepolarizations.

Gene expression analyses of TAS1R taste receptors relevant to the treatment of cardiometabolic disease

The sweet taste receptor (STR) is a G protein-coupled receptor (GPCR) responsible for mediating cellular responses to sweet stimuli. Early evidence suggests that elements of the STR signaling system are present beyond the tongue in metabolically active tissues, where it may act as an extraoral glucose sensor. This study aimed to delineate expression of the STR in extraoral tissues using publicly available RNA-sequencing repositories. Gene expression data was mined for all genes implicated in the structure and function of the STR, and control genes including highly expressed metabolic genes in relevant tissues, other GPCRs and effector G proteins with physiological roles in metabolism, and other GPCRs with expression exclusively outside the metabolic tissues. Since the physiological role of the STR in extraoral tissues is likely related to glucose sensing, expression was then examined in diseases related to glucose-sensing impairment such as type 2 diabetes. An aggregate co-expression network was then generated to precisely determine co-expression patterns among the STR genes in these tissues. We found that STR gene expression was negligible in human pancreatic and adipose tissues, and low in intestinal tissue. Genes encoding the STR did not show significant co-expression or connectivity with other functional genes in these tissues. In addition, STR expression was higher in mouse pancreatic and adipose tissues, and equivalent to human in intestinal tissue. Our results suggest that STR expression in mice is not representative of expression in humans, and the receptor is unlikely to be a promising extraoral target in human cardiometabolic disease.

Electrophysiological heterogeneity in large populations of rabbit ventricular cardiomyocytes

AIMS

Cardiac electrophysiological heterogeneity includes: (i) regional differences in action potential (AP) waveform, (ii) AP waveform differences in cells isolated from a single region, (iii) variability of the contribution of individual ion currents in cells with similar AP durations (APDs). The aim of this study is to assess intra-regional AP waveform differences, to quantify the contribution of specific ion channels to the APD via drug responses and to generate a population of mathematical models to investigate the mechanisms underlying heterogeneity in rabbit ventricular cells.

METHODS AND RESULTS

APD in ∼50 isolated cells from subregions of the LV free wall of rabbit hearts were measured using a voltage-sensitive dye. When stimulated at 2 Hz, average APD90 value in cells from the basal epicardial region was 254 ± 25 ms (mean ± standard deviation) in 17 hearts with a mean interquartile range (IQR) of 53 ± 17 ms. Endo-epicardial and apical-basal APD90 differences accounted for ∼10% of the IQR value. Highly variable changes in APD occurred after IK(r) or ICa(L) block that included a sub-population of cells (HR) with an exaggerated (hyper) response to IK(r) inhibition. A set of 4471 AP models matching the experimental APD90 distribution was generated from a larger population of models created by random variation of the maximum conductances (Gmax) of 8 key ion channels/exchangers/pumps. This set reproduced the pattern of cell-specific responses to ICa(L) and IK(r) block, including the HR sub-population. The models exhibited a wide range of Gmax values with constrained relationships linking ICa(L) with IK(r), ICl, INCX, and INaK.

CONCLUSION

Modelling the measured range of inter-cell APDs required a larger range of key Gmax values indicating that ventricular tissue has considerable inter-cell variation in channel/pump/exchanger activity. AP morphology is retained by relationships linking specific ionic conductances. These interrelationships are necessary for stable repolarization despite large inter-cell variation of individual conductances and this explains the variable sensitivity to ion channel block.

Sex-dependent transcription of cardiac electrophysiology and links to acetylation modifiers based on the GTEx database

Development of safer drugs based on epigenetic modifiers, e.g., histone deacetylase inhibitors (HDACi), requires better understanding of their effects on cardiac electrophysiology. Using RNAseq data from the genotype-tissue-expression database (GTEx), we created models that link the abundance of acetylation enzymes (HDAC/SIRT/HATs), and the gene expression of ion channels (IC) via select cardiac transcription factors (TFs) in male and female adult human hearts (left ventricle, LV). Gene expression data (transcripts per million, TPM) from GTEx donors (21–70 y.o.) were filtered, normalized and transformed to Euclidian space to allow quantitative comparisons in 84 female and 158 male LVs. Sex-specific partial least-square (PLS) regression models, linking gene expression data for HDAC/SIRT/HATs to TFs and to ICs gene expression, revealed tight co-regulation of cardiac ion channels by HDAC/SIRT/HATs, with stronger clustering in the male LV. Co-regulation of genes encoding excitatory and inhibitory processes in cardiac tissue by the acetylation modifiers may help explain their predominantly net-neutral effects on cardiac electrophysiology. ATP1A1, encoding for the Na/K pump, represented an outlier—with orthogonal regulation by the acetylation modifiers to most of the ICs. The HDAC/SIRT/HAT effects were mediated by strong (+) TF regulators of ICs, e.g., MEF2A and TBX5, in both sexes. Furthermore, for male hearts, PLS models revealed a stronger (+/-) mediatory role on ICs for NKX25 and TGF1B/KLF4, respectively, while RUNX1 exhibited larger (-) TF effects on ICs in females. Male-trained PLS models of HDAC/SIRT/HAT effects on ICs underestimated the effects on some ICs in females. Insights from the GTEx dataset about the co-expression and transcriptional co-regulation of acetylation-modifying enzymes, transcription factors and key cardiac ion channels in a sex-specific manner can help inform safer drug design.

Pathophysiological metabolic changes associated with disease modify the proarrhythmic risk profile of drugs with potential to prolong repolarisation

BACKGROUND AND PURPOSE

Hydroxychloroquine, chloroquine and azithromycin are three drugs that were proposed to treat COVID-19. While concern already existed around their proarrhythmic potential there is little data regarding how altered physiological states encountered in patients such as febrile state, electrolyte imbalances or acidosis might change their risk profiles.

EXPERIMENTAL APPROACH

Potency of hERG block was measured using high-throughput electrophysiology in the presence of variable environmental factors. These potencies informed simulations to predict population risk profiles. Effects on cardiac repolarisation were verified in human induced pluripotent stem cell-derived cardiomyocytes from multiple individuals.

KEY RESULTS

Chloroquine and hydroxychloroquine blocked hERG with IC₅₀ of 1.47±0.07 μM and 3.78±0.17 μM respectively, indicating proarrhythmic risk at concentrations effective against SARS-CoV-2 in vitro. Hypokalaemia and hypermagnesemia increased potency of chloroquine and hydroxychloroquine, indicating increased proarrhythmic risk. Acidosis significantly reduced potency of all drugs, whereas increased temperature decreased potency of chloroquine and hydroxychloroquine against hERG but increased potency for azithromycin. In silico simulations demonstrated that proarrhythmic risk was increased by female sex, hypokalaemia and heart failure, and identified specific genetic backgrounds associated with emergence of arrhythmia.

CONCLUSION AND IMPLICATIONS

Our study demonstrates how proarrhythmic risk can be exacerbated by metabolic changes and pre-existing disease. More broadly, the study acts as a blueprint for how high-throughput in vitro screening, combined with in silico simulations can help guide both preclinical screening and clinical management of patients in relation to drugs with potential to prolong repolarisation.

Computational models of atrial fibrillation: achievements, challenges, and perspectives for improving clinical care

Despite significant advances in its detection, understanding and management, atrial fibrillation (AF) remains a highly prevalent cardiac arrhythmia with a major impact on morbidity and mortality of millions of patients. AF results from complex, dynamic interactions between risk factors and comorbidities that induce diverse atrial remodelling processes. Atrial remodelling increases AF vulnerability and persistence, while promoting disease progression. The variability in presentation and wide range of mechanisms involved in initiation, maintenance and progression of AF, as well as its associated adverse outcomes, make the early identification of causal factors modifiable with therapeutic interventions challenging, likely contributing to suboptimal efficacy of current AF management. Computational modelling facilitates the multilevel integration of multiple datasets and offers new opportunities for mechanistic understanding, risk prediction and personalized therapy. Mathematical simulations of cardiac electrophysiology have been around for 60 years and are being increasingly used to improve our understanding of AF mechanisms and guide AF therapy. This narrative review focuses on the emerging and future applications of computational modelling in AF management. We summarize clinical challenges that may benefit from computational modelling, provide an overview of the different in silico approaches that are available together with their notable achievements, and discuss the major limitations that hinder the routine clinical application of these approaches. Finally, future perspectives are addressed. With the rapid progress in electronic technologies including computing, clinical applications of computational modelling are advancing rapidly. We expect that their application will progressively increase in prominence, especially if their added value can be demonstrated in clinical trials.

Ion current profiles in canine ventricular myocytes obtained by the "onion peeling" technique

The profiles of ion currents during the cardiac action potential can be visualized by the action potential voltage clamp technique. To obtain multiple ion current data from the same cell, the "onion peeling" technique, based on sequential pharmacological dissection of ion currents, has to be applied. Combination of the two methods allows recording of several ion current profiles from the same myocyte under largely physiological conditions. Using this approach, we have studied the densities and integrals of the major cardiac inward (I_Ca, I_NCX, I_Na-late) and outward (I_Kr, I_Ks, I_K1) currents in canine ventricular cells and studied the correlation between them. For this purpose, canine ventricular cardiomyocytes were chosen because their electrophysiological properties are similar to those of human ones. Significant positive correlation was observed between the density and integral of I_Ca and I_Kr, and positive correlation was found also between the integral of I_Ca and I_NCX. No further correlations were detected. The Ca²⁺-sensitivity of K⁺ currents was studied by comparing their parameters in the case of normal calcium homeostasis and following blockade of I_Ca. Out of the three K⁺ currents studied, only I_Ks was Ca²⁺-sensitive. The density and integral of I_Ks was significantly greater, while its time-to-peak value was shorter at normal Ca²⁺ cycling than following I_Ca blockade. No differences were detected for I_Kr or I_K1 in this regard. Present results indicate that the positive correlation between I_Ca and I_Kr prominently contribute to the balance between inward and outward fluxes during the action potential plateau in canine myocytes. The results also suggest that the profiles of cardiac ion currents have to be studied under physiological conditions, since their behavior may strongly be influenced by the intracellular Ca²⁺ homeostasis and the applied membrane potential protocol.

Network Pharmacology-Based Systematic Analysis of Molecular Mechanisms of Dingji Fumai Decoction for Ventricular Arrhythmia

Background

Dingji Fumai Decoction (DFD), a traditional herbal mixture, has been widely used to ventricular arrhythmia (VA) in clinical practice in China. However, research on the bioactive components and underlying mechanisms of DFD in VA is still scarce.

Methods

Components of DFD were collected from TCMSP, ETCM, and literature. The chemical structures of each component were obtained from PubChem. Next, SwissADME and SwissTargetPrediction were applied for compounds screening and targets prediction of DFD; meanwhile, targets of VA were collected from DrugBank and Online Mendelian Inheritance in Man (OMIM). Then, the H-C-T-D network and the protein-protein interaction (PPI) network were constructed based on the data obtained above. CytoNCA was utilized to filter hub genes and VarElect was used to analyze the relationship between genes and diseases. At last, Metascape was employed for systematic analysis on the potential targets of herbals against VA, and AutoDock was applied for molecular docking to verify the results.

Results

A total of 434 components were collected, 168 of which were qualified, and there were 28 shared targets between DFD and VA. Three function modules of DFD were found from the PPI network. Further systematic analysis of shared genes and function modules explained the potential mechanism of DFD in the treatment of VA; molecular docking has verified the interactions.

Conclusions

DFD could be employed for VA through mechanisms, including complex interactions between related components and targets, as predicted by network pharmacology and molecular docking. This work confirmed that DFD could apply to the treatment of VA and promoted the explanation of DFD for VA in the molecular mechanisms.