Screening and Identification of Cardioprotective Compounds From Wenxin Keli by Activity Index Approach and in vivo Zebrafish Model

Wenxin Keli (WXKL) is a widely used Chinese botanical drug for the treatment of arrhythmia, which is consisted of four herbs and amber. In the present study, we analyzed the chemical composition of WXKL using liquid chromatography coupled with high-resolution mass spectrometry (LC-HRMS) to tentatively identify 71 compounds. Through typical separate procession, the total extract of WXKL was divided into fractions for further bioassays. Cardiomyocytes and zebrafish larvae were applied for assessment. In vivo arrhythmia model in Cmlc2-GFP transgenic zebrafish was induced by terfenadine, which exhibited obvious reduction of heart rate and occurrence of atrioventricular block. Dynamic beating of heart was recorded by fluorescent microscope and sensitive camera to automatically recognize the rhythm of heartbeat in zebrafish larvae. By integrating the chemical information of WXKL and corresponding bioactivities of these fractions, activity index (AI) of each identified compound was calculated to screen potential active compounds. The results showed that dozens of compounds including ginsenoside Rg₁, ginsenoside Re, notoginsenoside R₁, lobetyolin, and lobetyolinin were contributed to cardioprotective effects of WXKL. The anti-arrhythmic activities of five compounds were further validated in larvae model and mature zebrafish by measuring electrocardiogram (ECG). Our findings provide a successful example for rapid discovery of bioactive compounds from traditional Chinese medicine (TCM) by activity index based approach coupled with in vivo zebrafish model.

Introduction

Natural products have played important roles in healthcare system throughout history and will continue to be served as huge and invaluable resource for the discovery of drug candidates. Traditional Chinese Medicine (TCM), widely used in Eastern Asian countries, has been regarded as an important part of natural products for the therapy of various diseases (Wang et al., 2012). The discovery of bioactive constituents from TCM is the key step in the modernization of TCM. As a consequence, developing high throughput methods with satisfied sensitivity for identifying active compounds from complex mixtures of TCM is in great demand. In past two decades, many efforts have been made for rapidly screening of compounds from TCM. Artemisinin (qinghaosu) with antimalarial effect is a successful and impressive example as the gift from TCM (Tu, 2011).

Arrhythmia occurs with abnormal beating of heart myocardium, generally represent disorder of ion channels or cardiomyopathy, can be classified with disorders of impulse formation or conduction. Arrhythmia has intricate pathogenesis, in general, most cardiovascular diseases such as heart failure always accompany with arrhythmia. Some types of arrhythmias are capable of triggering cardiac arrest and sudden death. Unfortunately, most of the antiarrhythmic drugs lack specificity and have numerous adverse effects (Page and Roden, 2005). TCMs have multitargets and synergy effect that benefit with those complex diseases (Li et al., 2011). Wenxin Keli (WXKL) is one of the widely used Chinese patent medicine for arrhythmia and heart failure, and is the first Chinese-developed anti-arrhythmic medicine approved by the China Food and Drug Administration (CFDA) and approved by Chinese Pharmacopeia (ChP), and consists of Codonopsis pilosula, Polygonatum sibiricum, Radix Notoginseng, Nardostachys jatamansi, Succinum. Researches in decades have proved that WXKL can suppress and prevent cardiac arrhythmias, including atrial and ventricular arrhythmias (Xing et al., 2013; He et al., 2016; Wang et al., 2016; Li et al., 2017), and inhibit multiple ion channels (Chen et al., 2013; Yang et al., 2017), especially the atrial-selective inhibition of sodium-channel current (Burashnikov et al., 2012), and effect on late Na current (Xue et al., 2013).

Arrhythmia is difficult to model in vitro, multiple elements may contribute to the final arrhythmia including genetic predisposition, extrinsic injury, environmental exposures, and stochastic processes. Cell-based model unable to fully reveal the pathological process of arrhythmia. On the other hand, using animal model for screening is costly. Zebrafish (Danio rerio) has been a model for biomedical research for decades and is ideal for phenotype-based screen in various ways. The generation time of adult fish is about 3 months and it is easy to maintain a great number of zebrafish with low cost and needn't much space. The embryogenesis can be finished in 24 h post-fertilization (hpf) and each pair of fish can produce more than hundred eggs and the mating is not depended on season. More attractively, their embryos are transparent and most organs including the heart, liver, intestine, and kidney develop in 96 hpf that can be clearly visualized (Barros et al., 2008). The larvae can be manipulated in well-plate and live with a little fluid and is permeable to small molecules (Zhang et al., 2003; Kari et al., 2007). These traits make it possible to establish an easy and high-throughput model to assay the effects of drug candidates on internal organs in the live organism. Zebrafish are easy to genetic manipulation to simulate human disease (Asnani and Peterson, 2014). As a developmental and genetic model, zebrafish has been used for anti-cancer compounds discovery, chemicals toxicity assessment, and so on. Zebrafish heart is highly comparable with human heart in structures, functions, signal pathways, and ion channels (Hu et al., 2000) and is particularly suitable for the study of the cardiovascular system. Here, we used terfenadine, an anti-histamine drug but also a potent hERG blocker and QT prolonger (Dhillon et al., 2013), and was reported that had pro-arrhythmic effects (Chaudhari et al., 2013), to induce the heart disturbance of zebrafish.

