Next Article in Journal
Wet-Chemical Synthesis of TiO2/PVDF Membrane for Energy Applications
Next Article in Special Issue
Comprehensive Evaluation of the Bioactive Composition and Neuroprotective and Antimicrobial Properties of Vacuum-Dried Broccoli (Brassica oleracea var. italica) Powder and Its Antioxidants
Previous Article in Journal
Surface Coordination of Pd/ZnIn2S4 toward Enhanced Photocatalytic Activity for Pyridine Denitrification
Previous Article in Special Issue
Natural Polyphenols for Treatment of Colorectal Cancer
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Fate and Intermediary Metabolism of Soyasapogenol in the Rat

1
Department of Food Science and Technology, School of Agriculture and Biology, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
2
Bor S. Luh Food Safety Research Center, Shanghai Jiao Tong University, 800 Dongchuan Road, Shanghai 200240, China
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(1), 284; https://doi.org/10.3390/molecules28010284
Submission received: 11 November 2022 / Revised: 20 December 2022 / Accepted: 27 December 2022 / Published: 29 December 2022
(This article belongs to the Special Issue Bioactive Compounds in Food Bioscience and Pharmacology)

Abstract

:
Research suggests that soyasaponins are poorly absorbed in the GI tract and that soyasaponin aglycones or soyasapogenols are absorbed faster and in greater amounts than the corresponding soyasaponins. Therefore, it is important to understand the bioavailability of these compounds for the potential development of functional foods containing their components. In this paper, to investigate the metabolic characteristics of soyasapogenols A and B, the pharmacokinetic parameters in rats were determined via oral and intravenous administration. The liver metabolites of soyasapogenols were identified using UPLC–/Q-TOF–MS/MS, and their metabolic pathways were also speculated. The results show that, after oral administration, there was a bimodal phenomenon in the absorption process. Tmax was about 2 h, and soyasapogenol was completely metabolized 24 h later. The bioavailability of soyasapogenol was superior, reaching more than 60%. There were sixteen metabolites of soyasapogenol A and fifteen metabolites of soyasapogenol B detected in rat bile. Both phase I and II metabolic transformation types of soyasapogenols, including oxidation, dehydrogenation, hydrolysis, dehydration, deoxidization, phosphorylation, sulfation, glucoaldehyde acidification, and conjugation with cysteine, were identified. In addition, soyasapogenol A could be converted into soyasapogenols B and E in the metabolic process. These results suggest that it is feasible to use soyasapogenols as functional ingredients in nutraceuticals or food formulations.

Graphical Abstract

1. Introduction

Soyasaponins are oleanene-type triterpenoid saponins, and they can be mainly categorized into two groups, A and B, according to their respective aglycones, soyasapogenol A and soyasapogenol B [1]. Soybeans and soy foods are the major sources of dietary soyasaponins consumed by humans. An increasing number of investigations have demonstrated the health-promoting activities of dietary soyasaponins, such as lowering serum lipids [2,3] and inhibiting the growth of cancer cells [4]. However, to date, detailed information on the fate or bioavailability of soyasaponins in the guts of animals and humans is generally lacking. Previous studies have suggested that soyasaponins are poorly absorbed early in the GI tract and that they are mainly metabolized in the colon, which means that the microbiota is the main agent responsible for their hydrolysis and the production of sapogenols [5,6]. In this respect, it is important to note that diverse sapogenols have demonstrated superior bioavailability and bioactivity compared to their precursor saponin due to their more favorable chemical properties caused by a lack of sugar chains [7,8]. Soyasapogenols have been reported to possess a wide range of health benefits, including anticancer, antivirus, and antioxidant activities, as well as hepatoprotective and cardiovascular protective effects [3,9,10]. However, evidence demonstrating that soyasapogenols are systemically bioavailable in the conjugated form following microbial metabolism in the colon is still limited [8]. Therefore, it is essential to elucidate the metabolism and bioavailability of soyasapogenols in order to understand their bioactivity and to properly design functional food ingredients.
Examining the kinetic dynamics of metabolism is an important method to reveal the physiological benefits and bioavailability of biologically active ingredients, such as saponins, after consumption [11,12]. Biologically active substances in the body are absorbed into the blood and need to reach a specific plasma concentration to be effective [8,12]. Therefore, plasma concentration data are essential for determining the pharmacokinetic (distribution and elimination) and pharmacodynamic (equilibration and sensitivity) characteristics of bioactive components.
The LC–MS technique combines high separation ability with the high selectivity of high-resolution mass spectrometry and the ability to provide information (molecular weight and structural information) on complex samples, and it is widely used to control quality standards for natural phytochemicals and/or their intestinal metabolites [13,14]. The composition of natural herbal medicines as well as the changes in their chemical composition before and after processing have been rapidly analyzed using UPLC–Q-TOF–MS/MS, making it an effective tool for identifying active ingredients, especially those of saponins and their derivatives, with improved sensitivity and accuracy [15,16,17].
The aims of the present work are to assess the bioavailability of dietary soyasaponins in rats by directly using soyasapogenols as a dietary component and to establish the concentrations and forms of the main metabolites following the consumption of soyasapogenol-rich diets. The findings of this study may provide further insights into the direct use of soyasapogenols as functional ingredients in food formulations.

2. Results and Discussion

2.1. Optimization of the UPLC–Q-TOF–MS/MS Method

