Skip to main content

ORIGINAL RESEARCH article

Front. Pharmacol., 26 November 2021
Sec. Ethnopharmacology

Interactions Between Ephedra sinica and Prunus armeniaca: From Stereoselectivity to Deamination as a Metabolic Detoxification Mechanism of Amygdalin

Updated
Yan Qin&#x;Yan Qin1Shanshan Wang&#x;Shanshan Wang2Qiuyu Wen,Qiuyu Wen1,2Quan Xia,Quan Xia2,3Sheng WangSheng Wang4Guanjun ChenGuanjun Chen4Jiayin SunJiayin Sun1Chenlin Shen,
Chenlin Shen1,5*Shuai Song,
Shuai Song2,3*
  • 1Inflammation and Immune Mediated Diseases Laboratory of Anhui Province, School of Pharmacy, Anhui Institute of Innovative Drugs, Institute for Liver Diseases of Anhui Medical University, Anhui Medical University, Hefei, China
  • 2Department of Pharmacy, The First Affiliated Hospital of Anhui Medical University, Hefei, China
  • 3The Grade 3 Pharmaceutical Chemistry Laboratory of State Administration of Traditional Chinese Medicine, Hefei, China
  • 4Center for Scientific Research of Anhui Medical University, Hefei, China
  • 5Hefei Kaifan Analytical Technology Co., Ltd., Hefei, China

Mahuang–Xingren (MX, Ephedra sinica Stapf-Prunus armeniaca L.) is a classic herb pair used in traditional Chinese medicine. This combined preparation reduces the toxicity of Xingren through the stereoselective metabolism of its main active ingredient amygdalin. However, whether stereoselectivity is important in the pharmacokinetic properties of amygdalin either in the traditional decoction or in the dispensing granules is unclear. Amygdalin is hydrolyzed to its metabolite, prunasin, which produces hydrogen cyanide by degradation of the cyano group. A comprehensive study of the metabolic pathway of amygdalin is essential to better understand the detoxification process. In this article, the potential detoxification pathway of MX is further discussed with regard to herb interactions. In this study, the pharmacokinetic parameters and metabolism of amygdalin and prunasin were investigated by comparing the traditional decoction and the dispensing granule preparations. In addition, several potential metabolites were characterized in an incubation system with rat liver microsomes or gut microbial enzymes. The combination of Xingren with Mahuang reduces exposure to D-amygdalin in vivo and contributes to its detoxification, a process that can be further facilitated in the traditional decoction. From the in vitro co-incubation model, 15 metabolites were identified and classified into cyanogenesis and non-cyanogenesis metabolic pathways, and of these, 10 metabolites were described for the first time. The level of detoxified metabolites in the MX traditional decoction was higher than that in the dispensing granules. The metabolism of amygdalin by the gut microbial enzymes occurred more rapidly than that by the rat liver microsomes. These results indicated that combined boiling both herbs during the preparation of the traditional decoction may induce several chemical changes that will influence drug metabolism in vivo. The gut microbiota may play a critical role in amygdalin metabolism. In conclusion, detoxification of MX may result 1) during the preparation of the decoction, in the boiling phase, and 2) from the metabolic pathways activated in vivo. Stereoselective pharmacokinetics and deamination metabolism have been proposed as the detoxification pathway underlying the compatibility of MX. Metabolic detoxification of amygdalin was quite different between the two combinations, which indicates that the MX decoctions should not be completely replaced by their dispensing granules.

Introduction

The combination of Ephedra sinica Stapf (Mahuang) and Prunus armeniaca L. (Xingren) is a classical herb combination used to treat cold-induced asthma in traditional Chinese medicine (TCM) (Ma et al., 2021). According to the Guidance for Corona Virus Disease 2019: Prevention, Control, Diagnosis and Management published in China, MX can be combined with other herbs, such as Glycyrrhiza uralensis Fisch. (Gancao), Astragalus mongholicus Bunge (Huangqi), and Atractylodes lancea (Thunb.) DC. (Cangzhu) has been shown to be effective against the new coronavirus and has significantly reduced the proportion of patients with mild symptoms transforming into more severe cases (Zhao et al., 2021a; Zhang et al., 2021). Ephedrines in Mahuang and amygdalin of Xingren have been considered as the main active constituents possessing anti-asthmatic (Wang et al., 2021), anti-inflammatory (Karimi et al., 2021), and antitussive properties (Wang et al., 2019). The main active ingredient of Xingren is amygdalin (3% in Xingren) (Zhang et al., 2020). Amygdalin is metabolized by β-glucosidase, found in plants or the animal gut, into glucose and prunasin. In turn, prunasin, a cyanogenic diglycoside, is metabolized into glucose and mandelonitrile, which are further hydrolyzed to benzaldehyde and hydrogen cyanide. Hydrogen cyanide is rapidly absorbed by the gastrointestinal tract and is distributed throughout the body by the blood circulation after absorption. Hydrogen cyanide is responsible for the antitussive properties of Xingren. However, the cyanogenesis metabolism pathway results in the release of hydrogen cyanide after oral administration of amygdalin and is a well-known mechanism for metabolic poisoning (Cressey and Reeve, 2019).

The essential and common characteristics of TCM preparations are defined by specific formulations that are developed in accordance with the compatibility theory. Different herbal preparations are combined to obtain multi-herb drug cocktails of TCM (Hu et al., 2005). The underlying principle of these formulations involves not only a simple mixture of botanical drugs but also processes able to reduce toxicity and increase treatment efficiency. The MX decoction is composed of Mahuang and Xingren, at a 1:1 dosage combination as described in the “Treatise on Cold Pathogenic and Miscellaneous Diseases.” MX was initially included in Shennong’s Classic of Materia Medica and ranked as a medium grade drug, which means it can present intermediate toxic effects (Dai et al., 2019). Some adverse reactions and even fatal poisoning have been reported in the literature (Blaheta et al., 2016; Han et al., 2018). In our previous study (Song et al., 2016), the contents and pharmacokinetics of alkaloids such as ephedrine and methylephedrine in MX were studied. Mahuang antagonizes the acute toxicity of Xingren and significantly improves the safety profile of the MX combination. Stereoselective pharmacokinetics and cytotoxicity data indicate that the antagonism properties of D-amygdalin metabolism may contribute to its detoxification. Whether the stereoselectivity of amygdalin pharmacokinetics depends on the traditional decoction or on the combination of the individual herbs is still unclear.

