• <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一区| 国产在线免费精品| 男女边摸边吃奶| a级毛片在线看网站| 精品一区二区免费观看| e午夜精品久久久久久久| 啦啦啦在线免费观看视频4| 亚洲欧美中文字幕日韩二区| 久久国产精品男人的天堂亚洲| 美女高潮到喷水免费观看| 青草久久国产| 青春草亚洲视频在线观看| 亚洲精品日本国产第一区| 涩涩av久久男人的天堂| videosex国产| 午夜老司机福利片| 国产99久久九九免费精品| 少妇人妻精品综合一区二区| 热99久久久久精品小说推荐| 国产黄频视频在线观看| 日韩成人av中文字幕在线观看| 亚洲精品成人av观看孕妇| 天天影视国产精品| a级毛片在线看网站| 高清av免费在线| 亚洲综合精品二区| 99久久人妻综合| 精品酒店卫生间| 欧美日韩视频精品一区| 亚洲 欧美一区二区三区| 免费看av在线观看网站| 在线精品无人区一区二区三| 日韩伦理黄色片| 久久99热这里只频精品6学生| a 毛片基地| 国产淫语在线视频| 男的添女的下面高潮视频| 国产男人的电影天堂91| 国产精品99久久99久久久不卡 | 久久性视频一级片| 黑人欧美特级aaaaaa片| 亚洲国产成人一精品久久久| 久久久久久久久久久免费av| 午夜老司机福利片| 99久久99久久久精品蜜桃| 嫩草影视91久久| 五月开心婷婷网| av福利片在线| av电影中文网址| 涩涩av久久男人的天堂| 中文字幕另类日韩欧美亚洲嫩草| 日韩中文字幕视频在线看片| 人妻人人澡人人爽人人| 精品国产一区二区三区四区第35| 亚洲精品视频女| 激情五月婷婷亚洲| 久久精品久久久久久久性| 亚洲,欧美精品.| av在线观看视频网站免费| 男女边吃奶边做爰视频| 亚洲成人一二三区av| 亚洲国产精品999| 老鸭窝网址在线观看| 熟女少妇亚洲综合色aaa.| 男人爽女人下面视频在线观看| 亚洲少妇的诱惑av| 永久免费av网站大全| 在线观看三级黄色| 侵犯人妻中文字幕一二三四区| 两个人免费观看高清视频| 久久精品国产亚洲av高清一级| 各种免费的搞黄视频| av不卡在线播放| 国产一区有黄有色的免费视频| av.在线天堂| 五月开心婷婷网| 婷婷色综合大香蕉| 国产成人免费无遮挡视频| 久久久久久久国产电影| 中文字幕另类日韩欧美亚洲嫩草| www.精华液| 少妇 在线观看| 国产成人午夜福利电影在线观看| 中文字幕人妻熟女乱码| 一级毛片我不卡| 在线观看人妻少妇| 电影成人av| 国产精品一区二区在线不卡| 狠狠精品人妻久久久久久综合| 岛国毛片在线播放| 久久久久国产一级毛片高清牌| 高清欧美精品videossex| 国产xxxxx性猛交| 交换朋友夫妻互换小说| 国产视频首页在线观看| 三上悠亚av全集在线观看| 国产精品.久久久| 亚洲伊人久久精品综合| 中文乱码字字幕精品一区二区三区| 日韩精品免费视频一区二区三区| 国产精品成人在线| 国产成人a∨麻豆精品| 十八禁高潮呻吟视频| 你懂的网址亚洲精品在线观看| 成年女人毛片免费观看观看9 | 亚洲精品国产区一区二| 夫妻性生交免费视频一级片| 七月丁香在线播放| 亚洲欧美中文字幕日韩二区| 久久久久精品人妻al黑| 18禁裸乳无遮挡动漫免费视频| 男人爽女人下面视频在线观看| 欧美黑人精品巨大| 免费观看a级毛片全部| 国产成人免费无遮挡视频| 亚洲av在线观看美女高潮| 香蕉丝袜av| 一级片'在线观看视频| 国产 精品1| 亚洲精品一二三| 亚洲天堂av无毛| 好男人视频免费观看在线| 波多野结衣av一区二区av| 国产成人精品无人区| 一级毛片电影观看| 九九爱精品视频在线观看| 国产高清不卡午夜福利| 自拍欧美九色日韩亚洲蝌蚪91| 久久久久精品性色| 啦啦啦 在线观看视频| 亚洲国产av影院在线观看| 咕卡用的链子| 1024香蕉在线观看| 操美女的视频在线观看| 国产精品国产三级国产专区5o| 捣出白浆h1v1| 亚洲av福利一区| 少妇 在线观看| 99国产精品免费福利视频| bbb黄色大片| 午夜福利影视在线免费观看| 18禁国产床啪视频网站| 菩萨蛮人人尽说江南好唐韦庄| 欧美国产精品va在线观看不卡| 久热这里只有精品99| 波多野结衣一区麻豆| 精品卡一卡二卡四卡免费| 巨乳人妻的诱惑在线观看| 亚洲激情五月婷婷啪啪| netflix在线观看网站| 国产精品久久久av美女十八| 制服诱惑二区| 99久国产av精品国产电影| 免费看不卡的av| 国产熟女午夜一区二区三区| 国产精品久久久久久精品电影小说| 亚洲男人天堂网一区| 亚洲美女视频黄频| 日韩 亚洲 欧美在线| 亚洲国产精品999| 亚洲国产av新网站| 亚洲国产最新在线播放| 中文欧美无线码| 男女国产视频网站| 亚洲av成人不卡在线观看播放网 | 90打野战视频偷拍视频| 国产精品欧美亚洲77777| 久久久久久免费高清国产稀缺| 久久久久久久久久久免费av| 久久ye,这里只有精品| 婷婷成人精品国产| 国产女主播在线喷水免费视频网站| 搡老乐熟女国产| 国产爽快片一区二区三区| 制服丝袜香蕉在线| 综合色丁香网| 日韩精品免费视频一区二区三区| 日韩不卡一区二区三区视频在线| 97精品久久久久久久久久精品| 99国产精品免费福利视频| 十分钟在线观看高清视频www| 97人妻天天添夜夜摸| 国产欧美亚洲国产| svipshipincom国产片| 美女视频免费永久观看网站| 亚洲欧美一区二区三区黑人| 亚洲,一卡二卡三卡| 搡老岳熟女国产| 国产精品成人在线| 国产精品久久久久成人av| 成人午夜精彩视频在线观看| 日韩精品有码人妻一区| 国产精品国产三级国产专区5o| 久久女婷五月综合色啪小说| 人人澡人人妻人| 十八禁高潮呻吟视频| 亚洲一区二区三区欧美精品| 一本大道久久a久久精品| 免费看不卡的av| 极品人妻少妇av视频| 午夜激情av网站| 岛国毛片在线播放| 久久午夜综合久久蜜桃| 精品国产一区二区三区四区第35| 99香蕉大伊视频| 亚洲精品日本国产第一区| 久久鲁丝午夜福利片| 亚洲成国产人片在线观看| 国产乱来视频区| 午夜福利在线免费观看网站| 人人澡人人妻人| 国产精品久久久av美女十八| 又大又黄又爽视频免费| 日韩 亚洲 欧美在线| 色94色欧美一区二区| 欧美中文综合在线视频| 亚洲成人免费av在线播放| 91国产中文字幕| 99九九在线精品视频| 亚洲欧美成人精品一区二区| 看免费av毛片| 国产精品香港三级国产av潘金莲 | 久久精品久久精品一区二区三区| 欧美黑人精品巨大| 日韩中文字幕视频在线看片| 乱人伦中国视频| 超色免费av| 无限看片的www在线观看| 免费在线观看完整版高清| 一本色道久久久久久精品综合| 另类精品久久| 在线看a的网站| 亚洲,一卡二卡三卡| 午夜免费观看性视频| 日韩精品有码人妻一区| 蜜桃国产av成人99| 18在线观看网站| 夜夜骑夜夜射夜夜干| 一本一本久久a久久精品综合妖精| 亚洲成国产人片在线观看| 国产成人欧美| 国产黄色免费在线视频| 欧美国产精品va在线观看不卡| 亚洲国产日韩一区二区| 国产高清国产精品国产三级| 午夜福利乱码中文字幕| 色吧在线观看| 国产欧美日韩一区二区三区在线| 欧美日韩成人在线一区二区| 日韩 亚洲 欧美在线| 伦理电影大哥的女人| 免费不卡黄色视频| 国产成人av激情在线播放| 国产精品无大码| 在线观看国产h片| 亚洲精品第二区| 日韩视频在线欧美| 国产女主播在线喷水免费视频网站| 大片电影免费在线观看免费| 少妇人妻精品综合一区二区| 欧美黄色片欧美黄色片| 日韩伦理黄色片| 免费在线观看视频国产中文字幕亚洲 | 伊人久久国产一区二区| 国产有黄有色有爽视频| 亚洲国产欧美在线一区| 秋霞在线观看毛片| 九草在线视频观看| 亚洲av男天堂| 丝袜脚勾引网站| 熟女av电影| 色吧在线观看| 亚洲少妇的诱惑av| bbb黄色大片| 亚洲精品在线美女| 免费高清在线观看视频在线观看| 一二三四在线观看免费中文在| 成人亚洲欧美一区二区av| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 国产激情久久老熟女| 亚洲少妇的诱惑av| 国产一区二区在线观看av| 精品国产超薄肉色丝袜足j| 伊人亚洲综合成人网| 亚洲欧美清纯卡通| 免费看不卡的av| av网站免费在线观看视频| 90打野战视频偷拍视频| 精品国产乱码久久久久久小说| 咕卡用的链子| 赤兔流量卡办理| 在线观看免费日韩欧美大片| 久久韩国三级中文字幕| 亚洲人成网站在线观看播放| 一区二区av电影网| 亚洲精品国产av成人精品| 日韩 欧美 亚洲 中文字幕| 亚洲一区二区三区欧美精品| 制服人妻中文乱码| 欧美亚洲日本最大视频资源| 女性生殖器流出的白浆| 爱豆传媒免费全集在线观看| 国产精品99久久99久久久不卡 | 日韩大片免费观看网站| 国产 精品1| 精品少妇黑人巨大在线播放| 七月丁香在线播放| 777久久人妻少妇嫩草av网站| 国产毛片在线视频| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 国产伦理片在线播放av一区| 精品一区二区三区四区五区乱码 | 天堂8中文在线网| 伊人久久国产一区二区| 国产精品偷伦视频观看了| 国产片特级美女逼逼视频| 99国产综合亚洲精品| 日日爽夜夜爽网站| 捣出白浆h1v1| 久久97久久精品| 亚洲男人天堂网一区| 成人影院久久| 亚洲精品一二三| 亚洲欧洲日产国产| 亚洲精品久久成人aⅴ小说| e午夜精品久久久久久久| 午夜精品国产一区二区电影| 久久99精品国语久久久| 97精品久久久久久久久久精品| 国产成人精品无人区| 又黄又粗又硬又大视频| 99国产综合亚洲精品| 午夜福利视频精品| 人人妻人人澡人人爽人人夜夜| 91老司机精品| 欧美黑人欧美精品刺激| 青春草亚洲视频在线观看| 国产午夜精品一二区理论片| 一级毛片电影观看| 亚洲欧美成人精品一区二区| 国产视频首页在线观看| 国产毛片在线视频| 校园人妻丝袜中文字幕| 人人妻人人爽人人添夜夜欢视频| 久久久久久久久免费视频了| 捣出白浆h1v1| 国精品久久久久久国模美| 深夜精品福利| av国产精品久久久久影院| 91精品三级在线观看| 丰满迷人的少妇在线观看| 国产一卡二卡三卡精品 | 天天影视国产精品| 大片电影免费在线观看免费| 两个人免费观看高清视频| 国产精品三级大全| 热99久久久久精品小说推荐| 两个人看的免费小视频| 亚洲久久久国产精品| 成人漫画全彩无遮挡| 1024香蕉在线观看| 人人妻人人添人人爽欧美一区卜| 国产探花极品一区二区| av国产精品久久久久影院| 久久鲁丝午夜福利片| 最近中文字幕高清免费大全6| 高清视频免费观看一区二区| 亚洲,一卡二卡三卡| 国产成人精品福利久久| 亚洲成人手机| 人体艺术视频欧美日本| 不卡av一区二区三区| 久久精品国产亚洲av高清一级| 国产精品熟女久久久久浪| 在线看a的网站| 久久久久精品性色| 汤姆久久久久久久影院中文字幕| 天堂俺去俺来也www色官网| 亚洲人成电影观看| 深夜精品福利| 亚洲av电影在线观看一区二区三区| 免费看不卡的av| 久久午夜综合久久蜜桃| 精品人妻熟女毛片av久久网站| 观看美女的网站| 在线精品无人区一区二区三| 少妇猛男粗大的猛烈进出视频| 人人妻人人澡人人爽人人夜夜| 好男人视频免费观看在线| 亚洲成人av在线免费| 天天躁日日躁夜夜躁夜夜| 老司机影院毛片| av国产久精品久网站免费入址| 国产av码专区亚洲av| 国产在视频线精品| 母亲3免费完整高清在线观看| 国产精品久久久久久精品电影小说| 母亲3免费完整高清在线观看| 一二三四在线观看免费中文在| 国产黄色免费在线视频| 亚洲av综合色区一区| 欧美精品一区二区免费开放| 欧美 日韩 精品 国产| 日韩精品有码人妻一区| 一个人免费看片子| 欧美少妇被猛烈插入视频| 国产一卡二卡三卡精品 | 久久精品aⅴ一区二区三区四区| 亚洲国产中文字幕在线视频| 亚洲成人国产一区在线观看 | 黄色 视频免费看| 免费观看a级毛片全部| 999久久久国产精品视频| 极品人妻少妇av视频| 无遮挡黄片免费观看| 久久久国产一区二区| 少妇被粗大的猛进出69影院| 久久久久久久精品精品| 色网站视频免费| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠久久av| 国产亚洲av高清不卡| 如日韩欧美国产精品一区二区三区| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 国产日韩欧美在线精品| 成人影院久久| 大香蕉久久网| 超碰成人久久| 最新的欧美精品一区二区| 在线观看www视频免费| 别揉我奶头~嗯~啊~动态视频 | 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 无遮挡黄片免费观看| 亚洲国产中文字幕在线视频| 深夜精品福利| netflix在线观看网站| 男女之事视频高清在线观看 | 熟妇人妻不卡中文字幕| 国产精品人妻久久久影院| 精品一品国产午夜福利视频| 青春草亚洲视频在线观看| 成年人免费黄色播放视频| 91精品伊人久久大香线蕉| 国产伦人伦偷精品视频| 国产精品一区二区精品视频观看| 亚洲天堂av无毛| 乱人伦中国视频| 自拍欧美九色日韩亚洲蝌蚪91| 久久久欧美国产精品| 天美传媒精品一区二区| 大陆偷拍与自拍| 老司机亚洲免费影院| 久久免费观看电影| av视频免费观看在线观看| 国产成人啪精品午夜网站| 欧美激情 高清一区二区三区| 国精品久久久久久国模美| 人体艺术视频欧美日本| 最近的中文字幕免费完整| 波多野结衣一区麻豆| 高清在线视频一区二区三区| 国产成人啪精品午夜网站| 亚洲,欧美,日韩| 美女福利国产在线| 国产成人精品久久二区二区91 | 交换朋友夫妻互换小说| 欧美亚洲日本最大视频资源| 十分钟在线观看高清视频www| 亚洲精品久久成人aⅴ小说| 综合色丁香网| 一级毛片黄色毛片免费观看视频| 最近的中文字幕免费完整| 国产熟女午夜一区二区三区| 国产精品麻豆人妻色哟哟久久| www.精华液| 精品少妇黑人巨大在线播放| 又黄又粗又硬又大视频| 超色免费av| 自线自在国产av| 午夜激情久久久久久久| 欧美激情 高清一区二区三区| 91精品伊人久久大香线蕉| 亚洲av日韩在线播放| 日韩中文字幕欧美一区二区 | 美国免费a级毛片| 男女边摸边吃奶| 亚洲精品国产av蜜桃| 在线观看三级黄色| 99国产精品免费福利视频| 人妻人人澡人人爽人人| 国产精品国产三级专区第一集| 熟妇人妻不卡中文字幕| 国产男女内射视频| 黄色视频不卡| 国产精品嫩草影院av在线观看| 国产极品天堂在线| 亚洲精品一区蜜桃| 国产精品 欧美亚洲| 19禁男女啪啪无遮挡网站| 不卡av一区二区三区| 久久精品亚洲av国产电影网| 人人妻,人人澡人人爽秒播 | 国产成人a∨麻豆精品| 国产女主播在线喷水免费视频网站| 亚洲精品国产区一区二| 各种免费的搞黄视频| 亚洲,欧美精品.| 亚洲国产最新在线播放| 免费高清在线观看视频在线观看| 极品少妇高潮喷水抽搐| 男女国产视频网站| 国产男女超爽视频在线观看| 亚洲欧美日韩另类电影网站| 日韩一卡2卡3卡4卡2021年| av视频免费观看在线观看| 美女福利国产在线| 日韩av不卡免费在线播放| 久久久久久人妻| 日本猛色少妇xxxxx猛交久久| 黄色视频不卡| a 毛片基地| 99热国产这里只有精品6| 久久狼人影院| av国产精品久久久久影院| 国产成人精品久久二区二区91 | 日韩大码丰满熟妇| 亚洲伊人色综图| 免费av中文字幕在线| 女人久久www免费人成看片| a级毛片黄视频| 少妇被粗大猛烈的视频| 久久精品aⅴ一区二区三区四区| 国产免费视频播放在线视频| 天天操日日干夜夜撸| 久久久久久久国产电影| 欧美国产精品一级二级三级| 国产精品一区二区在线观看99| 亚洲欧美中文字幕日韩二区| 日韩精品有码人妻一区| 考比视频在线观看| 欧美日韩国产mv在线观看视频| 亚洲一区二区三区欧美精品| av在线app专区| 中文字幕最新亚洲高清| 大陆偷拍与自拍| 日韩中文字幕欧美一区二区 | 国产老妇伦熟女老妇高清| 99精品久久久久人妻精品| 久久婷婷青草| 亚洲国产欧美网| 18在线观看网站| 日本91视频免费播放| 在线观看国产h片| 免费观看人在逋| 欧美亚洲日本最大视频资源| 2018国产大陆天天弄谢| 精品酒店卫生间| 久热爱精品视频在线9| 日本wwww免费看| 宅男免费午夜| 我要看黄色一级片免费的| 香蕉国产在线看| 精品视频人人做人人爽| 丝袜美足系列| av在线观看视频网站免费| 国产一区有黄有色的免费视频| 看免费av毛片| 亚洲国产欧美一区二区综合| 日韩制服丝袜自拍偷拍| 久久影院123| 欧美xxⅹ黑人| 亚洲国产成人一精品久久久| 一级毛片电影观看| 久久亚洲国产成人精品v| 亚洲欧美清纯卡通| 美女主播在线视频| 赤兔流量卡办理| 一本大道久久a久久精品| 观看av在线不卡| 国产午夜精品一二区理论片| 婷婷色综合www| 成人黄色视频免费在线看| 美女中出高潮动态图| 国产1区2区3区精品| 麻豆av在线久日| 女的被弄到高潮叫床怎么办| 亚洲欧洲日产国产| 国产av码专区亚洲av| 欧美 日韩 精品 国产| 免费av中文字幕在线| 你懂的网址亚洲精品在线观看| 晚上一个人看的免费电影| 久久久精品94久久精品| 欧美激情高清一区二区三区 | 又粗又硬又长又爽又黄的视频| 女人久久www免费人成看片| 亚洲欧美激情在线| 久久久久视频综合| 性高湖久久久久久久久免费观看| 麻豆精品久久久久久蜜桃| 中文精品一卡2卡3卡4更新| 女人久久www免费人成看片| 黄网站色视频无遮挡免费观看| 成人三级做爰电影| 久久国产精品大桥未久av| 国产一级毛片在线| 久久久久精品国产欧美久久久 | 青青草视频在线视频观看| 黄色一级大片看看| 日韩免费高清中文字幕av|