The choices of mobile-phase pairs for gradient elution and column separation were investigated. As a result, the acetonitrile–water pair was selected as the mobile-phase pair for the method. In addition, 0.1% formic acid was added to the mobile-phase pair to facilitate the protonation or deprotonation of the analytes for the mass spectrometric detection of the analytes in the positive ionization mode. It was also found that the Waters ACQUITY UPLC® BEH C18 column (2.1 mm i.d. × 100 mm, 1.7 μM) had greater separation efficiency and produced a larger signal-to-noise ratio and a better peak shape; therefore, it was adopted for the method. To optimize the MS/MS parameters for the compounds, both positive and negative ionization modes were applied to the analyses of soyasapogenol and its derived metabolite samples using QTOF–MS/MS, and comprehensive information about these components was obtained. The positive ion mode gave good ionization, and collision-induced dissociation energy was optimized in the range of 20–30 ev. For all the compounds, the protonated molecules [M + H]+ and their sodium adduct [M + Na]+ ions were selected based on their high abundance and stronger fragmentation regularity. To confirm the results and the structures of the metabolites, an ESI–QTOF–MS/MS analysis was performed. The peak results were identified based on their MS spectral data (accurate mass and fragmentation patterns), a comparison to standard compounds, and a search of public online databases (DNP, Reaxys, and SciFinder). The presence of a particular type of metabolite was proposed based on the occurrence of a characteristic molecular ion peak, and a further analysis of its fragmentation ions implied the structure of the metabolite.
At first, individual standard solutions of soyasapogenols A and B were directly injected. The recorded first-order mass spectra of the positive ions are displayed in Figure 1. The retention times of soyasapogenols A and B were 10.23 min and 10.72 min, respectively. The excimer ion peak [M + H]+ m/z of soyasapogenol A was 475.2774. The main fragment ions under the high-energy conditions were observed at m/z 185.1321, 283.2627, 333.2782, and 337.2521. The excimer ion peak [M + H]+ m/z of soyasapogenol B was 459.2299, and the main fragment ions were observed at m/z 163.1106, 283.2627, 401.2284, and 429.3716.

2.2. Comparison of Plasma Soyasapogenol Concentrations in Rats after Soyasapogenol Intake

2.2.1. Pharmacokinetics of Oral Administration

The pharmacokinetic parameters of soyasapogenols A and B in the plasma of the rats 0–24 h after oral administration are summarized in Table 1. As shown in Figure 2, following the oral administration of soyasapogenol A, there were rapid increases in the plasma concentrations to 29 ± 10 and 7.0 ± 6 μg/mL at the dosage levels of 10 mg/kg and 20 mg/kg by gavage at 1.0 h postadministration, respectively. The plasma concentrations reached the maxima (Cmax) of 8.59 mg /L at 2.0 h postadministration in the case of the 10 mg/kg dose and 28.99 mg /L at 1.0 h postadministration in the case of the 10 mg/kg dose, and then they decreased. The results also demonstrate that the area under the plasma or serum concentration–time curve from 0 to 24 h, AUC(0-∞), and Cmax were both linear and proportional to the oral dose. Likewise, after the oral administration of soyasapogenol B, it was also rapidly absorbed, with a Tmax of 2 h. At different concentrations, the Cmax values were 26.37 mg/L and 40.29 mg/L, and the half-life T1/2z was very similar, which demonstrates the linear kinetics absorption process of soyasapogenol B in the rats following oral administration. In addition, the plasma elimination half-life appeared to be independent of the dose. As shown, the area under the curve AUC(0-∞) gives a profile identical to that of soyasapogenol A in Figure 2 (left panel), and the results revealed that they were positively correlated (R = 0.63, p < 0.001).
Regarding the absorption and elimination of the orally administrated soyasapogenols, as indicated in Figure 2, in each dosage group, the kinetic profile revealed rapid absorption with a one-peak plasma concentration versus time course, followed by a different decreasing pattern consisting of a distribution phase and an elimination phase between the different dosage groups. This result is consistent with previous investigations [5], which demonstrated that the bioavailability of soyasapogenols is better than that of the corresponding soyasaponins and that the absorption of soyasapogenol B is better than that of soyasapogenol A. In this work, the results indicate that the time needed for the soyasapogenol A administered to the rats to reach peak concentration was 1 to 3 h, whereas the time needed for the soyasapogenol B administered to the rats to reach peak concentration was 8 h.
We previously reported the in vivo absorption characteristic properties of soyasapogenols A and B at different locations in four intestinal segments of mice—the duodenum, jejunum, ileum, and colon. The results demonstrated that, in terms of the apparent absorption coefficient (Papp) and the absorption rate constant (Ka), the rank order of the absorption of soyasapogenol from various administration sites was as follows: jejunum > duodenum > ileum > colon, and the absorption capacity of soyasapogenol A exposure significantly and concentration-dependently accumulated over time compared to soyasapogenol B [18]. The results suggest that some transporters may be involved in the regulation of soyasapogenol B absorption in addition to concentration gradient-dependent diffusion.

2.2.2. Pharmacokinetics of Intravenous Injection

Following the intravenous administration of a bolus dose of 3 mg/kg soyasapogenols to the rats, the elimination of soyasapogenol also followed a two-compartment or three-compartment model of elimination (Figure 3).
The plasma concentration–time curve of the soyasapogenol intravenous injection is shown in Figure 3. The results indicate that the absorption of aglycone through the dorsal foot vein injection in the rats was very fast, within 0–2 h; the blood concentration decreased rapidly, and then the speed of its clearance slowed down until most of it had disappeared in vivo after 24 h.
The main metabolic parameters of the soyasapogenols obtained after the intravenous administration are shown in Table 2. The soyasapogenols disappeared rapidly from the blood after the intravenous injections, reaching their peak at 0.03 h. The half-lives T1/2z of soyasapogenol A were 3.54 h and 3.48 h, and the half-lives T1/2z of soyasapogenol B were 5.14 h and 5.53 h, which indicates that the half-life of the soyasapogenols after the intravenous injections was fixed and independent of the concentration. The Cmax and AUC(0-∞) of the soyasapogenols were positively correlated with their concentrations. The metabolic regularity was similar to that of mangiferin aglycone [19], which undergoes extensive phase II metabolism and predominantly exists as glucuronidation and sulfation metabolites in vivo, and aglycone exhibits favorable pharmacokinetics properties, as partly indicated by its intermediate bioavailability.

2.3. Bioavailability of Soyasapogenols