In the present study, we simultaneously used cell-based and zebrafish-based model to assess the cardioprotective and anti-arrhythmia effect of WXKL and the separated fractions, and combined HRMS and chemometric analysis to identify bioactive compounds from WXKL. H9c2 cell damaged by H₂O₂ were conducted to evaluate the protective activity of fractions. In vivo arrhythmic model based on Cmlc2-GFP transgenic (Tg) zebrafish was applied and the heart rate and rhythm of larvae were measured to evaluate pharmacological effects. The cell viability and heart rate recovery of zebrafish were transformed as the bioactivity coefficient and correlated with the compounds constitute in fractions from WXKL to calculate active index of every compounds. The entire process is illustrated in Figure 1. With both in vitro and in vivo assessment and active index calculation, several compounds including ginsenoside Rg₁, ginsenoside Re, notoginsenoside R₁, lobetyolin, lobetyolinin were selected to validate the activity on larvae and mature fish.

FIGURE 1

Figure 1. Scheme of active index calculation based active compounds screening from WXKL.

Materials and Methods

Materials and Reagents

Wenxin Keli was obtained from Shanxi Buchang Pharmaceutical Co., Ltd (Shanxi, China). Ginsenoside Rg₁, Ginsenoside Re, Notoginsenoside R₁ were purchased from Winherb Medical Tech. Co., Ltd (Shanghai, China). Lobetyolin was obtained from Push Bio-Tech (Chengdu, China). Extract of codonopsis glycosides was obtained from Dasfbio (Nanjing, China).

HPLC-grade acetonitrile and methanol were purchased from Merck (Darmstadt, Germany). Formic acid (HPLC grade) was purchased from Roe Scientific (Newark, DE, USA). Ethanol was purchased from Zhejiang Changqing Chemical (Hangzhou, China). Deionized water was prepared with an Elga PURELAB flex system (ELGA LabWater, UK).

High-glucose Dulbecco's modified Eagle's medium, fetal bovine serum, trypsin-EDTA and antibiotics (100 U/ml penicillin G and 100 g/mL streptomycin) were obtained from Gibico BRL (Grand Island, NY, USA). Tert-Butyl hydroperoxide solution, Thiazolyl Blue Tereazolium Bromide, Terfenadine, DMSO, N-Phenylthiourea (PTU), and Tricaine were acquired from Sigma–Aldrich (St. Louis, MO, USA).

Apparatus

Tecan infinite M1000 system (Tecan, Zurich, Switzerland). AB TripleTOF 5600^plus System (AB SCIEX, Framingham, USA), coupled to a Waters ACQUITY UPLC™ system (Waters, MA, USA). Finnigan LCQ DecaXP^plus mass spectrometer equipped with an ESI source (Thermo, MA, USA) coupled to Agilent 1100 liquid chromatography (Agilent, Waldbronn, Germany). Agilent 1200 preparative performance liquid chromatography (Agilent, Waldbronn, Germany). Leica DMI 3000B Fluorescence Inversion Microscope System (Leica Microsystems Inc., USA), Andor Zyla 5.5 sCMOS Cameras (Oxford Instruments plc, Tubney Woods, Abingdon, UK). IX-100F Zebrafish system (iWorx Systems, Inc., USA).

Characterization of Major Compounds of WXKL by LC–HRMS

Chemical composition of WXKL was characterized by AB TripleTOF 5600^plus System coupled to a Waters Acquity UPLC system. The MS conditions: scan range m/z 100–2,000. Negative ion mode: source voltage was−4.5 kV, and the source temperature was 550°C. Positive ion mode: source voltage was +5.5 kV, and the source temperature was 600°C. The pressure of Gas 1 (Air) and Gas 2 (Air) were set to 50 psi. The pressure of Curtain Gas (N₂) was set to 35 psi. Maximum allowed error was set to ±5 ppm. Declustering potential (DP), 100 V; collision energy (CE), 10 V. For MS/MS acquisition mode, the parameters were almost the same except that the collision energy (CE) was set at 50 ± 20 V, ion release delay (IRD) at 67, ion release width (IRW) at 25. The acquisition parameters for Finnigan LCQ DecaXP^plus mass spectrometer were as follows: nebulizing gas, high purity nitrogen (N₂); collusion gas, high-purity helium (He); ion spray voltage:−3 kV; capillary temperature: 350°C; capillary voltage:−15 V; mass range: m/z 100–1,500. Chromatographic separation was carried out on an Agilent Zorbax SB-C18 analytical column (4.6 × 250 mm I.D., 5 μm; Agilent Technologies, USA). The mobile phase consisted of water (A) and acetonitrile (B) both containing 0.05% v/v formic acid. A gradient program was used as follows: 0–5 min, 10% B; 5–15 min, 10–25% B; 15–35 min, 25–35% B; 35–40 min, 35–40% B; 40–45 min, 40–70% B; 45–55 min, 70–95% B; 55–65 min: 95% B. The flow rate was 0.5 mL/min, the column temperature was 30°C, and the injection volume was 20 μL.

Cell Culture and Anti-Oxidation Assays

H9c2 cell were obtained from Cell Bank of the Chinese Academy of Science (Shanghai, China) and cultured in high glucose Dulbecco's modified Eagle's medium supplemented containing 10% fetal bovine serum (FBS) and antibiotics (100 units/mL penicillin and 100 μg/mL streptomycin). The cultures were maintained at 37°C in a humidified atmosphere of 5% CO₂. The anti-oxidation activity of each fraction was determined by tetrazolium based colorimetric assay (MTT assay). Briefly, cells (5 × 10⁴ cell/mL) were seeded to 96-well plates for 24 h and then treated with fractions of WXKL for another 24 h prior to 150 μM H₂O₂ exposure in fresh medium for 3 h, After that, 100 μL 0.5 mg/mL MTT in fresh medium replaced the former medium for 4 h at 37°C. Then, the medium was replaced by 100 μL DMSO and vibrated for 10 min. The cell viabilities of tested fractions were determined by measuring the optimal densities (ODs) of untreated cells (control), the cells exposed to H₂O₂ (model), and the cells pre-incubated with components (tested). The activities of the components were calculated using the following formula: Survival rate% = OD of tested/OD of control. Protection rate% = (OD of model – OD of tested)/(OD of model – OD of control).