The traditional decoction is the most common form of administration of TCM (Wang et al., 2018). Herbal materials are often extracted with water to form an aqueous extract (Wang et al., 2012). Dispensing granules, consisting of a single herb preparation in a single formula, have become increasingly popular in the TCM clinic. The granules not only share the benefits of a customized formulation but also allow convenient storage and quality control (Denesh et al., 2021). However, it is not clear whether stereoselective pharmacokinetics are relevant in MX dispensing granule preparation as there is still discussion on the metabolic pathways of amygdalin metabolism in vivo. Amygdalin is a cyanogenic glucoside naturally produced by plants (Kovacikova et al., 2019). Metabolic detoxification refers to the process of reducing exposure to harmful xenobiotics in vivo which could also contribute to its pharmacokinetic processes (Ren et al., 2020), including the first-pass effect and phase II conjugation. However, the metabolic detoxification pathways of amygdalin have not been well-clarified. Furthermore, there have also been no reports describing how Mahuang could influence the cyanogenic effects induced by amygdalin.

Ultrahigh-performance liquid chromatography coupled with high-resolution mass spectrometry (UPLC-HRMS) is a sensitive and potentially more robust technique than nuclear magnetic resonance and gas chromatography–mass spectrometry (GC-MS) for the identification of metabolites present in low concentrations and has better capacity to control the false discovery rate. A reliable and stable LC-MS/MS method is often used for the quantification of drugs and their metabolites in different matrices (Amer et al., 2017; Attwa et al., 2018a). In this study, a precise, accurate, and fast UPLC–tandem mass spectrometry (UPLC-MS/MS) method was used to measure the concentration of amygdalin and prunasin in rat plasma after oral administration. UPLC-HRMS was used to identify the metabolites of amygdalin in vitro and in vivo. We provide a comparative analysis of the stereoselective pharmacokinetics of amygdalin of the traditional decoction and dispensing granule formulation to clarify MX compatibility. Identification of metabolites and possible detoxification pathways of amygdalin in vivo was evaluated by incubating the preparation with rat liver microsomes and gut microbial enzymes. The possible effects of Mahuang on the detoxification pathway are also discussed later.

Materials and Methods

Chemicals and Materials

Amygdalin (99% purity, purchased from Sigma-Aldrich, Co.), prunasin (98% purity, purchased from Wuhan ChemFaces Biochemical Co., Ltd.), and benzaldehyde (98.5% purity, purchased from Shanghai Macklin Biochemical Co., Ltd.) were used as reference chemicals in the pharmacokinetics study. Diazepam (98% purity) and geniposide (97.5% purity), provided by the National Institutes for Food and Drug Control (Beijing, China), were selected as internal standards (IS) for studies of amygdalin pharmacokinetics and metabolism, respectively. Xingren (Gansu Province, batch number: 20170901) and Mahuang (Inner Mongolia Autonomous Region, batch number: 20180901) were purchased from Zisun Medicinal Material Co., Ltd (Guangzhou, China) and were authenticated by Professor Qunlin Zhang (School of Pharmacy, Anhui Medical University). The dispensing granules of Mahuang (1 g, equivalent to 5 g of the crude drug, batch number: 20170901) and Xingren (1 g, equivalent to 10 g of the crude drug, batch number: 20170901) were obtained from E-FONG Pharmaceutical Co. (Guangzhou, China). The inspection report of dispensing granules was provided by the E-FONG Pharmaceutical Co. (Supplementary Figure S1). Acetonitrile, methanol, and formic acid (HPLC grade) were purchased from Tedia Company, Inc (Fairfield, OH, United States). Ammonium acetate (HPLC grade) and SPE cartridges (ProElut, PLS 30 mg/ml) were provided by Dikma Technologies Inc (Lake Forest, CA, United States). Ultrapure water (18.2 MΩ) was purified by the Millipore system (Millipore Corp, Billerica, MA, United States). Ultrasonic cell disruption (JY92-II) was purified using a Scientz homogenizer. All reagents were analytical grades, unless otherwise stated. The BCA Protein Assay Kit for protein concentration quantification was purchased from Biosharp (Hefei, Anhui, China).

Preparation of Herb Extracts

According to an extraction method described in TCM, Mahuang (90 g) and Xingren (90 g) were soaked in a moderate amount of water for 30 min and then extracted by boiling in 10-fold volumes of water (1800 ml) for 60 min. After filtering with absorbent gauze, an 8-fold volume of water was subsequently added to the residues and boiled for 90 min. The solution was filtered with five-layer absorbent gauze. The two filtrates were combined and concentrated to obtain the MX decoction. The obtained solution was concentrated to 180 ml (equal to 1 g crude herb/ml) using a rotary evaporator. The Xingren decoction was prepared following the same procedure. Finally, a 180 ml solution (equal to 0.5 g crude herb/ml) of Xingren (90 g) was obtained. The dispensing granule solution was prepared by dissolving the formula granules directly in water. Mahuang and Xingren were mixed in a 1:1 ratio (w/w, crude drug). We determined the concentrations of two main compounds (D-amygdalin and neoamygdalin) present in the Xingren decoction, MX decoction, Xingren granules, and MX granules (Supplementary Figure S2), the contents of which are shown in Supplementary Table S1, and other related components were characterized in a follow-up study.

UPLC-MS/MS Conditions for Quantification

LC-MS/MS analysis was performed using a Shimadzu UPLC system (Shimadzu, Kyoto, Japan) coupled with a Triple Quada™ 5500 equipped with an electrospray ionization source (AB, Milwaukee, Wisconsin, United States). The UPLC system includes a DGU-20A5R online vacuum degasser, an LC-30AD pump, a SIL-30AC automatic sampler, a CBM-20A system controller, and a CTO-20AC column oven. Separation was carried out in a CDshell-RSP column (3.0 × 150 mm, 2.7 μm; Arlington, TX, United States), a novel chiral stationary phase composed of hydroxypropylated-β-cyclodextrin. The mobile phase at a flow rate of 0.4 ml/min consisted of phases A (5 mM ammonium acetate solution containing 0.1% formic acid) and B (methanol). Gradient separation was carried out by phase B of 12% during the first 2.5 min and then increasing from 12 to 95% during the following 4 min. Subsequently, phase B was linearly reduced to the initial ratio from 6.5 to 8 min. The temperatures of the column and the autosampler were fixed at 35°C and 4°C, respectively. A 5 µl sample was loaded into the system.