The bioavailability (F) of bioactive components in vivo has been frequently used to investigate their absorption, metabolism, tissue distribution, and bioactivity [2,8]. The absolute oral bioavailability was calculated from the percentage ratio of the AUCs derived from the plasma soyasapogenol concentrations after oral and I.V. administration of the different sources of soyasapogenols to the rats. As a result, for soyasapogenol A, the F values were 73.10% and 67.34% at the dosages of 10 mg/kg and 20 mg/kg, respectively, while for soyasapogenol B, the F values were 60.94% and 69.01% at the dosages of 25 mg/kg and 50 mg/kg, respectively (Table 1).
The high molecular weight and poor membrane permeability of saponins cause their low bioavailability, hence restricting their application as a drug candidate [20]. For example, the oral bioavailability of ginsenosides Ra3, Rb1, and Rd has been found to be 0.1–0.2%, whereas that of ginsenosides Re, Rg1, and notoginsenoside R1 has been found to be 0.2–0.6% [21]. An investigation has demonstrated that the bioavailability of soyasapogenols is better than that of the corresponding soyasaponins and that the bioavailability of group B soyasaponins is better than that of group A soyasaponins [5]. Previous research has also proved that sapogenols are present in human plasma following microbial metabolism in the colon; furthermore, soyasapogenol B has been found to be a major metabolite, and it has been found in the conjugated form [7]. It should be pointed out that a previous investigation has demonstrated that the amount (%) of soyasapogenol B that permeates into the basolateral side from the Caco-2 cell monolayer apical side is much higher than that of soyasapogenol A [5]. However, in the present study, the result revealed that the absorption of soyasapogenol B was not much higher than that of soyasapogenol A, and this is not in agreement with previous studies. This may be attributable to the different dosages used, as well as to the different experimental models; our data were obtained using in vivo experimental approaches, while the previously reported data were obtained using an in vitro Caco-2 cell line as a model system.
It was demonstrated that methylation, sulfate conjugation, S-cysteine conjugation, sulfate conjugation, and glycine conjugation were the main metabolic pathways in vivo for sulfate-related metabolites, indicating that the phase II reaction plays an important role in transforming soyasapogenols to less toxic metabolites, which has also been observed in other types of saponins [22,23,24].
As shown in Figure 1, the only structural difference between soyasapogenol A and soyasapogenol B is the presence of a hydroxyl group at position 21. Further research is needed to elucidate the mechanism by which this single hydroxyl group affects absorption.
It can be seen that the absorption of the soyasapogenols in the rats after intragastric administration follows a bimodal pattern that is similar to that of other sapogenols of ginsenosides. It is well-known that various endogenous and exogenous compounds undergo enterohepatic circulation and intestinal metabolism, which determine their fate in an organism. The absorption of saponin components across the gastrointestinal wall membrane may be affected by multiple factors, such as solubility, interactions with other dietary ingredients, molecular transformations, different cellular transporters, metabolism, and interactions with the gut microbiota, resulting in changes in their bioavailability.
For dietary soyasaponins and soyasapogenols that are excreted extensively into bile, an insight into the magnitude of enterohepatic circulation is of crucial importance, as it significantly affects pharmacokinetic parameters, such as plasma half-life and the area under the plasma concentration–time curve (AUC), as well as the estimates of bioavailability. Therefore, the findings of this study provide further insights into the direct use of soyasapogenols as functional ingredients in food formulations. However, their poor solubility and undesirable taste merit more studies in order to resolve these problems.

2.4. Identification of Liver Metabolites of Soyasapogenols

2.4.1. Metabolite Identification and Metabolic Characteristic Elucidation of Soyasapogenol A

UPLC–QTOF–MS was used to identify the liver metabolites of soyasapogenol A in the rats, and the extraction ion chromatograms are shown in Figure 4. According to the relative molecular weight value of each metabolite and the fragment ion information of the secondary mass spectrum, the possible structure of each metabolite was speculated according to the normal metabolic reaction in vivo. A total of 16 metabolites named M1–M16, principally glucuronides, were identified in the plasma of the rats after oral and intravenous administration. Their identities, summarized in Table 3, were determined using UPLC–QTOF–MS in the consecutive reaction monitoring mode.

2.4.2. Metabolite Identification and Metabolic Characteristic Elucidation of Soyasapogenol B

Moreover, using the same analysis conditions as those used for the analysis of the soyasapogenol A metabolites, as shown in Table 4 and Figure 4, a total of 15 of soyasapogenol B metabolites, named M1–M15, were identified unambiguously by comparing their retention times, mass accuracies, and fragmentation behaviors with data from the corresponding reference standards, or they were tentatively characterized based on MS data and reference to the available data in the literature [5,7].

2.5. Metabolic Pathway of Soyasapogenols In Vivo

In this work, 16 compounds were detected for soyasapogenol A metabolites, and 15 compounds were observed in soyasapogenol B metabolite samples. Summarizing all of the results presented here, the proposed metabolic pathways of the soyasapogenols in vivo are illustrated in Figure 5. The major biotransformation of the soyasapogenols involved oxidation, dehydrogenation, hydrolysis, dehydration, and deoxygenation reactions, as well as other phase I reactions, in the metabolic process; furthermore, phosphorylation, sulfation, gluconaldehyde acidification, and conjugation with cysteine, as well as other phase II reactions, also occurred in the metabolic process. According to the reactions, M1–M6 are one-step metabolites; M7–M9 are two-step metabolites of M1 conjugated with cysteine, cysteine sulfhydryl, and acetylcysteine; M10–M11 are two-step metabolites of M2; M12–M14 are two-step metabolites of M4; and M15–M16 are two-step metabolites of M6. In addition, soyasapogenol A can be metabolized into soyasapogenol B and soyasapogenol E under certain conditions.

3. Materials and Methods

3.1. Chemicals and Reagents

HPLC-grade methanol and analytical-grade ethanol were purchased from Shanghai Lingfeng Chemical Reagents Co., Ltd. (Shanghai, China).
Soyasapogenol A and B standards (98 %+), β-glucuronidase, and sulfatase were purchased from Sigma-Aldrich (Shanghai) Trading Co., Ltd.
Pentobarbital sodium (AR, 98 %) was bought from Nanjing Jiancheng Biological Engineering Research Institute (Nanjing, China). Sodium citrate as an anticoagulant (4%) and Krebs–Ringer buffer were bought from Beijing Leagene Biotech Co., Ltd. (Beijing, China). Saline (0.9%) was obtained from Shandong Kangning Pharmaceutical Co., Ltd. (Liaocheng, China). Formononetin (FMN) was purchased from Yuanye (Shanghai, China). Food-grade sodium carboxymethyl cellulose (CMC) was purchased from Shanghai Aladdin Biochemical Technology Co., Ltd. (Shanghai, China).
All other reagents were of analytical grade, and ultra-pure water was used throughout the experimentation.

3.2. Preparation of Soyasapogenol Standard Solution