Zebrafish Husbandry and Management

Heterozygotes and homozygote transgenic Cmlc2-GFP zebrafish expressing green fluorescent protein (GFP) exclusively in myocardium were provided by Zebrafish Resource Center, Zhejiang University School of Medicine (Hangzhou, China) and maintained according to established standard procedures. Two parent zebrafish were placed separately in a mating box equipped with a separator to protect the eggs from being eaten. Spawning was induced in the morning and embryos from each box were collected and rinsed with system fish water (containing 0.3% Instant Ocean Salt in deionized water with final pH 6.9–7.2, conductivity 450–550 μs/cm, and hardness of about 90 mg/L NaHCO₃). The embryos were maintained in the Petri dish with system fish water and transferred to the incubator and incubated at 28°C. This study was granted by the Institutional Animal Care and Use Committee of the Laboratory Animal Center, Zhejiang University. We followed the relevant guidelines from the Laboratory Animal Center of Zhejiang University.

Zebrafish Arrhythmia Model and Drug Incubation

In 24 hpf, larvae with fluorescence were picked under fluorescent microscope and membranes of these larvae were ruptured artificially. Larvae were distributed into a 24-well plate and 8–10 larvae in each well with system fish water added with 0.2 mM N-Phenylthiourea (PTU) and 6 nM methylene blue for treatment. Set groups by wells, including Control, Model, and Treat. Terfenadine was stocked in DMSO at 100 mM, fractions of WXKL was stocked in DMSO at 100 mg/mL and fish water was used to dilute the stock to appropriate concentration. The model group was only given terfenadine, and the treat groups were given terfenadine and corresponding fractions. In 48 hpf incubating, the previous medium was discarded, and added fractions and terfenadine working solution, according to the groups, and filled to 2 mL with fish water medium in each well. The final concentrations of terfenadine was 6 μM. The fractions were diluted to appropriate concentration, mostly 50 μg/mL and some were 25, 12.5, 6.25 μg/mL, depending on the toxicity refer to cell assay.

Heartbeat Recording

In 72 hpf, the beating of zebrafish heart was recorded under florescence with Leica DMI 3000B Fluorescence Inversion Microscope System (Leica). The readout speed of the sCMOS camera was set at 10 frames or 20 frames per second with 4 × 4 pixel binning. L5 filter cube (excitation wave length of 480 nm and emission wave length of 527 nm). One hundred continuous dynamic images were captured by Zyla 5.5 sCMOS Cameras (Andor), subsequently were recognized by Matlab. The area of heart in each picture was measured. The area change with time was supposed to exhibit the heart rhythm. The heartbeats were also recorded manually for accuracy. We calculated the heartbeat of ventricle uniformly.

Calculation of the Activity Indexes

The recovery rate (R_i) of the components were calculated using the following formula:

$\begin{array}{l}{R}_i=\frac{B_i-B_M}{B_C-B_M}\times 100\%\end{array}$

R_i: normalized heart rate recovery rate of fraction i; B_i: beats of larvae treated by fraction i; B_M: beats of larvae treated by terfenadine; B_C: beats of larvae treated in control group.

The peak area of each compound was normalized according to the following formula:

$\begin{array}{l}{A}_j=\frac{A_{i,j}}{{\displaystyle {\sum }_{i=1}^m}{A}_{i,j}}\end{array}$

A_j: normalized values of peak area of constituent j in fraction i; A_{i, j}: peak area of compound j in fraction i; m: the numbers of fractions obtained from whole extract.

The activity indexes of compounds were given by the following formula:

$\begin{array}{l}A{I}_j=\sum (R_i{\times A}_j)\end{array}$

AI_j: activity index of compound j.

Zebrafish ECG Measurement

The ECG of zebrafish was measured by IX-100F Zebrafish system (iWorx Systems, Inc., USA). Zebrafish was anesthetic at first and positioned on its back on a fish-bed. Use a paper to gently remove excess water and ensure that the fins are not crossing the belly of the fish. Place the fish-bed with the fish, head to the right, in the chamber and position it under Ag/AgCl surface electrodes. The two electrodes were placed axially along the center-line of the fish's belly and the forward electrode should be placed close to the gills. ECG was recorded by LabScribe v3 software (iWorx System Inc., USA).

Statistical Analysis

The data are expressed as mean ± standard deviation (SD). Parameter comparisons between groups were made with one-way ANOVA analysis of variance. GraphPad prism 7 software (GraphPad Software, USA) was used to carry out statistical analysis. P < 0.05 was considered statistically significant.

Results and Discussion

The Chemical Composition of Wenxin Keli Extract