Tandem MS was equipped with an electrospray ionization interface operating in the positive ion mode. For the optimization of each multiple reaction monitoring (MRM) parameter, other normal LC-MS/MS conditions were used, including ion source parameters (spray voltage: 5500 V; atomization temperature: 500°C; nebulizer pressure: 35 psi). The optimized quantitative ion pair (m/z) and qualitative ion pair selected for the MRM mode used the following optimized ion transitions for each precursor/product pair: amygdalin (457.3→163.3, 90 V, 22 eV), prunasin (313.0→163.1, 90 V, 14 eV), and diazepam (285.2→228.0, 140 V, 30 eV). AB-SCIEX Analyst (version 1.6.3) software was used for data acquisition and processing.

Preparation of Standard and Quality Control Samples

Standard solutions were obtained by mixing and serial dilution of stock solutions of amygdalin (914.8 μg/ml) and prunasin (1,181.4 μg/ml). The diazepam (IS, 10.0 ng/ml) was prepared by dilution of its stock solution (1.0 mg/ml). Acidic conditions are required to prevent its epimerization. Next, stock solutions and further dilutions were prepared in 50 and 10% aqueous methanol solution (v/v, containing 0.1% formic acid), respectively. Calibration curves and quality control (QC) samples were prepared by spiking 5 μl of working standard with 95 μl blank rat plasma. The calibration curve of amygdalin and prunasin was prepared in the concentration range between 1.0 and 128.1 ng/ml and 4.6 and 590.7 ng/ml, respectively. Three levels of QC samples were prepared containing amygdalin (2.2, 22.0, and 109.8 ng/ml) and prunasin (9.13, 91.3, and 456.5 ng/ml). All solutions were stored at 4°C and brought to room temperature prior to use.

Preparation of Rat Plasma Samples

For the pharmacokinetics study, 100 μl rat plasma was transferred into a 1.5-ml polypropylene centrifuge tube, followed by the addition of 400 μl of diazepam (IS, 10.0 ng/ml). The samples were mixed and loaded into an SPE cartridge which was preconditioned with methanol and deionized water. The cartridge was washed with 1 ml of 10% aqueous methanol solution containing 0.1% formic acid and vacuum-dried for 1 min. The analytes were then eluted with 1 ml of methanol containing 0.1% formic acid. The eluent was evaporated to dryness under a gentle stream of nitrogen gas. The residue was resuspended with 400 µl of 5 mM ammonium acetate containing 0.1% formic acid and centrifuged (15,000×g for 8 min) for quantitative analysis.

Analytical Method Validation

The method was validated for specificity, linearity, sensitivity, accuracy and precision, matrix effect, extraction recovery, and stability of the analytes in rat plasma (Commission, 2015).

Pharmacokinetics of Amygdalin in Traditional Decoction and Dispensing Granule

Male Sprague–Dawley rats, of specific-pathogen-free grade, were obtained from the Anhui Medical University Laboratory Animal Center (Certificate number SCXK 2017–001). All animal experiments were approved by the Experimental Animal Ethics Committee of Anhui Medical University (LLSC20170348) and complied with the requirements of the Animals Ethics Procedures and Guidelines of the People’s Republic of China. They were housed under controlled environmental conditions (ambient temperatures 23–25°C; 12-h light/12-h dark cycles, 45–55% relative humidity) and fed ad libitum for the first 2 weeks. Rats were fasted for 12 h with free access to water prior to and during the experiment.

The same dose of crude extract (3 g/kg) was administered orally to rats in these four groups, including the Xingren decoction (equivalent to 44.1 mg/kg D-amygdalin and 40.2 mg/kg neoamygdalin), Xingren dispensing granules (equivalent to 45.3 mg/kg D-amygdalin and 50.7 mg/kg neoamygdalin), MX decoction (equivalent to 38.4 mg/kg D-amygdalin and 14.5 mg/kg neoamygdalin), and MX dispensing granules (equivalent to 24.3 mg/kg D-amygdalin and 25.7 mg/kg neoamygdalin). A 500 μl sample of orbital venous blood was collected in heparinized 1.5-ml polythene tubes before the treatment and subsequently at 5 min, 0.25, 0.5, 0.75, 1, 1.5, 2, 3, 4, 5, 6, 9, and 12 h after dosing. Plasma samples were separated immediately from blood by centrifugation at 4,500×g for 10 min and stored at −80°C until analysis.

Extraction of Gut Microbial Enzyme

Rat fecal specimens (approximately 0.5 g each) were dissolved in 10-fold PBS (m/v) and centrifuged (3,500×g for 15 min at 4°C) to remove the supernatant. The precipitate was washed twice and resuspended in cold PBS mentioned earlier. The resulting suspension was sonicated (300W, 10 s over 20 s intervals) for 45 min in an ice bath to rupture bacterial cells and then centrifuged at 13,000×g for 30 min at 4°C to obtain a suspension containing microbial enzymes and to remove debris (Sun et al., 2020), and aliquots of the supernatant were subsequently stored in a refrigerator until use.

UPLC-HRMS Conditions for Metabolite Separation and Detection

A UHPLC Dionex Ultimate 3,000 (Thermo Scientific, San Jose, United States) equipped with a cooling autosampler and column oven was utilized. Separation was carried out on an HSS T3 column (2.1 × 100 mm, 1.8 µm, Waters, Ireland) with a column temperature maintained at 45°C. Following the gradient elution program with the same mobile phase composition was established: 0–3 min, 5% B; 3–5 min, 5–25% B; 5–11 min, 25% B; 11–13.8 min, 60% B; 13.9–16 min, 5% B. The injection volume was 2 µl, and the flow rate was 0.25 ml/min.