The stock solution of soyasapogenols A and B (IS) was prepared by dissolving 1.00 mg powder in 1.00 mL of methanol to a 1.00 mg/mL concentration. The working solution of IS was prepared via a 1/10 dilution of the stock solution in methanol to a concentration of 100 μg/mL, and it was stored in a refrigerator at 4 °C for later use.

3.3. Preparation and Application of Dosing Solutions

Soyasapogenol was prepared via the acid hydrolysis of the crude extract according to our previous methodology [25]. For the preparation of the soyasapogenol samples for oral and intravenous administration, 0.5% (w/v) sodium carboxymethyl cellulose (NaCMC) was dissolved in pure water as a vehicle; then, the soyasapogenols were added to the solution under continuous magnetic stirring to the desired concentration. For the intestine experiment, the in vivo intestinal perfusion consisted of soyasapogenol in Krebs–Ringer bicarbonate buffer (Beijing Leagene Biotech Co., Ltd., Beijing, China). Soyasapogenol was dissolved in 500 ml Krebs–Ringer bicarbonate buffer at two concentrations of 10 mg/L and 50 mg/L.

3.4. Animals

Twenty-four Sprague Dawley adult male rats (250 ± 20 g) were purchased from SPF Beijing Vital River Laboratory Animal Technology Co., Ltd. (Beijing, China). The eleven-week-old male Sprague Dawley rats were maintained in a ventilated room at 22 ± 3°C and 40–60% humidity with a 12 h light–dark cycle. The rats were fasted for 12 h and allowed to drink freely before the experiment.
This study was approved by the Experimental Animal Ethics Committee of the Shanghai Jiao Tong University, and it strictly followed the instructions of the Organizational Guidelines and Ethics Guidelines of the Experimental Animal Ethics Committee (No. IACUC, approval no. A2019086).

3.5. Metabolic Kinetics and Bile Metabolism Experiment of Soyasapogenols

3.5.1. Oral Administration

The experiments were performed according to Kamo et al [5]. The rats were divided into three groups (n = 8). All rats were fasted for 14–15 h prior to soyasapogenol A or B administration. All rats were orally administered soyasapogenol A or B via direct stomach intubation. The dosages of the test samples are summarized in Table 5.

3.5.2. Intravenous Injection Experiment

The rats’ dorsal foot veins were cannulated with a standard syringe needle to administer the soyasapogenol. As presented in Table 5, the dosages of soyasapogenol A were 10 mg/kg and 20 mg/kg, and the dosages of soyasapogenol B were 25 mg/kg and 50 mg/kg. Blood samples were collected from the tail vein just prior to soyasapogenol administration (0 h) and 0.03, 0.08, 0.25, 0.5, 1, 2, 4, 6, 8, and 24 h after soyasapogenol loading; blood soyasapogenol levels were measured immediately. The collected blood sample volume was approximately 200 μL, and the samples were subjected to the same sample pretreatment procedure as described below.

3.5.3. Perfusion of the Small Intestine

The preparation of the isolated, perfused small bowel was performed according to Guo et al. [19], with minor modifications, and the dosages of the test samples are presented in Table 5. Briefly, before the experiment, the rats were fasted for 12 h and supplied with distilled water. The rats were subjected to bile metabolism experiments, intraperitoneally injected with 1% sodium pentobarbital (40 mg/kg) as anesthesia, and placed on an electric heating pad after anesthesia administration to ensure a constant temperature of 37 °C; their backs were fixed to the operating table. The abdominal cavity was cut along the midline of the abdomen, with an incision of 3~4 cm, and the position of each intestinal segment was determined; the location of the bile was found, and a thin catheter was carefully inserted to collect the bile. Afterwards, the duodenum was incised at both ends, and the tubes were intubated and ligated with surgical sutures. The intestinal materials were slowly washed out with preheated normal saline at 37 °C. After the normal saline was drained, a constant flow pump was connected, and Krebs buffer was used first. The Ringer buffer (preheated at 37 °C in advance) was used to equilibrate each intestinal segment at a flow rate of 0.2 mL/min for 30 minutes, and then soyasapogenol perfusion fluid was used to start the experiment. Bile samples were collected in centrifuge tubes, and they were moistened during perfusion. A cotton cloth with physiological saline covered the abdominal cavity to keep it moist, and 1% sodium pentobarbital was given intermittently to maintain the anesthesia.
Blood samples were collected from the tail vein at 0 (pre-dose), 0.25, 0.5, 1, 2, 4, 6, 8, and 24 h after administration. The sample volume was 200 μL. The samples were immediately transferred to a centrifuge tube and treated with a sodium citrate anticoagulant. This fraction was then centrifuged at 10,000 rpm for 10 min, and the upper plasma was collected and placed in a refrigerator at −80 °C for later use.

3.6. Processing and Analysis of Plasma Samples

The soyasapogenol concentrations in the rat plasma solutions were determined according to the method published in [5]. At first, the plasma samples were thawed to 4 °C and transferred (20 to 100 μL) to an Eppendorf tube. After adding the same volume of 0.2 mol/L sodium acetate buffer to each sample, 3 µL aliquots of glucuronidase (Sigma; extracted from Helix pomatia, containing 95,600 units/mL glucuronidase and 1100 units/mL sulfatase) were added to analyze the glucuronyl- and sulfatide-conjugated forms of the soyasapogenols, and the mixture was incubated at 37 °C for 12 h. The samples were mixed thoroughly, and then they were centrifuged at 10,000 rpm for 15 min. The supernatant was blown with nitrogen gas until it was dry, and residual supernatant was dissolved with 1 mL methanol. Subsequently, after filtration, the samples were analyzed according to the method published in [5].

3.7. Processing and Analysis of Bile Samples

The bile samples were taken from the perfusion syringe before and after the perfusion experiment. Effluent perfusate samples were collected every 15 min for 120 min in tared 3 mL centrifuge tubes. Subsequently, 3 times the amount of acetonitrile was added to precipitate proteins. The supernatant was dried under nitrogen flow and re-dissolved in 400 μL acetonitrile with an ultrasonic bath for 5 min. Next, the samples were centrifuged for 10 min at 12,000× g to obtain the supernatant. The final step was to allow the samples to pass through a 0.22 μm millipore filter for an UPLC–QTOF–MS analysis. The processing of the blank group was the same as that of the experimental group.