The main compounds of WXKL include sugar, glycosides, lignans, polyynes, saponins, iridoid glycosides, detailed information is listed in Table 1. The negative ion model base peak LC-MS chromatograms of WXKL was showed as Figure 2. We collected 71 compounds information of MS/MS and identified 53 compounds primarily, including saponins, phenylpropanoids, polyacetylene, triterpenoid, and others. Twenty-seven compounds of them belong to Notoginseng, including notoginsenoside R₁, ginsenoside Rg₁, ginsenoside Re, ginsenoside Ra₃, ginsenoside Rb₁, notoginsenoside R₂, ginsenoside Rc, ginsenoside Rd and so on. Sixteen compounds belong to Codonopsis, including tangshenoside V, lobetyolinin, lobetyolin, atractylenolide III, gentisic acid β-D-glucoside, syringin, hexyl 6-O-β-D-glucopyranosyl-β-D-glucopyranoside, hexyl 2-O-β-D-glucopyranosyl-beta-D-glucopyranoside and others. Besides, 5 were identified from Polygonatum, and 2 were from Nard (Table 1).

TABLE 1

Table 1. Chromatographic and mass spectral data of the constituents of WXKL.

FIGURE 2

Figure 2. Negative ion model HPLC–MS base peak chromatograms of WXKL.

Evaluating Cardioprotective Effect of Components by Zebrafish Arrhythmia Model

We first performed standard isolation by preparative chromatography to obtain fractions, which were analyzed by Finnigan LCQ DecaXP^plus mass spectrometer. The mass spectrums of every fractions were shown in Supplementary Material.

Oxidative stress plays a key role in the pathogenesis of various diseases (Furukawa et al., 2004). These fractions were then evaluated the protective activity on H9c2 cell damaged by H₂O₂. The toxicity of all fractions were tested at first, and the safe concentration of most fraction were 50 μg/mL, and some were 25/12.5/6.25 μg/mL to insure no toxicity. The cell viability with these fractions treated were shown in Supplementary Material.

Zebrafish (D. rerio) has been an ideal model for drug screening (Goldsmith, 2004). There have been many applications of zebrafish as a high-throughput screening model and cardiotoxicity risk assessment of drug candidates (Wen et al., 2012; Zhu et al., 2014). Here we used Cmlc2-GFP Tg zebrafish as the base. This transgenic line expressing GFP exclusively in myocardium driven by promoter cmlc2 (cardiac myosin light chain 2 gene) (Huang et al., 2003). Terfenadine has been reported can induce QT prolongation in zebrafish and guinea pig (Milan et al., 2006; Lu et al., 2012; Chaudhari et al., 2013), which is associated with ventricular tachyarrhythmia (Gowda et al., 2004). Terfenadine causes QT prolongation in adult zebrafish, also demonstrate in zebrafish embryos (Langheinrich et al., 2003).

As shown in Figure 3, terfenadine had less effect on the structure of the heart (Figure 3A) but influenced the rhythm of beat obviously (Figure 3B). The rhythm of heartbeat was exhibited through the area change of heart analyzed by Matlab. It's obvious that the control groups performed fast and regular rhythm (about 180 beats/min), and in the model groups, the heart rate were down to 80–100 beats/min with irregular heartbeats after incubated with terfenadine for 24 h, and some showed typically atrioventricular block (Peal et al., 2011), while co-incubated with WXKL stabilized the rhythm (Figure 3B, Supplementary Videos). And WXKL and its fractions performed varying affection on heart rate. Fraction 1 and Fraction 2 and former part of Fraction 3 showed beneficial effect on heart rate, while the remaining parts of Fraction 3 lowered the rate even more (Figure 3C).

FIGURE 3

Figure 3. Cardioprotection of Wenxin Keli and fractions on zebrafish model. (A) Represent images of heart structures in zebrafish under bright-field and fluorescence. (B) Represent images of heartbeat rhythm of different groups exhibited by Matlab. (C) Normalized heart rate of zebrafish larvae influenced by WXKL fractions. n ≥ 8. ^##P < 0.01 vs. control, *P < 0.05, **P < 0.01 vs. model.

Screening Active Compounds by Activity Indexes Calculation and Ranking

Activity index (AI) of each compound was calculated according to mathematical formulae proposed in our previous study. It was assumed that the compounds with positive activity index might be active and has contribution to the activity of whole formula to some extent (Wang et al., 2014). The relative intensities of the identified compounds in each fraction were visually presented in a heatmap (Figure 4A). After multiply corresponding bioactivity coefficient (i.e., heart rate recovery rate) of each fraction, the heatmap was converted into a bio-active map, and the red and gray color represent good or bad effect, respectively. The calculated scores were exhibited as histogram on the right (Figure 4B). The detailed scores were listed in Supplementary Material. We plot compounds with the effect on cardiomyocytes and the heart rate of zebrafish (Figure 4C), the compounds in upper right region represented a better activity.

FIGURE 4

Figure 4. Identification of active compounds by coefficient ranking. (A) Heatmap for relative content of identified compounds in each WXKL fractions. (B) Bio-active map was converted from content map and corresponding bioactivity coefficient in each fraction, and active indexes of identified compounds were calculated and showed on the right. (C) Compounds were plot by cell protection and bio-active scores and mass spectrum of several compounds on top. X axis means the contribution of compounds to the cardiomyocytes protection. Y axis means the AI of compounds.

Validation of Active Compounds