A Q-Exactive plus hybrid quadrupole-orbitrap mass spectrometer (Thermo Scientific, San Jose, United States) with a heat electrospray ionization (HESI) probe was employed using a data-dependent MS/MS approach in a negative or positive ion mode. The mass conditions were optimized as follows: capillary temperature, 320°C; spray voltage, −3.5 kV; S-lens RF level, 50 V; auxiliary gas heater temperature, 200°C; sheath gas pressure, 40 psi; auxiliary gas pressure, 10 psi; scan range, 100–1,000 m/z; resolution, 17,500; solation window: 2.0 m/z; stepped normalized collision energy: 20, 40, 60; dynamic exclusion, 10 s; and apex trigger, 2–6 s. Nitrogen was used for spray stabilization, for high-energy collision dissociation, and as the damping gas in the C-trap. Q-Exactive 2.9, Xcalibur 4.1, Compound Discoverer 3.0, and MS Frontier 7.0 software (Thermo Fisher Scientific, San Jose, United States) were used for instrument control, data acquisition, and analysis.

Metabolic Transformations of Amygdalin in the Traditional Decoction and on Exposure to Rat Liver Microsomes and Gut Microbial Enzymes

In the TCM theory, the proper combination of herbs is considered essential to enhance therapeutic effects and reduce toxic side effects. Herb–herb interactions, inducing changes in the solubility or generation of a new compound, may occur during the decoction process. Liver and gut microbes play an important role in drug metabolism (Knudsen et al., 2021; Sodhi and Benet, 2021). Thus, chemical transformations or endogenous metabolism of amygdalin was investigated during decoction and incubation of the rat liver microsome (RLM) and gut microbial enzyme (RGME) preparations, respectively.

In this study, RLM was prepared according to the hypothermal differential centrifugation method, as described in the literature (Cheng et al., 2019). The liver tissue was homogenized with four volumes (w/v) of ice-cold PBS buffer. The homogenate was centrifuged at 9,000 g for 20 min at 4°C. Then the supernatant was centrifuged at 100,000×g for 60 min at 4°C, and the pink precipitate was used as the liver microsome fraction. The protein concentration of RLM and RGME was estimated using the BCA and subsequently aliquoted and stored in a refrigerator until use. The incubation solution for RLM consisted of RLM (0.5 mg/ml), 50 mM PBS (pH = 7.4), and the NADPH regeneration system (containing 1.5 mM NADP +, 5 mM glucose-6-phosphate, 5 mM MgCl2, and 1 unit/ml G6PDH). The reaction was started by adding 10 μl NADPH. The incubation preparation for RGME consisted of PBS, D-amygdalin, and RGME and was carried out at 37°C for 2 h, and the reaction was terminated by adding 3-fold ice-cold acetonitrile. Furthermore, 20 μM (42 μg/ml) D-amygdalin was added to an already prewarmed incubation system (800 μl).

All samples, including the incubation mixture of the traditional decoction and the dispensing granule solution, were filtered using a 0.22-µm nylon filter after a tenfold dilution with water, after which the samples were purified with the SPE protocol, as described before for the preparation rat plasma samples. The supernatant was obtained by centrifugation at 16,000×g for 10 min, and the residues were dried under gentle nitrogen steam. The residues were suspended with 400 μl initial mobile phase and then filtered through a syringe with a 0.22-µm nylon membrane. Principal component analysis (PCA) (Peterson et al., 2016) was performed to evaluate the changes in metabolite composition associated with amygdalin at a normalized concentration (40 μM, 84 μg/ml) of the decoction and dispensing granules.

Data Handing and Statistical Analysis

Pharmacokinetic parameters including half-life (t1/2), maximum plasma concentration (Cmax), the time for maximum concentration (Tmax), the area under the curve (AUC0–t), the volume of distribution (Vz), and the clearance rate (Hulin et al., 2019) were estimated by noncompartmental analysis using Phoenix WinNonlin software (version 8.3.1.5014, Pharsight, United States).

Generally, metabolites with the same or similar parent nucleus always have similar cleavage fragments in the ESI-MSn spectrum, so it is possible to identify metabolites and even the corresponding metabolic pathways. First, the MS1 (primary level of MS) and MS2 (secondary level of MS) spectrum data of amygdalin were collected to determine its MS/MS fragment behavior (precursor–product ion relationships). Next, Compound Discoverer 3.0 was employed to analyze the metabolites of amygdalin and generated extracted ion chromatograms (XICs) with retention time, accurate mass, and isotopic pattern. Considering a significant mass intensity, a suitable tolerance window for XIC extraction (±5 ppm), and established fragment ions, the structure of possible metabolites of amygdalin could be identified. Finally, the fragmentation behavior analysis predicted by Mass Frontier software was combined with the analysis of known metabolic pathways to assess whether the identification of metabolites was reasonable or not.

Variations in the chemical profile of amygdalin between the traditional decoction and dispensing granules at a normalized concentration of Xingren were explored by principal component analysis (PCA). All data are presented as mean ± SD. Statistical analysis was performed using the single-factor ANOVA. A p value less than 0.05 was considered statistically significant. Statistical analysis was performed using SPSS v.20 (IBM Corp., Armonk, NY, United States).

Results

Method Validation

The validated method showed good sensitivity, selectivity, and reproducibility. The XICs of blank rat plasma, blank plasma spiked with IS and analytes at LLOQ, and plasma samples after administration of traditional decoction and dispensing granules are shown in Supplementary Figure S3. No obvious interference was found in the retention time. Eight-point calibration plots exhibited satisfactory linearity (correlation coefficient >0.99) in the concentration range of 1.00–128.1 and 4.62–590.7 ng/ml for amygdalin and prunasin (Supplementary Table S2). Reproducible extraction recovery (79.6–90.3%, CV ≤ 7.01%) and the neglectable matrix effect (112.6–114.7%, CV ≤ 8.74%) were observed (Table 1). Acceptable intraday and intraday precision and accuracy were demonstrated (88.3–113.0%; CV ≤ 13.3%) for matrix QC samples (Supplementary Table S3). In addition, both amygdalin and prunasin were quite stable under current storage and processing conditions, with variability ranging from −5.7 to 14.7% (CV ≤ 14.3%) (Supplementary Table S4).

