- School of Chinese Materia Medica, Beijing University of Chinese Medicine, Beijing, China
Background: In traditional Mongolian or Tibetan medicine in China, Chebulae Fructus (CF) is widely used to process or combine with aconitums to decrease the severe toxicity of aconitums. Researches in this area have predominantly focused on tannins, with few research on other major CF components for cardiotoxicity mitigation. The present study aimed to clarify whether triterpenoids can attenuate the cardiotoxicity caused by mesaconitine (MA) and investigate the mechanism of cardiotoxicity attenuation.
Methods: Firstly, the pharmacophore model, molecular docking, and 3D-QSAR model were used to explore the mechanism of CF components in reducing the toxicity of MA mediated by the TRPV1 channel. Then three triterpenoids were selected to verify whether the triterpenoids had the effect of lowering the cardiotoxicity of MA using H9c2 cells combined with MTT, Hoechst 33258, and JC-1. Finally, Western blot, Fluo-3AM, and MTT assays combined with capsazepine were used to verify whether the triterpenoids reduced H9c2 cardiomyocyte toxicity induced by MA was related to the TRPV1 channel.
Results: Seven triterpenoids in CF have the potential to activate the TRPV1 channel. And they exhibited greater affinity for TRPV1 compared to other compounds and MA. However, their activity was relatively lower than that of MA. Cell experiments revealed that MA significantly reduced H9c2 cell viability, resulting in diminished mitochondrial membrane potential and nuclear pyknosis and damage. In contrast, the triterpenoids could improve the survival rate significantly and counteract the damage of MA to the cells. We found that MA, arjungenin (AR), and maslinic acid (MSA) except corosolic acid (CRA) upregulated the expression of TRPV1 protein. MA induced a significant influx of calcium, whereas all three triterpenoids alleviated this trend. Blocking the TRPV1 channel with capsazepine only increased the cell viability that had been simultaneously treated with MA, and AR, or MSA. However, there was no significant difference in the CRA groups treated with or without capsazepine.
Conclusion: The triterpenoids in CF can reduce the cardiotoxicity caused by MA. The MSA and AR function as TRPV1 agonists with comparatively reduced activity but a greater capacity to bind to TRPV1 receptors, thus antagonizing the excessive activation of TRPV1 by MA.
1 Introduction
According to traditional Mongolian medical theory, Chebulae Fructus (the fruits of Terminalia chebula Retz. or Terminalia chebula Retz, var. Tomentella Kurt., named Hezi in Chinese, CF) has a detoxifying function and may not only treat a variety of toxic diseases but also reconcile the drug properties and lessen the toxicity of poisonous medications (Nigam et al., 2020; Zhou et al., 2022; Ai et al., 2023). According to the records in the Tibetan medical classic “Jingzhu Materia Medica,” CF has a therapeutic effect on various ailments and is the king of all remedies (Dimar, 1986). It can be inferred from the phenomenon that among the commonly used 200 prescriptions in the Ministry standards of Tibetan medicine, 125 prescriptions all contain the CF (Song et al., 2023). In the use habits of ethnic minorities, it is more commonly used together with the aconitum herbs, such as Aconitum pendulum Busch, Aconiti kusnezoffii radix. In order to reduce the toxicity or preservation effect, CF decoction processing aconitums or compatibility aconitums with CF are popular methods in clinical practice (Li et al., 2019; Li et al., 2021), and they are often used in the form of pills or powders, such as Naru Sanwei Wan (Guo et al., 2023). However, the mechanism behind the attenuation toxicity remains unclear and requires further exploration.
Diester diterpenoid alkaloids (DDAs), represented by mesaconitine (C33H45NO11, CAS: 2752-64-9, MA), are the main toxic components of aconitums (Wei et al., 2019; Chen et al., 2023). The toxicity of MA is quite potent; the median lethal dose (LD50) of oral administration, subcutaneous injection, intraperitoneal injection, and intravenous injection in mice is 1.9 mg/kg, 0.20 mg/kg ∼ 0.38 mg/kg, 0.20 mg/kg ∼ 0.30 mg/kg, 0.068 mg/kg ∼ 0.13 mg/kg, respectively (Liu et al., 2016b; Yang et al., 2018; Chan et al., 2021). Cardiotoxicity is the main manifestation of aconitum herbs toxicity. It can lead to abnormal electrocardiogram changes, palpitations and arrhythmia, especially premature ventricular beats, which are the most common (Feng, 2019; Ye et al., 2021; Chen et al., 2023). For cells, it can lead to myocardial cell apoptosis, which is manifested in reducing the survival rate of rat H9c2 cardiomyocytes, causing cell nuclear condensation, damage and an increase in mitochondrial membrane potential (Li et al., 2020; Han et al., 2022c). There are multiple mechanisms of cardiotoxicity, including ion channels, protein expression, oxidative stress, energy metabolism damage, and so on (Sun et al., 2014; Zhou et al., 2021; Gao et al., 2022). Among them, ion channels are an important way to exert the efficacy and toxicity of diester alkaloids. When the function of ion channels is aberrant, the proportion of each channel in the myocardium will be destroyed, which will induce cardiac diseases such as arrhythmia (Fu et al., 2007; Yang et al., 2021).
The transient receptor potential vanilloid type 1 channel (TRPV1) is widely distributed in the heart (Peng and Li, 2010). It is a non-selective ion channel that dominates the influx of calcium ions. Activating the TRPV1 channel results in an influx of extracellular calcium and the release of intracellular calcium, causing an increase in intracellular calcium concentration (Sappington et al., 2009). Since calcium serves as the second messenger in many intracellular signaling pathways, TRPV1 can thus participate in various physiological processes of cells, including energy metabolism, cell proliferation, and apoptosis (Wu et al., 2015; Zhai et al., 2020). Researchers have discovered that TRPV1 channels are present in rat H9c2 cardiomyocytes and can be triggered in response to hypoxia/reoxygenation stimulation. The TRPV1 agonist capsaicin can dramatically worsen apoptosis, raise intracellular Ca2+ levels, increase mitochondrial superoxide release, and then lower mitochondrial membrane potential when mixed with it. However, the TRPV1 antagonist capsazepine can reverse this effect because it is a competitive small-molecule TRPV1 inhibitor and can block the activation of the TRPV1 channels by capsaicin or other agonists. Overactivation of TRPV1 channels in the heart induces apoptosis of cardiomyocytes, facilitating their death and contributing to ventricular tachycardia and fibrillation in individuals with myocardial infarction (Sun et al., 2014). This is a significant mechanism of cardiotoxicity and a notable pathway of cardiotoxicity as well. The TRPV1 channel is a significant pathway of cardiotoxicity induced by DDAs. The previous study found that the DDAs in aconitums caused intracellular calcium overload due to excessive activation of the TRPV1 channel in H9c2 cardiomyocytes. This led to damage to various organelles, resulting in a significant amount of apoptosis, which is one of the ways aconitum herbs produce cardiotoxicity (Han et al., 2022b; Han et al., 2022c).
CF is rich in tannins such as gallic acid and chebulic acid, as well as triterpenoids such as arjunolic acid, maslinic acid (MSA), corosolic acid (CRA), and chebulatriol. It also contains flavonoids, amino acids, and volatile compounds (Li et al., 2022). The components in CF might be involved in lowering toxicity, regardless of whether CF is compatibile with aconitum herbs or used as an excipient to process aconitums (Liu et al., 2016a; Zhi et al., 2020a). Academics have conducted studies and found that the chemical constituents interactions between aconitums and the auxiliary substance CF can lower the cardiotoxicity of aconitums (Liu et al., 2013). Network pharmacology was employed to investigate how CF reduces the cardiac toxicity of aconitums. A total of 15 components in CF were screened, including seven tannins, five triterpenoids, and three alkaloids. They reduce cardiotoxicity by involving numerous protein functions, signaling pathways, and biological processes (Li et al., 2018). Our research group has verified that tannins and phenolic acids from CF can significantly reduce the mitochondrial membrane potential elevation and intracellular Ca2+ overload caused by aconitine and MA in rat H9c2 cardiomyocytes through the TRPV1 channel, thereby improving mitochondrial function and inhibiting the cardiotoxicity of aconitine (Han et al., 2020; Han et al., 2022b; Han et al., 2022c). When we were committed to the study of tannins, we accidentally found that triterpenoids in CF, while not as abundant as tannins, have a significant impact on antioxidant and cardioprotective properties. These effects were achieved through antioxidation, inhibition of intracellular calcium overload in cardiomyocytes, and prevention of oxidative stress damage in cardiomyocytes (Sumitra et al., 2001; Zou et al., 2017; Wang, 2019; An et al., 2020). Coincidentally, this aligns with the mechanism by which DDAs in aconitums exert cardiotoxicity. Furthermore, research has shown that the triterpenoids in CF have the potential to activate TRPV1, which may play a role in attenuating cardiotoxicity (Mi et al., 2023). Consequently, further research is vital to identify whether the triterpenoids present in CF possess the ability to alleviate the cardiotoxicity of MA and if the attenuation mechanism involves the TRPV1 channel.
In the present study, the pharmacophore model based on the common feature structure of ligands (TRPV1 agonists) was used to search for drug molecules that are identical or similar to the pharmacophore model screened in the database (Acharya et al., 2011; Qiao and Zhang, 2014). Then the binding ability of biological macromolecules and small molecules was evaluated according to the scoring function of molecular docking, and the molecules with reasonable binding modes and high prediction scores were selected (Yang, 2010). Subsequently, the 3D-QSAR model was used to investigate the structure-activity relationship and perform activity prediction afterward (Hansch et al., 1962), thereby screening out the components of CF that attenuate the cardiotoxicity of MA through the TRPV1 channel. The attenuation effect was then verified by cell experiments, and the chosen attenuation concentration was used to further verify whether the effect was related to the TRPV1 channel. This provided experimental support for the theoretical underpinnings and molecular mechanism of CF in reducing cardiotoxicity induced by aconitums.
2 Materials and methods
2.1 Virtual screening based on pharmacophore model
2.1.1 Preparation of ligands
The pharmacophore model is particularly useful when the receptor structure is unknown or the mechanism of action is unclear. Based on the common features of a series of small molecule ligands, the pharmacophore model is generated, and then the compound database is searched for active compounds. Pharmacophore model has become an important tool in drug discovery, particularly for fast and effective screening of new bioactive molecules.
By searching the small molecule database PubChem (https://pubchem.ncbi.nlm.nih.gov/), TCMSP (http://tcmspw.com/tcmsp.php), TCMID (http://www.Megabionet.org/tcmid/) and relevant literature, we found 60 components in CF (Supplementary Table S1). And then, we searched for 10 TRPV1 agonists and 10 TRPV1 antagonists from the literature (Supplementary Table S2) and downloaded small molecule ligands (Gavva, 2008; Gunthorpe and Szallasi, 2008; Burgess and Williams, 2010; Li et al., 2011; Szolcsányi and Sándor, 2012; Sugimoto et al., 2013; Yin, 2019; Choi et al., 2020; Duarte et al., 2020; Mazeto et al., 2020; More et al., 2020; Nie, 2020). We also had to randomly download 90 compounds from the PubChem database that did not have well-defined TRPV1 agonist activity for later verification. All small molecular structures were reconstructed with energy minimization in the MMFF94 forcefield through Chem 3D (version 15.1.0) software. The obtained conformations were used for subsequent analysis.
2.1.2 Generation of pharmacophore model
Firstly, seven representative components were selected from the downloaded components with TRPV1 agonists to form a training set (Figure 1), including piperine, capsaicin, 6-gingerol, 6-shogaol, zingerone, resiniferatoxin (RTX), and DA-5018. And then the HipHop approach of the Common Feature Pharmacophore Generation module in Discovery Studio (version 4.0) was used to automatically generate the pharmacophore models based on the common structural features of the small molecules in the training set. The Principal value of the active compounds was set to 2, indicating that the small molecules are active, and the MaxOmitFeat value was set to 1, indicating that the number of characteristic elements in each molecule that allow not match to the pharmacophore model is 1. Then the Hydrophobe, Donor, Acceptor, Ionizable Positive and Ring Aromatic were selected as the characteristic elements of pharmacophore in the HipHop module. The best mode was adopted for superposition, and the Maximum Conformation was set to 255. Only the model whose energy difference with the lowest conformation was less than 20 kcal/mol was saved. Eventually, we would get the 10 pharmacophore models with the higher score.
2.1.3 Pharmacophore model validation
Using the Ligand Profiler module in Discovery Studio, the 10 generated pharmacophores were evaluated. The molecules in the test set consisted of 10 TRPV1 agonists (Supplementary Table S2) and 90 components with no reported TRPV1 agonist activity. According to the Rank value and the matching degree of the test set molecules, one pharmacophore was selected for subsequent screening, and the final selected pharmacophore was characterized by a high score, a high matching degree with the active ingredient, and a low matching degree with the inactive ingredient.
2.1.4 Virtual screening
The Search 3D Database module in Discovery Studio was used to screen components with potentially activated TRPV1 channel activity based on the generated pharmacophore model.
2.2 Virtual screening based on molecular docking
2.2.1 Preparations for molecular docking
Molecular docking is a method used to study the interaction between small molecules and target proteins. It can efficiently and accurately predict drug activity and selectivity, thus avoiding blind screening and reducing the cost of screening active compounds.
In the molecular docking experiment, the structure of MA also needed to be downloaded from PubChem (https://pubchem.ncbi.nlm.nih.gov/). We still use the small molecule ligands in Supplementary Table S1. Then we handled them with energy minimization and picked root and rotatable keys via Autodock Vina software (version 1.5.6), which was also used for the subsequent molecular docking calculations.
The two TRPV1 protein receptors that we selected for molecular docking were downloaded from the Protein Data Bank (PDB) database (http://www.ncbi.nlm.nih.gow/protein/). They are respectively the structure of TRPV1 in complex with DkTx and RTX (PDB ID code: 5IRX; Resolution: 2.95Å), and the structure of TRPV1 in complex with capsazepine (PDB ID code: 5IS0; Resolution: 3.43Å) (Gao et al., 2016).
To conduct virtual screening, firstly the active binding sites of molecular targets were determined by Discovery Studio. RTX is a TRPV1 agonist that interacts with residues: ARG557, LEU553, ALA566, ASN551, THR550, ILE569, MET547, ALA546, PHE591, LEU669, ALA665, TYR511, PHE543, LEU515, PHE507, ILE514, VAL518 of 5IRX (Supplementary Figures S1A, B). Capsazepine is a TRPV1 antagonist which interacts with residues: GLU570, LEU553, TYR554, TYR511, SER512, THR550, MET547, ALA665, PHE543, LEU662 of 5IS0 (Supplementary Figures S1C, D). The comparison shows that the TRPV1 agonist and the TRPV1 antagonist act on the same region of the TRPV1 protein, while the residues that act are not the same.
We needed to pretreat the protein receptors, including dehydration, hydrogenation, etc. The active sites of proteins also needs to be determined. What we obtained from the PDB database were the protein-ligand complexes. By analyzing the spatial position relationship between the receptor and the original ligand at the docking site, the information of the docking site can be obtained (Feinstein and Brylinski, 2015). The coordinates of RTX in 5IRX were (X = 93.69, Y = 114.04, Z = 129.88), and generated a square docking box with a side length of 21.275 centered on this site. The coordinates of capsazepine in 5IS0 were (X = 108.115, Y = 94.986, Z = 102.379), and generated a square docking box with a side length of 21.141 centered on this site. In addition, we researched the necessary literature to corroborate the location (Li and Yang, 2018).
2.2.2 Molecular docking
Molecular docking was carried out on Autodock Vina. During the procedure, the semi-flexible docking method was used, that is, the protein receptor was rigid and the partial bond of the small molecules ligand could rotate. The TRPV1 agonists were placed into the indicated cave of the 5IRX receptor, and the TRPV1 antagonists were placed in the specified box in the 5IRX receptor. And then the components of CF (Supplementary Table S1) and the MA were docked to the two receptors severally and evaluated their binding ability.
2.3 Virtual screening based on 3D-QSAR model
2.3.1 Construction of 3D-QSAR model
The 3D-QSAR model is a pharmacophore model with activity prediction ability based on a series of compounds with clear activity values for specific biological targets. The model was used to predict the related activities of new compounds that were not experimentally determined or even not synthesized. It is highly efficient, especially in the screening of compounds in large databases has an irreplaceable advantage.
Twelve TRPV1 agonists with the same EC50 data source were selected from the BindingDB database (https://www.bindingdb.org/). The corresponding small molecule ligands were downloaded from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) and the pEC50 data were added to the small molecule ligands (Supplementary Table S3). Chem 3D software was used to minimize the MMFF94 force field energy of all small molecule ligands. The first seven small molecules in the table were selected to form the training set, while all 12 TRPV1 agonists in Supplementary Table S3 will be used as the test set.
Using the Create 3D QSAR Model module in Discovery Studio, a 3D-QSAR model with activity prediction ability was constructed based on the common structural features and pEC50 values of small molecules selected in the training set, and the constructed model was verified by the test set. Representative components and a 3D-QSAR model were selected to generate visual contour maps of the electrostatic field coefficient and stereo field coefficient.
2.3.2 Prediction of activity based on 3D-QSAR model
The Calculate Molecular Properties module in Discovery Studio was used to predict the activity of prepared triterpenoids and MA using the already constructed 3D-QSAR model, and output the results.
2.4 Verification of triterpenoids alleviating the H9c2 cardiomyocyte toxicity caused by MA
2.4.1 Preparations of materials and cell culture
Rat cardiomyocyte H9c2 cells were purchased from the Cell Resource Center, IBMS, CAMS/PUMC (CRC/PUMC, Beijing, China). The H9c2 cells were cultured in high-glucose dulbecco’s modified eagle medium (DMEM) medium (VivaCell Biosciences, China) supplemented with 10% fetal bovine serum (FBS) (VivaCell Biosciences, China) and 1% penicillin-streptomycin (Lablead, China), and then incubated at 37°C with 5% CO2.
Corosolic acid (CRA, >98%, DSTDX000501) Maslinic acid (MSA, >98%, DST220526-039) Arjungenin (AR, >98%, DST210106-442) Mesaconitine (MA, >98%, DSTDX002501) were purchased from DeSiTe (Chengdu, China). Capsazepine (>98%, D15O9F2190) was obtained from Yuanye Biochemical (Shanghai, China). All test samples were accurately weighed, and dissolved by dimethyl sulfoxide (DMSO) (Sigama-Aldrich, United States) to obtain concentrated solutions, and then preserved at 4°C. They were diluted to the corresponding concentration with medium prior to use and allowed the DMSO concentration less than 0.5%.
2.4.2 Cell viability of triterpenoids on H9c2 cardiomyocyte
Cell viability was assessed using the 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) (Lablead, Beijing, China) assay. Initially, the MTT assay was utilized to identify a safe and non-toxic concentration range of triterpenoids to prevent the interference of triterpenoids with the survival rate of H9c2 cells and the subsequent experimental results. H9c2 cells were cultured with a medium containing various concentrations of triterpenoids for 24 h. The following grouping and dosages were used: 1) control group: H9c2 cells were cultured without any treatment; 2) CRA group: H9c2 cells were treated with different concentrations of CRA (8, 16, 32, 64, 128 μmol/L); 3) MSA group: H9c2 cells were treated with different concentrations of MSA (8, 16, 32, 64, 128 μmol/L); 4) AR group: H9c2 cells were treated with different concentrations of AR (8, 16, 32, 64, 128 μmol/L).
Upon incubation for 24 h, MTT (0.5 mg/mL) was added to each well, followed by a 4-h incubation period inside a light-protected incubator. Subsequently, DMSO was added and the optical density (OD) was measured at 490 nm wavelength. The cell survival was calculated using the relevant Formula 1.
2.4.3 Cell viability of triterpenoids on MA-induced H9c2 cardiomyocyte
We conducted the MTT assay to evaluate the ability of triterpenoids from CF to attenuate cardiac toxicity induced by MA. The best-attenuated concentration was identified and used in subsequent experiments. The MTT assay was performed as described in Section 2.4.2. They were distributed into four groups as follows: 1) control group: H9c2 cells were cultured for 24 h without any treatment; 2) MA group: H9c2 cells were treated with 300 μmol/L of MA for 24 h; 3) MA+CRA group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of CRA (2, 4, 8, 16, 32 μmol/L) for 24 h; 4) MA+MSA group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of MSA (2, 4, 8, 16, 32, 64 μmol/L) for 24 h; 5) MA+AR group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of AR (2, 4, 8, 16, 32, 64, 128 μmol/L) for 24 h.
2.4.4 Nucleus morphology of triterpenoids on MA-induced H9c2 cardiomyocyte
Cell apoptosis was analyzed using Hoechst 33258 staining (Lablead, Beijing, China). H9c2 cells were treated with a medium containing MA or different concentrations of triterpenoids for 24 h. The grouping and dosages were as follows: 1) Control group: H9c2 cells were cultured without any treatment; 2) MA group: H9c2 cells were treated with 300 μmol/L of MA; 3) MA+CRA group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of CRA (28, 32, 36 μmol/L); 4) MA+MSA group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of MSA (12, 16, 20 μmol/L); 5) MA+AR group: H9c2 cells were treated with 300 μmol/L of MA and different concentrations of AR (60, 64, 68 μmol/L).
Following the treatment, the cells were stained with the Hoechst 33258 for 20 min in the absence of light. After then observed under a fluorescence microscope.
2.4.5 Mitochondrial membrane potential of triterpenoids on MA-induced H9c2 cardiomyocyte
We assessed changes in mitochondrial membrane potential in H9c2 cells by using JC-1 dye staining (Solarbio, Beijing, China) to confirm the mechanism by which triterpenoids from CF can reduce cardiac damage induced by mesaconitine. H9c2 cells were treated with a medium containing MA or triterpenoids for 24 h. The grouping and dosages were as follows: 1) control group: H9c2 cells were cultured without any treatment; 2) MA group: H9c2 cells were treated with 300 μmol/L of MA; 3) MA+CRA group: H9c2 cells were treated with 300 μmol/L of MA and 32 μmol/L of CRA; 4) MA+MSA group: H9c2 cells were treated with 300 μmol/L of MA and 20 μmol/L of MSA; 5) MA+AR group: H9c2 cells were treated with 300 μmol/L of MA and 64 μmol/L of AR.
Subsequently, the H9c2 cells were treated with JC-1 staining solution. And then incubated for 20 min in a 5% CO2 incubator at 37°C. Apoptotic cells were directly visualized with inverted fluorescence microscopy.
2.5 Verification of triterpenoids reducing H9c2 cardiomyocyte toxicity induced by MA based on TRPV1 channel
2.5.1 TRPV1 expression in the triterpenoids or MA-treated H9c2 cardiomyocyte
The Western blot assay was performed to validate the results of the previous pharmacophore screening and molecular docking, which screened the triterpenoids with activation activity on the TRPV1 channel. This laid the groundwork for the subsequent confirmation of whether the effect of triterpenoids alleviating cardiotoxicity caused by MA was related to TRPV1.
H9c2 cells were treated with a medium containing MA or triterpenoids for 24 h. The grouping and dosages were as follows: 1) control group: H9c2 cells were cultured without any treatment; 2) MA group: H9c2 cells were treated with 300 μmol/L of MA; 3) CRA group: H9c2 cells were treated with 32 μmol/L of CRA; 4) MSA group: H9c2 cells were treated with 20 μmol/L of MSA; 5) AR group: H9c2 cells were treated with 64 μmol/L of AR.
Following treatment, the H9c2 cells were lysed with the total cell protein extraction reagent and then centrifuged at 13,000 g for 5 min at 4°C. The resulting supernatant collected was the total protein solution. The levels of proteins in the supernatant were determined using a BCA kit (Baiqiandu, Wuhan, China). Afterward, the total protein was separated by SDS-polyacrylamide gel electrophoresis (SDS-PAGE, Baiqiandu, Wuhan, China) and transferred onto polyvinylidene difluoride membranes (PVDF, Millipore, United States). The membranes were treated with the various primary antibodies (GAPDH, 1:6,000; TRPV1, 1:1,000) overnight after being blocked with 5% non-fat dry milk at 37°C for 2 h. After washing the membrane with TBST (Baiqiandu, Wuhan, China), the membranes were incubated with secondary antibodies for 30 min at room temperature. Finally, an ECL chemiluminescent (Baiqiandu, Wuhan, China) working buffer was treated with the protein bands. After that, the bands’ quantitative grayscale analysis was done using the IPWIN60 program.
2.5.2 Intracellular calcium ions content of triterpenoids on MA-induced H9c2 cardiomyocyte
The TRPV1 channel is a calcium ion channel, and it plays a crucial role in the flow of calcium ions into cells. The effect of this channel on calcium ions in cardiomyocytes was detected by Fluo-3 AM fluorescent probe method. The H9c2 cells were treated with a medium containing MA or various doses of triterpenoids for 24 h. Groups allocation and dose were the same as in Section 2.4.5. Following the treatment, the Fluo-3 AM (Solarbio, Beijing, China) (5 μmol/L) was used to stain the cells. The cells were then rinsed with HBBS after being incubated for 60 min at 5% CO2 at 37°C, and the fluorescent intensity was assessed using a fluorescence microscope.
2.5.3 Cell viability of triterpenoids on MA-induced H9c2 cardiomyocyte after TRPV1 channel blocked by the capsazepine
To further explore whether the attenuated effect is related to the TRPV1 channel, we observed the effect of triterpenoids on the survival rate of H9c2 cells treated with MA before and after the TRPV1 channel blockade by MTT assay. H9c2 cells were treated with a medium containing MA or triterpenoids for 24 h. The grouping and dosages were as follows: 1) control group: H9c2 cells were cultured without any treatment; 2) MA group: H9c2 cells were treated with 300 μmol/L of MA; 3) triterpenoid + MA group: H9c2 cells were treated with 300 μmol/L of MA and one of the triterpenoids (32 μmol/L of CRA; 20 μmol/L of MSA; 64 μmol/L of AR); 4) Capsazepine + triterpenoid + MA group: H9c2 cells were pretreated with capsazepine (10 μmol/L) for 30 min before 300 μmol/L of MA and one of the triterpenoids (32 μmol/L of CRA; 20 μmol/L of MSA; 64 μmol/L of AR).
2.6 Statistical analysis
All of the above cell experiments were independently repeated three times. The GraphPad Prism (version 8.0.1) was utilized to create all statistical calculations and graphs. Two-tailed Student’s t-test was implemented to compare the two groups. ANOVA was used to compare three or more groups against one another. The mean and standard deviation (SD) were employed to express the results. When p < 0.05, the differences were deemed statistically significant.
3 Results
3.1 Pharmacophore model screening
3.1.1 Pharmacophore model establishment and validation
Based on the commonly shared features amongst seven agonists of TRPV1 selected from the training set, ten pharmacophoric models were constructed and their pharmacophoric characteristics along with Rank values were reported in Table 1. All the pharmacophores contained four features, namely, A: ACCEPTOR (Hydrogen bond acceptor) and H: HYDROPHOBIC had two each. The Rank values of these ten pharmacophoric models were greater than 60, signifying the prominence of the molecular common features while constructing the pharmacophores. Additionally, the differences among the top scores were extremely small. Hence, it was imperative to validate and screen the generated 10 pharmacophore models using a test set (Figure 2). The ligands that depict higher responses to the pharmacophore models have a warm tone and those with lower responses have a cold tone. As seen, pharmacophore model 02 portrayed high responses against the active ingredients present in the test set and relatively lower responses to the inactive ingredients, thereby meeting the study requirements. Hence, pharmacophore model 02 was selected for virtual screening.
FIGURE 2. Validation of 10 pharmacophore models. (The bottom 10 color blocks in each column correspond to the TRPV1 agonists, and the rest 90 correspond to the components that have not been verified for activity. The ligands that depict higher responses to the pharmacophore models have a warm tone while those with lower responses have a cold tone.)
3.1.2 Virtual screening based on pharmacophore model 02
The 60 components in CF were virtually screened using the pharmacophore model 02 selected in Section 3.1.1, and eight components were finally selected (Table 2), comprising seven triterpenoids and one alkaloid. Furthermore, the matching diagram of the pharmacophore and some potential active components were shown in Figure 3. An observation from the figure depicted that the triterpenoids exhibited similar structural features as that of the pharmacophore model.
FIGURE 3. Pharmacophore model 02 and mapping with the selected triterpenoids. (The green spheres represent hydrophobic features, while blue spheres represent hydrogen bond acceptor features.)
3.2 Molecular docking analysis
3.2.1 Molecular docking validation
When assessing the stability of the molecular docking model, re-docking the co-crystallized ligand into the empty binding pocket of the protein was employed. The calculated pose was then compared to the bioactive conformation from the crystal structure and the root mean square distance (RMSD) was calculated. Typically, RMSD values below 2 Å are considered acceptable (Temml and Kutil, 2021). Upon re-docking, the structural and spatial positions were fairly consistent with the native ligand. The RMSD (5IRX) = 1.2826 Å < 2 Å, RSMD (5IS0) = 1.7092 Å < 2 Å. Lower Score indicates higher affinity towards the TRPV1 receptor and greater propensity for interaction. The outcome (Table 3) exhibits that most of the other small molecule ligands displays high affinity and are in close proximity to the aforementioned ligands. This affirmed the high reliability of the molecular docking method utilized in this study. Meanwhile, The docking data also provides reference values for the subsequent molecular docking simultaneously.
TABLE 3. Docking scores of TRPV1 agonists and TRPV1 antagonists with 5IRX and 5IS0 respectively (Kcal/mol).
3.2.2 Virtual screening based on molecular docking
The docking outcomes of CF components with TRPV1 receptors were presented in Supplementary Table S4. All of them and MA docking results were scrutinized and subsequently compared (Table 4; Figure 4). The results indicates that pharmacophore-selected triterpenoids exhibits higher affinity towards TRPV1 receptor 5IRX than MA (**p < 0.01) as well as other CF components (##p < 0.01). Furthermore, these constituents demonstrates a higher affinity for the 5IRX receptor than the 5IS0 receptor (∆∆∆p < 0.001).
TABLE 4. Docking scores of mesaconitine and the selected triterpenoids by the pharmacophore with 5IRX and 5IS0 (Kcal/mol).
FIGURE 4. Docking score analysis. Molecular docking results of MA (A0), CF all components (A1), triterpenoids selected by the pharmacophore model 02 (A2) with 5IRX; molecular docking results of MA (B0), CF all components (B1), triterpenoids selected by the pharmacophore model 02 (B2) with 5IS0 (Compare A2 to A0, **p < 0.01; compare A2 to A1, ##p < 0.01; compare B2 to A2, △△△p < 0.001).
3.3 Activity prediction of triterpenoids and MA based on 3D-QSAR model
3.3.1 3D-QSAR model construction and validation
The 3D-QSAR model is a pharmacophore model with activity prediction ability based on a series of compounds with clear activity values for specific biological targets. The model constructed in this section was based on the electrostatic field and steric field characteristics of TRPV1 agonists, which can be used to predict the activity of unknown compounds, thus providing a basis for subsequent experimental verification and mechanism exploration (Figures 5A, B). In the electrostatic map, the red area indicates that high negative charges enhance the compound’s activity, while the blue area suggests high positive charges improve the compound activity. In the stereo field map, the yellow region denotes that expanding this area’s volume does not enhance the compound activity, while the green area indicates that increasing the region’s volume improves the compound activity.
FIGURE 5. 3D-QSAR model electrostatic and stereo field coefficient contour maps. (A) the high map of electrostatic field coefficient. (B) the high map of stereo field coefficient. (C) The fitting diagram of the prediction value of the 3D-QSAR model and the experimental value based on the training set. (D) The fitting diagram of the prediction value of the 3D-QSAR model and the experimental value based on the test set.
The training set’s predicted activity values were compared to experimental activity values (Figure 5C), with the R2 value reflecting the accuracy of the curve. The curve’s R2 value of 1, indicating accurate prediction and closely matching experimental values.
To verify the accuracy of the 3D-QSAR model, the fitted curve based on the prediction and experimental values of the test set was created (Figure 5D). The R2 value of 0.526 indicates that the model has some confidence, with the predicted values reflecting the experimental values to a reasonable extent. The model has the potential to predict the activity order of triterpenoids and MA’s potential as a TRPV1 agonist.
3.3.2 Results of activity prediction based on the 3D-QSAR model
The developed 3D-QSAR model was utilized for predicting the activity of triterpenoids and MA, as shown in Table 5. All of the eight components in the table exhibited activity, among which MA had the strongest predicted activity with a predicted pEC50 of −0.446712. The triterpenoids exhibited relatively weak activity, with CRA, MSA, and AR exhibiting the weakest activity among them.
3.4 Verification of triterpenoids alleviating the H9c2 cardiomyocyte toxicity caused by MA
3.4.1 Effect of triterpenoids on H9c2 cell viability
The structures of the three triterpenoids were shown in Figure 6A. As shown in Figures 6B–D, the cell survival rate decreased when the CRA was in the range of 64 μmol/L to 128 μmol/L and MSA at 128 μmol/L (##p < 0.01, ###p < 0.001). In contrast, AR not only exhibited safety and non-toxicity towards H9c2 cells in the concentration range of all 128 μmol/L, but also promoted cell growth (#p < 0.05), and the proliferation reached its maximum effect at 64 μmol/L of AR (##p < 0.01). For subsequent experimental verification, concentration ranges of 0 ∼ 32 μmol/L, 0 ∼64 μmol/L, and 0 ∼ 128 μmol/L were selected for CRA, MSA, and AR, respectively.
FIGURE 6. Triterpenoids attenuated cardiomyocyte toxicity caused by MA. (A) The structure of selected triterpenoids. Effect of CRA (B), MSA (C), and AR (D) on cell viability by MTT assay. The attenuating effect of CRA (E), MSA (F), and AR (G) on MA caused myocardial cytotoxicity by MTT assay. Effect of triterpenoids or MA on cell nucleus apoptosis in H9c2 cells by Hoechst 33258 staining, (H) fluorography (×100); (I) quantification is expressed as the Mean fluorescence indensity (MFI). Effect of triterpenoids or MA on mitochondrial membrane potential in H9c2 cells by JC-1 dye staining, (J) fluorography (×100); (K) quantification is expressed as the MFI. Data were presented as mean ± SD (n = 3. Compared with normal control group, #p < 0.05, ##p < 0.01; compared with MA group, *p < 0.05, **p < 0.01, ***p < 0.01).
3.4.2 Effect of triterpenoids on MA-induced cell viability in H9c2 cardiomyocyte
To confirm the triterpenoids can lessen the H9c2 cardiomyocyte toxicity brought on by MA, the survival rate of H9c2 cells was examined using the MTT assay (Figures 6E–G). MA (300 μmol/L) significantly decreased the survival rate of H9c2 cells as compared to the control group (##p < 0.01). All three triterpenoids have the potential to reduce the toxicity of MA on cardiomyocytes and increase cell survival (32 μmol/L of CRA, 16 μmol/L of MSA, 64 μmol/L of AR, **p < 0.01). And the optimized attenuated concentration were selected for subsequent experiments.
3.4.3 Effect of triterpenoids on MA-induced nuclear apoptosis in H9c2 cardiomyocyte
The cells nuclear change was obtained by Hoechst 33258 staining (Figures 6H, I). Hoechst 33258 specifically stains the nucleus, and when the cell becomes apoptotic, its nuclei is colored blue and fragmented. The control group’s nucleus displayed light blue fluorescence following staining and the nucleus shape were unaltered in the Hoechst 33258 staining results. Compared to the control group, the nucleus’ apoptotic features became visible in the MA group. It was obviously that the nucleus’ chromatin was compressed, and some fragments could be observed. The average fluorescence intensity was much higher than that of the control group on the entire, and it displayed dense and strong bright blue fluorescence (###p < 0.001). In comparison to the MA group, the addition of a small amount of CRA, MSA, and AR greatly enhanced cell apoptosis, dramatically enhanced nuclear chromatin condensation, and lowered fluorescence intensity (***p < 0.001, **p < 0.01).
3.4.4 Effect of triterpenoids on MA-induced mitochondrial membrane potential in H9c2 cardiomyocyte
The cells mitochondrial membrane potential change was measured by JC-1 staining (Figures 6J, K). During the JC-1 staining experiment, healthy cells produce red fluorescence due to their higher mitochondrial membrane potential. Conversely, apoptotic cells produce green fluorescence due to their lower membrane potential. A decrease in mitochondrial membrane potential is an early sign of apoptosis, and this decrease can be detected easily by the transition from red fluorescence to green fluorescence when cellular membrane potential drops. In the control group, orange-red fluorescence was predominantly observed, with only a small amount of green fluorescence indicating higher mitochondrial membrane potential. The MA group showed decreased mitochondrial membrane potential, as indicated by greater green fluorescence and weaker orange-red fluorescence (###p < 0.001). In comparison to the MA group, the inclusion of triterpenoids significantly increased mitochondrial membrane potential (***p < 0.001, *p < 0.01, *p < 0.05).
Through this part of the experiment, we selected the concentrations of CRA, MSA, and AR as 32 μmol/L, 20 μmol/L, and 64 μmol/L for the subsequent exploration of the mechanism of triterpenoids reducing H9c2 cardiomyocyte toxicity induced by MA.
3.5 Verification of triterpenoids reducing MA-induced H9c2 cardiomyocyte toxicity via TRPV1 channel
3.5.1 Effect of triterpenoids or MA on TRPV1 expression in H9c2 cardiomyocyte
As depicted in Figures 7A, B, MA, MSA, and AR significantly increased the expression of the TRPV1 protein in H9c2 cells (**p < 0.01, ***p < 0.001), although MSA and AR had a weaker effect when used in attenuated toxicity concentrations. CRA specifically failed to induce TRPV1 protein expression.
FIGURE 7. Triterpenoids attenuated cardiomyocyte toxicity caused by MA through TRPV1 channel. Effect of triterpenoids or MA on expression of TRPV1 in H9c2 cells detected by Western blot, (A) Western blot result, (B) the semiquantitative result of (A). Effect of triterpenoids or MA on Ca2+ concentration in H9c2 cells by Fluo-3 AM, (C) fluorography (×100); (D) quantification is expressed as the MFI. (E) Effect of triterpenoids or MA on cell survival before and after TRPV1 channel blockade by MTT assay and capsazepine. Data were presented as mean ± SD (n = 3. Compared with normal control group, #p < 0.05, ##p < 0.01; compared with MA group, *p < 0.05, **p < 0.01, ***p < 0.001; compared with the group of the TRPV1 channels that were not blocked, △p < 0.05, △△△p < 0.001).
3.5.2 Effect of triterpenoids on MA-induced intracellular calcium ions content in H9c2 cardiomyocyte
The changes in calcium ion concentrations in H9c2 cells were estimated. As illustrated in Figures 7C, D, the intracellular Ca2+ level was lower in the control group, and only a trace of green fluorescence was observed. However, MA may greatly boost the green fluorescence intensity of H9c2 cells as well as the intracellular Ca2+ content (###p < 0.001). The addition of CRA, MSA, and AR significantly reduced both the green signal and intracellular Ca2+ levels compared to the MA group (*p < 0.05, **p < 0.01).
3.5.3 Effect of triterpenoids on MA-induced H9c2 cardiomyocyte toxicity after TRPV1 channel blocked by the capsazepine
MA dramatically decreased the H9c2 cells’ survival rate as compared to the control group (#p < 0.01) (Figure 7E). Compared to the group without capsazepine, the AR and MSA groups that had previously received capsazepine were able to further improve the survival rate of the H9c2 cells (∆p < 0.05, ∆∆∆p < 0.001), but the CRA group was not significantly different regardless of whether they had received capsazepine in advance treatment.
4 Discussion
Aconitum herbs have a long history of medication, such as Aconitx kusnezoffll radx., and Aconitum pendulum Busch. They have good anti-inflammatory and analgesic properties, but their toxicity is also extremely strong (Singhuber et al., 2009). It is necessary to study the measures and mechanism of toxicity reduction of aconitums, which is of great significance for safe drug use in clinical practice. In Mongolian and Tibetan medicine theory, CF is often used to process or combine aconitums to reduce their toxicity (Liu et al., 2017; Liet al., 2022b). Studies have shown that CF has the effect of reducing the toxicity of aconitums (Liu et al., 2015a; Zhang et al., 2017; Yuan et al., 2022), and the mechanism of the attenuation toxicity was not the same as the conventional heating and boiling to convert the toxic DDAs into less toxic and more powerful monoester-diterpenoid alkaloids (MDAs) or non-esterified diterpene alkaloids (NDAs) (Liu et al., 2013), involving multiple pathways in vitro and in vivo, such as material basis, internal absorption and metabolism, myocardial electrophysiology, etc. (Zhang, 2019; Tan et al., 2023). On the basis of the previous studies, we aimed to clarify the mechanism of reducing cardiotoxicity from the aspect of TRPV1 channel that DDAs induced cardiotoxicity in this study.
The pharmacophore model, molecular docking and 3D-QSAR model were used for virtual screening in this study. To screen for components in CF that may activate the TRPV1 channel, a pharmacophore model based on common features of TRPV1 agonists was constructed. However, the pharmacophore model can only identify components with activation potential; it cannot determine their relative activity or ability to bind to the TRPV1 protein. Furthermore, molecular docking was used to quantify each component’s binding capacity to the TRPV1 protein receptor, which served as a confirmation and secondary screening of the pharmacophore results. However, further clarification of their activities is required. Consequently, a 3D-QSAR model based on TRPV1 agonists was developed to predict the activity of the screened chemical components. The pharmacophore model revealed seven triterpenoids in CF with TRPV1 agonistic effects. Subsequently, the molecular docking and 3D-QSAR model results indicated that the triterpenoids of CF, which were screened by pharmacophore, had a higher affinity to bind to the TRPV1 agonist receptor than MA. However, their relative activity was lower than that of MA, with CRA, MSA, and AR exhibiting the lowest predicted activity. Thus, we initially integrated the results of the virtual screening section and speculated that the triterpenoids in CF function as TRPV1 agonists. These compounds have a higher affinity for TRPV1 receptors and can competitively antagonize the binding of MA to TRPV1 receptors, whereas having lower activity. This prevents the over-activation of TRPV1 channels and exerts an attenuated effect.
The rat H9c2 cardiomyocytes were then used to validate the results of the virtual screen. The administration concentration of MA was determined according to the previous pre-experiments. The concentrations of triterpenoids were selected step by step during the experiment. First, the triterpenoids were screened for safe and non-toxic ranges. Then the best attenuation was screened in the attenuated verification experiment section by degrees. Obviously, MA damaged the viability of H9c2 cardiomyocytes, seriously reduced the survival rate. It also caused nuclear damage and pyknosis, and decreased mitochondrial membrane potential. These are landmark events in the early stage of apoptosis and also prove the cardiotoxicity of MA. Moreover, MA promoted the expression of TRPV1 protein in cardiomyocytes, which in turn significantly increased the content of calcium ions in cardiomyocytes. This was consistent with the previous study that the MA had myocardial cytotoxicity, and the generation of toxicity was related to the excessive activation of TRPV1 channels. It was also established that the triterpenoids could reduce the myocardial cytotoxicity of MA to some extent, increase the mitochondrial membrane potential and relieve the damage to the nucleus caused by MA. In addition, MSA and AR can similarly stimulated TRPV1 protein expression in H9c2 cells, although their capacities were inferior to that of MA. In addition, the H9c2 cells survival rate increased when TRPV1 channels were blocked by capsazepine, indicating that the attenuated effects of MSA and AR were related to TRPV1 channels and acted as agonists of TRPV1 channels. Not only MSA and AR but also CRA, could attenuate this tendency for large calcium influx caused by MA and maintain homeostasis of intracellular calcium. These tests provided additional evidence that the TRPV1 channel played a role in reducing cardiotoxicity caused by MA, as the TRPV1 channel is a crucial channel for calcium ions to flow into cells.
So far, efforts have been devoted to the study of the mechanism of CF reducing the toxicity of aconitums, with the expectation of providing theoretical support for the safe application of aconitums. The mechanisms can be summarized in four categories: 1) The acidic components in the CF can inhibit the dissolution of DDAs and detoxify the aconitums (Liu et al., 2013; Yang et al., 2013; Xin et al., 2014). 2) Tannins in CF complexed with DDAs in aconitums to form insoluble substances, which slows down the dissolution rate of alkaloids in gastric and intestinal fluids and plays a certain role in reducing toxicity (Liu et al., 2015b). 3) CF can play a role in reducing toxicity by affecting the absorption, distribution and metabolic process of aconitums in the body. For absorption: the pharmacokinetics of alkaloids can be changed after the processing of aconitums or the compatibility with CF, and the absorption level and absorption rate of alkaloids are reduced, thus reducing the toxicity (Zhi et al., 2020b; Liu et al., 2023). For distribution: CF can reduce the distribution of DDAs in various organs and tissues, and the tannins can bind to DDAs and slowly distribute and eliminate them in the body to play a role in detoxification (Wang et al., 2002). Tannic acid in can reduce the plasma protein binding rate and bioavailability of DDAs (Zhang, 2019). For metabolism: CF can change the absorption and metabolism sites of toxic components, reduce its absorption in the intestinal wall, and increase intestinal bacteria and liver metabolism. This can prevent the absorption of toxic components from being too fast, and promote the metabolism, reduce the concentration of drugs in the blood and prevent poisoning (Celimuge et al., 2022). 4) The components in CF have myocardial protective effect, which can effectively inhibit the cardiotoxicity of aconitums. The use of CF can reduce the myocardial tissue damage caused by the single administration of aconitums in rats (Zhang et al., 2017). Scholars have proved that CF had a good intervention effect on arrhythmia caused by aconitine by Langendorff isolated heart perfusion technique (Liang et al., 2017). CF can resolve the toxicity of Aconitums and has a protective effect on the heart (Pan et al., 2011).
CF alleviates the cardiotoxicity of aconitums mainly due to the active ingredients in CF, such as tannins and triterpenoids, protect the heart from the serious toxicity of aconitums through various ways in vivo. Researchers have found that CF could regulate the contraction and relaxation of myocardial cells by multi-component, multi-target, multi-channel, further protect the normal function of the heart, and reduce the cardiac toxicity caused by aconitiums by using the method of network pharmacology (Li et al., 2018). Metabolomics technology was used to investigate the effect of CF on endogenous metabolites in rats. The results showed that glutamine and lysophosphatidylcholine may be the key to the detoxification and regulation efficacy of CF, and also have a certain role in improving cardiac function and acute myocardial infarction (Bilige et al., 2022). The combination of ellagic acid, liquiritin and aconitine can significantly increase the expression of cardiac metabolic enzyme CYP2J3, promote the metabolism of arachidonic acid to produce epoxyeicosatrienoic acids, and then reduce the cardiotoxicity (Li et al., 2023). The effective components of CF can inhibit the release of lactate dehydrogenase and creatine kinase isoenzyme and promote the release of superoxide dismutase to achieve the protective effect on aconitine-injured cardiomyocytes (Duolan, 2022). The compatibility of CF with aconitums can regulate the levels of FKBP1B and RyR2 genes in sarcoplasmic reticulum by affecting TRP gene, especially TRPM8, so as to protect myocardium from damage and play a role in reducing toxicity (Meng et al., 2023). Our research group found that Caowu could induce cardiotoxicity by activating the TRPV1 channel, while CF also played an attenuated role through the TRPV1 channel. The attenuated effect may be related to the synergistic effect of the acidic components in CF and the alkaloids in Caowu on the TRPV1 channel (Han et al., 2022c). Subsequently, we used the uniform design method to identify the primary components responsible in CF for the detoxification effect, and then gallic acid and ellagic acid were screened and the two a synergistic effect (Han et al., 2022a). We also demonstrated that chebulagic acid in CF has a protective effect on myocardial mitochondria, and can inhibit aconitine-induced Ca2+ overload in myocardial cells, thereby reducing the occurrence of arrhythmia (Han et al., 2020). In addition, we also found that the TRPV1 channel have an important role in exerting their attenuated effects, the gallic acid in CF may achieve the effect of detoxification by rapid desensitization of TRPV1 channel on the cells (Han et al., 2022b; Mi et al., 2023).
In summary, the research on the active ingredients reducing cardiotoxicity caused by aconitums mainly focused on tannins and phenolic acids. These components play a mitigating cardiotoxicity role by affecting the expression of relevant genes, enzyme activity, and ligand-receptor interactions to maintain a steady state within the cardiomyocyte and protect the organelles to maintain normal cellular function. In addition to tannins, the triterpenoids in CF are also noteworthy, and most of them have cardiac and myocardial protective effects. MSA has a protective impact against cardiomyocyte damage by exerting both an anti-oxidative action and an ability to inhibit calcium overload in cardiomyocytes (Liu et al., 2008). The cardioprotective effect of arjunolic acid may lie in its protective potential against damage caused by myocardial necrosis (Sumitra et al., 2001). CRA not only upregulates antioxidant levels in cardiomyocytes and hence inhibits apoptosis and cardiomyocyte stress, but also effectively inhibits the development of fibrosis after myocardial infarction, thereby improving cardiac function (Wang, 2019; An et al., 2020). They are in line with how DDAs in Aconitum herbs exert cardiotoxicity. It was also discovered that the triterpenoids in CF showed a certain potential for TRPV1 activation, which may play an attenuated role (Mi et al., 2023).
Combined with the virtual screening and the attenuated validation experiment, the results preliminarily demonstrated that the triterpenoids in CF played a certain role in CF reducing aconitums cardiotoxicity. MSA and AR most likely acted as TPPV1 agonists with less potency that acted on the TRPV1 channel in cardiomyocytes like MA. Because of their high binding capacity, they can bind better to TRPV1 channels but exert weaker agonism, thus preventing MA from excessively stimulating the TRPV1 channel, inhibiting the large influx of calcium ions, keeping cells viable and achieving the purpose of attenuating cardiotoxicity. Besides, it was discovered that CRA failed to increase the expression of TRPV1 protein on its own, and its attenuating effect did not change before and after capsazepine blocked the TRPV1 channel. Nevertheless, CRA was still able to alleviate the calcium ion influx caused by MA. In this context, it is hypothesized that CRA did not achieve an attenuated effect through the TRPV1 channel, but it was likely through other channels related to calcium ions, such as the MAPK pathway (Schaub et al., 2006; Ai et al., 2016), and we will also explore this in depth next.
The toxicity and efficacy of aconitums are closely related to the TRPV1 channel, so reducing their cardiotoxicity through this channel while maintaining good efficacy is very meaningful. In this study, based on the theory of compatibility or processing in Traditional Chinese Medicine Science, the mechanism of reducing the toxicity of MA by triterpenoids of CF was elucidated from the perspective of antagonistic or synergistic effects of active components of drugs. At present, we have only explored the mechanism through computer simulation and at the cellular level. Next, we will further investigate the mechanism of attenuation at the animal level. Of course, other pathways also play a role in the reduction process, and the link and leading role between them are unclear. The DDAs toxicity is produced by the TRPV1 channel, and its analgesic anti-inflammatory effect is also related to TRPV1 channel. We only have evidence that triterpenoids can reduce MA toxicity at the cardiomyocyte level, but we are not clear about their treatment effect. So how can reduce the toxicity and at the same time show good treatment effects needs to further explored.
5 Conclusion
The triterpenoids in CF can reduce the H9c2 cardiomyocyte toxicity caused by MA. Through excessive TRPV1 channel activation, MA causes significant cardiomyocyte death, while MSA and AR also exert a detoxifying effect on cells via TRPV1 channels. Due to their higher affinity with TRPV1 protein receptors, they can competitively antagonize the binding of MA to TRPV1 channels, but at the same time, they have lower TRPV1 channel activation activity than MA, thus preventing MA from overactivating TRPV1 channels, thereby protecting cardiomyocytes and reducing toxicity. CRA may play an attenuated role through other pathways that affect intracellular calcium homeostasis. This paper provides ideas and credence to investigate the mechanism of CF alleviating the cardiotoxicity caused by aconitums in vitro experiments.
Abbreviation
AR, Arjungenin; CF, Chebulae Fructus; CRA, Corosolic acid; DDAs, Diesterditerpenoid alkaloids; DMEM, Dulbecco’s modified eagle medium; DMSO, Dimethyl sulfoxide; EC50, Concentration for 50% of maximal effect; FBS, Fetal bovine serum; LD50, Median lethal dose; MA, Mesaconitine; MDAs, Monoester-diterpenoid alkaloids; MFI, Mean fluorescence indensity; MSA, Maslinic acid; MTT, 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; NDAs, Non-esterified diterpene alkaloids; OD, Optical density; RMSD, Root-mean-square deviation; RTX, Resiniferatoxin; SD, Standard deviation; TRPV1, Transient receptor potential vanilloid type 1 channel; WB, Western blot.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics statement
Ethical approval was not required for the studies on animals in accordance with the local legislation and institutional requirements because only commercially available established cell lines were used.
Author contributions
LS: Data curation, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. SM: Software, Validation, Writing–review and editing. YZ: Formal Analysis, Methodology, Project administration, Validation, Visualization, Writing–review and editing. JnW: Conceptualization, Formal Analysis, Methodology, Validation, Writing–review and editing. ZL: Conceptualization, Data curation, Methodology, Writing–review and editing. HoW: Conceptualization, Data curation, Writing–review and editing. WL: Conceptualization, Data curation, Writing–review and editing. JaW: Conceptualization, Data curation, Writing–review and editing. WZ: Conceptualization, Data curation, Resources, Writing–review and editing. HD: Funding acquisition, Resources, Supervision, Writing–review and editing.
Funding
The authors declare financial support was received for the research, authorship, and/or publication of this article. National Natural Science Foundation of China (No. 82130113), and National Natural Science Foundation of China (No. 82274206).
Acknowledgments
This work is grateful to the above fundings.
Conflict of interest
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.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2024.1367682/full#supplementary-material
References
Acharya, C., Coop, A., Polli, J. E., and Mackerell, A. D. (2011). Recent advances in ligand-based drug design: relevance and utility of the conformationally sampled pharmacophore approach. Curr. Comput. Aided Drug Des. 7 (1), 10–22. doi:10.2174/157340911793743547
Ai, L., Yang, F., Hu, W., Guo, L., Liu, W., Xue, X., et al. (2023). Hepatotoxic components effect of Chebulae Fructus and associated molecular mechanism by integrated transcriptome and molecular docking. Mol. Basel, Switz. 28 (8), 3427. doi:10.3390/molecules28083427
Ai, X., Yan, J., Carrillo, E., and Ding, W. (2016). The stress-response MAP kinase signaling in cardiac arrhythmias. Rev. Physiol. Biochem. Pharmacol. 172, 77–100. doi:10.1007/112_2016_8
An, C., Qi, Z., and Su, X. (2020). Research advance in Terminalia antioxidants. Chem. Life 40 (06), 925–932. doi:10.13488/j.smhx.20200109
Bilige, W., Li, L., Chen, H., Meng, X., and Bolor, R. (2022). Effect of Terminalia Chebula Retz on serum metabolites in rats. Yunnan J. Traditional Chin. Med. Materia Medica 43 (08), 85–88. doi:10.16254/j.cnki.53-1120/r.2022.08.008
Burgess, G., and Williams, D. (2010). The discovery and development of analgesics: new mechanisms, new modalities. J. Clin. Invest. 120 (11), 3753–3759. doi:10.1172/jci43195
Celimuge, S., Xu, L., Song, L., and Tu, Y. (2022). Study on the mechanism of processing Radix Aconiti with hezi decoction based on UPLC-Q-TOF/MS technology and sequential metabolism method. Mod. Traditional Chin. Med. Materia Medica-World Sci. Technol. 24 (08), 2995–3013. doi:10.11842/wst.20210817002
Chan, Y. T., Wang, N., and Feng, Y. (2021). The toxicology and detoxification of Aconitum: traditional and modern views. Chin. Med. 16 (1), 61. doi:10.1186/s13020-021-00472-9
Chen, Q., Deng, X., Zhang, K., Kang, Y., Jiao, M., Zhang, J., et al. (2023a). Changes to PUFA-PPAR pathway during mesaconitine induced myocardial coagulative necrosis. Food Chem. Toxicol. 177, 113831. doi:10.1016/j.fct.2023.113831
Chen, Y., Tang, Y., Yang, L., Xiong, A., and Wang, Z. (2023b). Advances on toxicity and metabolism of alkaloids from the "toxic" traditional Chinese medicines. Acta Pharm. Sin. 2023, 1–22. doi:10.16438/j.0513-4870.2023-0241
Choi, J. K., Cho, W., Lee, J. H., Choi, G., and Park, M. (2020). A TRPV1 antagonist, PAC-14028 does not increase the risk of tumorigenesis in chemically induced mouse skin carcinogenesis. Regul. Toxicol. Pharmacol. 112, 104613. doi:10.1016/j.yrtph.2020.104613
Duarte, Y., Cáceres, J., Sepúlveda, R. V., Arriagada, D., Olivares, P., Díaz-Franulic, I., et al. (2020). Novel TRPV1 channel agonists with faster and more potent analgesic properties than capsaicin. Front. Pharmacol. 11, 1040. doi:10.3389/fphar.2020.01040
Duolan (2022). Protective effect of Terminalia chebula extract on H9C2 myocardialcell injury induced by aconitine and its mechanism. Inner Mongolia Medical University, 12–15.
Feinstein, W. P., and Brylinski, M. (2015). Calculating an optimal box size for ligand docking and virtual screening against experimental and predicted binding pockets. J. Cheminformatics 7 (1), 18. doi:10.1186/s13321-015-0067-5
Feng, A. (2019). Observation on the treatment of arrhythmia caused by severe acute aconitine poisoning. World Latest Med. Inf. 19 (94), 131–134. doi:10.19613/j.cnki.1671-3141.2019.94.082
Fu, M., Wu, M., Wang, J. F., Qiao, Y. J., and Wang, Z. (2007). Disruption of the intracellular Ca2+ homeostasis in the cardiac excitation-contraction coupling is a crucial mechanism of arrhythmic toxicity in aconitine-induced cardiomyocytes. Biochem. Biophys. Res. Commun. 354 (4), 929–936. doi:10.1016/j.bbrc.2007.01.082
Gao, Y., Cao, E., Julius, D., and Cheng, Y. (2016). TRPV1 structures in nanodiscs reveal mechanisms of ligand and lipid action. Nature 534 (7607), 347–351. doi:10.1038/nature17964
Gao, Y., Fan, H., Nie, A., Yang, K., Xing, H., Gao, Z., et al. (2022). Aconitine: a review of its pharmacokinetics, pharmacology, toxicology and detoxification. J. Ethnopharmacol. 293, 115270. doi:10.1016/j.jep.2022.115270
Gavva, N. R. (2008). Body-temperature maintenance as the predominant function of the vanilloid receptor TRPV1. Trends Pharmacol. Sci. 29 (11), 550–557. doi:10.1016/j.tips.2008.08.003
Gunthorpe, M. J., and Szallasi, A. (2008). Peripheral TRPV1 receptors as targets for drug development: new molecules and mechanisms. Curr. Pharm. Des. 14 (1), 32–41. doi:10.2174/138161208783330754
Guo, J., Xue, J., He, Z., Jia, H., and Yang, X. (2023). The mechanism by which Naru 3 pill protects against intervertebral disc cartilage endplate degeneration based on network pharmacology and experimental verification. J. Orthop. Surg. Res. 18 (1), 552. doi:10.1186/s13018-023-04014-x
Han, S., Bao, L., Han, X., Tang, Y., Liu, Z., Wang, H., et al. (2022a). Effects of characteristic components in Aconitii Kusnezoffii Radix processed with Chebulae Fructus on the cytotoxicity of H9c2 cardiomyocytes based on uniform design experiment. China Pharm. 25 (01). doi:10.19962/j.cnki.issn1008-049X.2022.01.005
Han, S., Bao, L., Li, W., Liu, K., Tang, Y., Han, X., et al. (2022b). Gallic acid inhibits mesaconitine-activated TRPV1-channel-induced cardiotoxicity. Evid. Based Complement. Altern. Med. 2022, 5731372. doi:10.1155/2022/5731372
Han, S., Bao, L., Liu, K., Han, X., Tang, Y., Liu, Z., et al. (2022c). Mechanism of Aconiti kusnezoffii radix processed with Chebulae Fructus against H9c2 cardiomyocyte toxicity based on TRPV1 channel. Chin. J. Exp. Traditional Med. Formulae 28 (01), 173–181. doi:10.13422/j.cnki.syfjx.20211748
Han, S., Liu, S., Zhang, Q., Zhi, M., Liu, K., Tang, Y., et al. (2020). Protective effect and mechanism of chebulagic acid on H9C2 cells injury induced by aconitine. J. Beijing Univ. Traditional Chin. Med. 43 (08), 645–652. doi:10.3969/j.issn.1006-2157.2020.08.006
Hansch, C. M., Peyton, P., Fujita, T., and Muir, R. M. (1962). Correlation of biological activity of phenoxyacetic acids with hammett substituent constants and partition coefficients. Nature 194 (4824), 178–180. doi:10.1038/194178b0
Li, C. Y., Zhou, Z., Xu, T., Wang, N. Y., Tang, C., Tan, X. Y., et al. (2022a). Aconitum pendulum and Aconitum flavum: a narrative review on traditional uses, phytochemistry, bioactivities and processing methods. J. Ethnopharmacol. 292, 115216. doi:10.1016/j.jep.2022.115216
Li, H., Liu, Y., Yang, H., and Liu, Y. (2022b). Research progress on chemical constituents, pharmacological mechanism, quality control and processing of Terminalia chebula Retz. Nat. Prod. Res. Dev. 34, 2130–2141. doi:10.16333/j.1001-6880.2022.12.017
Li, H., Song, L., Gao, Y., Chen, T., Hou, H., Peng, B., et al. (2023). Detoxification mechanism of combined Chebulae Fructus, Glycyrrhizae Radix et Rhizoma and prepared Aconiti Kusnezoffii Radix Cocta in regulating cardiac metabolic enzyme CYP2J3. Chin. J. Exp. Traditional Med. Formulae 29 (17), 88–95. doi:10.13422/j.cnki.syfjx.202202326
Li, H., Zhang, G., Chen, T., Gao, Y., Song, L., Hou, H., et al. (2021). Research progress on mechanism of reducing toxicity and increasing efficiency of compatibility of Mongolian medicinal materials Terminalia chebula and aconitum kusnezoffii. Word Chin. Med. 16 (24), 3701–3707. doi:10.3969/j.issn.1673-7202.2021.24.025
Li, H. Q., Xu, J. Y., Fan, X. H., and Wu, S. S. (2019). Optimization of the traditional processing method for precision detoxification of CaoWu through biomimetic linking kinetics and human toxicokinetics of aconitine as toxic target marker. J. Ethnopharmacol. 242, 112053. doi:10.1016/j.jep.2019.112053
Li, J., Liang, H., Cai, S., Li, Z., and Tu, Y. (2018). Mechanism of detoxification of Chebulae Fructus against Aconiti kusnezoffii radix toxicity based on network pharmacology. Acta Pharm. Sin. 53 (10), 1670–1679. doi:10.16438/j.0513-4870.2018-0506
Li, K., li, C., and Liu, X. (2011). Progress on the capsaicin receptor-TRPV1:Biological actions and as atherapeutic target for novel drugs. J. Shenyang Pharm. Univ. 28 (11), 917–927. doi:10.14066/j.cnki.cn21-1349/r.2011.11.014
Li, M., Xie, X., Chen, H., Xiong, Q., Tong, R., Peng, C., et al. (2020). Aconitine induces cardiotoxicity through regulation of calcium signaling pathway in zebrafish embryos and in H9c2 cells. J. Appl. Toxicol. 40 (6), 780–793. doi:10.1002/jat.3943
Li, X., and Yang, X. (2018). Study on the molecular targeting docking virtual screening of potential active components in tong-xie-yao-fang. J. Nanjing Normal Univ. Sci. Ed. 41 (03), 85–94. doi:10.3969/j.issn.1001-4616.2018.03.014
Liang, H., Li, J., Yang, Y., Tu, Y., and Song, L. (2017). Traditional Mongolian Medicine Chebulae Fructus protect arhythmiainduced by aconitine in rats. J. Med. Pharm. Chin. Minorities 23 (11), 75–77. doi:10.16041/j.cnki.cn15-1175.2017.11.036
Liu, J., Li, Y., Gong, X., Fang, W., and Sun, H. (2008). The protective effect of maslinic acid on the cardiac myocytes with ischemia in vitro. Chin. J. New Drugs 17 (09), 743–747. doi:10.3321/j.issn:1003-3734.2008.09.007
Liu, S., Li, F., Hou, Y., Yang, C., Tan, P., Li, X., et al. (2013). Influence of tannins from Chebulae Fructus on aconitum alkaloids of aconitum kusnezoffi processed with Chebulae fructus-principal of aconitum processed with Chebulae Fructus I. Chin. J. Exp. Traditional Med. Formulae 19 (05), 158–160. doi:10.13422/j.cnki.syfjx.2013.05.068
Liu, S., Li, F., Li, Y., Li, W., Xu, J., and Du, H. (2017). A review of traditional and current methods used to potentially reduce toxicity of Aconitum roots in Traditional Chinese Medicine. J. Ethnopharmacol. 207, 237–250. doi:10.1016/j.jep.2017.06.038
Liu, S., Li, W., Li, Y., Xu, j., Liu, X., Li, F., et al. (2016a). Lnfluence of Fructus Chebulae water solution on quality of Mongolian medicine Aconiti kusnezoffi radix praeparata. J. Beijing Univ. Traditional Chin. Med. 39 (11). doi:10.3969/j.issn.1006-2157.2016.11.012
Liu, S., li, Y., Li, W., Xu, J., Li, F., and Du, H. (2016b). Advances in studies on toxicity and modern toxicology of species in Aconitum L. Chin. Traditional Herb. Drugs 47 (22), 4095–4102. doi:10.7501/j.issn.0253-2670.2016.22.027
Liu, S., Liu, X., Li, Y., Yang, C., Liu, Q., Li, F., et al. (2015a). Experimental study on acute toxicity of the Aconiti kusnezo ffii radix and lts processed products. Chin. J. Pharmacovigil. 12 (09), 513–517. doi:10.19803/j.1672-8629.2015.09.001
Liu, S., Liu, X., Lin, S., Yang, C., Li, F., and Du, H. (2015b). Hydrolysis of simulation processing products of Aconitum kusnezoffii Processed with Terminalia chebula in simulated gastric fluid and intestinal fluid. China Pharm. 26 (13), 1752–1754. doi:10.6039/j.issn.1001-0408.2015.13.07
Liu, X., Tao, H., Tian, R., Huang, W., Zhang, T., Liu, Y., et al. (2023). Hezi inhibits Tiebangchui-induced cardiotoxicity and preserves its anti-rheumatoid arthritis effects by regulating the pharmacokinetics of aconitine and deoxyaconitine. J. Ethnopharmacol. 302, 115915. doi:10.1016/j.jep.2022.115915
Mazeto, T. K., Picada, J. N., Correa Á, P., Rebelo, I. N., Ribeiro, M. T., Gomez, M. V., et al. (2020). Antinociceptive and genotoxic assessments of the antagonist TRPV1 receptor SB-366791 on morphine-induced tolerance in mice. Naunyn Schmiedeb. Arch. Pharmacol. 393 (3), 481–490. doi:10.1007/s00210-019-01748-6
Meng, X., Huang, X., Dong, Y., Xu, X., Liu, S., Wu, Z., et al. (2023). Influence and mechanism of toxicity-efficacy relationship between the combination of Terminalia chebula Retz. and Aconitum pendulum Bush. On pentobarbital sodium-induced H9c2 myocardial cells. Chin. J. Mod. Appl. Pharm. 40 (08). doi:10.13748/j.cnki.issn1007-7693.20222649
Mi, S., Han, S., Liu, K., Bao, L., Zhang, F., Wang, H., et al. (2023). Mechanism of aconitum kusnezoffii processed with Terminalia chebula in reducing toxicity mediatedby TRPV1 based on pharmacophore and molecular docking. World Chin. Med. 18 (12), 1645–1652. doi:10.3969/j.issn.1673-7202.2023.12.001
More, K. N., Lee, Y. J., Kim, K. J., Suh, Y. G., Son, Y. J., and Chang, D. J. (2020). Effect of TRPV1 antagonist SC0030, a potent painkiller, on RANKL-mediated osteoclast differentiation involved in bone resorption. Bull. Korean Chem. Soc. 41 (4), 488–491. doi:10.1002/bkcs.11992
Nie, C. (2020). Chemical modification and analgesic activity of TRPV1 antagonist N-(4-tert-butylphenyl)-4 -(3-chloropyridine-2-yl)piperazine-1-formamide. Henan University, 22.
Nigam, M., Mishra, A. P., Adhikari-Devkota, A., Dirar, A. I., Hassan, M. M., Adhikari, A., et al. (2020). Fruits of Terminalia chebula Retz.: a review on traditional uses, bioactive chemical constituents and pharmacological activities. Phytother. Res. 34 (10), 2518–2533. doi:10.1002/ptr.6702
Pan, Y., Zhang, S., Li, F., Hou, J., and Zhou, X. (2011). Protective effects of decoction of Terminalia on myocardial cells damaged by Aconitine. J. Med. Pharm. Chin. Minorities 17 (01), 43–46. doi:10.16041/j.cnki.cn15-1175.2011.01.041
Peng, J., and Li, Y. (2010). The vanilloid receptor TRPV1: role in cardiovascular and gastrointestinal protection. Eur. J. Pharmacol. 627 (1-3), 1–7. doi:10.1016/j.ejphar.2009.10.053
Qiao, L., and Zhang, Y. (2014). Application of CADD on multi-target drug R&D in natural products. China J. Chin. Materia Medica 39 (11), 1951–1955. doi:10.4268/cjcmm20141101
Sappington, R. M., Sidorova, T., Long, D. J., and Calkins, D. J. (2009). TRPV1: contribution to retinal ganglion cell apoptosis and increased intracellular Ca2+ with exposure to hydrostatic pressure. Invest. Ophthalmol. Vis. Sci. 50 (2), 717–728. doi:10.1167/iovs.08-2321
Schaub, M. C., Hefti, M. A., and Zaugg, M. (2006). Integration of calcium with the signaling network in cardiac myocytes. J. Mol. Cell Cardiol. 41 (2), 183–214. doi:10.1016/j.yjmcc.2006.04.005
Singhuber, J., Zhu, M., Prinz, S., and Kopp, B. (2009). Aconitum in traditional Chinese medicine: a valuable drug or an unpredictable risk? J. Ethnopharmacol. 126 (1), 18–30. doi:10.1016/j.jep.2009.07.031
Song, Z., Chen, T., Wang, X., Wang, D., Shen, C., Wang, S., et al. (2023). Standardization of chebulagic acid as quality marker of Tibetan medicine Terminalia chebula. Analysis Test. Technol. Instrum. 29 (2), 140–148. doi:10.16495/j.1006-3757.2023.02.002
Sugimoto, Y., Kojima, Y., Inayoshi, A., Inoue, K., Miura-Kusaka, H., Mori, K., et al. (2013). K-685, a TRPV1 antagonist, blocks PKC-sensitized TRPV1 activation and improves the inflammatory pain in a rat complete Freund's adjuvant model. J. Pharmacol. Sci. 123 (3), 256–266. doi:10.1254/jphs.13088fp
Sumitra, M., Manikandan, P., Kumar, D. A., Arutselvan, N., Balakrishna, K., Manohar, B. M., et al. (2001). Experimental myocardial necrosis in rats: role of arjunolic acid on platelet aggregation, coagulation and antioxidant status. Mol. Cell Biochem. 224 (1-2), 135–142. doi:10.1023/a:1011927812753
Sun, G. B., Sun, H., Meng, X. B., Hu, J., Zhang, Q., Liu, B., et al. (2014a). Aconitine-induced Ca2+ overload causes arrhythmia and triggers apoptosis through p38 MAPK signaling pathway in rats. Toxicol. Appl. Pharmacol. 279 (1), 8–22. doi:10.1016/j.taap.2014.05.005
Sun, Z., Han, J., Zhao, W., Zhang, Y., Wang, S., Ye, L., et al. (2014b). TRPV1 activation exacerbates hypoxia/reoxygenation-induced apoptosis in H9C2 cells via calcium overload and mitochondrial dysfunction. Int. J. Mol. Sci. 15 (10), 18362–18380. doi:10.3390/ijms151018362
Szolcsányi, J., and Sándor, Z. (2012). Multisteric TRPV1 nocisensor: a target for analgesics. Trends Pharmacol. Sci. 33 (12), 646–655. doi:10.1016/j.tips.2012.09.002
Tan, X., He, Q., Pei, Z., Liu, Y., Feng, Z., Li, C., et al. (2023). Rapid visual characterization of alkaloid changes in traditional processing of Tibetan medicine Aconitum pendulum by high-performance thin-layer chromatography coupled with desorption electrospray ionization mass spectrometry imaging. Front. Pharmacol. 14, 1104473. doi:10.3389/fphar.2023.1104473
Temml, V., and Kutil, Z. (2021). Structure-based molecular modeling in SAR analysis and lead optimization. Comput. Struct. Biotechnol. J. 19, 1431–1444. doi:10.1016/j.csbj.2021.02.018
Wang, M., Zhang, S., and Zhai, H. (2002). The influence of terminalia on the toxicokinetics of aconitum decoction. J. Inn. Mong. Med. Univ. (04), 219–222. doi:10.16343/j.cnki.issn.2095-512x.2002.04.001
Wang, Z. (2019). Effects of corosolic acid on post-myocardial infarction cardiac remodeling in mice. Wuhan University, 27–30.
Wei, X., Qiu, Z., Chen, J., Sun, R., Huang, L., and Lai, C. (2019). Research advancement in mechanisms of processing and compatibility for detoxication of Aconitums. China J. Chin. Materia Medica 44 (17), 3695–3704. doi:10.19540/j.cnki.cjcmm.20190629.301
Wu, Z., He, E. Y., Scott, G. I., and Ren, J. (2015). α,β-Unsaturated aldehyde pollutant acrolein suppresses cardiomyocyte contractile function: role of TRPV1 and oxidative stress. Environ. Toxicol. 30 (6), 638–647. doi:10.1002/tox.21941
Xin, Y., Wang, S., Zhang, W., Lv, J., Zhang, Y., and Zhang, Z. (2014). Chemical detoxification of fruit of medicine Terminalia on Root of Kusnezoffii Monkshood by mass spectrometry. Chin. J. Exp. Traditional Med. Formulae 20 (23), 51–56. doi:10.13422/j.cnki.syfjx.2014230051
Yang, C., Li, F., hOU, y., Liu, S., Tan, P., Li, X., et al. (2013). Comparision of contents of aconitum alkaloids between the decoctions of aconitum processed with myrobalan and co-decoctions of aconitum and myrobalan - principal of aconitum processed with myrobalan Ⅱ. Chin. J. Exp. Traditional Med. Formulae 19 (04), 130–132. doi:10.13422/j.cnki.syfjx.2013.04.051
Yang, C., Zeng, X., Cheng, Z., Zhu, J., and Fu, Y. (2021). Aconitine induces TRPV2-mediated Ca(2+) influx through the p38 MAPK signal and promotes cardiomyocyte apoptosis. Evid. Based Complement. Altern. Med. 2021, 9567056. doi:10.1155/2021/9567056
Yang, M., Ji, X., and Zuo, Z. (2018). Relationships between the toxicities of radix Aconiti lateralis preparata (fuzi) and the toxicokinetics of its main diester-diterpenoid alkaloids. Toxins (Basel) 10 (10), 391. doi:10.3390/toxins10100391
Yang, S. Y. (2010). Pharmacophore modeling and applications in drug discovery: challenges and recent advances. Drug Discov. Today 15 (11-12), 444–450. doi:10.1016/j.drudis.2010.03.013
Ye, Q., Liu, H., Fang, C., Liu, Y., Liu, X., Liu, J., et al. (2021). Cardiotoxicity evaluation and comparison of diterpene alkaloids on zebrafish. Drug Chem. Toxicol. 44 (3), 294–301. doi:10.1080/01480545.2019.1586916
Yin, Y. (2019). Structural mechanism underlying TRPV1 activation by pungent compounds in gingers. Qingdao University, 18–19.
Yuan, X., Gao, Y., Wang, Q., Wang, L., Zhang, C., Gu, C., et al. (2022). Study on the pharmacological effects of the compatibility of Mongolian medicine Aconiti kusnezoffii radix and Chebulae Fructus on the kidney and liver toxicity caused by Aconiti kusnezoffii radix. J. Inn. Mong. Med. Univ. 44 (06), 581–584+590. doi:10.16343/j.cnki.issn.2095-512x.2022.06.006
Zhai, K., Liskova, A., Kubatka, P., and Büsselberg, D. (2020). Calcium entry through TRPV1: a potential target for the regulation of proliferation and apoptosis in cancerous and healthy cells. Int. J. Mol. Sci. 21 (11), 4177. doi:10.3390/ijms21114177
Zhang, X. (2019). The mechanism of detoxification and retention effect of Terminalia Chebula on Radix Aconiti Kusnezoffii based on pharmacokinetics. Inner Mongolia Medical University, 45–46.
Zhang, X., Cui, Y., Miao, X., Liu, D., Ma, Z., and Li, G. (2017). Study on the protective effect of Mongolian medicine Chebulae Fructus on cardiac toxicity induced by Aconiti Kusnezoffii Radix in rats. J. Chin. Med. Mater. 40 (11), 2693–2696. doi:10.13863/j.issn1001-4454.2017.11.044
Zhi, M., Gu, X., Han, S., Liu, K., Liu, Z., Tang, Y., et al. (2020a). Chemical variation in aconti kusnezoffii radix before and after processing based on UPLC-orbitrap-MS. China J. Chin. Materia Medica 45 (05), 1082–1089. doi:10.19540/j.cnki.cjcmm.20191221.301
Zhi, M., Liu, K., Han, S., Xu, J., Li, W., Li, F., et al. (2020b). Influence of different dosage forms on pharmacokinetics of 6 alkaloids in raw Aconiti kusnezoffii radix (Caowu) and Chebulae Fructus- (hezi-) processed Caowu by UPLC-MS/MS. Biomed. Res. Int. 2020, 1942849. doi:10.1155/2020/1942849
Zhou, B., Qian, Z., Li, Q., Gao, Y., and Li, M. (2022). Assessment of pulmonary infectious disease treatment with Mongolian medicine formulae based on data mining, network pharmacology and molecular docking. Chin. Herb. Med. 14 (3), 432–448. doi:10.1016/j.chmed.2022.07.001
Zhou, W., Liu, H., Qiu, L. Z., Yue, L. X., Zhang, G. J., Deng, H. F., et al. (2021). Cardiac efficacy and toxicity of aconitine: a new frontier for the ancient poison. Med. Res. Rev. 41 (3), 1798–1811. doi:10.1002/med.21777
Keywords: Chebulae Fructus triterpenoids, mesaconitine, TRPV1, cardiotoxicity, pharmacophore model, molecular docking, 3D-QSAR model
Citation: Song L, Mi S, Zhao Y, Liu Z, Wang J, Wang H, Li W, Wang J, Zu W and Du H (2024) Integrated virtual screening and in vitro studies for exploring the mechanism of triterpenoids in Chebulae Fructus alleviating mesaconitine-induced cardiotoxicity via TRPV1 channel. Front. Pharmacol. 15:1367682. doi: 10.3389/fphar.2024.1367682
Received: 09 January 2024; Accepted: 19 February 2024;
Published: 04 March 2024.
Edited by:
Qianfeng Gong, Jiangxi University of Traditional Chinese Medicine, ChinaReviewed by:
Zhipeng Qu, University of Adelaide, AustraliaGang Fan, Chengdu University of Traditional Chinese Medicine, China
Copyright © 2024 Song, Mi, Zhao, Liu, Wang, Wang, Li, Wang, Zu and Du. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong Du, duhong@vip.163.com