According to the scores, Ginsenoside Rg₁, Ginsenoside Re, Notoginsenoside R₁, Lobetyolin, Lobetyolinin were selected to validate activity, considering the available. Lobetyolinin was prepared and enriched by ourselves from commercial codonopsis glycosides. Their toxicity was confirmed before. We increased the dosage of terfenadine and shortened the incubation time as an acute injury model to improve significance when validating the active of pure compounds by reason of the pure compounds were not strong enough to exhibit activity in original method. After pre-treated with compounds (50 μM) for 24 h, the larvae were treated with 15 μM terfenadine for 2 h and recorded heartbeat under fluorescent. As the consequence, the heart rate of larvae was recovered in varying degree (Figure 5A). Ginsenoside Rg₁ and lobetyolinin exhibited better activities. Meanwhile, ECG of adult zebrafish treated with compound was measured. The heart rate of normal zebrafish was around 100 beats/min, and after treated with 25 μM terfenadine for 1 h, the heart rate was down and occurred irregular rhythm (Figure 5B). Lobetyolinin pre-treated for 6 h recovered the heart rate and represent electrocardiograms were showed as Figure 5C.

FIGURE 5

Figure 5. (A) Normalized heart rate of zebrafish larvae treated with compounds. n = 7. ^#P < 0.05 vs. control, ^##P < 0.01 vs. control, *P < 0.05, **P < 0.01 vs. model. (B) The schematic of ECG measurement, and the result of different group. n = 3. (C) Represent electrocardiograms of adult zebrafish treated in different groups.

Discussion

A few researches described the chemical components of WXKL. Wang et al. established a database for the chemical components of the five herbs in WXKL for active compounds predication (Wang et al., 2017), however, all the compounds were acquired from database refer to the herb not the real composition of the patent drug, and the chemical components probably change during the manufactory process. We analyzed the extract of WXKL directly with LC-HRMS at beginning, but the analytic method we established still has limitation. Actually, the compounds in Fraction 2 weren't separated clearly and seemed have low mass spectrum response, which make this portion of fractions have similar composition. Compounds identified from Nard and Succinum were rare, maybe for the reason that the main constitutes of Nard and Succinum are volatile oil, which are more appropriate analyzed by gas chromatography–mass spectrometry (GC–MS). According to report, the extracts of Nard significantly blocked I_Na and I_Ito of rat ventricular myocytes (Liu et al., 2009). The active compounds we predicted especially on the top are minor composition were difficult to get standard substance for bioactive assays except ginsenoside Re, ginsenoside Rg₁ and notoginsenoside R₁, which limited the further validation. We attempted to isolate substances such as lobetyolinin from extracts of C. pilosula. It has reported that ginsenoside Re has negative effect on cardiac contractility and autorhythmicity (Peng et al., 2012), ginsenoside Rg₁ prolonged ventricular refractoriness and repolarization (Wu et al., 1995), notoginsenoside R₁ has protective effects on cardiovascular system (Li et al., 2014). Related activity of lobetyolin has few reports.

Drug-induced model is a common approach, verapamil and terfenadine were applied to develop a zebrafish heart failure model (Zhu et al., 2018). QT prolonging is a typical characteristic of arrhythmia which also can be induced by cisapride and astemizole besides terfenadine (Langheinrich et al., 2003). It has to be considered that drug treated by oral may cause unstable effect, so design rational approaches of drug treatment is necessary. Chaudhari et al. performed parenteral administration of terfenadine with different doses and recorded ECG to assess drug-induced QTc prolongation in zebrafish. Those with the doses above 1 mg/kg were observed some proarrhythmic effects such as Ventricular Premature Contractions, Ventricular Tachycardia, Atrio-Ventricular (AV) Block, and Torsade de pointes (TdP) (Chaudhari et al., 2013). However, the cardiotoxic of terfenadine is possible not associated with QT prolongation and the occurrence of TdPs, but with marked widening of the QRS complex and other cardiac arrhythmias (Hondeghem et al., 2011). It's reported that terfenadine caused non-TdP like VT/VF by slowing of conduction via blockade of I_Na (Lu et al., 2012). In addition, transgenic zebrafish lines are also feasible to avoid the unstable results of drug induce method. Several mutants were identified that exhibited arrhythmias (Milan and Macrae, 2008) like the bradycardic line slo mo, with variable degrees of sinoatrial or atrioventricular heart block (Baker et al., 1997). Besides, some mutants exhibited recessive lethal phenotypes included mutants such as tremblor (Langenbacher et al., 2005), island beat (Rottbauer et al., 2001), and reggae (Hassel et al., 2008). The Tg line that express fluorescent proteins (e.g., Cmlc2-GFP) was beneficial for optical measurement from the other perspective.

ECG (electrocardiogram) abnormalities is the critical characteristic of arrhythmia. Several ECG measurement equipments for zebrafish have been developed, and mostly of them consist of electrode or micropipette and electrical filter (Milan et al., 2006; Chaudhari et al., 2013; Dhillon et al., 2013). Detection of 3 dpf larva is also possible (Chi et al., 2008). However, the invasive injury and anesthesia could cause damage to the individual. As for high throughput screening, a fast and stable method of ECG measurement is required, but existing devices seem not compatibility. Computerized recognition of ECGs has become a well-established practice, assisting to classify long-term ECG recordings, which suggests new approaches like Machine Learning are able to recognize and classify the rhythm signal. For instance, automatic classification of single-lead ECG signals with Deep learning (also known as unsupervised feature learning or representation learning) was established (Singh et al., 2018). A new semi-supervised approach based on deep learning and active learning for classification of electrocardiogram signals is proposed (Sayantan et al., 2018). Though several algorithms have focused on automatically classifying heartbeats in ECGs, the scalability failure to handle large intra-class variations wherein the robustness of many existing ECG classification techniques remains limited. We have acquired plenty of dynamic images of different conditions of heartbeat and attempt to establish the relationship between waveform and define characteristic to classify the phenotype. Rhythm classify method based on image processing will be a non-invasive measurement of heart regulation.