TABLE 1
www.frontiersin.org

TABLE 1. Recovery and matrix effect of amygdalin and prunasin in rat plasma using solid-phase extraction (n = 6).

Bioequivalence of Dispensing Granule Was not Consistent With the Traditional Decoction

The mean plasma concentration–time profiles (n = 5) are shown in Figure 1. Plasma concentration versus time was simulated and followed the first-order rate law. Therefore, dose-normalized Cmax and AUC0-t (Cmax/D and AUC0-t/D) were used to assess drug exposure based on the different concentrations of amygdalin in the decoction or the dispensing granule. Prunasin, presenting a cleaved glucose molecule, was the primary deglycosylated metabolite of amygdalin. Both amygdalin (including D-amygdalin and neoamygdalin) and prunasin (including D-prunasin and sambunigrin) (Supplementary Figure S2) were used to estimate the pharmacokinetic properties in the following analysis.

FIGURE 1
www.frontiersin.org

FIGURE 1. Mean plasma concentration profiles over time for (1) D-amygdalin, (2) neoamygdalin, (3) D-prunasin, and (4) sambunigrin in rat (mean ± SD, n = 5) following oral administration of Xingren, Mahuang–Xingren decoction, and dispensing granule preparations.

No significant differences in Tmax or t1/2 of amygdalin and prunasin (except the shorter Tmax of prunasin in the dispensing granule preparation) were observed in the MX combination. However, the Cmax/D and AUC0-t/D of amygdalin (except for the dispensing granule) and prunasin were all improved with the addition of Mahuang (p < 0.05 or p < 0.01) (Tables 2, 3). Higher plasma exposure (presented as Cmax/and AUC) of amygdalin and prunasin was achieved in the MX decoction (p < 0.05 or p < 0.01). In addition, the clearance (p < 0.01) and the biodistribution of amygdalin (p < 0.05) were markedly lower than for those of the dispensing granule preparation (Supplementary Tables S5, S6). In summary, Mahuang increased the bioavailability of total amygdalin, which was consistent with our previous report. Meanwhile, reduced drug clearance and biodistribution were only confirmed for the decoction. The data suggested that the traditional decoction was superior in terms of its lower toxicity and also showed better therapeutic efficacy than the MX dispensing granule preparation.

TABLE 2
www.frontiersin.org

TABLE 2. Pharmacokinetic parameters of D-amygdalin, neoamygdalin, D-prunasin, and sambunigrin in rats after oral administration of Xingren and MX decoction.

TABLE 3
www.frontiersin.org

TABLE 3. Pharmacokinetic parameters of D-amygdalin, neoamygdalin, D-prunasin, and sambunigrin in rats after oral administration of Xingren and MX dispensing granules.

D-amygdalin appears to be better absorbed than neoamygdalin due to its higher Cmax/D and AUC0-t/D (Supplementary Table S7). Previous studies have also indicated that antagonism of D-amygdalin (i.e., amygdalin) metabolism may be critical for detoxification in the MX combination. The ratio of D-amygdalin to neoamygdalin was considered an exposure index of D-amygdalin in vivo. The variations of this ratio were compared for the MX combination and the single herbs. As shown in Supplementary Table S7, the ratio of AUC0-t/D in the MX decoction decreased from 7.19 to 5.72 (p < 0.05) (by approximately 20.4%). On the contrary, no significant decrease was observed for the dispensing granule preparations. These findings indicated that Mahuang decreased the risk of exposure relative to D-amygdalin and reduced its toxicity compared to traditional decoction. In general, better bioavailability of D-amygdalin, neoamygdalin, D-prunasin, or sambunigrin in the decoction was observed than the dispensing granules, which indicated that in vivo drug disposition of D-amygdalin differed significantly for the two combinations (Supplementary Tables S5, S6). MX decoctions should not be completely replaced by the dispensing granule formulation in terms of toxicity.

Identification of Metabolites and the Metabolic Detoxification Pathway of Amygdalin

By comparing the retention time, exact mass and possible MS/MS fracture behavior of metabolites, a total of 16 different compounds were identified or tentatively characterized (Table 4). The fragmentation patterns of the 16 compounds are shown in Supplementary Figure S4. The retention time of amygdalin amide was 7.14/7.24 min, and the MS spectrum showed [M + HCOO]- ion at m/z 520.16608 (C20H27NO11, 0.58 ppm). According to the representative product ion (m/z 474.16170, 323.09837, 179.05611), it could be inferred that the nitrile group was introduced into H2O. Therefore, it was presumed to be the amidated product of amygdalin epimers.

TABLE 4
www.frontiersin.org

TABLE 4. Identification of amygdalin-related metabolites in gut microbial enzyme incubation system by HR-MS.

The drug incubation test indicated that amygdalin can be deglycosylated by β-glucosidase to yield prunasin, mandelonitrile, benzaldehyde, and laetrile, which is recognized as a classical metabolism pathway in many studies. Furthermore, amidated, carboxylic, hydroxylated, and methylated metabolites were also identified (Figure 2; Table 4). The in vitro half-life (t1/2) indicated that the metabolism of amygdalin in RGME was superior to that obtained in the RLM preparation (Supplementary Table S8). Similarly, the CLint obtained in vitro indicated that amygdalin achieved a higher clearance rate in RGME than in RLM, which also demonstrated that the gut microbiota may perform the primary functions in amygdalin metabolism.

FIGURE 2
www.frontiersin.org

FIGURE 2. Proposed major metabolic pathways for amygdalin occurring in the incubation system and in the decoction preparation.

As shown in Figure 2, the cyanogenesis and non-cyanogenesis metabolism pathways depended on whether the carbon and nitrogen atoms presented a triple bond. Amination and carboxylation processes were demonstrated to be the detoxification pathways of amygdalin-induced cyanide toxicity. Compounds derived from the traditional decoction and dispensing granules were identified or tentatively characterized according to the 16 compounds emerging from the MS/MS analysis (Supplementary Figure S5). The chemical transformation of amygdalin was compared for the traditional decoction and dispensing granule preparations. As shown in Figure 3, obvious metabolic transformations occurred in the MX decoction preparation compared to the dispensing granule preparation. Overall, these results indicated that the detoxification of Xingren might occur during the boiling of the preparation of the herbal mixture and by metabolism in vivo.

