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

    Extraction process in Chuantieling (CTL) based on the quality by design (QbD)concept

    2022-07-02 07:46:26LILIZiruiGUOZhihu
    Digital Chinese Medicine 2022年2期

    LI Y,LI Zirui,GUO Zhihu

    a. College of Pharmacy,Hunan University of Chinese Medicine,Changsha,Hunan 410208,China

    b. College of Traditional Chinese Medicine,Hunan University of Chinese Medicine,Changsha,Hunan 410208,China

    ABSTRACT Objective This study aimed to design and optimize the water extraction process for Chuantieling (喘貼靈,CTL) based on the concept of quality by design (QbD).Methods The single-factor experiments were used to select the best experimental points of CTL water extraction. On this basis,the transfer rate of ephedrine hydrochloride and sinapine thiocyanate,and the yield of the extract were evaluated as the evaluation indicators.The liquid-solid ratio,extraction time,and pH value were selected as the main factors to carry out the Box-Behnken design (BBD). A mathematical model of the critical process parameters(CPPs) and critical quality attributes (CQAs) was established,the interaction between CQAs and CPPs was examined,and the CTL extraction process design space was established and optimized,which guaranteed the stability of the process. The high performance liquid chromatography (HPLC) was used to analyze the main active compounds in the CTL extract.Results Through single-factor experiments,the best experimental parameters were found to be a liquid-solid ratio of 8∶1,extraction time of 90 min,pH value of 4,and extraction times of three. The experimental variance analysis results of the BBD showed that the P value of the regression model was less than 0.05,and the lack-of-fit value was greater than 0.01,indicating that the constructed model had good predictive ability. The operating space of the CPPs of the CTL water extraction process was combined with the actual production situation. In the production situation,the best extraction process was nine times of water addition,72 min of extraction time,and 4.5 of pH value. The HPLC results showed that the peak areas of ephedrine hydrochloride,sinapine thiocyanate,tetrahydropalmatine,methyl eugenol,cinnamaldehyde,and 6-gingerol in the CTL extract accounted for 0.94%,14.32%,0.78%,31.23%,0.34%,and 0.44% of the total peak area,respectively.Conclusion The water extraction process design space of CTL based on QbD was conducive to actual production operations,ensuring the stability of the process.

    Keywords Chuantieling (CTL) Quality by design (QbD) Box-Behnken design (BBD) Critical quality attributes (CQAs) Critical process parameters (CPPs) Mahuang (Ephedrae Herba) Jiezi (Sinapis Semen) Ephedrine hydrochloride Sinapine thiocyanate

    1 Introduction

    Bronchial asthma is a chronic inflammatory disease of the airway involving multiple cells (such as eosinophils,mast cells,T lymphocytes,neutrophils,and airway epithelial cells) and cellular components[1]. Its clinical manifestations include wheezing,shortness of breath,tightness of chest,coughing,and other respiratory symptoms accompanied with reversible expiratory airflow limitation. The economic burden caused by the treatment of bronchial asthma exceeds the total economic burden caused by tuberculosis and human immunodeficiency virus (HIV) infection[2,3]. It causes great physical and psychological pain to patients,and affects daily life. Therefore,asthma is considered as a serious health problem and has attracted global attention[4]. Traditional Chinese medicine (TCM) has the unique advantage of eliminating the cause,controlling the symptoms,reducing adverse drug reactions,and also reducing the recurrence of asthma[5]. Acupoint application is a characteristic therapy in TCM. It is a complex intervention method that is affected by factors such as drugs,acupoints,time,and stimulation methods. By stimulating local acupoints,it stimulates the body's internal meridian energy,enhances immune function,and reduces asthma recurrence[6-8].Chuantieling (喘貼靈,CTL) is a hospital preparation of the Second Affiliated Hospital of Hunan University of Chinese Medicine. It is composed of eight herbs:Mahuang (Ephedrae Herba),Jiezi (Sinapis Semen),Xixin(Asari Radix et Rhizoma),Dingxiang (Caryophylli Flos),Rougui (Cinnamomi Cortex),Yanhusuo (Corydalis Rhizoma),Gansui (Kansui Radix),and Shengjiang (Zingiberis Rhizoma Recens). It has been used clinically for many years,and clinical research have shown that CTL plaster exerts a substantial effect on bronchial asthma[9,10].However,there are still problems in the current CTL extraction process such as unclear effective substance content,lack of reasonable extraction process parameters,and weak quality standards. In the prescription,Mahuang (Ephedrae Herba) and Jiezi (Sinapis Semen)are the principal drugs. Ephedrine hydrochloride and sinapine thiocyanate,the main compounds of Mahuang(Ephedrae Herba) and Jiezi (Sinapis Semen),exert definite and strong curative effects on bronchial asthma[9,10].Therefore,in the optimization of the water extraction process,the ephedrine hydrochloride and sinapine thiocyanate content,and extract yield have been proposed to be used as evaluation indicators. To ensure the quality and clinical efficacy of the patch,a single-factor experiment combined with Response Surface Methodology (RSM) was used to optimize the CTL extraction process that can maximize the extraction and reduce the loss of active compounds (ephedrine hydrochloride and sinapine thiocyanate)[11].

    Quality by design (QbD) is a scientific,risk-based,comprehensive,and systematic research method for improving product quality by emphasizing the design. The QbD system is composed of a quality target product profile (QTPP),critical material attribute (CMA),critical process parameters (CPPs),and critical quality attributes(CQAs). Through the establishment of mathematical models to characterize the relationships among CMA,CPPs,and CQAs with mutual influence,the corresponding design space and control strategy were adopted to ensure the stability of the QTPP (Figure 1). Recently,QbD has promoted the transformation of the pharmaceutical production model and received widespread attention[12].It is now believed that quality is not only purely dependent on inspection and production processes,but is also conferred by design[13]. QbD has been widely used in the research and development of chemical and biological medicine production technology[14]; it has also been widely utilized in the research of TCM production technology[15,16]. This study focused on QbD by applying experimental design,design space,and other methods to study the extraction process of CTL,reduce the risk of preparation production and post-marketing,and provide an experimental basis for the late-stage research of CTL and development of new drugs.

    Figure 1 Application of QbD in products

    2 Materials

    2.1 Instruments and equipments

    A 1200 high-performance liquid chromatography (HPLC)column (Agilent Technologies) with a Hibar C18bonded silica column (150 mm × 4.6 mm,5 μm); an Agilent ChemStation chromatography workstation (Agilent Technologies); a DK-98-11A electric heating constant temperature water bath (Tianjin Test Instrument Co.,Ltd.),a CP114 electronic balance (Ohaus Instrument Co.,Ltd.),and a KQ-500DE CNC ultrasonic cleaner (Kunshan Ultrasonic Instrument Co.,Ltd.) were used in this study.

    2.2 Reagents and drugs

    Methanol and acetonitrile (chromatographically pure,American Tiandi Company,USA),phosphoric acid (analytically pure,Sinopharm Chemical Reagent Co.,Ltd.,China),and ultrapure water (Option RT ultra AN ultrapure water system,Gemany) were used in this study.Cinnamaldehyde,6-gingerol,tetrahydropalmatine,sinapine thiocyanate,and methyl eugenol (purity ≥ 98.0%)were purchased from Chengdu Aifa Biological Technology Co.,Ltd. (China). Ephedrine hydrochloride (purity ≥98.0%) was purchased from the Chinese Medicine Food and Drug Control Institute (China). The ingredients Mahuang (Ephedrae Herba),Jiezi (Sinapis Semen),Xixin(Asari Radix et Rhizoma),Dingxiang (Caryophylli Flos),Rougui (Cinnamomi Cortex),Yanhusuo (Corydalis Rhizoma),Gansui (Kansui Radix),and Shengjiang (Zingiberis Rhizoma Recens) were all purchased from Hunan Zhenxing Chinese Medicine Co.,Ltd. (China),in compliance with the requirements of an appendix ofChinese Pharmacopoeia(2020 edition).

    3 Methods and results

    3.1 QTPP and CQAs establishment of CTL extraction process

    The appropriate dosage form design,route of administration,and packaging form of preparation were determined by investigating relevant research literature and combining it with previous experiments[17,18]. The QTPP and CQAs of the CTL extract were established,as shown in Table 1 and 2.

    Table 1 QTPP of CTL extract

    Table 2 CQAs identification of CTL extract

    3.2 CTL test solution preparation and composition identification

    3.2.1 CTL test solution preparation The prescribed amount of CTL was weighed including 30 g Jiezi (Sinapis Semen),20 g Mahuang (Ephedrae Herba),6 g Xixin (Asari Radix Et Rhizoma),20 g Dingxiang (Caryophylli Flos),20 g Rougui (Cinnamomi Cortex),20 g Yanhusuo (Corydalis Rhizoma),20 g Gansui (Kansui Radix),and 20 g Shengjiang (Zingiberis Rhizoma Recens). The liquidsolid ratio was 10∶1,the pH value was adjusted to 4,and the extraction time was 120 min. The extraction was performed in triplicate,filtered,the filtrate combined,concentrated,filtered through a microporous membrane,and the filtrate was ready for use.

    3.2.2 Reference solution preparation In total,10.8 mg ephedrine hydrochloride,10.0 mg sinapine thiocyanate,11.8 mg tetrahydropalmatine,8.8 mg cinnamaldehyde,9.1 mg 6-gingerol,and 10.0 mg methyl eugenol reference substance were precisely weighed and dissolved in an appropriate amount of methanol. The volume was made up to 10 mL in a volumetric flask with methanol,and 1 mL was transferred to 25 mL volumetric flasks and mixed evenly. A mixed reference solution was prepared.The mass concentrations of ephedrine hydrochloride,sinapine thiocyanate,tetrahydropalmatine,cinnamaldehyde,6-gingerol,and methyl eugenol were 43.2,40.0,47.2,35.2,36.4,and 40.0 μg/mL,respectively.

    3.2.3 High Perfomance Liquid Chromatography (HPLC)conditions Hypersil BDS C18column (250 mm × 4.6 mm,5 μm); mobile phase: acetonitrile (A) - 0.1% phosphoricacid (B); elution gradient: 0 - 10 min,4% A; 10 - 20 min,4% - 12% A; 20 - 40 min,12% - 25% A; 40 - 60 min,25% - 35% A; 60 - 80 min,35% - 50% A; 80 - 90 min,50% - 90% A; 90 - 100 min,90% - 4% A; flow rate:1.0 mL/min; column temperature: 30 °C; detection wavelength: 210 nm,310 nm; injection volume: 10 μL.

    3.2.4 Linear relationship investigation The appropriate amounts of ephedrine hydrochloride and sinapine thiocyanate reference samples were taken and accurately weighed. Methanol was used to prepare mixed solutions containing 0.26 mg/mL ephedrine hydrochloride and 0.38 mg/mL sinapine thiocyanate. 1,1.5,2,3,and 5 mL of the above-mentioned reference solution were placed in a 10 mL measuring flask,diluted to the mark with methanol,shaken well,and analyzed. Linear regression was performed on the mass concentration of ephedrine hydrochloride and sinapine thiocyanate by the peak area. The regression equations of ephedrine hydrochloride and sinapine thiocyanate were Y = 20 895X - 28.348 (R = 0.999 6)and Y = 19 454X - 226.960 (R = 0.999 7),rspectively. The linear relationship between the acid salt concentration was 0.038 - 0.190 mg/mL,and the linear relationship between ephedrine hydrochloride was 0.026 - 0.130 mg/mL.

    3.2.5 Precision inspection The same test solution of CTL extract was injected six times for analysis (conditions according to section 3.2.3) and the chromatograms were recorded. Using ephedrine hydrochloride as an example,the Relative Standard Deviation (RSD) value of the relative retention time of the six main chromatographic peaks was 1.02%,the RSD value of the relative peak area was 1.24%,and the similarity value of each chromatogram was 0.998,indicating that the precision of HPLC was high.

    3.2.6 Stability investigation The same test solution of the CTL extract was used,and samples were injected at 0,2,4,6,8,12,and 24 h (conditions according to section 3.2.3) for analysis. The chromatograms were recorded at 210 and 310 nm. Using ephedrine hydrochloride as an example,the RSD value of the relative retention time of the six main chromatographic peaks was 0.73%,the RSD value of the relative peak area was 2.51%,and the similarity value of each chromatogram was 0.998,indicating that the test solution remained stable within 24 h,and the stability was good.

    3.2.7 Repeatability investigation The precision inspection test (according to section 3.2.5) was repeated to determine repeatability. Using ephedrine hydrochloride as an example,the RSD value of the relative retention time of the six main chromatographic peaks was 0.93%,the RSD value of the relative peak area was 0.63%,and the similarity value of each chromatogram was 0.998,indicating that the repeatability of this method was good. Relative peak area was 2.51%,and the similarity value of each chromatogram was 0.998,indicating that the test solution remained stable within 24 h,and the stability was good.

    3.2.8 Transfer rate calculations for ephedrine hydrochloride and sinapine thiocyanate 10 μL each of the reference substance solution was prepared (preparation according to section 3.2.2) and measured under the chromatographic conditions (according to section 3.2.3). The peak areas of ephedrine hydrochloride and sinapine thiocyanate were recorded and the ephedrine hydrochloride content (M0) and sinapine thiocyanate content(N0) in ephedra decoction pieces,ephedrine hydrochloride content (M) and sinapine thiocyanate content (N) in CTL extract were calculated. The transfer rate of ephedrine hydrochloride in CTL extract was calculated using formula (1),and the transfer rate of sinapine thiocyanate in CTL extract was calculated using formula (2).

    3.2.9 HPLC analysis of the main compounds of CTL extract The test solution (according to section 3.2.1) and the mixed control solution (according to section 3.2.2)were injected under the chromatographic conditions (according to section 3.2.3),and the chromatogram was recorded. Ephedrine hydrochloride and sinapine thiocyanate were the active compounds in the principle drugs Mahuang (Ephedrae Herba) and Jiezi (Sinapis Semen) in the prescription of CTL; tetrahydropalmatine was the sedative compound in Yanhusuo (Corydalis Rhizoma);methyl eugenol was the active compound in Xixin (Asari Radix et Rhizoma) and Dingxiang (Caryophylli Flos); cinnamaldehyde was the active compound in Rougui (Cinnamomi Cortex),and 6-gingerol was the active compound in Shengjiang (Zingiberis Rhizoma Recens). Gansui (Kansui Radix) was used as an envoy drug in the prescription,and the main active compound had no evident effect on the treatment of bronchial asthma. Therefore,we selected these six active compounds as detection indicators.

    As shown in Figure 2,peaks one to six in the chromatogram of the mixed control solution were ephedrine hydrochloride,sinapine thiocyanate,tetrahydropalmatine,methyl eugenol,cinnamaldehyde,and 6-gingerol. Compared with that seen in the reference chromatogram,it could be seen that the CTL extract contained ephedrine hydrochloride,sinapine thiocyanate,tetrahydropalmatine,methyl eugenol,cinnamaldehyde,and 6-gingerol.The peak area of these components accounted for 0.94%,14.32%,0.78%,31.23%,0.34%,and 0.44% of the total peak area of the CTL extract.

    Figure 2 HPLC chromatogram of mixed reference substance and CTL extract

    3.3 Determination of extract yield

    25.00 mL of the CTL water extract preparation refered to section 3.2.1 was precisely pipetted and placed in an evaporating dish that had dried to a constant weight,evaporated to dryness in a water bath,dried at 105 °C for 3 h,and transferred to a desiccator after cooling for 30 min,weighed,and calculated according to formula (3).

    3.4 Single-factor investigation

    3.4.1 Liquid-solid ratio Five doses of CTL were weighed for the single-factor experiment of the liquid-solid ratios of 6∶1,8∶1,10∶1,12∶1,and 14∶1,pH value was adjusted to 4,the extraction time was 120 min,extraction was done three times,and then filtered. The filtrate was combined,500 mL of which was concentrated to approximately 100 mL,the content of ephedrine hydrochloride and mustard thiocyanate,and their extraction rates were determined. The comprehensive score was calculated as(ephedrine hydrochloride transfer rate × 0.4)/maximum ephedrine hydrochloride transfer rate + (sinapine thiocyanate transfer rate × 0.4)/maximum sinapine thiocyanate transfer rate + (extract yield × 0.2)/maximum extract transfer rate.

    When the liquid-solid ratio was approximately 8∶1,the comprehensive score peaked (Figure 3A). When the liquid-solid ratio was greater than 8∶1,a substantial increase in the liquid-solid ratio could exert a certain impact on the pH value of the extract. This would change the content of ephedrine hydrochloride,sinapine thiocyanate,and other substances in CTL that were pH-sensitive; therefore,8∶1 was selected as the best liquid-solid ratio.

    3.4.2 Extraction time Five doses of CTL were weighed for the single-factor experiment of extraction time. The extraction duration was 60,90,120,150,and 180 min,the liquid-solid ratio was 10∶1,pH value was adjusted to 4,extracted three times,and filtered. The filtrate was combined,500 mL of which was concentrated to approximately 100 mL,and the content of ephedrine hydrochloride and mustard thiocyanate and the extraction rate were determined. The comprehensive score was calculated as(ephedrine hydrochloride transfer rate × 0.4)/maximum ephedrine hydrochloride transfer rate + (sinapine thiocyanate transfer rate × 0.4)/maximum sinapine thiocyanate transfer rate + (extract yield × 0.2)/maximum extract transfer rate.

    The comprehensive score peaked when the extraction time was approximately 90 min (Figure 3B). Because ephedrine hydrochloride had a small molecular weight and a certain volatility,and sinapine thiocyanate would be degraded with the extension of the extraction time,resulting in a downward trend in the comprehensive score of extraction time greater than 90 min,90 min was selected as the best extraction time.

    3.4.3 pH value Five doses of CTL were weighed for the single-factor experiment of pH value. Hydrochloric acid was added to adjust the pH value to 1,2,3,4 and 5,the ratio of liquid-solid was 8∶1,extracted three times and filtered. The filtrate was combined,500 mL of which was concentrated to approximately 100 mL,and the content of ephedrine hydrochloride and mustard thiocyanate and the extraction rate were determined. The comprehensive score was calculated as (ephedrine hydrochloride transfer rate × 0.4)/maximum ephedrine hydrochloride transfer rate + (sinapine thiocyanate transfer rate ×0.4)/maximum sinapine thiocyanate transfer rate + (extract yield × 0.2)/maximum extract transfer rate.

    When the pH value was approximately 4,the comprehensive score reached its peak (Figure 3C).

    3.4.4 Times of extraction screening Under the conditions of a extraction time of 90 min,a liquid-solid ratio of 8∶1,and a pH value of 4,the effects of different of decoction extraction times (1,2,3,4,and 5) on the comprehensive score were investigated.

    The results showed that the extraction rate increased with the increased extraction times,but when the extraction times was greater than three,the comprehensive score decreased (Figure 3D). Excessive extraction times would increase the cost and affect extraction efficiency.Therefore,three times of extration were selected for the follow-up experiments.

    Figure 3 Results of single-factor experiments

    3.5 Response surface optimization of extraction process

    3.5.1 Box-Behnken Design (BBD) and RSM According to the single-factor experiment,the best experimental point for water extraction of CTL ingredients was selected. Taking the ephedrine hydrochloride and sinapine thiocyanate content,and the comprehensive score of the extract yield as the response value,a mathematical model was established to optimize the prescription. The factor levels are presented in Table 3.

    Table 3 Factor levels of BBD

    3.5.2 Model fitting and analysis of variance The selected variables were studied at three different levels,and experiments were performed according to BBD (Table 4).

    Table 4 Box-Behnken design and results

    3.5.3 BBD and RSM The quadratic polynomial regression equation referring to the three independent variables was as follows: Y = - 117.19 + 29.00A + 201.90B -7.66C + 0.33AB - 1.01AC - 7.36BC - 1.10A2- 42.04B2+1.44C2,R = 0.945 1,indicating that the model fitted well.Table 5 shows that the model was statistically significant(P< 0.05),and described the experimental results well.ThePvalue of the model's lack-of-fit was 0.6987,indicating that there was no significant difference in the lack-offit term. The probability that the actual value did not fit the predicted value was low,the model was significant,and the fitting equation had a good predictability for the experimental results. It can be seen that among the three factors investigated in the extraction process,the liquidsolid ratio (A) was a very significant item,the extraction time (B) and pH value (C) were significant,and B2(P =0.006 9) indicated the relationship of the factor and the response surface value was not a simple linear relationship; the BC interaction was significant,indicating that the interaction between the extraction time (B) and pH value (C) was evident. The response surface three dimensions (3D) diagram is shown in Figure 4.

    3.6 Establishment and verification of design space

    A single-factor experiment was used to select extraction times,and liquid-solid ratio was to be as the main investigation factor. BBD was then used to establish the mathematical models of CQAs and CPPs,and the MODDE software (Umetrics) self-edited program was used for calculations. Based on the literature[19],a comprehensive score of the yield of the CTL extract and the transfer rate of active compounds through water extraction greater than 80 was ideal,and the comprehensive score optimization standard was set to 80.

    The design space for the CTL water extraction process was obtained. The operating space range was8.9 -10.0 times the water addition,extraction time was 72.0 - 103.8 min,and the pH value was 4.47 - 5.00(Figure 5).

    Table 5 Response surface analysis of variance and significance test

    Figure 4 Contour plots and 3D response surface plots of the comprehensive evaluation

    Figure 5 Design spatial result

    The verification results (Table 6) showed that the design space operation of the process parameters could ensure the stable quality of the CTL water extraction process and effectively reduce the variability in product quality. The design space was irregular. To facilitate the operation and subsequent industrialization,the recommended operating space range was nine times the water addition,the extraction time was 72 min,and the pH value was 4.5.

    Table 6 Verification points and results (±s,n = 3)

    Table 6 Verification points and results (±s,n = 3)

    P1,P2,and P3 represent three random experimental points in the design space.

    Experimental point Liquidsolid ratio Extraction time(min)pH Comprehensive score P1 9.0 80 4.6 93 P2 9.6 94 5.0 85 P3 10.0 100 4.5 87

    4 Discussion

    Mahuang (Ephedrae Herba),Jiezi (Sinapis Semen),Xixin(Asari Radix et Rhizoma),Dingxiang (Caryophylli Flos),Rougui (Cinnamomi Cortex),and Shengjiang (Zingiberis Rhizoma Recens) in the prescription of CTL are rich in volatile oil components that can be extracted by water,and they exert certain curative effects on bronchial asthma[20]. The HPLC method established in this study was simple and effective,with accurate results and good precision. It can simultaneously determine the six effective components in CTL. With the improvement in people's awareness regarding drug safety,it is extremely important to strengthen the quality control of the extraction process of TCM ingredients.

    Based on the influencing factors,different pH values were selected as the investigating factors. When the pH value was 4.5 - 5.0,the acid promoted alkaloid leaching,improved alkaloid stability,and freed the organic acid to facilitate the removal of acid-insoluble impurities[21]. The results of this experiment showed that if the extraction time was too short or too long,the comprehensive score would be reduced,and both ephedrine hydrochloride and sinapine thiocyanate showed a downward trend,suggesting that high-temperature conditions might cause the destruction of active compounds,thereby affecting the test results. Therefore,in the preparation process of the CTL patch,it is particularly necessary to pay attention to the quality control in the extraction process and select the appropriate pH value,liquid-solid ratio,and extraction time to avoid reducing of the efficacy and the transfer of effective components due to unreasonable extraction processes.

    In this study,the extraction process of CTL as an example. By combining a single-factor test and BBD and introducing QbD into the experimental design method of CTL extraction,a quantitative mathematical model and process of CQAs and CPPs were constructed[22]. This method effectively improved the transfer rate and extract yield of ephedrine hydrochloride and sinapine thiocyanate,reduced the loss during the extraction process,and demonstrated a stable process and good reproducibility.The design space provided an experimental basis for future research on CTL and the development of new drugs.Combining the operating space obtained in the experiment with the actual production situation,the best extraction process was nine times of water addition,72 min of extraction time,and 4.5 of pH value.

    5 Conclusion

    Compared with the results obtained by using fixed process parameters,the establishment of a design space made process operations more flexible and effectively reduced the variability of product quality. The application of QbD in the pharmaceutical industry can effectively improve production efficiency and is of great importance to the development of the pharmaceutical field.

    Fundings

    Hunan University of Traditional Chinese Medicine Graduate Training Quality Engineering Project(2019CX57),and the First-class Discipline Project on Chinese Pharmacology of Hunan University of Chinese Medicine (201803).

    Competing interests

    The authors declare no conflict of interest.

    麻豆成人av在线观看| 国产亚洲精品久久久久5区| 亚洲精品美女久久久久99蜜臀| 蜜桃在线观看..| 无遮挡黄片免费观看| 婷婷成人精品国产| 又紧又爽又黄一区二区| 国产精品1区2区在线观看. | 国产精品一区二区精品视频观看| 久久精品国产综合久久久| 午夜福利视频精品| 国产成人av激情在线播放| 精品国产乱子伦一区二区三区| 亚洲欧洲精品一区二区精品久久久| 自拍欧美九色日韩亚洲蝌蚪91| 91麻豆av在线| 欧美日韩视频精品一区| 亚洲精品中文字幕一二三四区 | 在线观看免费午夜福利视频| 交换朋友夫妻互换小说| 久久99一区二区三区| 亚洲一码二码三码区别大吗| 狠狠婷婷综合久久久久久88av| 精品国产一区二区三区久久久樱花| 成人免费观看视频高清| 日韩三级视频一区二区三区| 国产一卡二卡三卡精品| 日韩欧美一区二区三区在线观看 | 欧美国产精品va在线观看不卡| 成年人午夜在线观看视频| 成年版毛片免费区| a级毛片黄视频| 亚洲五月色婷婷综合| 国产日韩欧美视频二区| 亚洲专区国产一区二区| 久久精品亚洲熟妇少妇任你| 成年人午夜在线观看视频| 麻豆成人av在线观看| 国产视频一区二区在线看| 飞空精品影院首页| 欧美成人午夜精品| 亚洲精品成人av观看孕妇| 国产精品亚洲一级av第二区| 国产精品久久久久久精品古装| 国产老妇伦熟女老妇高清| 在线观看免费视频日本深夜| 9热在线视频观看99| 亚洲一区中文字幕在线| 久热爱精品视频在线9| 国产视频一区二区在线看| 建设人人有责人人尽责人人享有的| 亚洲成人免费电影在线观看| 侵犯人妻中文字幕一二三四区| 亚洲成人免费电影在线观看| 欧美日韩亚洲高清精品| 天天躁夜夜躁狠狠躁躁| 男女午夜视频在线观看| 男女午夜视频在线观看| 啦啦啦免费观看视频1| av视频免费观看在线观看| 人人妻,人人澡人人爽秒播| 久久久久久久国产电影| 国产在线一区二区三区精| 久久久国产欧美日韩av| 日韩视频在线欧美| 一区二区日韩欧美中文字幕| 高清视频免费观看一区二区| 中文字幕人妻丝袜一区二区| 天堂动漫精品| 99国产综合亚洲精品| 国产免费福利视频在线观看| 天天操日日干夜夜撸| 亚洲精品自拍成人| 日本av免费视频播放| 免费观看a级毛片全部| 女人精品久久久久毛片| 国产成人一区二区三区免费视频网站| 成人国语在线视频| 大片电影免费在线观看免费| 两个人看的免费小视频| 人人妻人人添人人爽欧美一区卜| 黄色怎么调成土黄色| 一区二区三区国产精品乱码| 在线观看免费高清a一片| 不卡av一区二区三区| 热99国产精品久久久久久7| 在线观看66精品国产| 少妇 在线观看| 蜜桃国产av成人99| 亚洲精品在线观看二区| 搡老岳熟女国产| 一区二区日韩欧美中文字幕| 中文字幕人妻丝袜一区二区| 日本精品一区二区三区蜜桃| 男女床上黄色一级片免费看| 日韩熟女老妇一区二区性免费视频| 性色av乱码一区二区三区2| 99riav亚洲国产免费| 欧美日韩一级在线毛片| 男人操女人黄网站| 国产成人免费观看mmmm| 黑人猛操日本美女一级片| 国产精品 欧美亚洲| 精品国内亚洲2022精品成人 | 叶爱在线成人免费视频播放| 五月开心婷婷网| 99国产精品一区二区蜜桃av | 国产福利在线免费观看视频| 亚洲精品自拍成人| 十分钟在线观看高清视频www| 91大片在线观看| 搡老乐熟女国产| 免费一级毛片在线播放高清视频 | 午夜视频精品福利| 精品福利观看| 美女视频免费永久观看网站| 亚洲成人免费av在线播放| 久久久精品免费免费高清| 老汉色av国产亚洲站长工具| 狠狠狠狠99中文字幕| 91国产中文字幕| 无人区码免费观看不卡 | 国产高清激情床上av| 99国产精品免费福利视频| 肉色欧美久久久久久久蜜桃| 国产成人影院久久av| 亚洲精品乱久久久久久| 我的亚洲天堂| 精品久久久精品久久久| 曰老女人黄片| www.自偷自拍.com| 国产欧美日韩一区二区三| 捣出白浆h1v1| 国产av精品麻豆| 欧美精品一区二区免费开放| 欧美日韩亚洲国产一区二区在线观看 | 久久精品国产综合久久久| 免费日韩欧美在线观看| 欧美日韩亚洲高清精品| 日韩欧美一区二区三区在线观看 | 日韩欧美国产一区二区入口| 国产精品久久久人人做人人爽| 成在线人永久免费视频| 国产黄色免费在线视频| 国产一区二区在线观看av| 9色porny在线观看| 一区二区三区乱码不卡18| 在线观看www视频免费| 日韩视频在线欧美| videosex国产| 中国美女看黄片| 极品人妻少妇av视频| 欧美久久黑人一区二区| 我要看黄色一级片免费的| 国产精品久久久久久精品古装| 日韩 欧美 亚洲 中文字幕| 啦啦啦免费观看视频1| 十八禁网站网址无遮挡| 久久国产精品大桥未久av| 男女边摸边吃奶| 无人区码免费观看不卡 | 国产在线观看jvid| 在线永久观看黄色视频| 高清黄色对白视频在线免费看| bbb黄色大片| 色尼玛亚洲综合影院| 亚洲精品av麻豆狂野| 一边摸一边抽搐一进一出视频| 老汉色∧v一级毛片| 色在线成人网| 欧美激情极品国产一区二区三区| 久久久久久久久久久久大奶| 99香蕉大伊视频| 91麻豆av在线| 丝袜美足系列| 欧美精品亚洲一区二区| 丰满饥渴人妻一区二区三| 亚洲国产看品久久| 国产三级黄色录像| 亚洲av国产av综合av卡| 亚洲精华国产精华精| 69精品国产乱码久久久| 最近最新中文字幕大全电影3 | 国产不卡av网站在线观看| 国产一区二区三区在线臀色熟女 | 色老头精品视频在线观看| 国产精品久久久久久精品古装| 国产又色又爽无遮挡免费看| 成人亚洲精品一区在线观看| 久久午夜亚洲精品久久| 精品久久蜜臀av无| 午夜日韩欧美国产| 999久久久精品免费观看国产| 国产野战对白在线观看| 他把我摸到了高潮在线观看 | 高清在线国产一区| 夜夜骑夜夜射夜夜干| 丝袜在线中文字幕| 亚洲成人国产一区在线观看| 日韩视频一区二区在线观看| 日韩人妻精品一区2区三区| 狂野欧美激情性xxxx| 淫妇啪啪啪对白视频| 国产精品久久久久成人av| 免费看a级黄色片| 久久精品熟女亚洲av麻豆精品| 蜜桃在线观看..| 精品国产乱码久久久久久男人| 欧美精品一区二区免费开放| 夜夜夜夜夜久久久久| 久久热在线av| 国产精品久久久久久精品电影小说| 国产亚洲精品第一综合不卡| 国产欧美日韩精品亚洲av| 日韩欧美免费精品| 18禁黄网站禁片午夜丰满| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| avwww免费| 久久天堂一区二区三区四区| 国产一区二区激情短视频| tocl精华| 国产精品一区二区在线不卡| 一边摸一边抽搐一进一出视频| 国产日韩欧美在线精品| 91成年电影在线观看| 欧美亚洲 丝袜 人妻 在线| 免费av中文字幕在线| 无限看片的www在线观看| 黄色 视频免费看| 精品乱码久久久久久99久播| 18在线观看网站| 国产亚洲午夜精品一区二区久久| 成人精品一区二区免费| 成人av一区二区三区在线看| 久久天堂一区二区三区四区| 男女边摸边吃奶| 欧美国产精品一级二级三级| 女同久久另类99精品国产91| 一区福利在线观看| 中文字幕另类日韩欧美亚洲嫩草| 人人妻人人澡人人看| 精品国内亚洲2022精品成人 | 国产精品一区二区免费欧美| 黄色成人免费大全| tube8黄色片| 午夜福利在线免费观看网站| 免费在线观看黄色视频的| 国产精品一区二区免费欧美| 亚洲成人免费电影在线观看| 亚洲成人国产一区在线观看| 又大又爽又粗| 国产日韩一区二区三区精品不卡| 日韩一区二区三区影片| 另类亚洲欧美激情| 亚洲精华国产精华精| 国产精品久久电影中文字幕 | 免费在线观看日本一区| 丝袜在线中文字幕| 捣出白浆h1v1| a在线观看视频网站| 欧美另类亚洲清纯唯美| 激情在线观看视频在线高清 | 久久国产精品男人的天堂亚洲| 国产成人精品无人区| www.精华液| 热99国产精品久久久久久7| 亚洲欧美色中文字幕在线| 免费观看人在逋| 欧美成人免费av一区二区三区 | www.熟女人妻精品国产| 美女国产高潮福利片在线看| 国产淫语在线视频| 国产精品久久久久久人妻精品电影 | 亚洲精品久久午夜乱码| 在线观看免费视频日本深夜| 在线永久观看黄色视频| 久久精品91无色码中文字幕| 亚洲国产欧美在线一区| 丰满少妇做爰视频| 狠狠精品人妻久久久久久综合| 亚洲成av片中文字幕在线观看| 亚洲人成电影观看| 黑人欧美特级aaaaaa片| 麻豆成人av在线观看| 午夜福利一区二区在线看| 狠狠狠狠99中文字幕| 久久久国产一区二区| 欧美精品啪啪一区二区三区| 丁香欧美五月| 国产三级黄色录像| 国产欧美日韩综合在线一区二区| 成人18禁在线播放| 老司机深夜福利视频在线观看| 三上悠亚av全集在线观看| 日韩熟女老妇一区二区性免费视频| 99久久国产精品久久久| 精品第一国产精品| 女同久久另类99精品国产91| 久久天堂一区二区三区四区| 中文欧美无线码| 精品少妇久久久久久888优播| 在线亚洲精品国产二区图片欧美| 欧美日本中文国产一区发布| 一级片免费观看大全| 99久久精品国产亚洲精品| netflix在线观看网站| 精品亚洲乱码少妇综合久久| 国产免费现黄频在线看| 国产一区二区在线观看av| 国产在视频线精品| 久久国产精品男人的天堂亚洲| 久久久久久久精品吃奶| 精品久久蜜臀av无| 午夜激情av网站| 蜜桃在线观看..| 久久午夜亚洲精品久久| 窝窝影院91人妻| 精品人妻1区二区| 叶爱在线成人免费视频播放| 色婷婷av一区二区三区视频| 国产成+人综合+亚洲专区| 天天影视国产精品| 亚洲av国产av综合av卡| 亚洲精品国产色婷婷电影| 在线亚洲精品国产二区图片欧美| 欧美在线黄色| 国产在线免费精品| bbb黄色大片| 久久青草综合色| 丝袜在线中文字幕| 久久国产精品大桥未久av| 欧美成人免费av一区二区三区 | 两人在一起打扑克的视频| 我要看黄色一级片免费的| 亚洲精品久久成人aⅴ小说| 丝袜人妻中文字幕| 啦啦啦视频在线资源免费观看| 99久久人妻综合| 少妇裸体淫交视频免费看高清 | 午夜福利免费观看在线| 女警被强在线播放| 亚洲自偷自拍图片 自拍| 日韩欧美一区二区三区在线观看 | 女性生殖器流出的白浆| 最新在线观看一区二区三区| 美国免费a级毛片| 国产精品国产av在线观看| 日本av免费视频播放| 精品国产国语对白av| 久久国产精品男人的天堂亚洲| 国产黄色免费在线视频| 黄片播放在线免费| 999精品在线视频| 午夜福利在线免费观看网站| 2018国产大陆天天弄谢| 激情视频va一区二区三区| 国产欧美亚洲国产| 国产精品秋霞免费鲁丝片| 欧美 亚洲 国产 日韩一| av在线播放免费不卡| 久久精品91无色码中文字幕| 91老司机精品| 成人国产一区最新在线观看| 国内毛片毛片毛片毛片毛片| 美女扒开内裤让男人捅视频| 精品国产亚洲在线| 国产精品.久久久| 亚洲avbb在线观看| 国产一区二区三区视频了| 欧美激情高清一区二区三区| 国产成人啪精品午夜网站| 午夜视频精品福利| 最近最新中文字幕大全免费视频| 久久热在线av| 日韩欧美国产一区二区入口| 在线看a的网站| 天堂8中文在线网| 亚洲精品在线观看二区| a级毛片在线看网站| 国产成人精品久久二区二区免费| 精品视频人人做人人爽| 久久九九热精品免费| 国产亚洲精品一区二区www | 国产主播在线观看一区二区| 在线观看免费日韩欧美大片| 色播在线永久视频| 一区福利在线观看| 久久精品亚洲av国产电影网| 国产1区2区3区精品| 欧美av亚洲av综合av国产av| 成人手机av| 日韩免费高清中文字幕av| 免费在线观看完整版高清| 一区二区av电影网| 97人妻天天添夜夜摸| 色94色欧美一区二区| 午夜激情久久久久久久| 免费在线观看影片大全网站| 亚洲av电影在线进入| 成人18禁在线播放| 欧美亚洲日本最大视频资源| 51午夜福利影视在线观看| 最新美女视频免费是黄的| 97在线人人人人妻| 这个男人来自地球电影免费观看| 国产av国产精品国产| av欧美777| 99国产综合亚洲精品| 后天国语完整版免费观看| 精品一区二区三区四区五区乱码| 欧美日韩福利视频一区二区| 国产精品自产拍在线观看55亚洲 | 国产在线精品亚洲第一网站| 亚洲综合色网址| 久久久久网色| 久久久水蜜桃国产精品网| 亚洲国产中文字幕在线视频| 久久久久久人人人人人| 电影成人av| 啦啦啦视频在线资源免费观看| 考比视频在线观看| 国产成人免费观看mmmm| 国产亚洲一区二区精品| 男女之事视频高清在线观看| 50天的宝宝边吃奶边哭怎么回事| 99re在线观看精品视频| 国产精品秋霞免费鲁丝片| 色婷婷久久久亚洲欧美| av片东京热男人的天堂| av一本久久久久| 一个人免费看片子| 美女主播在线视频| 国产精品98久久久久久宅男小说| 欧美日韩中文字幕国产精品一区二区三区 | 精品人妻在线不人妻| 欧美日韩一级在线毛片| 亚洲色图av天堂| 女人久久www免费人成看片| 国产亚洲午夜精品一区二区久久| 大型黄色视频在线免费观看| 在线观看免费日韩欧美大片| 搡老熟女国产l中国老女人| 欧美老熟妇乱子伦牲交| 成年女人毛片免费观看观看9 | 黄色怎么调成土黄色| 精品一品国产午夜福利视频| 国产熟女午夜一区二区三区| 精品福利观看| 精品国产乱码久久久久久男人| 国产成人免费无遮挡视频| kizo精华| 久久久国产一区二区| 91成人精品电影| 美女主播在线视频| 亚洲欧美日韩另类电影网站| 亚洲九九香蕉| 国精品久久久久久国模美| 1024香蕉在线观看| 大码成人一级视频| 大型黄色视频在线免费观看| a在线观看视频网站| 亚洲色图av天堂| 精品国产一区二区三区久久久樱花| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 久久精品国产亚洲av高清一级| 久久青草综合色| 18禁国产床啪视频网站| 可以免费在线观看a视频的电影网站| 嫩草影视91久久| 这个男人来自地球电影免费观看| 又黄又粗又硬又大视频| 国产精品.久久久| 淫妇啪啪啪对白视频| 日韩欧美一区视频在线观看| 日韩中文字幕视频在线看片| 亚洲专区中文字幕在线| 亚洲人成电影观看| 美女国产高潮福利片在线看| 国产免费现黄频在线看| 亚洲欧美日韩高清在线视频 | 国产成人欧美在线观看 | tube8黄色片| 热99re8久久精品国产| 日韩一区二区三区影片| 免费观看人在逋| 丝袜美足系列| 一二三四社区在线视频社区8| 成人18禁高潮啪啪吃奶动态图| 巨乳人妻的诱惑在线观看| 日本wwww免费看| 久久久久久人人人人人| 成人手机av| 在线观看66精品国产| 精品一区二区三卡| 亚洲av日韩在线播放| 黑人操中国人逼视频| 国产午夜精品久久久久久| 国产在线精品亚洲第一网站| 欧美午夜高清在线| 黑人猛操日本美女一级片| 天天影视国产精品| e午夜精品久久久久久久| 999久久久国产精品视频| 久久久久国内视频| 国产精品熟女久久久久浪| 欧美激情高清一区二区三区| 亚洲中文日韩欧美视频| 久久久久精品人妻al黑| 久久午夜综合久久蜜桃| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 亚洲伊人色综图| 女警被强在线播放| 亚洲国产欧美日韩在线播放| 又大又爽又粗| 欧美日韩av久久| 51午夜福利影视在线观看| 黑人操中国人逼视频| 日韩熟女老妇一区二区性免费视频| 亚洲 欧美一区二区三区| 欧美黑人欧美精品刺激| 91老司机精品| 亚洲第一青青草原| 9色porny在线观看| 男女无遮挡免费网站观看| 婷婷丁香在线五月| 国产成人系列免费观看| 精品一区二区三区四区五区乱码| 久久久久久久国产电影| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩黄片免| 大香蕉久久成人网| 亚洲中文字幕日韩| 汤姆久久久久久久影院中文字幕| 变态另类成人亚洲欧美熟女 | 99精品久久久久人妻精品| 热re99久久精品国产66热6| 青草久久国产| 亚洲精品中文字幕在线视频| 久久九九热精品免费| 成人特级黄色片久久久久久久 | 中文字幕最新亚洲高清| 国产精品国产av在线观看| 黑人猛操日本美女一级片| 欧美黑人精品巨大| 在线播放国产精品三级| 另类亚洲欧美激情| 中文字幕最新亚洲高清| 欧美黄色片欧美黄色片| 50天的宝宝边吃奶边哭怎么回事| 99国产精品免费福利视频| 啦啦啦在线免费观看视频4| 国产一区二区激情短视频| 精品卡一卡二卡四卡免费| 欧美日韩亚洲综合一区二区三区_| 亚洲av片天天在线观看| 久久国产精品大桥未久av| 一区二区三区乱码不卡18| 久久中文字幕一级| 我要看黄色一级片免费的| 18禁美女被吸乳视频| 久久久久网色| 久久天躁狠狠躁夜夜2o2o| 国产视频一区二区在线看| 欧美精品一区二区免费开放| 亚洲精品国产区一区二| 99精品久久久久人妻精品| 久久久久久久国产电影| aaaaa片日本免费| 搡老岳熟女国产| 2018国产大陆天天弄谢| 91麻豆精品激情在线观看国产 | 国产男女超爽视频在线观看| 一本综合久久免费| 久久精品国产亚洲av高清一级| 国产精品国产高清国产av | 亚洲av欧美aⅴ国产| 少妇 在线观看| 国产欧美日韩一区二区精品| 久久人人97超碰香蕉20202| 桃花免费在线播放| 精品国产一区二区三区四区第35| 2018国产大陆天天弄谢| 久久久欧美国产精品| 一级毛片女人18水好多| 成人免费观看视频高清| 国产亚洲精品一区二区www | 又大又爽又粗| 久久久久久久大尺度免费视频| 久久这里只有精品19| 少妇精品久久久久久久| svipshipincom国产片| 亚洲第一青青草原| 亚洲色图综合在线观看| 天天添夜夜摸| a级毛片黄视频| a在线观看视频网站| 国产高清视频在线播放一区| 操美女的视频在线观看| 99re在线观看精品视频| 欧美大码av| 视频区欧美日本亚洲| av电影中文网址| 久久精品国产亚洲av高清一级| 麻豆成人av在线观看| 在线观看免费午夜福利视频| 欧美 日韩 精品 国产| 国产成人av教育| 日韩制服丝袜自拍偷拍| 大码成人一级视频| 国产精品久久久人人做人人爽| 精品国产乱子伦一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 |