3.8. Instrumentation and LC–MS/MS Analytical Conditions

The UPLC–QTOF–MS/MS analysis was performed on an ACQUITY UPLC I-Class (Waters, Milford, MA, USA) system, which was coupled with tandem quadrupole (Xevo TQ-S) mass spectrometry (UPLC–Xevo TQ-S–MS/MS).
Gradient UPLC elution was carried out on the column BECH C18 (100 mm × 2.1 mm, 1.7 μm, Waters, Milford, MA, USA) with deionized water containing 0.1% formic acid as mobile-phase A and acetonitrile (LC grade) containing 0.1% formic acid as mobile-phase B. The column temperature was kept at 45 °C, and the flow rate was set to 0.4 mL/min. The gradient elution program was finally optimized as follows: 5–20% B from 0 to 3 min, 20–100% B from 3 to 10 min, 100% B from 12 to 15 min, and 100–95% B from 15 to 19 min. The samples were kept at 20 °C in the auto-sampler manager, and the injection volume was 1.0 μL.

3.9. Mass Spectrometry Analysis Using Xevo G2-XS QTOF Instrument

A mass spectrometry analysis of the soyasapogenol metabolites was completed using a Xevo G2-XS QToF (Waters Corporation, USA) instrument. All mass spectrometry analyses were performed in the positive ionization mode, with mass spectra acquired over the 50 to 1000 m/z range. The following parameter settings were used: a capillary voltage of 2 kV; a cone voltage of 40 V; a turbo spray temperature of 450 °C; and a collision energy (CE) of 6 eV. Nitrogen was used as the nebulizer and the auxiliary gas, the atomizing gas flow rate was set to 900 L/h, and the taper hole blowback gas was set to 50 L/h.

3.10. Data Processing

Curve-fitting simulations were performed using MicroCal Origin software, version 6.1 (MicroCal Software Inc., Northampton, MA, USA).
The metabolites of the soyasapogenols in the rat blood and bile samples were identified using the MetaboLynx™ (v4.1) program (Waters Corp., Milford, MA, USA) based on the UPLC–QTOF–MS chromatograms. In post-acquisition analyses, the MetaboLynx™ (v4.1) program could employ an extensive list of potential biotransformation reactions (e.g., oxidation, phosphorylation, sulfation, and conjugation) of the prototype molecule, in combination with the elemental compositions, to generate a series of extracted ion chromatograms (XICs). These mass spectral peaks, using a mass accuracy of two decimal places, were then compared between the control and experimental samples to eliminate the chromatographic peaks in the samples that also appeared in the blank matrix. Eventually, the metabolites that were derived from the major compounds in vitro and the absorbed prototype constituents were predicted using MetaboLynx software (Waters). The molecular formulas and structures of the metabolites investigated were elucidated based on the information of their precursor ions, fragmentation ions, and chromatographic retention times.
The bioavailability parameters were calculated by comparing AUCoral postadministration with AUCiv post-soyasapogenol I.V. injection:
%Bioavailability = [(AUCoral) × Doseiv/(AUCiv) × Doseoral] × 100
The calculation formula is as follows: F = oral administration AUC(0-∞)/intravenous injection AUC(0-∞).
Data are presented as the percentage change in LV thrombus minus cavity gray scale (%). Data were analyzed by using SigmaStat version 2.03 statistical software (SPSS Inc). The results are expressed as mean ± SEM. Differences between groups were analyzed using the one-way ANOVA Student–Newman–Keuls method; a value of p < 0.05 was considered statistically significant.

4. Conclusions

In conclusion, the present research provides useful information on the in vivo pharmacokinetics and bioavailability of soyasapogenols in rats, helps in the design of new nutraceuticals by providing kinetics-based data, and offers a better understanding of the fate and intermediary metabolism of dietary soyasapogenols. Further investigations are needed to identify the conditions that can enhance their solubility and to clarify their physiologic activities in the digestive organs.

Author Contributions

C.P., investigation, data curation and formal analysis; Y.Y., writing—original draft preparation; D.Z., funding acquisition, project administration, supervision and writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number: 31972069.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Shanghai, and approved by the Ethics Committee of SHANGHAI JIAO TONG UNIVERSITY (protocol code: A2019086, approval date: 25 June 2020).

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The present study was financially supported by a grant from the National Natural Science Foundation of China (31972069). The authors are grateful to Xi Zhang and Dongyue Qi for their assistance during the preparation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available.