FIGURE 3
www.frontiersin.org

FIGURE 3. The PCA plots (A) and heatmap (B) of 16 compounds in Xingren decoction, MX (Mahuang–Xingren) decoction, Xingren dispensing granules, and MX dispensing granules.

Discussion

Optimization of UPLC-MS/MS Condition

Amygdalin and prunasin, glycosides with high polarity and inefficient protonation, showed weak signal intensity in the MS/MS electrospray ionization mode. Although the bioanalysis of amygdalin has been reported in several pharmacokinetic studies, satisfactory chiral separation has not been achieved of the two epimers to date (Picmanova et al., 2015). Herein, LC-MS/MS optimization was performed by selecting optimal mobile-phase modifiers, to achieve the chiral separation of the epimers, and by adjusting the SPE conditions in the sample preparation.

The addition of a hydrogen cation [M + H]+ or of a sodium ion [M + Na]+ was observed in precursor ion scans, but no significant or stable product ions were detected with increasing collision energy. The adduct ion of ammonium [M + NH4]+ was shown to achieve high abundance and reproducible precursor and product ions. Cyclodextrin-based chiral stationary phases were also obtained successfully for the chromatography separation of chiral molecules with aromatic substituents (Li et al., 2017). The hydroxylpropyl group of the derivatized cyclodextrins extended beyond the mouth of the cyclodextrin cavity. Considering the superiority of liquid chromatographic enantiomer separation, CDShell-RSP (beta-cyclodextrin, hydroxylpropyl), the best column stationary phases available for superficially porous particle silica, was adopted to separate the two epimers of amygdalin and prunasin (Barhate et al., 2017). As shown in Supplementary Figure S3, the retention times of D-amygdalin, neoamygdalin, D-prunasin, and sambunigrin were 3.52, 3.20, 4.22, and 5.17 min, respectively, with similar baseline resolution. Protein precipitation and liquid–liquid extraction were performed according to published data, but the epimers could not be separated efficiently due to a higher proportion of organic solvents. Nitrogen steam evaporation has also been shown to be time-consuming (Zhao et al., 2021b). SPE, or conditioning, loading of the sample, washing, and elution (CLWE), is an efficient protocol for sample preparation in LC-MS/MS quantification (Bury et al., 2018). In this study, the composition of the washing solvent was optimized for reproducible recovery and acceptable matrix interference using gradient elution. The data indicated that satisfactory extraction of amygdalin and prunasin was achieved using a 10% methanol aqueous solution containing 0.1% formic acid.

MX Compatibility Superiority: Decoction Maintain a Balance Between Therapeutic Effectiveness and Safety

In our previous urinary excretion study, the XIC of MRM chromatograms of amygdalin displayed two pairs of peaks indicating different retention times, suggesting that there could still be an unknown metabolite that shares the same chemical structure as amygdalin. Amygdalin may be aminated, methylated, hydroxylated, and so on in vivo (Attwa et al., 2018b). In this study, a powerful LC-HRMS method was employed to screen and identify the potential metabolites of amygdalin. The deamination metabolism pathway was first proposed as the detoxification pathway based on the identification of amygdalin amide and other metabolites. Interestingly, the pathway was also identified during the preparation of the MX decoction. These results suggested that the traditional MX decoction was superior to its single-herb application in maintaining a balance between therapeutic effectiveness and safety. Conversely, MX reduced the relative exposure risk of D-amygdalin and increased the bioavailability of amygdalin. In contrast, amygdalin might be partially detoxified by MX decoction and by in vivo metabolism. Herb pairings embody the balanced wisdom of Chinese herbal compatibility (Wang et al., 2020), which is also confirmed by the MX compatibility of the stereoselectivity and deamination of amygdalin.

Conclusion

Lower exposure of D-amygdalin and a higher abundance of detoxified metabolites were present in the MX decoction. Based on our findings, stereoselective pharmacokinetics and deamination metabolism are the detoxification pathways underlying the compatibility of the MX preparation. Thus, MX decoctions should not be completely replaced by dispensing granules from the perspective of toxicity. However, the detoxification mechanism of amygdalin was not well-presented in this study. The deamination pathway of amygdalin in vivo and its regulation by Mahuang remain unclear and merit further investigation. Therefore, transport mechanisms are being explored using the Caco-2 cell model.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding authors.

Ethics Statement

The animal study was reviewed and approved by the Experimental Animal Ethics Committee of Anhui Medical University.

Author Contributions

SS and CS conceived the idea and designed the study. YQ conducted the experiment, analyzed the data, and wrote the manuscript. SSW provided technical support, performance data analysis, and revised the manuscript. QW, QX, JS, SW, and GC conducted in the experiments. All authors read and approved the final version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (No. 81703800), the First Affiliated Hospital of Anhui Medical University for the Cultivation Program of Youth Science Foundation (No. 2792), Grants for Scientific Research of BSKY (No. XJ201625) from Anhui Medical University, the fund of Anhui medical university doctoral start research (No. 0601069201) and Innovation and Entrepreneurship Training Program for Undergraduates of Anhui Medical University (No.S202010366011).

Conflict of Interest

CS was employed by Hefei Kaifan Analytical Technology Co., Ltd.

The remaining 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.

Acknowledgments

We thank the Center for Scientific Research of Anhui Medical University for its assistance.

Supplementary Material

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

References

Amer, S. M., Kadi, A. A., Darwish, H. W., and Attwa, M. W. (2017). LC-MS/MS Method for the Quantification of Masitinib in RLMs Matrix and Rat Urine: Application to Metabolic Stability and Excretion Rate. Chem. Cent. J. 11 (1), 136. doi:10.1186/s13065-017-0365-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Attwa, M. W., Kadi, A. A., Alrabiah, H., and Darwish, H. W. (2018a). LC-MS/MS Reveals the Formation of Iminium and Quinone Methide Reactive Intermediates in Entrectinib Metabolism: In Vivo and In Vitro Metabolic Investigation. J. Pharm. Biomed. Anal. 160, 19–30. doi:10.1016/j.jpba.2018.07.032