In conclusion, we identified 71 compounds from extract of WXKL by LC-HRMS, and firstly utilized a transgenic zebrafish cmlc2-GFP induced by terfenadine as an animal model for screening active compounds from WXKL. After recording heartbeat that affected by fractions under a fluorescent microscope, a convenient image process was applied to exhibit the rhythm of heartbeat. Subsequently, we integrated chemometric analysis with bio-activity in vivo model of corresponding fractions and calculated active index of identified compounds. Ginsenoside Rg₁, ginsenoside Re, notoginsenoside R₁, lobetyolin, and lobetyolinin were selected to validate activity. Measurement of ECG of adult zebrafish also performed as a complement. Our results suggest that integrate bio-assay and substantial analysis to perform active index calculation improve the efficiency of active compounds discovering from TCM, and this approach is possible to be applied for the research of complex diseases.

Author Contributions

HL and XC designed and performed the experimental work. BZ, KQ, and YS provided the WXKL patent drug and related herb and extract. MB guided the theory of cardiac electrophysiology and pharmacology. All authors proofread the paper and provided feedback.

Funding

This study was supported by the National Key Scientific and Technological Project of China (grant 2017ZX09301012), and National Natural Science Foundation of China (No. 81774151, No. 81822047), the Fundamental Research Funds for the Central Universities (2016FZA7016).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.01288/full#supplementary-material

Supplementary Table 1. Bioactive-coefficient ranking of identified compounds.

Supplementary Figure 1. LC-MS chromatograms of fractions of WXKL. Some fractions were not enough for analysis so were not showed here.

Supplementary Figure 2. Anti-oxidation activity of WXKL fractions. (A) Protection rate of all fractions of WXKL, n = 3. (B) Dose-response of selected fractions that has protective activity.

Supplementary Figure 3. Dose effect of Terfenadine on heart rate of zebrafish. **P < 0.01 vs. control. n = 8.

Video 1. Side view of heart of zebrafish larva in normal group.

Video 2. Vertical view of heart of zebrafish larva in normal group.

Video 3. Side view of heart of zebrafish larva in terfenadine group.

Video 4. Vertical view of heart of zebrafish larva in terfenadine group.

Video 5. Side view of heart of zebrafish larva in WXKL group.

Video 6. Vertical view of heart of zebrafish larva in WXKL group.

References

Asnani, A., and Peterson, R. T. (2014). The zebrafish as a tool to identify novel therapies for human cardiovascular disease. Dis. Model. Mech. 7, 763–767. doi: 10.1242/dmm.016170

PubMed Abstract | CrossRef Full Text | Google Scholar

Baker, K., Warren, K. S., Yellen, G., and Fishman, M. C. (1997). Defective “pacemaker” current (Ih) in a zebrafish mutant with a slow heart rate. Proc. Natl. Acad. Sci. U.S.A. 94, 4554–4559. doi: 10.1073/pnas.94.9.4554

PubMed Abstract | CrossRef Full Text | Google Scholar

Barros, T. P., Alderton, W. K., Reynolds, H. M., Roach, A. G., and Berghmans, S. (2008). Zebrafish: an emerging technology for in vivo pharmacological assessment to identify potential safety liabilities in early drug discovery. Br. J. Pharmacol. 154, 1400–1413. doi: 10.1038/bjp.2008.249

PubMed Abstract | CrossRef Full Text | Google Scholar

Burashnikov, A., Petroski, A., Hu, D., Barajas-Martinez, H., and Antzelevitch, C. (2012). Atrial-selective inhibition of sodium-channel current by Wenxin Keli is effective in suppressing atrial fibrillation. Heart Rhythm 9, 125–131. doi: 10.1016/j.hrthm.2011.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaudhari, G. H., Chennubhotla, K. S., Chatti, K., and Kulkarni, P. (2013). Optimization of the adult zebrafish ECG method for assessment of drug-induced QTc prolongation. J. Pharmacol. Toxicol. Methods 67, 115–120. doi: 10.1016/j.vascn.2013.01.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y., Li, Y., Guo, L., Chen, W., Zhao, M., Gao, Y., et al. (2013). Effects of wenxin keli on the action potential and L-type calcium current in rats with transverse aortic constriction-induced heart failure. Evid. Based Complement. Alternat. Med. 2013:572078. doi: 10.1155/2013/572078

PubMed Abstract | CrossRef Full Text | Google Scholar

Chi, N. C., Shaw, R. M., Jungblut, B., Huisken, J., Ferrer, T., Arnaout, R., et al. (2008). Genetic and physiologic dissection of the vertebrate cardiac conduction system. PLoS Biol. 6:e109. doi: 10.1371/journal.pbio.0060109

PubMed Abstract | CrossRef Full Text | Google Scholar

Dhillon, S. S., Doro, E., Magyary, I., Egginton, S., Sik, A., and Muller, F. (2013). Optimisation of embryonic and larval ECG measurement in zebrafish for quantifying the effect of QT prolonging drugs. PLoS ONE 8:e60552. doi: 10.1371/journal.pone.0060552

PubMed Abstract | CrossRef Full Text | Google Scholar

Furukawa, S., Fujita, T., Shimabukuro, M., Iwaki, M., Yamada, Y., Nakajima, Y., et al. (2004). Increased oxidative stress in obesity and its impact on metabolic syndrome. J. Clin. Invest. 114, 1752–1761. doi: 10.1172/J.C.I.21625

PubMed Abstract | CrossRef Full Text | Google Scholar