References

  1. Zhang, W.; Popovich, D.G. Chemical and Biological Characterization of Oleanane Triterpenoids from Soy. Molecules 2009, 14, 2959–2975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhao, D.Y. Challenges associated with elucidating the mechanisms of the hypocholesterolaemic activity of saponins. J. Funct. Foods 2016, 23, 52–65. [Google Scholar] [CrossRef]
  3. Guang, C.; Chen, J.; Sang, S.; Cheng, S. Biological Functionality of Soyasaponins and Soyasapogenols. J. Agric. Food Chem. 2014, 62, 8247–8255. [Google Scholar] [CrossRef] [PubMed]
  4. Zhou, Y.F.; Farooqi, A.A.; Xu, B.J. Comprehensive review on signaling pathways of dietary saponins in cancer cells suppression. Crit. Rev. Food Sci. Nutr. 2021, 1–26. [Google Scholar] [CrossRef] [PubMed]
  5. Kamo, S.; Suzuki, S.; Sato, T. Comparison of bioavailability (I) between soyasaponins and soyasapogenols, and (II) between group A and B soyasaponins. Nutrition 2014, 30, 596–601. [Google Scholar] [CrossRef] [PubMed]
  6. Navarro del Hierro, J.; Herrera, T.; Fornari, T.; Reglero, G.; Martin, D. The gastrointestinal behavior of saponins and its significance for their bioavailability and bioactivities. J. Funct. Foods 2018, 40, 484–497. [Google Scholar] [CrossRef]
  7. Neacsu, M.; Raikos, V.; Benavides-Paz, Y.; Duncan, S.H.; Duncan, G.J.; Christie, J.S.; Johnstone, A.M.; Russell, W.R. Sapogenol is a Major Microbial Metabolite in Human Plasma Associated with High Protein Soy-Based Diets: The Relevance for Functional Food Formulations. Foods 2020, 9, 422. [Google Scholar] [CrossRef] [Green Version]
  8. Gao, S.; Basu, S.; Yang, Z.; Deb, A.; Hu, M. Bioavailability Challenges Associated with Development of Saponins As Therapeutic and Chemopreventive Agents. Curr. Drug Targets 2012, 13, 1885–1899. [Google Scholar] [CrossRef]
  9. Kang, J.; Badger, T.M.; Ronis, M.J.J.; Wu, X.L. Non-isoflavone Phytochemicals in Soy and Their Health Effects. J. Agric. Food Chem. 2010, 58, 8119–8133. [Google Scholar] [CrossRef]
  10. Zhang, W.; Popovich, D.G. Effect of Soyasapogenol A and Soyasapogenol B Concentrated Extracts on Hep-G2 Cell Proliferation and Apoptosis. J. Agric. Food Chem. 2008, 56, 2603–2608. [Google Scholar] [CrossRef]
  11. Oleszek, M.; Oleszek, W. Saponins in Food; Springer: Berlin/Heidelberg, Germany, 2020. [Google Scholar]
  12. Furtado, N.A.J.C.; Pirson, L.; Edelberg, H.; Miranda, L.M.; Loira-Pastoriza, C.; Preat, V.; Larondelle, Y.; André, C.M. Pentacyclic Triterpene Bioavailability: An Overview of In Vitro and In Vivo Studies. Molecules 2017, 22, 400. [Google Scholar] [CrossRef] [Green Version]
  13. Chen, Y.H.; Bi, J.H.; Xie, M.; Zhang, H.; Shi, Z.Q.; Guo, H.; Yin, H.B.; Zhang, J.N.; Xin, G.Z.; Song, H.P. Classification-based strategies to simplify complex traditional Chinese medicine (TCM) researches through liquid chromatography-mass spectrometry in the last decade (2011–2020): Theory, technical route and difficulty. J. Chromatogr. A 2021, 1651, 462307. [Google Scholar] [CrossRef]
  14. Mao, Q.; Xu, J.D.; Kong, M.; Shen, H.; Zhu, H.; Zhou, S.S.; Li, S.L. LC-MS-based Metabolomics in Traditional Chinese Medicines Research: Personal Experiences. Chin. Herb. Med. 2017, 9, 14–21. [Google Scholar] [CrossRef]
  15. Gao, X.; Li, Y.; Meng, M.; Wang, P.; Feng, Y.; Jia, J.; Qin, X. Exploration of chemical composition and absorption characteristics of Chaigui granules based on UHPLC-Q-orbitrap-MS/MS. J. Pharm. Biomed. Anal. 2020, 187, 113293. [Google Scholar] [CrossRef]
  16. Dong, G.M.; Yu, H.; Pan, L.B.; Ma, S.R.; Xu, H.; Zhang, Z.W.; Han, P.; Fu, J.; Yang, X.Y.; Keranmu, A.; et al. Biotransformation of Timosaponin BII into Seven Characteristic Metabolites by the Gut Microbiota. Molecules 2021, 26, 3861. [Google Scholar] [CrossRef]
  17. Farag, M.A.; Porzel, A.; Wessjohann, L.A. Comparative metabolite profiling and fingerprinting of medicinal licorice roots using a multiplex approach of GC-MS, LC-MS and 1D NMR techniques. Phytochemistry 2012, 76, 60–72. [Google Scholar] [CrossRef]
  18. Pan, C.H.; Qi, D.Y.; Zhao, D.Y. Determination of Soyasapogenols Absorption Characteristics in Rat Gut Using in Situ Single-pass Intestinal Perfusion Model. Sci. Technol. Food Ind. (Chin.) 2021, 42, 336–342. [Google Scholar]
  19. Guo, X.; Cheng, M.; Hu, P.; Shi, Z.; Chen, S.; Liu, H.; Shi, H.; Xu, Z.; Tian, X.; Huang, C. Absorption, Metabolism, and Pharmacokinetics Profiles of Norathyriol, an Aglycone of Mangiferin, in Rats by HPLC-MS/MS. J. Agric. Food Chem. 2018, 66, 12227–12235. [Google Scholar] [CrossRef]
  20. Singh, D.; Chaudhuri, P.K. Structural characteristics, bioavailability and cardioprotective potential of saponins. Integr. Med. Res. 2018, 7, 33–43. [Google Scholar] [CrossRef]
  21. Yu, K.; Chen, F.; Li, C. Absorption, Disposition, and Pharmacokinetics of Saponins from Chinese Medicinal Herbs: What Do We Know and What Do We Need to Know More? Curr. Drug Metab. 2012, 13, 577–598. [Google Scholar] [CrossRef]
  22. Xu, T.; Li, S.; Sun, Y.; Pi, Z.; Liu, S.; Song, F.; Liu, Z. Systematically characterize the absorbed effective substances of Wutou Decoction and their metabolic pathways in rat plasma using UHPLC-Q-TOF-MS combined with a target network pharmacological analysis. J. Pharm. Biomed. Anal. 2017, 141, 95–107. [Google Scholar] [CrossRef] [PubMed]
  23. Ghosh, J.; Das, J.; Manna, P.; Sil, P.C. Arjunolic acid, a triterpenoid saponin, prevents acetaminophen (APAP)-induced liver and hepatocyte injury via the inhibition of APAP bioactivation and JNK-mediated mitochondrial protection. Free. Radic. Biol. Med. 2010, 48, 535–553. [Google Scholar] [CrossRef] [PubMed]
  24. He, C.; Zhou, D.; Li, J.; Han, H.; Ji, G.; Yang, L.; Wang, Z. Identification of 20(S)-protopanaxatriol metabolites in rats by ultra-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight tandem mass spectrometry and nuclear magnetic resonance spectroscopy. J. Pharm. Biomed. Anal. 2014, 88, 497–508. [Google Scholar] [CrossRef] [PubMed]
  25. Pan, C.H.; Qi, D.Y.; Zhao, D.Y. Isolation and purification of soyasapogenol A and B starting from soy hypocotyls: High-speed counter-current chromatography versus preparative high-performance liquid chromatography. J. Liq. Chromatogr. Relat. Technol. 2021, 44, 103–113. [Google Scholar] [CrossRef]