CrossRef Full Text | Google Scholar

Attwa, M. W., Kadi, A. A., Darwish, H. W., Amer, S. M., and Alrabiah, H. (2018b). A Reliable and Stable Method for the Determination of Foretinib in Human Plasma by LC-MS/MS: Application to Metabolic Stability Investigation and Excretion Rate. Eur. J. Mass. Spectrom. (Chichester) 24 (4), 344–351. doi:10.1177/1469066718768327

CrossRef Full Text | Google Scholar

Barhate, C. L., Regalado, E. L., Contrella, N. D., Lee, J., Jo, J., Makarov, A. A., et al. (2017). Ultrafast Chiral Chromatography as the Second Dimension in Two-Dimensional Liquid Chromatography Experiments. Anal. Chem. 89 (6), 3545–3553. doi:10.1021/acs.analchem.6b04834

PubMed Abstract | CrossRef Full Text | Google Scholar

Blaheta, R. A., Nelson, K., Haferkamp, A., and Juengel, E. (2016). Amygdalin, Quackery or Cure? Phytomedicine 23 (4), 367–376. doi:10.1016/j.phymed.2016.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Bury, D., Belov, V. N., Qi, Y., Hayen, H., Volmer, D. A., Brüning, T., et al. (2018). Determination of Urinary Metabolites of the Emerging UV Filter Octocrylene by Online-SPE-LC-MS/MS. Anal. Chem. 90 (1), 944–951. doi:10.1021/acs.analchem.7b03996

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, Y., Tang, S., Chen, A., Zhang, Y., Liu, M., and Wang, X. (2019). Evaluation of the Inhibition Risk of Shikonin on Human and Rat UDP-Glucuronosyltransferases (UGT) through the Cocktail Approach. Toxicol. Lett. 312, 214–221. doi:10.1016/j.toxlet.2019.05.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Commission, N. P. (2015). Pharmacopoeia of the People's Republic of China. Beijing: China Medical Science and Technology Press.

Google Scholar

Cressey, P., and Reeve, J. (2019). Metabolism of Cyanogenic Glycosides: A Review. Food Chem. Toxicol. 125, 225–232. doi:10.1016/j.fct.2019.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, Y., Li, Q., Tong, J., Verpoorte, R., Zhao, S. J., Qin, X. M., et al. (2019). Quality Marker Identification Based on Standard Decoction of Differently Processed Materials of Ephedrae Herba. J. Ethnopharmacol. 237, 47–54. doi:10.1016/j.jep.2019.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Denesh, D., Carbonell, J., Kane, J. S., Gracie, D., and Selinger, C. P. (2021). Patients with Inflammatory Bowel Disease (IBD) Prefer Oral Tablets over Other Modes of Medicine Administration. Expert Rev. Gastroenterol. Hepatol. 15 (9), 1091–1096. doi:10.1080/17474124.2021.1898944

PubMed Abstract | CrossRef Full Text | Google Scholar

Han, H. Y., Huh, J. I., Han, S. R., Kang, M. G., Yoon, S., Han, J. S., et al. (2018). Assessing the Safety of an Ephedrae Herba Aqueous Extract in Rats: A Repeat Dose Toxicity Study. Regul. Toxicol. Pharmacol. 94, 144–151. doi:10.1016/j.yrtph.2018.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Z., Yang, X., Ho, P. C., Chan, S. Y., Heng, P. W., Chan, E., et al. (2005). Herb-drug Interactions: a Literature Review. Drugs 65 (9), 1239–1282. doi:10.2165/00003495-200565090-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

Hulin, A., Stocco, J., and Bouattour, M. (2019). Clinical Pharmacokinetics and Pharmacodynamics of Transarterial Chemoembolization and Targeted Therapies in Hepatocellular Carcinoma. Clin. Pharmacokinet. 58 (8), 983–1014. doi:10.1007/s40262-019-00740-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Karimi, Z., Firouzi, M., Dadmehr, M., Javad-Mousavi, S. A., Bagheriani, N., and Sadeghpour, O. (2021). Almond as a Nutraceutical and Therapeutic Agent in Persian Medicine and Modern Phytotherapy: A Narrative Review. Phytother. Res. 35 (6), 2997–3012. doi:10.1002/ptr.7006

PubMed Abstract | CrossRef Full Text | Google Scholar

Knudsen, C., Neyrinck, A. M., Leyrolle, Q., Baldin, P., Leclercq, S., Rodriguez, J., et al. (2021). Hepatoprotective Effects of Indole, a Gut Microbial Metabolite, in Leptin-Deficient Obese Mice. J. Nutr. 151 (6), 1507–1516. doi:10.1093/jn/nxab032

CrossRef Full Text | Google Scholar

Kovacikova, E., Kovacik, A., Halenar, M., Tokarova, K., Chrastinova, L., Ondruska, L., et al. (2019). Potential Toxicity of Cyanogenic Glycoside Amygdalin and Bitter Apricot Seed in Rabbits-Health Status Evaluation. J. Anim. Physiol. Anim. Nutr. (Berl) 103 (2), 695–703. doi:10.1111/jpn.13055

CrossRef Full Text | Google Scholar

Li, L., Cheng, B., Zhou, R., Cao, Z., Zeng, C., and Li, L. (2017). Preparation and Evaluation of a Novel N-Benzyl-Phenethylamino-β-Cyclodextrin-Bonded Chiral Stationary Phase for HPLC. Talanta 174, 179–191. doi:10.1016/j.talanta.2017.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, J.-x., Xiao, X., Zhou, K.-f., Huang, G., Ao, B., Zhang, Y., et al. (2021). Herb Pair of Ephedrae Herba-Armeniacae Semen Amarum Alleviates Airway Injury in Asthmatic Rats. J. Ethnopharmacology 269, 113745. doi:10.1016/j.jep.2020.113745

CrossRef Full Text | Google Scholar