Goldsmith, P. (2004). Zebrafish as a pharmacological tool: the how, why and when. Curr. Opin. Pharmacol. 4, 504–512. doi: 10.1016/j.coph.2004.04.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Gowda, R. M., Khan, I. A., Wilbur, S. L., Vasavada, B. C., and Sacchi, T. J. (2004). Torsade de pointes: the clinical considerations. Int. J. Cardiol. 96, 1–6. doi: 10.1016/j.ijcard.2003.04.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Hassel, D., Scholz, E. P., Trano, N., Friedrich, O., Just, S., Meder, B., et al. (2008). Deficient zebrafish ether-a-go-go-related gene channel gating causes short-QT syndrome in zebrafish reggae mutants. Circulation 117, 866–875. doi: 10.1161/CIRCULATIONAHA.107.752220

PubMed Abstract | CrossRef Full Text

He, M., Lv, Z., Yang, Z. W., Huang, J. L., and Liu, F. (2016). Efficacy and safety of Chinese herbal medicine Wenxin Keli for ventricular premature be ats: a systematic review. Complement. Ther. Med. 29, 181–189. doi: 10.1016/j.ctim.2016.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Hondeghem, L. M., Dujardin, K., Hoffmann, P., Dumotier, B., and De Clerck, F. (2011). Drug-induced QTC prolongation dangerously underestimates proarrhythmic potential: lessons from terfenadine. J. Cardiovasc. Pharmacol. 57, 589–597. doi: 10.1097/FJC.0b013e3182135e91

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, N., Sedmera, D., Yost, H. J., and Clark, E. B. (2000). Structure and function of the developing zebrafish heart. Anat. Rec. 260, 148-157. doi: 10.1002/1097-0185(20001001)260:2<148::AID-AR50>3.0.CO;2-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, C. J., Tu, C. T., Hsiao, C. D., Hsieh, F. J., and Tsai, H. J. (2003). Germ-line transmission of a myocardium-specific GFP transgene reveals critical regulatory elements in the cardiac myosin light chain 2 promoter of zebrafish. Dev. Dyn. 228, 30–40. doi: 10.1002/dvdy.10356

PubMed Abstract | CrossRef Full Text | Google Scholar

Kari, G., Rodeck, U., and Dicker, A. P. (2007). Zebrafish: an emerging model system for human disease and drug discovery. Clin. Pharmacol. Ther. 82, 70–80. doi: 10.1038/sj.clpt.6100223

PubMed Abstract | CrossRef Full Text | Google Scholar

Langenbacher, A. D., Dong, Y., Shu, X., Choi, J., Nicoll, D. A., Goldhaber, J. I., et al. (2005). Mutation in sodium-calcium exchanger 1 (NCX1) causes cardiac fibrillation in zebrafish. Proc. Natl. Acad. Sci. U.S.A. 102, 17699–17704. doi: 10.1073/pnas.0502679102

PubMed Abstract | CrossRef Full Text | Google Scholar

Langheinrich, U., Vacun, G., and Wagner, T. (2003). Zebrafish embryos express an orthologue of HERG and are sensitive toward a range of QT-prolonging drugs inducing severe arrhythmia. Toxicol. Appl. Pharmacol. 193, 370–382. doi: 10.1016/j.taap.2003.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Li, Q., Liu, Y. Y., Wang, M. X., Pan, C. S., Yan, L., et al. (2014). Protective effects of Notoginsenoside R1 on intestinal ischemia-reperfusion injury in rats. Am. J. Physiol. Gastrointest. Liver Physiol. 306, G111–G122. doi: 10.1152/ajpgi.00123.2013

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, M., Qiu, R., Tian, G., Zhang, X., Li, C., Chen, S., et al. (2017). Wenxin Keli for ventricular premature complexes with heart failure: a systematic review and meta-analysis of randomized clinical trials. Complement. Ther. Med. 33, 85–93. doi: 10.1016/j.ctim.2017.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Zhang, B., and Zhang, N. (2011). Network target for screening synergistic drug combinations with application to traditional Chinese medicine. BMC Syst. Biol. 5(Suppl. 1):S10. doi: 10.1186/1752-0509-5-S1-S10

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Guo, J., Zhang, P., Li, J., and Li, C. (2009). The effects of nardostachys chinensis batal extract on the sodium current and transient outward potassium current of rat ventricular myocytes. Chinese J. Card. Pacing Electrophysiol. 23, 533–535. doi: 10.13333/j.cnki.cjcpe.2009.06.002

CrossRef Full Text

Lu, H. R., Hermans, A. N., and Gallacher, D. J. (2012). Does terfenadine-induced ventricular tachycardia/fibrillation directly relate to its QT prolongation and Torsades de Pointes? Br. J. Pharmacol. 166, 1490–1502. doi: 10.1111/j.1476-5381.2012.01880.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Milan, D. J., Jones, I. L., Ellinor, P. T., and Macrae, C. A. (2006). In vivo recording of adult zebrafish electrocardiogram and assessment of drug-induced QT prolongation. Am. J. Physiol. Heart Circ. Physiol. 291, H269–H273. doi: 10.1152/ajpheart.00960.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Milan, D. J., and Macrae, C. A. (2008). Zebrafish genetic models for arrhythmia. Prog. Biophys. Mol. Biol. 98, 301–308. doi: 10.1016/j.pbiomolbio.2009.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Page, R. L., and Roden, D. M. (2005). Drug therapy for atrial fibrillation: where do we go from here? Nat. Rev. Drug Discov. 4, 899–910. doi: 10.1038/nrd1876

PubMed Abstract | CrossRef Full Text | Google Scholar

