• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    An integrated strategy for comprehensive characterization of metabolites and metabolic profiles of bufadienolides from Venenum Bufonis in rats

    2022-04-07 12:46:50WenLongWeiHoJvLiWenZhiYngHuQuZhenWeiLiChngLingYoJinJunHouWnYingWuDeAnGuo
    Journal of Pharmaceutical Analysis 2022年1期

    Wen-Long Wei,Ho-Jv Li,Wen-Zhi Yng,Hu Qu,Zhen-Wei Li,b,Chng-Ling Yo,Jin-Jun Hou,Wn-Ying Wu,*,De-An Guo,b,**

    aShanghai Research Center for Modernization of Traditional Chinese Medicine,National Engineering Laboratory for TCM Standardization Technology,Shanghai Institute of Materia Medica,Chinese Academy of Sciences,Shanghai,201203,China

    bUniversity of Chinese Academy of Sciences,Beijing,100049,China

    Keywords:

    Metabolic profiles

    Extension-mass defect filter

    Multidimensional data acquiring

    Metabolic network prediction

    Bufadienolides of Venenum Bufonis

    A BS T R A C T

    Comprehensive characterization of metabolites and metabolic profiles in plasma has considerable significance in determining the efficacy and safety of traditional Chinese medicine(TCM)in vivo.However,this process is usually hindered by the insufficient characteristic fragments of metabolites,ubiquitous matrix interference,and complicated screening and identification procedures for metabolites.In this study,an effective strategy was established to systematically characterize the metabolites,deduce the metabolic pathways,and describe the metabolic profiles of bufadienolides isolated from Venenum Bufonis in vivo.The strategy was divided into five steps.First,the blank and test plasma samples were injected into an ultra-high performance liquid chromatography/linear trap quadrupole-orbitrap-mass spectrometry(MS)system in the full scan mode continuously five times to screen for valid matrix compounds and metabolites.Second,an extension-mass defect filter model was established to obtain the targeted precursor ions of the list of bufadienolide metabolites,which reduced approximately 39% of the interfering ions.Third,an acquisition model was developed and used to trigger more tandem MS(MS/MS)fragments of precursor ions based on the targeted ion list.The acquisition mode enhanced the acquisition capability by approximately four times than that of the regular data-dependent acquisition mode.Fourth,the acquired data were imported into Compound Discoverer software for identification of metabolites with metabolic network prediction.The main in vivo metabolic pathways of bufadienolides were elucidated.A total of 147 metabolites were characterized,and the main biotransformation reactions of bufadienolides were hydroxylation,dihydroxylation,and isomerization.Finally,the main prototype bufadienolides in plasma at different time points were determined using LC-MS/MS,and the metabolic profiles were clearly identified.This strategy could be widely used to elucidate the metabolic profiles of TCM preparations or Chinese patent medicines in vivo and provide critical data for rational drug use.

    1.Introduction

    Metabolism is undoubtedly significantly correlated with the efficacy and safety of drugs,especially those used in traditional Chinese medicine(TCM)[1-3].However,the diversity,complexity,and low content of active agents make global characterization of metabolites and elucidation of their metabolic profiles in vivo very challenging[4,5].The challenge of elucidating these substances is usually associated with the following disadvantages:1)the low content of prototype compounds in plasma,2)the trace amount of metabolites in vivo[6,7],3)the diversity of metabolic pathways[8],and 4)the restricted generation characteristic fragments with coelution of diverse metabolites and matrix interference.

    Recently,hyphenated techniques such as liquid chromatography/mass spectrometry(LC-MS)have become an effective platform for the rapid characterization of metabolites and elucidation of metabolic profiles in vivo owing to their powerful separation capacity,remarkable sensitivity,and high throughput[9].In particular,ultra-high performance LC(UHPLC)coupled with highresolution MS(HR-MS)has often been used to characterize the metabolites in biological samples based on accurate mass and MS/MS information[10-12].UHPLC coupled with triple quadrupoletandem MS has been used to quantitate the concentration of analytes in plasma for metabolic profile elucidation based on the multiple reaction monitoring(MRM)mode[13-15].

    For chemical characterization based on HR-MS,two models were used to acquire MS data and characteristic fragments.The first was data-independent acquisition(DIA),which provided highthroughput fragmentation information for the structural elucidation of chemical components[16-18].However,the deconvolution and attribution of diagnostic fragments were laborious and tedious[16,17].Compared with DIA,the second mode,data-dependent acquisition(DDA)combined with diverse selection criteria such as the precursor ion list(PIL)[19],exclusion list[20],neutral loss[21],product ion[22],and mass defect filter(MDF)[23],may be a superior option to improve sensitivity and selectivity in obtaining MS/MS or MSndata[24-26].

    However,DDA scans usually inevitably miss fractional PIL information,which does not trigger MS/MS fragmentation.Recently,gas-phase fractionation with staggered mass range and timestaggered PIL models have been developed to notably enhance MS/MS coverage to improve the identification of chemical components compared to that achieved with the conventional DDA mode[27].Furthermore,a sequentially step-wise targeted MS/MS method was used to acquire the MS/MS spectra of all the plasma metabolites for targeted metabolomics studies[28].

    The time-staggered PIL was also developed to improve the performance of MS/MS acquisition in untargeted metabolomics[29].Regarding the post-processing of MS data,a number of remarkable process models,such as the diagnostic product ion network[30],mass defect filtering[31],mass spectral tree substructure recognition and statistical analysis[32],mass spectral tree similarity filtering[33],and diagnostic ion-guided network bridging[34],could dramatically simplify the MS spectral elucidation and make the characterization convenient.In addition,sectional commercial software such as UNIFI[35,36],Mass Frontier[37],and Compound Discoverer[38]could be used to match the measured MS and MS/MS information with theoretical data obtained from the prediction and in-house libraries.This is achieved using the scoring algorithms and structure matching to obtain the metabolic pathway prediction in an automatic mode for the identification of metabolites.For the elucidation of metabolic profiles,MRM with QqQ-MS enabled the detection of multiple analytes simultaneously covering an extensive concentration range in biological samples[39,40].Furthermore,the MRM model exhibited high sensitivity and specificity,and has been widely used for pharmacokinetic and metabolic profile studies[41,42].

    Venenum Bufonis(VB),the dry secretions of Bufo gargarizans Cantor or Bufo melanostictus Schneider as stipulated in the Chinese Pharmacopoeia(2015 edition),has been used to treat sunstroke,fainting,and acute filthy disease such as abdominal pain or vomiting and diarrhea.Bufadienolides are the main active constituents that exhibit diverse antitumor activities against numerous cancer cells[43-45],but their strong cardiotoxicity and neurotoxicity limit their clinical applications[46,47].The efficacy and toxicity of these compounds are closely related to their metabolism,so the comprehensive characterization of their metabolites and metabolic profiles is important for drug evaluation in vivo.However,no effective strategy has been developed to systematically characterize the metabolites of bufadienolides,and few studies have integrally evaluated the metabolic profiles of VB.

    In this study,an integral and convenient five-step strategy was established to study the metabolic profile of VB.First,test and blank plasma samples were detected using UHPLC/linear ion trap(LTQ)-Orbitrap MS in the full scan mode,and the data were used to screen the valid metabolites whereas the blank was deducted using MZmine2.Second,an extension-MDF(E-MDF)model was established to screen the targeted PIL.Third,the multidimensional data acquisition model was optimized and applied to trigger more MS/MS fragmentation of metabolites.Fourth,the obtained data were imported into Compound Discoverer software to identify metabolites.Finally,the main prototypes were simultaneously quantitated for the entire metabolic profile study of the VB sample.

    2.Experimental

    2.1.Chemicals and materials

    HPLC-grade acetonitrile and methanol were purchased from Merck KGaA(Merck,Darmstadt,Germany).Formic acid(FA)was purchased from Sigma-Aldrich(St.Louis,MO,USA),and deionized water was prepared using a Millipore Alpha-Q water purification system(18.2 MΩ?cm at 25°C;Millipore,Bedford,MA,USA).VB was purchased from Sanyitang Chinese Medicine Decoction Slice Co.,Ltd.(Bozhou,China).Authentication of the drug materials was confirmed using thin layer chromatography and HPLC assays according to the Chinese Pharmacopoeia(2015edition).A voucher specimen(CS20161101)of the plant was deposited at the Shanghai Institute of Materia Medica,Chinese Academy of Sciences(Shanghai,China).The following 14 bufadienolides reference standards were isolated from the VB sample:pseudobufarenogin(CS01),arenobufagin(CS03),hellebrigenin(CS04),19-oxocinobufotalin(CS11),telocinobufagin(CS14),bufotalin (CS15),resibufagin(CS16),19-oxocinobufagin(CS17),cinobufotalin(CS19),bufalin(CS21),resibufogenin(CS22),cinobufagin(CS23),desacetylcinobufagin(CS30),and gamabufotalin(CS36)with purity>98% as determined using HPLC-ultraviolet(UV)spectroscopy.Seven bufadienolides reference standards,namely,desacetylbufotalin,cinobufaginol,19-oxobufalin,1beta-hydroxybufalin,bufarenogin,marinobufagenin,and bufotalinin,were isolated from the VB sample with purities>90%,as determined using HPLC-UV spectroscopy.BF211 was synthesized in our laboratory for use as the internal standard(IS)and the purity was>98%,as determined using HPLC-UV spectroscopy.The structures of the various compounds are shown in Table S1.

    2.2.Blood collection and sample preparation

    Male Wistar rats(200±20 g)purchased from SLAC Lab Animal Center(Shanghai,China)were used in the various experiments according to the protocol approved by the Review Committee of Animal Care and Use at Shanghai Institute of Materia Medica(Shanghai,China).Wistar rats were housed in a room with controlled temperature(20-24°C),relative humidity(40%-70%),and unidirectional airflow,under a 12 h light/dark cycle.The rats were fed food and water ad libitum for 7 days before the experiments and allowed to acclimate to the environment.The 0.5% dimethyl sulfoxide was used to dissolve bufadienolides preliminarily,and 5 mg/mL sodium carboxymethylcellulose was used for dispersion of bufadienolides.The rats were intragastrically administered the crude extract at a dose of 200 mg/kg.Blood samples were collected in heparinized tubes on ice at 0,0.083,0.25,0.5,0.75,1,2,4,6,8,10,and 24 h after intragastric administration and then immediately centrifuged at 3,500 r/min for 10 min.The supernatants were transferred into new tubes and stored at-20°C until analysis.

    Sample preparation for qualitative analysis:1 mL of ethyl acetate and 100μL of plasma sample were added to a centrifuge tube.The mixture was vortexed for 3 min,centrifuged for 10 min at 3,500 r/min,and then 900μL of the supernatant was transferred into a new tube and evaporated to dryness under a nitrogen gas(N2)stream.Subsequently,100μL of a 20% acetonitrile-water(V/V)solution was used to redissolve the residue.Finally,the solution was centrifuged at 14,000 r/min for 10 min after vortexing for 3 min,and 10μL of supernatant was injected into the UHPLC for analysis.

    Sample preparation for quantitative analysis:the 14 bufadienolides reference standards were dissolved in methanol to a concentration of 1 mg/mL as stock solutions,which were each further diluted with 50% methanol-water(1:1,V/V)to obtain seven standard working solutions of varying concentrations and low,medium,and high quality control(QC)samples.For the standard curve,each standard solution(5 μL),blank plasma(95 μL),and IS working solution(10μL)were added to a 1.5 mL centrifuge tube,vortexed,and then 1 mL of ethyl acetate was added for extraction.The mixture was vortexed and centrifuged at 3,500 r/min for 10 min,and then the supernatant(900μL)was transferred into a 1.5 mL tube and evaporated to dryness at 40°C.The residue was redissolved in 100μL of 20% acetonitrile,and the supernatant was transferred to a sample vial for analysis after centrifugation.For the experimental samples,100 μL of plasma samples,10 μL of IS working solution,and 1 mL of ethyl acetate were added to a 1.5 mL centrifuge tube,followed by the same steps described above.

    2.3.UHPLC/LTQ-Orbitrap MS

    An Ultimate?3000 UHPLC system(Thermo Fisher Scientific,San Jose,CA,USA)equipped with an ACQUITY UPLC?HSS T3 column(2.1 mm×100 mm,1.8μm)was used for the separation.The column temperature was 30°C and the flow rate was 0.5 mL/min.The elution was performed using a binary mobile phase consisting of 0.1% FA-water(A)and acetonitrile(B)with the following gradient program:0-2 min,0-5%B;2-4 min,5%-17%B;4-6 min,17%-28%B;6-9 min,28%B;9-14 min,28%-35%B;14-20 min,35%-39%B;20-30 min,39%-95%B;and 30-35 min,95%B.

    An LTQ-Orbitrap Velos Pro hybrid mass spectrometer(Thermo Fisher Scientific,San Jose,CA,USA)connected to the UHPLC system via an electrospray ionization source was used to obtain HR-MS data in the positive mode.The parameters used were as follows:spray voltage,3 kV(Fig.S1);capillary temperature,320°C;source heater temperature,300°C;sheath gas(N2),40 a.u.;and auxiliary gas(N2),10 a.u..In scan event 1,the conditions were as follows:analyzer,Fourier transform MS;mass range,m/z 100-1500;scan type,full;and resolution,30,000.In event 2,the conditions were as follows:activation type,collision-induced dissociation;default charge state,1;isolation width,m/z 2;normalized collision energy,30 V(Fig.S2);activation,0.25;and activation time,10 ms.

    2.4.Generation of target PIL and data acquisition

    The plasma samples were injected into the UHPLC/LTQ-Orbitrap mass spectrometer in the full-scan mode with five replicate injections.The raw data format was transferred to mzXML using MSconvert and imported into MZmine2 for the chromatographic peak extraction,alignment,and data analysis using default settings.The chromatographic peak that appeared three times and the relative standard deviation(RSD)of the peak area that was not>30% were considered the stable precursor ion.The PIL of the blank plasma samples was obtained using the same process.

    To acquire the targeted PIL,the PIL obtained after deducting the blank was further screened using an E-MDF.Then,the plasma samples were injected into the UHPLC/LTQ-Orbitrap mass spectrometer to trigger the MS/MS fragments with the targeted PIL.The data acquisition method of the MS/MS process and the multidimensional data acquisition method were selected.

    2.5.Elucidation of bufadienolide metabolic pathways

    The Compound Discoverer software is an efficient tool for the identification of metabolites.A five-step metabolic network prediction method was used to elucidate the metabolic pathways and the workflow is shown in Fig.1.First,the structure files of bufadienolides(in-house library)were imported into the software to generate the theoretical metabolite library.The theoretical pathways of the bufadienolides were dehydration,desaturation,hydration,oxidation,oxidative deamination to alcohol,oxidative deamination to ketone,reduction,conversion of thiourea to urea in a phase I metabolic reaction,acetylation,arginine conjugation,cysteine conjugation 1 (C3H5NO2S),cysteine conjugation 2(C3H7NO2S),glucoside conjugation,glucuronide conjugation,glutamine conjugation,glycine conjugation,GSH conjugation 1(C10H15N3O6S),GSH conjugation 2 (C10H17N3O6S),ornithine conjugation,palmitoyl conjugation,stearyl conjugation,sulfation,and taurine conjugation in phase II metabolic reactions.Then,the peak alignment and peak extraction were implemented after importing the MS/MS data into the Compound Discoverer.Third,a theoretical metabolite library and database supplied by Thermo with a mass tolerance of 5 ppm and a threshold of S/N≥3 was screened for possible metabolites.Finally,the match degree between the predicted and theoretical metabolites was evaluated by comparing the acquired MS/MS fragments and predicted characteristic fragments,and the metabolite list was ranked according to Fragment Ion Search(FISh)scores.

    Fig.1.Workflow of Compound Discoverer software.MS:mass spectrometry;FISh:Fragment Ion Search.

    Metabolic network prediction was useful for identifying the metabolites.The MS/MS spectra of bufalin and its hydroxylated metabolite(compound B)are presented as an example(Fig.2).Fig.2 shows that the characteristic peaks were obviously different between bufalin and compound B,with a difference of 16 Da.Thus,compound B was likely the hydroxylated metabolite of bufalin.

    Fig.2.Tandem mass spectrometry(MS/MS)spectra of bufalin and compound B.

    2.6.Metabolic profile study using UHPLC-MS/MS

    A UHPLC-MS/MS system consisting of an Agilent UHPLC system(Agilent Technologies,Palo Alto,CA,USA)equipped with a column manager,sample manager,binary pump,and Applied Biosystems 4000 QTRAP?LC-MS/MS system(Toronto,Canada)was used in the metabolic profile study of bufadienolides.The ACQUITY UPLC?HSS T3 column(2.1 mm×100 mm,1.8μm)was used for the separation,and the eluent was composed of 0.1% FA-water(V/V,A)and acetonitrile(B).The gradient elution program was as follows:0-2 min,5%B;2-4 min,5%-17%B;4-6 min,17%-28%B;6-9 min,28%B;9-14 min,28%-35%B;14-20 min,35%-39%B;20-21 min,39%-95%B;and 21-25 min 95%B.The flow rate was 0.4 mL/min.Other parameters were as follows:column temperature,40°C;injection volume,10μL;curtain gas,35 psi;ion spray voltage,5500 V;ion spray temperature,500°C;heater gas,50 psi;nebulizer gas,50 psi;and collision gas,medium.The AB Sciex Analyst 1.6.3 software was used for instrument control,data acquisition,and raw data processing.The pharmacokinetic parameters were analyzed using the Drug and Statistics(DAS,version 2.0,Anhui,China)software and the concentration-time curve of the bufadienolides was constructed using GraphPad Prism 5.01.

    2.7.UHPLC-MS/MS method validation

    The method was developed in accordance with the US Food and Drug Administration guidelines.The specificity,accuracy,precision,recovery,matrix effect,dilution effect,and stability were evaluated as described in the following sections.

    2.7.1.Specificity

    The specificity of the following bufadienolides was evaluated at the lower limit of quantitation(LLOQ):pseudobufarenogin(CS01),arenobufagin(CS03),hellebrigenin(CS04),19-oxocinobufotalin(CS11),telocinobufagin (CS14),bufotalin (CS15),resibufagin(CS16),19-oxocinobufagin(CS17),cinobufotalin(CS19),bufalin(CS21), resibufogenin (CS22), cinobufagin (CS23),desacetylcinobufagin(CS30),and gamabufotalin(CS36).For the LLOQ,the peak area of the blank was expected to be<20% of that of the LLOQ at the retention time of the analytes,and the peak area of the IS was expected to be<5%.

    2.7.2.Accuracy and precision

    The QCs including the LLOQ and low-,medium-,and highconcentration(3,25,and 160 ng/mL)QC samples were analyzed over 3 consecutive days(n=6)to determine the intra-and inter-day precisions.The accuracy was calculated as the relative error(RE;RE=absolute error/true value),which was expected to be within±15%.The precision was calculated as the RSD,which was expected to be within±20% for the LLOQ and within±15% for other QCs.

    2.7.3.Recovery,matrix effect,and dilution effect

    The recovery was investigated by calculating the ratio of the mean peak areas of the extracted QC samples to that of the unextracted sample,and the RSD was expected to be within±15%.The matrix effect was evaluated by comparing the peak areas of analytes resolved in extracted blank plasma with those in which blank plasma was substituted with water at three concentration levels.The matrix effect was considered acceptable if the RSD was not more than 15%.The concentration of CS03 was above the linear range in plasma so the dilution effect was evaluated.

    2.7.4.Stability

    The stability of 14 bufadienolides was evaluated by evaluating the QC samples following storage under different conditions at low(3 ng/mL)and high(160 ng/mL)concentration levels(n=6).The autosampler and freeze-thaw stability(short-term and long-term freeze-thaw stability)were evaluated by analyzing the samples following storage at an ambient temperature for 4 h,4°C for 24 h in an autosampler,freeze-thaw three times with every freeze-thaw cycle not less than 12 h.The analytes were considered stable when the bias(%)between the measured and theoretical QC values was not more than±15%.

    3.Results and discussion

    3.1.Construction of E-MDF for screening target PIL

    MDF is a powerful data processing method for removing a mass of interfering ions.To eliminate interfering ions in the PIL,an EMDF was developed using the expanded metabolite library(Fig.3,Table 1)for screening the targeted PIL.The E-MDF was established in four steps.First,140 bufadienolides(bufogenin without the side chain),which are the primary bioactive constituents of VB,were screened as prototype compounds.The mass(m/z)of these bufadienolides was previously characterized(Table S2)[48].Second,the 140 prototype bufadienolides were extended with 44 regularmetabolic pathways(Table S3),and an in-house library containing 6160 metabolites was obtained.Then,the E-MDF was established with a mass(m/z)of 6160 metabolites.Finally,the 4954 precursor ions acquired from plasma samples after deducting the blank using the full scan mode were further screened using E-MDF.A total of 3040 precursor ions were selected as the targeted precursor ions(Table S4).This model obviously reduced approximately 39% of the interfering ions.

    Fig.3.Extension-mass defect filter.Green,blue,and orange represent 140 prototype bufadienolides,6160 predicted metabolites,and mass-to-charge ratio(m/z)acquired in full scan model after blank deduction,respectively.

    Table 1Inequality of extension-mass defect filter.

    3.2.Optimization of multidimensional data acquisition

    To trigger more MS/MS fragmentation of precursor ions,the different data acquisition modes were compared and evaluated,and the results are presented in Fig.4.As shown in Fig.4A,the mass segmentation acquisition mode indicated that the target precursor ion list was divided into 3-6 mass sections in the range of m/z 331-784,and the amount of precursor ions that triggered MS/MS fragmentation increased slightly.In Fig.4B,the target PIL was divided into three individual ion lists based on the staggered time points,which obviously improved the acquisition ability of MS/MS fragments.There were 671 precursor ions that triggered MS/MS fragmentation.To obtain additional MS/MS information,the mass segmentation acquisition mode and the time segmentation acquisition mode were combined,and as many as 1000 precursor ions were triggered MS/MS fragmentation.Considering the efficiency and capability,the mass range,which was divided into four sections,was selected as mass segmentation to combine with time segmentation.The results showed that the multidimensional data acquisition enhanced the acquisition capability by approximately four times than that of the regular DDA mode,and increased the single segmentation pattern by approximately1.5-3 times(Fig.4C).

    Fig.4.(A)Mass segmentation mode,(B)time segmentation mode,and(C)amount of triggered MS/MS fragmentation in different data acquisition modes.

    3.3.Characterization of metabolites and metabolic pathways

    All raw data were imported into the Compound Discoverer for characterization of the metabolites with metabolic network prediction.After merging the identified metabolites with the same retention time and deleting those with low FISh scores,147 components were identified,and the following 15 prototype bufadienolides:arenobufagin,bufalin,bufarenogin,cinobufagin,bufotalin,cinobufotalin,telocinobufagin,resibufogenin,cinobufaginol,desacetylbufotalin, desacetylcinobufagin, hellebrigenin, marinobufagin,resibufagin,and pseudobufarenogin were detected in the rat plasma(Table S5).Then,the metabolic pathways of the seven main bufadienolides were elucidated(Figs.S3-S9).Taking cinobufotalin as an example,deacetylation,isomerization,acetylation,and oxidation reactions occurred in the first step.The deacetylated cinobufotalin metabolite was further metabolized through isomerization,oxidation,hydrolysis,and reoxidation(Fig.S6).In summary,the main identified metabolic pathways of the bufadienolides were hydroxylation,dihydroxylation,and isomerization.

    3.4.UHPLC-MS/MS method validation

    Fig.S10 shows the specificity of the 14 bufadienolides and number(1-3)shows the 14 reference standards in the blank plasma samples with and without the standards at LLOQ concentration and plasma samples from rats intragastrically administered the crude VB extract.The results indicated that bufadienolides exhibited satisfactory specificity and no obvious endogenous interference or carryover effect at the retention time.The linear ranges of CS01,CS03,CS04,CS17,CS30,and CS36 were from 2.0 to 200 ng/mL and those of CS11,CS14,CS15,CS16,CS19,CS21,CS22,and CS23 were 1.0-200 ng/mL,indicating that they possessed good linearity(correlation coefficient,r>0.995,Table 2).The results of the stability,precision,accuracy,recovery,and matrix effect analyses of the bufadienolides,which are shown in Tables S6-S8,respectively,were acceptable.The dilution effect of CS03 was evaluated because its concentration in the plasma was over the linearity range,and no obvious dilution effect was observed(the value of dilution effect was 101.9%).

    Table 2Linear equation and correlation coefficient(r)values of bufadienolides.

    3.5.Metabolic profiles of VB

    The MRM mode was used for the quantitative analysis,and the ion pairs of the 14 bufadienolides are listed in Table 3.Fig.5 shows the concentration-time curves of 10 bufadienolides and the pharmacokinetic parameters of the main bufadienolides,except for CS15,CS16,CS17,and CS19,are listed in Table 4.The results showed that the absorption and the elimination of CS21 and CS22 were rapid(time to reach the maximum concentration(tmax)was(0.396±0.292)h and(0.333±0.204)h;half-life(t1/2)values were(1.775±0.821)h and(1.195±0.477)h,respectively),and the result was consistent with that of previous studies[49].CS15,CS16,CS17,and CS19 were not detected in rat plasma samples over the sampling period(their concentrations were lower than that of the LLOQ).CS03 displayed the highest maximum concentration(Cmax;(466.5±116.920)ng/mL)and the highest area under curve from time 0 to t(AUC(0-t);(1487.316 ± 407.361)ng/mL·h)in rat plasma,probably because it was one of the main components of the bufadienolides and had the highest concentration in the crude extract.Moreover,CS30 showed the longest,t1/2of CS30 was(10.849±5.235)h,which might have a close relationship with the metabolism of CS30,which was one of the metabolites of CS23.In addition,almost all the main components of the bufadienolides reached their maximum concentrations in 1 h except for CS04,which had a tmaxof(1.438±0.657)h.

    Table 3Ion pairs of main bufadienolides in multiple reaction monitoring modea.

    Table 4Pharmacokinetic parameters of 10 bufadienolides(n=4).

    Fig.5.Concentration-time curves of 10 bufadienolides.

    4.Conclusion

    Elucidation of the metabolic profiles of TCM preparations could clarify the in vivo metabolism and provide significant data for their clinical application.To this end,an efficient strategy was developed to characterize the metabolites of VB and describe their metabolic profiles in vivo.The strategy is summarized into three primary points:1)E-MDF screening,multidimensional data acquisition,and metabolic network prediction were used to identify the metabolites in rat plasma.In total,147 metabolites were tentatively characterized,and the main metabolic pathways were hydroxylation,dihydroxylation,and isomerization.2)A method for the simultaneous quantitation of 14 prototype bufadienolides in vivo was established and validated.3)Concentration-time curves of 10 bufadienolides were constructed after intragastric administration of 200 mg/kg crude VB extract to rats.The results showed that CS03 exhibited the highest concentration in rat plasma and the elimination of bufadienolides was rapid.This strategy could be used to elucidate the dynamic change of VB in vivo and provide critical data for drug development and application.

    CRediT author statement

    Wen-Long Wei:Conceptualization,Methodology,Writing-Reviewing and Editing;Hao-Jv Li:Data curation,Visualization,Formal analysis,Writing-Original draft preparation;Wen-Zhi Yang:Conceptualization;Hua Qu:Supervision;Zhen-Wei Li:Validation,Formal analysis;Chang-Liang Yao:Resources,Writing-Reviewing and Editing;Jin-Jun Hou:Software,Data curation;Wan-Ying Wu:Project administration;De-An Guo:Funding acquisition.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China(Grant Nos.:81530095 and 81673591),Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.:XDA12020348),National Standardization of Traditional Chinese Medicine Project(Grant No.:ZYBZH-K-LN-01),and Science and Technology Commission Foundation of Shanghai(Grant No.:15DZ0502800).

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jpha.2021.02.003.

    国产淫语在线视频| 偷拍熟女少妇极品色| videossex国产| 久久精品国产亚洲网站| 欧美变态另类bdsm刘玥| 国产高潮美女av| 亚洲伊人久久精品综合| 啦啦啦视频在线资源免费观看| 少妇人妻精品综合一区二区| av专区在线播放| 黄色日韩在线| 亚洲国产av新网站| 男男h啪啪无遮挡| 免费看日本二区| av视频免费观看在线观看| 国产精品久久久久久久电影| 亚洲综合色惰| 最后的刺客免费高清国语| 如何舔出高潮| 久久鲁丝午夜福利片| 亚洲精品,欧美精品| 欧美国产精品一级二级三级 | 丝袜脚勾引网站| 乱码一卡2卡4卡精品| 国产黄片美女视频| 日本wwww免费看| 一区二区三区四区激情视频| 在线亚洲精品国产二区图片欧美 | 国产色婷婷99| 亚洲av日韩在线播放| 最近最新中文字幕免费大全7| 午夜福利视频精品| 老熟女久久久| 久久国产精品大桥未久av | 91久久精品国产一区二区成人| 久久99热6这里只有精品| 久久婷婷青草| 国产乱人视频| 久久ye,这里只有精品| 熟女人妻精品中文字幕| 亚洲人成网站在线观看播放| 午夜日本视频在线| 国产亚洲最大av| 最近2019中文字幕mv第一页| 午夜免费男女啪啪视频观看| 男的添女的下面高潮视频| 91aial.com中文字幕在线观看| 国产真实伦视频高清在线观看| 国产真实伦视频高清在线观看| 国产亚洲最大av| 国产成人aa在线观看| 日韩电影二区| 亚洲激情五月婷婷啪啪| 久久久精品94久久精品| 免费在线观看成人毛片| 久久av网站| 国产亚洲5aaaaa淫片| 亚洲精品国产色婷婷电影| 日韩 亚洲 欧美在线| 午夜免费观看性视频| 91久久精品电影网| 欧美少妇被猛烈插入视频| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验| 国产亚洲5aaaaa淫片| 国产精品精品国产色婷婷| 国产乱来视频区| kizo精华| 日本色播在线视频| 美女福利国产在线 | freevideosex欧美| 久久人人爽av亚洲精品天堂 | 亚洲伊人久久精品综合| 国语对白做爰xxxⅹ性视频网站| 日日啪夜夜爽| 插逼视频在线观看| 日韩大片免费观看网站| 女性被躁到高潮视频| 成人综合一区亚洲| 成人毛片60女人毛片免费| 久久人人爽人人爽人人片va| 国产 精品1| 亚洲av中文字字幕乱码综合| 国产黄频视频在线观看| 高清不卡的av网站| 精品亚洲成国产av| 永久免费av网站大全| 国产精品免费大片| 国产精品一二三区在线看| 亚洲国产精品999| 黄片无遮挡物在线观看| 啦啦啦视频在线资源免费观看| 国产成人一区二区在线| 看非洲黑人一级黄片| 91精品国产国语对白视频| 精品亚洲乱码少妇综合久久| 亚洲av日韩在线播放| 日韩av免费高清视频| 丝袜脚勾引网站| 午夜免费男女啪啪视频观看| 观看av在线不卡| 女人十人毛片免费观看3o分钟| 欧美极品一区二区三区四区| 99视频精品全部免费 在线| 中国三级夫妇交换| 人妻少妇偷人精品九色| 亚洲,欧美,日韩| 美女主播在线视频| 丰满少妇做爰视频| 久久国产精品男人的天堂亚洲 | 国产精品国产三级国产专区5o| 国产成人一区二区在线| 国产伦精品一区二区三区视频9| 久久午夜福利片| 国产精品麻豆人妻色哟哟久久| 亚洲av日韩在线播放| 成人特级av手机在线观看| 麻豆精品久久久久久蜜桃| 日韩国内少妇激情av| 26uuu在线亚洲综合色| 美女内射精品一级片tv| 联通29元200g的流量卡| 国产极品天堂在线| 大片免费播放器 马上看| 久久久久久九九精品二区国产| 久久久久视频综合| 国产美女午夜福利| 亚洲国产日韩一区二区| 十分钟在线观看高清视频www | 最近2019中文字幕mv第一页| 看免费成人av毛片| 国产一区亚洲一区在线观看| 国产精品一区www在线观看| 人人妻人人看人人澡| 黄片无遮挡物在线观看| 亚洲国产最新在线播放| 免费人成在线观看视频色| 亚洲av国产av综合av卡| 99久久综合免费| 亚洲人成网站高清观看| 久久综合国产亚洲精品| 人妻系列 视频| 日本黄色日本黄色录像| 亚洲伊人久久精品综合| 狂野欧美白嫩少妇大欣赏| 国产毛片在线视频| 亚洲综合色惰| 免费黄色在线免费观看| 精品酒店卫生间| 久久久久久伊人网av| 中文乱码字字幕精品一区二区三区| 亚洲人成网站在线观看播放| 欧美日韩精品成人综合77777| 18禁在线无遮挡免费观看视频| 亚洲人成网站高清观看| 欧美精品一区二区免费开放| 午夜福利影视在线免费观看| 欧美高清成人免费视频www| av免费观看日本| 日韩欧美精品免费久久| 久久久精品免费免费高清| 街头女战士在线观看网站| 久久 成人 亚洲| 又黄又爽又刺激的免费视频.| av网站免费在线观看视频| 偷拍熟女少妇极品色| 久久精品国产a三级三级三级| 亚洲欧美精品自产自拍| 视频中文字幕在线观看| 亚洲国产毛片av蜜桃av| 亚洲精品久久久久久婷婷小说| 亚洲色图综合在线观看| 精品一品国产午夜福利视频| 欧美日韩视频精品一区| 成人18禁高潮啪啪吃奶动态图 | 精品人妻一区二区三区麻豆| 国产一区二区三区综合在线观看 | 大片电影免费在线观看免费| 性色avwww在线观看| 亚洲精品国产av蜜桃| 日韩不卡一区二区三区视频在线| 丰满迷人的少妇在线观看| 精品人妻熟女av久视频| 美女中出高潮动态图| 久久人妻熟女aⅴ| 五月玫瑰六月丁香| 婷婷色av中文字幕| 寂寞人妻少妇视频99o| 成人漫画全彩无遮挡| 老熟女久久久| 自拍偷自拍亚洲精品老妇| 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 亚洲美女视频黄频| 国产精品国产三级国产专区5o| 亚洲成色77777| 日本猛色少妇xxxxx猛交久久| 一级av片app| 成年美女黄网站色视频大全免费 | 各种免费的搞黄视频| 99热全是精品| h日本视频在线播放| 国产精品三级大全| 成年av动漫网址| 国产高潮美女av| 久久6这里有精品| 大又大粗又爽又黄少妇毛片口| 精品人妻视频免费看| 插阴视频在线观看视频| 欧美激情国产日韩精品一区| 777米奇影视久久| 亚洲精品中文字幕在线视频 | 亚洲精品乱码久久久久久按摩| 观看美女的网站| 国产 一区 欧美 日韩| 18禁在线无遮挡免费观看视频| 午夜免费观看性视频| 99热国产这里只有精品6| 精品一区二区三区视频在线| 日韩国内少妇激情av| 国产成人精品婷婷| 伊人久久国产一区二区| 欧美日韩在线观看h| 免费观看的影片在线观看| 国产成人一区二区在线| 麻豆国产97在线/欧美| 亚洲国产成人一精品久久久| 亚洲精品视频女| 国产精品久久久久久av不卡| 欧美人与善性xxx| 爱豆传媒免费全集在线观看| 免费看日本二区| 两个人的视频大全免费| 欧美xxxx黑人xx丫x性爽| 另类亚洲欧美激情| 久久久成人免费电影| 麻豆乱淫一区二区| 在线天堂最新版资源| 搡老乐熟女国产| 老师上课跳d突然被开到最大视频| 国产极品天堂在线| 国产综合精华液| 纯流量卡能插随身wifi吗| 男女免费视频国产| 黑人高潮一二区| 秋霞在线观看毛片| 人妻夜夜爽99麻豆av| 亚洲,一卡二卡三卡| 成人午夜精彩视频在线观看| 一本—道久久a久久精品蜜桃钙片| 日本黄色片子视频| 欧美精品人与动牲交sv欧美| 国产精品久久久久成人av| 99re6热这里在线精品视频| 国产男人的电影天堂91| 国产在线一区二区三区精| av在线观看视频网站免费| 中文字幕精品免费在线观看视频 | 能在线免费看毛片的网站| 中文资源天堂在线| 精品久久国产蜜桃| 国产成人a区在线观看| 26uuu在线亚洲综合色| 免费黄色在线免费观看| 国产免费一级a男人的天堂| 卡戴珊不雅视频在线播放| 欧美变态另类bdsm刘玥| 九色成人免费人妻av| 高清毛片免费看| 免费人妻精品一区二区三区视频| 欧美精品亚洲一区二区| 大话2 男鬼变身卡| 日韩免费高清中文字幕av| 国产成人freesex在线| 九色成人免费人妻av| 日日啪夜夜爽| 亚洲美女视频黄频| 午夜精品国产一区二区电影| 中文欧美无线码| 久久久久网色| 久久人妻熟女aⅴ| 在线观看免费高清a一片| 欧美高清成人免费视频www| 亚洲国产欧美在线一区| 乱码一卡2卡4卡精品| 观看av在线不卡| 亚洲精品成人av观看孕妇| 午夜福利在线观看免费完整高清在| 久久99热这里只频精品6学生| 欧美变态另类bdsm刘玥| 多毛熟女@视频| videos熟女内射| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| 我的老师免费观看完整版| 国产免费一区二区三区四区乱码| 大香蕉久久网| 18禁在线播放成人免费| 中国美白少妇内射xxxbb| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 亚洲成人中文字幕在线播放| 秋霞在线观看毛片| 婷婷色麻豆天堂久久| 男人爽女人下面视频在线观看| 久久国产精品大桥未久av | 亚洲av不卡在线观看| 亚洲av中文av极速乱| 国内精品宾馆在线| 熟女电影av网| 夜夜骑夜夜射夜夜干| 看十八女毛片水多多多| 亚洲aⅴ乱码一区二区在线播放| 国产毛片在线视频| 免费大片18禁| 人人妻人人爽人人添夜夜欢视频 | av国产精品久久久久影院| 熟女av电影| 欧美另类一区| 久久这里有精品视频免费| 一本色道久久久久久精品综合| 99热这里只有精品一区| 最新中文字幕久久久久| 简卡轻食公司| 精品少妇黑人巨大在线播放| 国产精品成人在线| 久久人人爽av亚洲精品天堂 | 我要看日韩黄色一级片| 人妻 亚洲 视频| 久久久久视频综合| 亚洲第一av免费看| 大码成人一级视频| 国产免费又黄又爽又色| 国产精品人妻久久久影院| 日日啪夜夜爽| av免费在线看不卡| 国产日韩欧美亚洲二区| 国产成人免费观看mmmm| 乱码一卡2卡4卡精品| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区三区在线 | 天天躁夜夜躁狠狠久久av| 久久久久网色| 亚洲国产欧美一区二区综合| 水蜜桃什么品种好| 1024视频免费在线观看| 日韩一区二区三区影片| 两性夫妻黄色片| 亚洲熟女毛片儿| 日韩av在线免费看完整版不卡| 亚洲欧洲日产国产| avwww免费| 久久久国产一区二区| 亚洲第一av免费看| 美女扒开内裤让男人捅视频| h视频一区二区三区| 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 日韩中文字幕视频在线看片| 久久狼人影院| 久久久精品94久久精品| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 热99久久久久精品小说推荐| 国产片内射在线| 国产午夜精品一二区理论片| 国产视频首页在线观看| 欧美激情 高清一区二区三区| 成人影院久久| 欧美人与善性xxx| 久久热在线av| 91精品三级在线观看| 天堂8中文在线网| 亚洲国产精品999| 精品第一国产精品| 91精品三级在线观看| 校园人妻丝袜中文字幕| 女性生殖器流出的白浆| 精品久久久精品久久久| 国产成人精品久久久久久| 国产精品久久久久成人av| 中文字幕人妻丝袜制服| 在线 av 中文字幕| 日韩人妻精品一区2区三区| 精品福利观看| 十八禁人妻一区二区| 国产亚洲av高清不卡| 又大又爽又粗| 国产日韩欧美亚洲二区| 777米奇影视久久| 男人添女人高潮全过程视频| 丰满人妻熟妇乱又伦精品不卡| 欧美变态另类bdsm刘玥| 久久热在线av| 黄色一级大片看看| 精品卡一卡二卡四卡免费| 成在线人永久免费视频| 国产精品久久久久成人av| 大话2 男鬼变身卡| 亚洲欧美中文字幕日韩二区| 1024视频免费在线观看| 免费看av在线观看网站| 日韩av不卡免费在线播放| 国产成人免费无遮挡视频| 亚洲久久久国产精品| 国产精品欧美亚洲77777| 成人国产一区最新在线观看 | 亚洲第一av免费看| 丰满人妻熟妇乱又伦精品不卡| 久久天堂一区二区三区四区| 午夜免费鲁丝| 免费高清在线观看日韩| 老司机影院成人| 美女大奶头黄色视频| 狂野欧美激情性bbbbbb| 如日韩欧美国产精品一区二区三区| 999精品在线视频| 美女主播在线视频| 精品人妻一区二区三区麻豆| 一级片'在线观看视频| 考比视频在线观看| 精品一区二区三区av网在线观看 | 欧美亚洲日本最大视频资源| 日本wwww免费看| 久久久国产欧美日韩av| 男女边吃奶边做爰视频| 亚洲成av片中文字幕在线观看| 日日爽夜夜爽网站| 久久人妻福利社区极品人妻图片 | 成人手机av| 国产一级毛片在线| 亚洲国产av新网站| 日韩 欧美 亚洲 中文字幕| 高清av免费在线| 蜜桃在线观看..| 色视频在线一区二区三区| 看免费成人av毛片| 一本大道久久a久久精品| 亚洲黑人精品在线| 免费看av在线观看网站| 日韩一卡2卡3卡4卡2021年| 欧美精品高潮呻吟av久久| 日本猛色少妇xxxxx猛交久久| 亚洲精品一二三| 天堂8中文在线网| 国产亚洲欧美在线一区二区| av电影中文网址| 欧美精品一区二区大全| 亚洲中文字幕日韩| 18禁黄网站禁片午夜丰满| 国产亚洲精品第一综合不卡| 欧美中文综合在线视频| 精品少妇黑人巨大在线播放| 天天躁日日躁夜夜躁夜夜| 国产成人精品久久二区二区91| 亚洲欧美日韩另类电影网站| 亚洲熟女精品中文字幕| 国产亚洲精品第一综合不卡| 欧美少妇被猛烈插入视频| 亚洲久久久国产精品| 丝袜美足系列| 交换朋友夫妻互换小说| 丰满饥渴人妻一区二区三| 精品熟女少妇八av免费久了| 涩涩av久久男人的天堂| 男女边吃奶边做爰视频| 日韩,欧美,国产一区二区三区| 国产成人精品久久二区二区免费| 久久99一区二区三区| 国产亚洲av高清不卡| 别揉我奶头~嗯~啊~动态视频 | 国产av精品麻豆| 久久久久网色| 国产一区二区三区av在线| 亚洲国产精品国产精品| 国产激情久久老熟女| 午夜福利,免费看| 亚洲中文日韩欧美视频| 最近最新中文字幕大全免费视频 | 热99久久久久精品小说推荐| 国产一区二区 视频在线| 女性生殖器流出的白浆| 成年人午夜在线观看视频| 超色免费av| 韩国精品一区二区三区| 国产精品免费大片| 日本黄色日本黄色录像| 高清欧美精品videossex| 1024视频免费在线观看| 男人舔女人的私密视频| 久久久国产一区二区| 深夜精品福利| 午夜福利影视在线免费观看| 丝袜人妻中文字幕| 老熟女久久久| svipshipincom国产片| 一二三四社区在线视频社区8| 好男人视频免费观看在线| 精品福利永久在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲成国产人片在线观看| 午夜免费男女啪啪视频观看| 少妇猛男粗大的猛烈进出视频| 天堂中文最新版在线下载| 亚洲成人免费电影在线观看 | 久久亚洲精品不卡| av国产精品久久久久影院| 精品国产乱码久久久久久男人| 国产黄色视频一区二区在线观看| 最黄视频免费看| 久久国产精品男人的天堂亚洲| 国产精品一区二区免费欧美 | 最新在线观看一区二区三区 | 可以免费在线观看a视频的电影网站| 午夜福利一区二区在线看| 亚洲,一卡二卡三卡| 大话2 男鬼变身卡| 亚洲五月色婷婷综合| 日韩av不卡免费在线播放| 国产精品免费视频内射| 欧美xxⅹ黑人| 国产免费一区二区三区四区乱码| 美女高潮到喷水免费观看| 超碰成人久久| 亚洲色图 男人天堂 中文字幕| 一个人免费看片子| 91老司机精品| 男女之事视频高清在线观看 | 精品免费久久久久久久清纯 | 老司机午夜十八禁免费视频| 国产精品国产av在线观看| 国产1区2区3区精品| 天天操日日干夜夜撸| 久久久久久人人人人人| 美女扒开内裤让男人捅视频| 大香蕉久久网| 日韩 亚洲 欧美在线| 女警被强在线播放| 成人午夜精彩视频在线观看| 大片免费播放器 马上看| 天天躁夜夜躁狠狠久久av| 好男人电影高清在线观看| 美女扒开内裤让男人捅视频| 五月开心婷婷网| 国产精品一区二区精品视频观看| 又大又黄又爽视频免费| 国产av一区二区精品久久| 美女福利国产在线| 1024视频免费在线观看| 亚洲少妇的诱惑av| 两个人看的免费小视频| 性色av乱码一区二区三区2| 国产精品一二三区在线看| 成人国语在线视频| 午夜福利,免费看| 免费在线观看视频国产中文字幕亚洲 | 免费看十八禁软件| 十八禁高潮呻吟视频| 久久影院123| 欧美激情高清一区二区三区| 国产精品久久久久久精品古装| 亚洲国产欧美在线一区| 九草在线视频观看| 国产黄色视频一区二区在线观看| 亚洲伊人久久精品综合| 免费在线观看完整版高清| 免费在线观看视频国产中文字幕亚洲 | 视频区图区小说| 精品人妻在线不人妻| 王馨瑶露胸无遮挡在线观看| 十八禁人妻一区二区| 免费黄频网站在线观看国产| 青青草视频在线视频观看| 成年美女黄网站色视频大全免费| 亚洲国产欧美日韩在线播放| 国产激情久久老熟女| 亚洲成色77777| 亚洲熟女毛片儿| 成年女人毛片免费观看观看9 | 男女高潮啪啪啪动态图| 最近中文字幕2019免费版| 妹子高潮喷水视频| 日韩av免费高清视频| 久久女婷五月综合色啪小说| 国产精品国产三级国产专区5o| 各种免费的搞黄视频| 成年人免费黄色播放视频| 国产精品一区二区精品视频观看| 免费不卡黄色视频| 亚洲男人天堂网一区| 亚洲精品日本国产第一区| 女警被强在线播放| 亚洲精品一二三| 国产精品熟女久久久久浪| 国产极品粉嫩免费观看在线| 蜜桃在线观看..| 日本欧美国产在线视频| 日本一区二区免费在线视频| 蜜桃在线观看..| 亚洲人成电影免费在线| 国产精品国产av在线观看| 高清欧美精品videossex| 观看av在线不卡| 亚洲人成网站在线观看播放| 夜夜骑夜夜射夜夜干| 欧美日韩一级在线毛片| 精品人妻在线不人妻| 一本久久精品| 久9热在线精品视频| 欧美变态另类bdsm刘玥| 我要看黄色一级片免费的| 美女主播在线视频| 亚洲视频免费观看视频| 老司机深夜福利视频在线观看 | 啦啦啦中文免费视频观看日本|