Peterson, G., Nilsson, D., Peterson, S., Dedering, A., Trygg, J., Wallman, T., et al. (2016). Changes in Dorsal Neck Muscle Function in Individuals with Chronic Whiplash-Associated Disorders: A Real-Time Ultrasound Case-Control Study. Ultrasound Med. Biol. 42 (5), 1090–1102. doi:10.1016/j.ultrasmedbio.2015.12.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Pičmanová, M., Neilson, E. H., Motawia, M. S., Olsen, C. E., Agerbirk, N., Gray, C. J., et al. (2015). A Recycling Pathway for Cyanogenic Glycosides Evidenced by the Comparative Metabolic Profiling in Three Cyanogenic Plant Species. Biochem. J. 469 (3), 375–389. doi:10.1042/BJ20150390

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, J. L., Dong, H., Han, Y., Yang, L., Zhang, A. H., Sun, H., et al. (2020). Network Pharmacology Combined with Metabolomics Approach to Investigate the Protective Role and Detoxification Mechanism of Yunnan Baiyao Formulation. Phytomedicine 77, 153266. doi:10.1016/j.phymed.2020.153266

PubMed Abstract | CrossRef Full Text | Google Scholar

Sodhi, J. K., and Benet, L. Z. (2021). Successful and Unsuccessful Prediction of Human Hepatic Clearance for Lead Optimization. J. Med. Chem. 64 (7), 3546–3559. doi:10.1021/acs.jmedchem.0c01930

CrossRef Full Text | Google Scholar

Song, S., Ma, Q., Tang, Q., Chen, F., Xing, X., Guo, Y., et al. (2016). Stereoselective Metabolism of Amygdalin-Based Study of Detoxification of Semen Armeniacae Amarum in the Herba Ephedrae-Semen Armeniacae Amarum Herb Pair. J. Ethnopharmacol. 179, 356–366. doi:10.1016/j.jep.2015.12.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, Y., Zhang, J., Zhao, A., Li, W., Feng, Q., and Wang, R. (2020). Effects of Intestinal flora on the Pharmacokinetics and Pharmacodynamics of Aspirin in High-Altitude Hypoxia. PLoS One 15 (3), e0230197. doi:10.1371/journal.pone.0230197

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Wong, Y. K., and Liao, F. (2018). What Has Traditional Chinese Medicine Delivered for Modern Medicine? Expert Rev. Mol. Med. 20, e4. doi:10.1017/erm.2018.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, N., Du, N., Peng, Y., Yang, K., Shu, Z., Chang, K., et al. (2020). Network Patterns of Herbal Combinations in Traditional Chinese Clinical Prescriptions. Front. Pharmacol. 11, 590824. doi:10.3389/fphar.2020.590824

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Zhang, D., Tang, D., Sun, K., Peng, J., Zhu, W., et al. (2021). Amygdalin Inhibits TGFβ1-Induced Activation of Hepatic Stellate Cells (HSCs) In Vitro and CCl4-Induced Hepatic Fibrosis in Rats In Vivo. Int. Immunopharmacol. 90, 107151. doi:10.1016/j.intimp.2020.107151

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, S., Hu, Y., Tan, W., Wu, X., Chen, R., Cao, J., et al. (2012). Compatibility Art of Traditional Chinese Medicine: from the Perspective of Herb Pairs. J. Ethnopharmacol. 143 (2), 412–423. doi:10.1016/j.jep.2012.07.033

CrossRef Full Text | Google Scholar

Wang, Z., Fang, K., Wang, G., Guan, X., Pang, Z., Guo, Y., et al. (2019). Protective Effect of Amygdalin on Epithelial-Mesenchymal Transformation in Experimental Chronic Obstructive Pulmonary Disease Mice. Phytother. Res. 33 (3), 808–817. doi:10.1002/ptr.6274

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, F., Huang, J., Liu, W., Wang, C. R., Liu, Y. F., Tu, D. Z., et al. (2021). Inhibition of Drug-Metabolizing Enzymes by Qingfei Paidu Decoction: Implication of Herb-Drug Interactions in COVID-19 Pharmacotherapy. Food Chem. Toxicol. 149, 111998. doi:10.1016/j.fct.2021.111998

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Guo, C., Lv, J., Hou, S., Zhang, Y., Gao, J., et al. (2020). Aqueous Chlorination of Ephedrine: Kinetic, Reaction Mechanism and Toxicity Assessment. Sci. Total Environ. 740, 140146. doi:10.1016/j.scitotenv.2020.140146

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, C., Wei, M., Zheng, Y., Tao, W., Lv, Q., Wang, Q., et al. (2021a). The Analyses of Chemical Components from Oldenlandia Hedyotidea (DC.) Hand.-Mazz and Anticancer Effects In Vitro. Front. Pharmacol. 12, 624296. doi:10.3389/fphar.2021.624296

CrossRef Full Text | Google Scholar

Zhao, J., Tian, S., Lu, D., Yang, J., Zeng, H., Zhang, F., et al. (2021b). Systems Pharmacological Study Illustrates the Immune Regulation, Anti-infection, Anti-inflammation, and Multi-Organ protection Mechanism of Qing-Fei-Pai-Du Decoction in the Treatment of COVID-19. Phytomedicine 85, 153315. doi:10.1016/j.phymed.2020.153315

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: amygdalin, LC-MS/MS, pharmacokinetic, stereoselectivity, compatibility, metabolism, detoxification

Citation: Qin Y, Wang S, Wen Q, Xia Q, Wang S, Chen G, Sun J, Shen C and Song S (2021) Interactions Between Ephedra sinica and Prunus armeniaca: From Stereoselectivity to Deamination as a Metabolic Detoxification Mechanism of Amygdalin. Front. Pharmacol. 12:744624. doi: 10.3389/fphar.2021.744624

Received: 20 July 2021; Accepted: 18 October 2021;
Published: 26 November 2021.

Edited by:

Ruiwen Zhang, University of Houston, United States

Reviewed by:

Guang-Bo Ge, Shanghai University of Traditional Chinese Medicine, China
Huaiyou Wang, Henan University, China
Mohamed Attwa, King Saud University, Saudi Arabia

Copyright © 2021 Qin, Wang, Wen, Xia, Wang, Chen, Sun, Shen and Song. 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: Chenlin Shen, shchlin520@163.com; Shuai Song, vae0558@163.com

These authors have contributed equally to this work and share first authorship

Disclaimer: 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.