Peal, D. S., Mills, R. W., Lynch, S. N., Mosley, J. M., Lim, E., Ellinor, P. T., et al. (2011). Novel chemical suppressors of long QT syndrome identified by an in vivo functional screen. Circulation 123, 23–30. doi: 10.1161/CIRCULATIONAHA.110.003731

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, L., Sun, S., Xie, L. H., Wicks, S. M., and Xie, J. T. (2012). Ginsenoside Re: pharmacological effects on cardiovascular system. Cardiovasc. Ther. 30, e183–188. doi: 10.1111/j.1755-5922.2011.00271.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Rottbauer, W., Baker, K., Wo, Z. G., Mohideen, M. A., Cantiello, H. F., and Fishman, M. C. (2001). Growth and function of the embryonic heart depend upon the cardiac-specific L-type calcium channel alpha1 subunit. Dev. Cell 1, 265–275. doi: 10.1016/S1534-5807(01)00023-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sayantan, G., Kien, P. T., and Kadambari, K. V. (2018). Classification of ECG beats using deep belief network and active learning. Med. Biol. Eng. Comput. 56, 1887–1898. doi: 10.1007/s11517-018-1815-2

CrossRef Full Text | Google Scholar

Singh, R. S., Saini, B. S., and Sunkaria, R. K. (2018). Detection of coronary artery disease by reduced features and extreme learning machine. Clujul Med 91, 166–175. doi: 10.15386/cjmed-882

PubMed Abstract | CrossRef Full Text | Google Scholar

Tu, Y. (2011). The discovery of artemisinin (qinghaosu) and gifts from Chinese medicine. Nat. Med. 17, 1217–1220. doi: 10.1038/nm.2471

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Chen, P., Jiang, W., Wu, L., Chen, L., Fan, X., et al. (2014). Identification of the effective constituents for anti-inflammatory activity of Ju-Zhi-Jiang-Tang, an ancient traditional Chinese medicine formula. J. Chromatogr. A 1348, 105–124. doi: 10.1016/j.chroma.2014.04.084

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, T., Lu, M., Du, Q., Yao, X., Zhang, P., Chen, X., et al. (2017). An integrated anti-arrhythmic target network of a Chinese medicine compound, Wenxin Keli, revealed by combined machine learning and molecular pathway analysis. Mol. Biosyst. 13, 1018–1030. doi: 10.1039/C7MB00003K

CrossRef Full Text | Google Scholar

Wang, X., Wang, Y., Feng, X., Lu, Y., Zhang, Y., Wang, W., et al. (2016). Systematic review and meta-analysis of randomized controlled trials on Wenxin keli. Drug Des. Devel. Ther. 10, 3725–3736. doi: 10.2147/DDDT.S112333

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Fan, X., Qu, H., Gao, X., and Cheng, Y. (2012). Strategies and techniques for multi-component drug design from medicinal herbs and traditional Chinese medicine. Curr. Top. Med. Chem. 12, 1356–1362. doi: 10.2174/156802612801319034

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, D., Liu, A., Chen, F., Yang, J., and Dai, R. (2012). Validation of visualized transgenic zebrafish as a high throughput model to assay bradycardia related cardio toxicity risk candidates. J. Appl. Toxicol. 32, 834–842. doi: 10.1002/jat.2755

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, W., Zhang, X. M., Liu, P. M., Li, J. M., and Wang, J. F. (1995). Effects of Panax notoginseng saponin Rg1 on cardiac electrophysiological properties and ventricular fibrillation threshold in dogs. Zhongguo Yao Li Xue Bao 16, 459–463.

PubMed Abstract | Google Scholar

Xing, Y., Gao, Y., Chen, J., Zhu, H., Wu, A., Yang, Q., et al. (2013). Wenxin-Keli regulates the calcium/calmodulin-dependent protein kinase II signal transduction pathway and inhibits cardiac arrhythmia in rats with myocardial infarction. Evid. Based Complement. Alternat. Med. 2013:464508. doi: 10.1155/2013/464508

PubMed Abstract | CrossRef Full Text | Google Scholar

Xue, X., Guo, D., Sun, H., Wang, D., Li, J., Liu, T., et al. (2013). Wenxin Keli suppresses ventricular triggered arrhythmias via selective inhibition of late sodium current. Pacing Clin. Electrophysiol. 36, 732–740. doi: 10.1111/pace.12109

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, X., Chen, Y., Li, Y., Ren, X., Xing, Y., and Shang, H. (2017). Effects of Wenxin Keli on cardiac hypertrophy and arrhythmia via regulation of the calcium/calmodulin dependent kinase II signaling pathway. Biomed Res. Int. 2017:1569235. doi: 10.1155/2017/1569235

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Willett, C., and Fremgen, T. (2003). Zebrafish: an animal model for toxicological studies. Curr. Protoc. Toxicol. Chapter 1, Unit 1.7. doi: 10.1002/0471140856.tx0107s17

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, J. J., Xu, Y. Q., He, J. H., Yu, H. P., Huang, C. J., Gao, J. M., et al. (2014). Human cardiotoxic drugs delivered by soaking and microinjection induce cardiovascular toxicity in zebrafish. J. Appl. Toxicol. 34, 139–148. doi: 10.1002/jat.2843

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, X. Y., Wu, S. Q., Guo, S. Y., Yang, H., Xia, B., Li, P., et al. (2018). A zebrafish heart failure model for assessing therapeutic agents. Zebrafish 15, 243–253. doi: 10.1089/zeb.2017.1546

PubMed Abstract | CrossRef Full Text | Google Scholar