Figure 1. ESI–QTOF–MS/MS spectra of soyasapogenols A and B in the positive ion mode.
Figure 1. ESI–QTOF–MS/MS spectra of soyasapogenols A and B in the positive ion mode.
Molecules 28 00284 g001
Figure 2. The concentration–time curve of soyasapogenol after intravenous administration in rat plasma (n = 6).
Figure 2. The concentration–time curve of soyasapogenol after intravenous administration in rat plasma (n = 6).
Molecules 28 00284 g002
Figure 3. Plasma concentration–time curve of soyasapogenol after different dosages of oral administration. Left panel, soyasapogenol A: 10 and 20 mg/kg body weight (BW). Right panel, soyasapogenol B: 25 and 50 mg/kg BW. Values are mean ± SE; n = 6.
Figure 3. Plasma concentration–time curve of soyasapogenol after different dosages of oral administration. Left panel, soyasapogenol A: 10 and 20 mg/kg body weight (BW). Right panel, soyasapogenol B: 25 and 50 mg/kg BW. Values are mean ± SE; n = 6.
Molecules 28 00284 g003
Figure 4. Extracted ion chromatograms of soyasapogenol metabolites in bile.
Figure 4. Extracted ion chromatograms of soyasapogenol metabolites in bile.
Molecules 28 00284 g004
Figure 5. Proposed metabolism of soyasapogenols following oral and intravenous administration in rats.
Figure 5. Proposed metabolism of soyasapogenols following oral and intravenous administration in rats.
Molecules 28 00284 g005
Table 1. Summary of pharmacokinetic parameter values (mean, n = 6) for soyasapogenols A and B after oral (P.O.) dose administration.
Table 1. Summary of pharmacokinetic parameter values (mean, n = 6) for soyasapogenols A and B after oral (P.O.) dose administration.
ParameterUnitSoyasapogenol ASoyasapogenol B
10 mg/kg20 mg/kg25 mg/kg50 mg/kg
AUC(0-t)mg/L*h35.19888.61996.029157.218
AUC(0-∞)mg/L*h50.70095.507110.513178.341
T1/2zh2.5521.9662.2882.266
Tmaxh2122
CLz/FL/h/kg0.3940.2090.1810.112
Cmaxmg/L8.5928.9926.3740.29
Ratio of B/A in test sample - 2.52.5
Ratio of B/A in plasma level (AUC) - 2.171.87
Bioavailability (F,%) 73.1067.3460.9469.01
Cmax: peak concentration; Tmax: peak time; T1/2z: biological half-life; AUC: area under concentration–time curve; CLz/F: clearance; Ratio of B/A: ratio of group B soyasapogenol B to soyasapogenol A in test samples or in plasma at the same concentration.
Table 2. Summary of pharmacokinetic parameter values (mean, n = 6) for soyasapogenols A and B after intravenous (I.V.) dose administration.
Table 2. Summary of pharmacokinetic parameter values (mean, n = 6) for soyasapogenols A and B after intravenous (I.V.) dose administration.
ParameterUnitSoyasapogenol ASoyasapogenol B
10 mg/kg20 mg/kg25 mg/kg50 mg/kg
AUC(0-t)mg/L*h56.029114.47376.214154.392
AUC(0-∞)mg/L*h69.355141.82181.347258.42
T1/2zh3.5353.4755.1385.528
Tmaxh0.030.030.030.03
CLz/FL/h/kg0.2880.1410.2460.077
Cmaxmg/L31.7101.6543.2666.85
Table 3. Summary of the Mass Spectral Data of Soyasapogenol A and Its Metabolites Detected in Rat Bile.
Table 3. Summary of the Mass Spectral Data of Soyasapogenol A and Its Metabolites Detected in Rat Bile.
The Serial NumberMetabolitesMolecular Formula[M + H]+ m/zRetention Time (min)Fragment Ion InformationMetabolic Pathways
M0Soyasapogenol AC30H50O4475.380310.28185.1321, 283.2627, 333.2782, 337.2521Prototype
M2Soyasapogenol A–H2OC30H48O3457.36789.92166.2149, 187.1435, 403.2801Dehydration
M3Soyasapogenol A–OC30H50O3459.383110.36190.2349, 320.2433, 408.2959Deoxidation
M4Soyasapogenol A + H2O + C5H7NO3SC35H59NO8S654.40059.24117.0907, 292.1592, 462.2684, 534.3287Hydrolysis + with acetylcysteine
M5Soyasapogenol A + OC30H50O5491.37429.73129.1607, 255.3044, 273.2602, 347.3641Oxidation
M6Soyasapogenol A + H2C30H52O4477.393111.04273.1810, 315.2651, 359.2912Reduction
M7Soyasapogenol A + C5H7NO3SC35H57NO7S636.39467.29193.1584, 253.1942, 460.2505, 478.2607Binding to acetylcysteine
Soyasapogenol A–H2O + C3H5NOSC33H53NO4S560.37849.34375.3043, 393.3139, 439.2520Dehydration + binding to cysteine
M8Soyasapogenol A–H2O + C3H5NO2SC33H53NO5S576.37327.63119.0851, 285.2203, 306.1489, 522.2855Dehydration + with cysteine sulfhydryl group
M9Soyasapogenol A–H2O + C5H7NO3SC35H55NO6S618.38517.24255.2099, 357.2779, 423.3239, 462.2661Dehydration + binding with acetylcysteine
M10Soyasapogenol A–O–H2O + C3H5NOSC33H53NO3S544.38449.73273.2208, 485.2630Deoxidation + dehydration + binding to cysteine
M11Soyasapogenol A–O + C5H7NO3SC35H57NO6S620.40078.42141.1271, 340.1576, 365.2829Deoxidation + binding to acetylcysteine
M12Soyasapogenol A + O + C5H8N2O3SC35H58N2O8S667.39809.92337.2521, 459.2476, 565.3347, 609.3608Oxidation + binding to cysteinylglycine mercapto
M13Soyasapogenol A + O + C3H5NO2SC33H55NO7S610.37549.2689.0595, 133.0857, 247.2050, 333.2428Oxidation + binding to cysteine sulfhydryl groups
M14Soyasapogenol A + O + SO3C30H50O8S571.32988.73209.1163, 277.2159, 335.2600, 427.3197Oxidation + sulfation
M15Soyasapogenol A + H2 + C6H8O6C36H60O10653.42539.68449.2132, 491.2973, 535.3234, 579.3498Reduction + gluconaldehyde acidification
M16Soyasapogenol A + H2 + HPO3C30H53O7P557.36069.20124.9994, 291.2311, 339.1717Reduction + phosphorylation
Table 4. Summary of the Mass Spectral Data of Soyasapogenol B and Its Metabolites Detected in Rat Bile.
Table 4. Summary of the Mass Spectral Data of Soyasapogenol B and Its Metabolites Detected in Rat Bile.
The Serial NumberMetabolitesMolecular Formula[M + H]+ m/zRetention Time (min)Fragment Ion InformationMetabolic Pathways
M0Soyasapogenol BC30H50O3459.383810.94163.1106, 283.2627, 401.2284, 429.3716Prototype
M1Soyasapogenol B + H2OC30H52O4477.392611.04185.1252, 342.2854, 370.3192Hydrolysis
M2Soyasapogenol B–H2O + C3H5NOSC33H53NO3S544.38489.81199.1472, 239.2363, 515.3103Dehydration + binding to cysteine
M3Soyasapogenol B + H2O + HPO3C30H53O7P557.36019.19139.1112, 253.1939, 413.3035, 507.2914Hydrolysis + phosphorylation
M4-1Soyasapogenol B–H2 + C3H5NOSC33H53NO4S560.37599.0689.0596, 149.0956, 411.2872, 465.3320Dehydrogenation + binding to cysteine
M4-2Soyasapogenol B–H2O + C3H5NO2SC33H53NO4S560.37789.34145.1008, 199.1474, 427.3192, 490.2986Dehydration + with cysteine sulfhydryl group
M5-1Soyasapogenol B + O-H2 + C3H5NOSC33H53NO5S576.37347.84133.0857, 273.2206, 481.2611, 525.2875Oxidation + dehydrogenation + binding to cysteine
M5-2Soyasapogenol B–H2 + C3H5NO2SC33H53NO5S576.37348.41215.1785, 339.2672, 409.3090, 488.2865Dehydrogenation + binding to the thiol group of cysteine
M6-1Soyasapogenol B–H2 + O2 + C3H5NOSC33H53NO6S592.36906.21121.0644, 179.1064, 351.2526, 427.2812Dehydrogenation + deoxidation + binding to cysteine
M6-2Soyasapogenol B + O–H2 + C3H5NO2SC33H53NO6S592.36928.29152.0374, 337.2518, 462.2663, 480.2770Oxidation + dehydrogenation + binding to cysteine sulfhydryl groups
M7Soyasapogenol B–O + C5H7NO3SC35H57NO5S604.40469.33149.0956, 231.1738, 275.1996, 526.2961Deoxidation + binding to acetylcysteine
M8Soyasapogenol B–H2 + C5H7NO3SC35H59NO5S606.42177.76103.0749, 276.1252, 339.2673, 464.2819Deoxidation + reduction + binding to acetylcysteine
M9Soyasapogenol B–O + H2 + C5H7NO3SC35H55NO6S618.38437.23163.1113, 243.2094, 393.2793, 502.2581Dehydrogenation + combining with cysteinylglycine mercapto
M10Soyasapogenol B + C5H7NO3SC35H57NO6S620.39977.66261.1839, 357.2776Binding to acetylcysteine
M11Soyasapogenol B–H2 + C5H8N2O3SC35H56N2O6S633.39409.0489.0594, 287.1460, 407.2023, 520.3314Dehydrogenation + binding to acetylcysteine
M12Soyasapogenol B + O–H2 + C5H7NO3SC35H55NO7S634.37946.43215.1424, 355.2622, 480.2767, 503.3029Oxidation + dehydrogenation + gluconaldehyde acidification
M13Soyasapogenol B + O + C5H7NO3SC35H57NO7S636.38997.29133.0853, 210.0425, 289.2152, 462.2657Oxidation + gluconaldehyde acidification
M14Soyasapogenol B–H2 + O2 + C5H7NO3SC35H55NO8S650.37386.8697.0643, 205.1216, 353.2458, 522.2850Dehydrogenation + deoxidation + gluconaldehyde acidification
M15Soyasapogenol B + H2O + C6H8O6C36H60O10653.42809.68221.2259, 361.1606, 518.3208, 546.3546Hydrolysis + gluconaldehyde acidification
Table 5. Dose and composition of the tested samples in each treatment group.
Table 5. Dose and composition of the tested samples in each treatment group.
TreatmentsBulk SolutionSoyasapogenol CompositionLow DosageHigh Dosage
Dose (μ mol/kg BW)Dose (mg/kg BW)Dose (μ mol/kg BW)Dose (mg/kg BW)
Oral administration (P.O.)0.5% (w/v) sodium carboxymethyl cellulose solution used as a vehicleA21104220
B552510950
Intravenous injection (I.V.)0.5% (w/v) sodium carboxymethyl cellulose solution used as a vehicleA21104220
B552510950
Perfusion of the small intestineKrebs–Ringer bicarbonate bufferA(B)10 mg/L(50 mg/L)
BW, body weight.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pan, C.; Yan, Y.; Zhao, D. The Fate and Intermediary Metabolism of Soyasapogenol in the Rat. Molecules 2023, 28, 284. https://doi.org/10.3390/molecules28010284

AMA Style

Pan C, Yan Y, Zhao D. The Fate and Intermediary Metabolism of Soyasapogenol in the Rat. Molecules. 2023; 28(1):284. https://doi.org/10.3390/molecules28010284

Chicago/Turabian Style

Pan, Chenghui, Yonggang Yan, and Dayun Zhao. 2023. "The Fate and Intermediary Metabolism of Soyasapogenol in the Rat" Molecules 28, no. 1: 284. https://doi.org/10.3390/molecules28010284

APA Style

Pan, C., Yan, Y., & Zhao, D. (2023). The Fate and Intermediary Metabolism of Soyasapogenol in the Rat. Molecules, 28(1), 284. https://doi.org/10.3390/molecules28010284

Article Metrics

Back to TopTop