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

    基于存在狀態(tài)定量計算模型研究丹參注射液的醇沉精制機制

    2022-10-21 06:54:02李存玉林亞娟李明明彭國平支興蕾
    中草藥 2022年20期
    關(guān)鍵詞:兒茶素鈉冪函數(shù)

    李存玉,林亞娟,李明明,章 蓮,李 碩,彭國平,支興蕾*

    基于存在狀態(tài)定量計算模型研究丹參注射液的醇沉精制機制

    李存玉1, 2, 3,林亞娟4,李明明5,章 蓮1,李 碩1,彭國平1, 2, 3,支興蕾1*

    1. 南京中醫(yī)藥大學藥學院,江蘇 南京 210023 2. 江蘇省經(jīng)典名方工程研究中心,江蘇 南京 210023 3. 江蘇省中藥資源產(chǎn)業(yè)化過程協(xié)同創(chuàng)新中心,江蘇 南京 210023 4. 南京中醫(yī)藥大學校醫(yī)院,江蘇 南京 210023 5. 江蘇神龍藥業(yè)有限公司,江蘇 東臺 224200

    基于成分存在狀態(tài)定量計算模型,分析丹參注射液(Danshen Injection,DI)的醇沉工藝中酚酸類成分的傳遞規(guī)律。以DI中水溶性酚酸類成分丹參素鈉、原兒茶醛、迷迭香酸、丹酚酸B為檢測指標,基于納濾傳質(zhì)數(shù)學模型,收集膜通量、截留率數(shù)據(jù),擬合溶質(zhì)傳質(zhì)系數(shù),以分子態(tài)單體成分為參照,構(gòu)建成分存在狀態(tài)定量計算模型,進而分析中間體密度1.05、1.15、1.25 g/mL時,指標性成分醇沉至75%、85%的轉(zhuǎn)移率與成分存在狀態(tài)及相對分子質(zhì)量的相關(guān)性,以期解明醇沉精制機制。存在狀態(tài)定量計算冪函數(shù)方程的相關(guān)性系數(shù)均大于0.99,DI中間體的密度為1.05~1.25 g/mL時,丹參素鈉、原兒茶醛、迷迭香酸、丹酚酸B分子態(tài)比例分別為1.14%~5.25%、66.70%~80.25%、10.51%~19.11%、1.13%~11.44%,對密度-分子態(tài)比例-醇沉轉(zhuǎn)移率進行相關(guān)性分析,醇沉體積分數(shù)由75%增加至85%,轉(zhuǎn)移率下降比例排序為丹參素鈉>丹酚酸B>迷迭香酸>原兒茶醛,成分存在狀態(tài)為影響醇沉轉(zhuǎn)移率的主導因素,且與成分相對分子質(zhì)量存在一定相關(guān)性。構(gòu)建了DI中間體中酚酸類成分分子狀態(tài)定量計算方法,初步闡明醇沉精制機制,有助于制劑醇沉工藝的標準化控制。

    存在狀態(tài);丹參;丹參注射液;醇沉;丹參素鈉;原兒茶醛;迷迭香酸;丹酚酸B;納濾;冪函數(shù)方程

    中藥制藥過程工藝參數(shù)的標準化控制有利于提升制劑質(zhì)量的均一性[1-2],2015年國家中醫(yī)藥管理局發(fā)布中藥標準化行動,旨在促進中藥產(chǎn)業(yè)提升標準化水平,切實提高中藥產(chǎn)品質(zhì)量。醇沉是中藥制劑生產(chǎn)中常用的精制技術(shù),根據(jù)成分溶解度的差異去除雜質(zhì),保留有效成分。在實際生產(chǎn)過程中,通過數(shù)字化控制乙醇加入速度、溶液攪拌轉(zhuǎn)速保障制劑批次間的一致性,但是醇沉引起質(zhì)量均一性的差異仍是制藥行業(yè)難以解決的技術(shù)難題[3-5]。

    中藥成分的溶解度與存在狀態(tài)直接相關(guān),以離子態(tài)、締合態(tài)、復合態(tài)等形式存在的中藥成分,在醇沉過程中損失明顯。在明確有效成分存在狀態(tài)的前提下,精準化調(diào)節(jié)醇沉參數(shù),將有助于提升制劑質(zhì)量的均一性。目前,用于分析成分結(jié)構(gòu)的紅外光譜、掃描電鏡、差示熱分析等分析技術(shù)難以獲取原溶液環(huán)境下的成分真實存在狀態(tài)[6-8]。納濾截留相對分子質(zhì)量為100~1000,以電荷、孔道位阻為主要分離效應,適用于荷電中藥成分的常溫化精制,前期研究發(fā)現(xiàn)溶質(zhì)存在狀態(tài)決定納濾分離行為,且分離過程中傳質(zhì)系數(shù)與溶質(zhì)質(zhì)量濃度呈冪值相關(guān)[9],基于成分存在狀態(tài)與限域環(huán)境傳質(zhì)行為的相關(guān)性,構(gòu)建存在狀態(tài)定量計算模型,探索中藥制劑醇沉工藝機制。丹參素鈉、原兒茶醛、迷迭香酸和丹酚酸B是丹參注射液(Danshen Injection,DI)中的指標性成分,相對分子質(zhì)量、pKa呈現(xiàn)出一定的梯度性,且醇沉過程中的轉(zhuǎn)移率存在明顯差異,參考DI制法中醇沉環(huán)節(jié)“丹參提取液密度1.15~1.28 g/mL時,分別醇沉至含醇量為75%、85%”的工藝參數(shù),以4種成分丹參素鈉、原兒茶醛、迷迭香酸和丹酚酸B為考察指標,分析“密度-成分狀態(tài)-成分轉(zhuǎn)移率”的相關(guān)性,以期為探索中藥制劑醇沉技術(shù)的內(nèi)在作用機制提供研究基礎(chǔ)。

    1 儀器與材料

    1.1 儀器

    Agilent 1260型高效液相色譜儀,VWD檢測器,美國Agilent公司;10-13型密度計,規(guī)格1.0~1.3 g/L,衡水耀華器械儀表廠;TNZ-1型納濾分離設(shè)備、Fog-machine型中壓泵,南京拓鉒醫(yī)藥科技有限公司;Buchi R-300型旋轉(zhuǎn)蒸發(fā)儀,瑞士Buchi公司;MS105型十萬分之一電子天平、AL204型萬分之一電子天平,瑞士Mettler Toledo集團;PB-10 pH計,德國Sartorius公司;NDJ-5S型黏度計,上海力辰邦西儀器科技有限公司;Millipore Direct-Q5型純水/超純水一體機,美國Millipore公司。

    1.2 材料

    1812型復合聚酰胺納濾膜,截留相對分子質(zhì)量800,星達(泰州)膜科技有限公司;丹參飲片(批號180201)購自安徽省金芙蓉中藥飲片廠,經(jīng)南京中醫(yī)藥大學藥學院嚴輝教授鑒定為唇形科鼠尾草屬植物丹參Bge.的干燥根和根莖;對照品丹參素鈉(批號110855-201915,質(zhì)量分數(shù)97.8%)、原兒茶醛(批號110810-201007,質(zhì)量分數(shù)98.2%)、迷迭香酸(批號111871-201102,質(zhì)量分數(shù)98.6%)、丹酚酸B(批號111562-201615,質(zhì)量分數(shù)96.2%),均購自中國食品藥品檢定研究院;提取物丹參素鈉(批號FY1784,質(zhì)量分數(shù)≥98%)、原兒茶醛(批號FY1784,質(zhì)量分數(shù)≥98%)、迷迭香酸(批號FY1143,質(zhì)量分數(shù)≥98%)和丹酚酸B(批號FY1167,質(zhì)量分數(shù)≥95%),均購自南通飛宇生物科技有限公司;甲醇、乙腈、三氟乙酸為色譜純;乙醇、鹽酸為分析純;水為超純水。

    2 方法與結(jié)果

    2.1 溶液的制備

    2.1.1 DI中間體 參照DI制法,取丹參9.0 kg加水煎煮3次,第1次2 h,第2、3次各1.5 h,合并煎液,濾液減壓濃縮至所考察的密度,即得。

    2.1.2 混合對照品溶液 精密稱取干燥至恒定質(zhì)量的丹參素鈉、原兒茶醛、迷迭香酸和丹酚酸B對照品適量,置于5 mL棕色量瓶中,加10%甲醇溶液制成含丹參素鈉0.420 mg/mL、原兒茶醛0.152 mg/mL、迷迭香酸0.196 mg/mL和丹酚酸B 1.105 mg/mL的混合對照品溶液。

    2.1.3 提取物溶液 參照不同密度時DI中間體中指標性成分質(zhì)量濃度和黏度,分別精密稱取干燥至恒定質(zhì)量的丹參素鈉、原兒茶醛、迷迭香酸和丹酚酸B提取物,超純水超聲溶解,加入蔗糖調(diào)節(jié)溶液黏度,即得4種成分的提取物溶液。

    2.2 樣品檢測

    2.2.1 色譜條件 Agilent 1260型高效液相色譜儀;Thermo Scientific C18色譜柱(250 mm×4.6 mm,5 μm);流動相為乙腈-0.1%甲酸水溶液,梯度洗脫:0~10 min,10%~20%乙腈;10~40 min,20%~45%乙腈;體積流量1.0 mL/min;進樣量10 μL;柱溫30 ℃;檢測波長280 nm。色譜圖見圖1。

    2.2.2 線性關(guān)系考察 分別精密吸取混合對照品0.1、0.2、0.5、1.0、2.0 mL置于5 mL棕色量瓶中,10%甲醇稀釋至刻度,高效液相色譜儀檢測,以峰面積為縱坐標(),指標性成分質(zhì)量濃度為橫坐標(),繪制標準曲線,進行線性回歸,得回歸方程:丹參素鈉=1119+517.67,2=0.999 2;原兒茶醛=1 420.1-459.4,2=0.999 5;迷迭香酸=5252+4 551.3,2=0.999 2;丹酚酸B=2 204.7-8 596.4,2=0.999 6;結(jié)果表明丹參素鈉在8.4~168.0 μg/mL、原兒茶醛在3.04~60.80 μg/mL、迷迭香酸在3.92~78.40 μg/mL、丹酚酸B在22.1~442.0 μg/mL線性關(guān)系良好。

    1-丹參素鈉 2-原兒茶醛 3-迷迭香酸 4-丹酚酸B

    2.2.3 精密度試驗 精密吸取混合對照品溶液10 μL,重復進樣6次,按“2.2.1”項下色譜條件檢測,色譜峰峰面積積分值的RSD分別為丹參素鈉1.06%、原兒茶醛0.61%、迷迭香酸0.86%、丹酚酸B 0.72%。

    2.2.4 穩(wěn)定性試驗 精密吸取按“2.1.1”項下方法處理的DI中間體溶液10 μL,分別于0、1、2、4、8、12、24 h進樣,按“2.2.1”項下色譜條件檢測,色譜峰峰面積積分值的RSD分別為丹參素鈉1.10%、原兒茶醛0.96%、迷迭香酸1.02%、丹酚酸B 1.65%,結(jié)果表明DI中間體在24 h內(nèi)穩(wěn)定。

    2.2.5 重復性試驗 按“2.1.1”項下方法平行制備5份DI中間體溶液,按“2.2.1”項下色譜條件檢測,各成分質(zhì)量濃度的RSD分別為丹參素鈉1.52%、原兒茶醛1.16%、迷迭香酸1.44%、丹酚酸B 2.05%。

    2.2.6 加樣回收率試驗 精密吸取DI中間體溶液6份,分別精密加入混合對照品適量,按“2.2.1”項下色譜條件進行定量測定,計算得平均回收率分別為丹參素鈉101.06%、原兒茶醛102.50%、迷迭香酸98.85%、丹酚酸B 97.21%,其RSD分別為1.64%、1.32%、1.95%、2.18%。

    2.3 密度對成分轉(zhuǎn)移率的影響

    為了探索密度對醇沉過程中指標性成分轉(zhuǎn)移率的影響規(guī)律,取DI中間體,分別采用減壓濃縮(80 ℃)至密度1.05、1.15、1.25 g/mL,分別醇沉至含醇量為75%、85%,靜置12 h,取上清液檢測并分析丹參素鈉、原兒茶醛、迷迭香酸和丹酚酸B的轉(zhuǎn)移率。

    轉(zhuǎn)移率=上清上清/(原原)

    上清、原分別為醇沉上清液和原液中指標性成分質(zhì)量濃度,上清、原分別為醇沉上清液和原液體積

    醇沉至含醇量75%和85%時,指標性成分轉(zhuǎn)移率見表1。隨著乙醇體積分數(shù)和DI中間體密度的升高,丹參中4種酚酸類成分的轉(zhuǎn)移率表現(xiàn)為原兒茶醛>迷迭香酸>丹酚酸B>丹參素鈉。根據(jù)指標性成分的解離常數(shù)(原兒茶醛7.56>迷迭香酸4.01>丹參素鈉3.82>丹酚酸B 2.77),相對分子質(zhì)量(丹酚酸B>迷迭香酸>丹參素鈉>原兒茶醛),結(jié)合DI中間體溶液pH 4.85~5.16,其中原兒茶醛主要以分子態(tài)形式存在,丹參素鈉、迷迭香酸、丹酚酸B主要以離子態(tài)形式存在,醇沉引起的成分轉(zhuǎn)移率變化存在狀態(tài)均的相關(guān)性表現(xiàn)出離子態(tài)轉(zhuǎn)移率低于分子態(tài),且與相對分子質(zhì)量存在相關(guān)性。

    表1 DI中間體密度對醇沉工藝中酚酸類成分轉(zhuǎn)移率的影響

    Table 1 Effect ofDI intermediatedensity on transfer rate of phenolic acids in alcohol precipitation process

    密度/(g·mL?1)75%醇沉轉(zhuǎn)移率/%85%醇沉轉(zhuǎn)移率/% 丹參素鈉原兒茶醛迷迭香酸丹酚酸B丹參素鈉原兒茶醛迷迭香酸丹酚酸B 1.0572.2594.4186.7673.5046.2390.7577.1652.08 1.1553.8891.0581.9067.7139.8186.7068.8747.62 1.2542.3081.6174.2756.2727.0672.2556.2231.60

    2.4 存在狀態(tài)的定量擬合

    將1812型復合聚酰胺納濾膜裝入納濾分離設(shè)備中,通過管路依次連接中壓泵、壓力表、納濾分離設(shè)備、截留閥。取DI中間體0.45 μm微孔濾膜濾過,濾液采用超純水稀釋2、4、8倍,根據(jù)成分的一級解離常數(shù)pKa(丹參素鈉3.82、原兒茶醛7.56、丹酚酸B 2.77,迷迭香酸4.01[10]),調(diào)節(jié)單體成分溶液pH低于pKa 2個酸度值,使得成分分子態(tài)比例≥99%。分別將提取物溶液置于儲液罐中,通過管路連接儲液罐和納濾液端,并進行納濾分離,調(diào)節(jié)中壓泵轉(zhuǎn)速和截留閥,控制跨膜壓力差為0.6、0.8、1.0、1.2、1.4、1.6 MPa,收集相應溶質(zhì)質(zhì)量濃度和跨膜壓力差下的膜通量(v),待儲液罐中溶質(zhì)質(zhì)量濃度相對穩(wěn)定時,取樣平衡液并收集納濾液,HPLC檢測平衡液及納濾液中溶質(zhì)質(zhì)量濃度,根據(jù)式(1)計算截留率(o)。

    o=1-納濾/平衡(1)

    根據(jù)溶質(zhì)傳質(zhì)系數(shù)()與o、v的相關(guān)性,擬合ln[(1-o)v/o]與的線性方程[11],如式(2)所示。收集并計算溶質(zhì)初始質(zhì)量濃度(0)下的、o、v,根據(jù)式(3)擬合與0的冪函數(shù)方程。

    ln[(1-o)v/o]=ln(/)+v/(2)

    1/為斜率,ln(/)為截距,為膜傳質(zhì)性能常數(shù),為膜厚度

    =0(3)

    為相關(guān)系數(shù),為冪值

    以單體成分構(gòu)建的分子態(tài)溶質(zhì)質(zhì)量濃度與傳質(zhì)系數(shù)的冪函數(shù)方程為參照,對系列密度的DI中間體進行納濾分離,擬合與0的冪函數(shù)方程。對比單體成分與DI中間體相同值時,根據(jù)式(3)計算中間體目標成分相較于單體成分分離環(huán)境下的質(zhì)量濃度(丹),根據(jù)式(4)計算成分分子態(tài)比例。

    分子態(tài)比例=1-丹/0(4)

    2.4.1 丹參素鈉的分子態(tài)定量擬合 如圖2所示,在密度為1.05 g/mL時,丹參素鈉單體及提取物的ln[(1-o)v/o]和v線性關(guān)系良好,相關(guān)性系數(shù)均大于0.98,其中單體成分的相關(guān)性數(shù)據(jù)點中,質(zhì)量濃度為48.60 μg/mL相對分散,在丹參素鈉的單一溶液中,溶質(zhì)質(zhì)量濃度對傳質(zhì)系數(shù)影響強于中間體。根據(jù)計算得到的傳質(zhì)系數(shù),擬合與溶質(zhì)質(zhì)量濃度的冪函數(shù)方程,結(jié)果見表2,相關(guān)性系數(shù)均大于0.99,丹參素鈉以分子態(tài)形式存在的比例為5.25%。

    如圖3所示,密度升高至1.15 g/mL,丹參素鈉質(zhì)量濃度對傳質(zhì)系數(shù)的影響程度下降,單體成分的ln[(1-o)v/o]和v相關(guān)性數(shù)據(jù)更為集中,丹參素鈉及中間體的ln[(1-o)v/o]和v線性關(guān)系、質(zhì)量濃度與傳質(zhì)系數(shù)的冪函數(shù)方程,結(jié)果見表3,相關(guān)性系數(shù)均大于0.99。計算得到丹參素鈉在DI中間體中分子態(tài)比例為1.40%。

    圖2 密度1.05 g?mL?1時丹參素鈉ln[(1-Ro)Jv/Ro]對JV的相關(guān)性

    表2 密度1.05 g?mL?1時丹參素鈉傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 2 Mass transfer coefficient, power function equation and molecular state percent of sodium danshensu at density of 1.05 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 48.6016.0510.20k=10.509 0 C00.109 5k=5.039 1 C00.181 25.25 24.3014.539.01R2=0.998 2R2=0.997 8 12.6512.977.90 6.8311.837.18

    圖3 密度1.15 g?mL?1時丹參素鈉ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表3 密度1.15 g?mL?1時丹參素鈉傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 3 Mass transfer coefficient, power function equation and molecular state percent of sodium danshensu at density of 1.15 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 79.2011.349.48k=9.465 4 C00.041 9k=7.384 4 C00.057 71.40 39.6011.109.14R2=0.990 8R2=0.995 4 19.8010.708.81 9.9010.428.40

    密度升高至1.25 g/mL,ln[(1-o)v/o]和v線性方程成立(圖4),相關(guān)性系數(shù)均大于0.98,且溶質(zhì)質(zhì)量濃度與傳質(zhì)系數(shù)的冪函數(shù)方程的相關(guān)性系數(shù)也均大于0.99,丹參素鈉以分子態(tài)形式存在的比例為1.14%,結(jié)果見表4。

    隨著DI中間體密度由1.05 g/mL升高至1.25 g/mL,同時存在丹酚酸類成分的氧化分解,丹參素鈉的質(zhì)量濃度也隨之升高48.6~126.6 μg/mL,以分子態(tài)形式存在的比例為5.25%~1.14%,結(jié)合表1中醇沉引起的轉(zhuǎn)移率變化,以離子態(tài)、締合態(tài)等形式存在的丹參素鈉隨著乙醇質(zhì)量濃度的增加,溶解度下降,引起轉(zhuǎn)移率下降。

    2.4.2 原兒茶醛的分子態(tài)定量擬合 如圖5所示,在密度為1.05 g/mL時,原兒茶醛在醇沉過程中的轉(zhuǎn)移率優(yōu)于其他3種指標性成分,與其存在狀態(tài)和相對分子質(zhì)量等因素相關(guān),結(jié)果見表5;在DI中間體密度為1.05 g/mL時,ln[(1-o)v/o]和v線性關(guān)系良好(2>0.98)。相較于原兒茶醛單體溶液,在不同溶質(zhì)質(zhì)量濃度環(huán)境下,中間體納濾分離時的傳質(zhì)系數(shù)擬合數(shù)據(jù)相似。從表5中可以看出,中間體中原兒茶醛傳質(zhì)系數(shù)變化程度低于單體溶液。結(jié)合表1中原兒茶醛醇沉轉(zhuǎn)移率均大于90%,說明以分子態(tài)(80.25%)形式存在的原兒茶醛易溶于高體積分數(shù)乙醇中,且相對分子質(zhì)量低不易被醇沉物包裹。

    隨著密度增加至1.15 g/mL,丹參酚酸類成分因熱處理水解,引起原兒茶醛質(zhì)量濃度的升高,在進行l(wèi)n[(1-o)v/o]和v線性擬合時,原兒茶醛在單體及中間體溶液中表現(xiàn)出相似的分離行為(圖6),且傳質(zhì)系數(shù)也相對接近(表6),在傳質(zhì)系數(shù)線性方程和冪函數(shù)方程均成立的前提下(2>0.98),醇沉至75%及85%時,原兒茶醛轉(zhuǎn)移率下降均低于5%,且此時分子態(tài)比例為78.50%,說明原兒茶醛在醇沉過程中轉(zhuǎn)移率與其分子態(tài)比例息息相關(guān)。

    圖4 密度1.25 g?mL?1時丹參素鈉ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表4 密度1.25 g?mL?1時丹參素鈉傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 4 Mass transfer coefficient, power function equation and molecular state percent of sodium danshensu at density of 1.25 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 126.606.635.52k=5.445 0 C00.040 6k=3.286 2 C00.106 71.14 63.306.445.07R2=0.999 7R2=0.990 5 31.656.274.81 15.836.094.39

    圖5 密度1.05 g?mL?1時原兒茶醛ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表5 密度1.05 g?mL?1時原兒茶醛傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 5 Mass transfer coefficient, power function equation and molecular state percent of protocatechuic aldehyde at density of 1.05 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 18.4216.0815.27k=8.948 5 C00.206 1k=11.847 C00.086 180.25 9.2114.3314.25R2=0.991 2R2=0.996 1 4.6112.4713.57 2.3010.4612.72

    在密度為1.25 g/mL的溶液環(huán)境下,原兒茶醛單體溶液的ln[(1-o)v/o]和V相關(guān)性數(shù)據(jù)點較中間體分散(圖7),但相關(guān)性系數(shù)均大于0.98。隨著原兒茶醛質(zhì)量濃度的升高,傳質(zhì)系數(shù)呈現(xiàn)上升趨勢,符合限域空間的溶解-擴散理論[12],結(jié)果見表7。以單體成分為參照,擬合出中間體中原兒茶醛的分子態(tài)比例為66.70%。

    DI中間體密度由1.05 g/mL升高至1.25 g/mL,分子態(tài)形式存在原兒茶醛比例為80.25%~66.70%,結(jié)合表1中醇沉轉(zhuǎn)移率變化規(guī)律,原兒茶醛的醇沉機制是以相對分子質(zhì)量為前提,分子狀態(tài)為主導的精制過程。

    圖6 密度1.15 g?mL?1時原兒茶醛ln[(1-Ro)Jv/Ro]對JV的相關(guān)性

    表6 密度1.15 g?mL?1時原兒茶醛傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 6 Mass transfer coefficient, power function equation and molecular state percent of protocatechuic aldehyde at density of 1.15 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 30.5412.3611.72k=6.874 0 C00.171 7k=8.837 0 C00.080 978.50 15.2710.9310.95R2=0.998 1R2=0.992 6 7.649.8410.37 3.828.619.90

    圖7 密度1.25 g?mL?1時原兒茶醛ln[(1-Ro)Jv/Ro]對JV的相關(guān)性

    表7 密度1.25 g?mL?1時原兒茶醛傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 7 Mass transfer coefficient, power function equation and molecular state percent of protocatechuic aldehyde at density of 1.25 g?mL?1

    質(zhì)量濃度/(μg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 73.208.488.18k=4.179 7 C00.167 7k=4.586 4 C00.133 466.70 36.607.747.35R2=0.990 2R2=0.997 3 18.306.906.76 9.155.986.18

    2.4.3 迷迭香酸的分子態(tài)定量擬合 結(jié)果見圖8,密度為1.05 g/mL的迷迭香酸單體和DI中間體,進行l(wèi)n[(1-o)v/o]和v進行線性擬合,相關(guān)性系數(shù)均大于0.98。對比圖8中的散點分布特征,中間體較單體成分相對分散,在相同黏度環(huán)境下,中間體溶液的膜通量明顯偏低。以迷迭香酸單體溶液為參照,擬合溶質(zhì)質(zhì)量濃度與傳質(zhì)系數(shù)的冪函數(shù)方程,結(jié)果見表8,相關(guān)性系數(shù)均大于0.99,DI中間體在質(zhì)量濃度為0.580 mg/mL時,迷迭香酸以分子態(tài)形式存在的比例為19.11%。

    密度增加至1.15 g/mL,膜通量均呈現(xiàn)下降趨勢,ln[(1-o)v/o]和v相關(guān)性系數(shù)均大于0.95(圖9),DI中間體在密度1.15 g/mL時,迷迭香酸在不同質(zhì)量濃度下的相關(guān)性數(shù)據(jù)點存在一定差異,從表9中也可看出,隨著中間體溶液的稀釋,迷迭香酸的傳質(zhì)系數(shù)相對穩(wěn)定。迷迭香酸在中間體質(zhì)量濃度為0.63 mg/mL時,迷迭香酸以分子態(tài)形式存在的比例為14.60%。

    在密度為1.25 g/mL的溶液環(huán)境下,在濃差極化效應的作用下,迷迭香酸單體及DI中間體的膜通量和跨膜壓力差的線性關(guān)系降低(圖10),膜通量進一步的下降,但ln[(1-o)v/o]和v相關(guān)性系數(shù)均大于0.92,且傳質(zhì)系數(shù)與溶質(zhì)質(zhì)量濃度的冪函數(shù)方程均成立(表10),以分子態(tài)形式存在的迷迭香酸比例為10.51%。

    對表1中醇沉過程中迷迭香酸的轉(zhuǎn)移率進行相關(guān)性分析,隨著DI中間體密度的增加,迷迭香酸分子態(tài)比例為19.11%~10.51%,分子態(tài)比例的降低,75%醇沉轉(zhuǎn)移率由86.76%下降至74.27%,85%醇沉轉(zhuǎn)移率由77.16%下降至56.22%。

    圖8 密度1.05 g?mL?1時迷迭香酸ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表8 密度1.05 g?mL?1時迷迭香酸傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 8 Mass transfer coefficient, power function equation and molecular state percent of rosmarinic acid at density of 1.05 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 0.58026.8820.33k=29.570 0 C00.167 6k=21.171 0 C00.068 719.11 0.29024.3919.57R2=0.992 3R2=0.994 8 0.14521.0518.48 0.07319.0517.64

    圖9 密度1.15 g?mL?1時迷迭香酸ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表9 密度1.15 g?mL?1時迷迭香酸傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 9 Mass transfer coefficient, power function equation and molecular state percent of rosmarinic acid at density of 1.15 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 0.6316.2914.84k=16.627 0 C00.046 2k=15.043 0 C00.027 414.60 0.3215.7214.60R2=0.996 6R2=0.997 5 0.1615.3114.29 0.0814.7714.03

    圖10 密度1.25 g?mL?1時迷迭香酸ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表10 密度1.25 g?mL?1時迷迭香酸傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 10 Mass transfer coefficient, power function equation and molecular state percent of rosmarinic acid at density of 1.25 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 0.728.767.44k=8.986 0 C00.078 1k=7.523 0 C00.033 310.51 0.368.317.26R2=0.998 6R2=0.993 3 0.187.837.13 0.097.466.93

    2.4.4 丹酚酸B的分子態(tài)定量擬合 結(jié)果見圖11,對密度為1.05 g/mL溶液環(huán)境下的丹酚酸B ln[(1-o)v/o]和v進行線性關(guān)系計算,相關(guān)性系數(shù)均大于0.98,從圖11中可以看出,單體較中間體溶液的數(shù)據(jù)點分散,說明質(zhì)量濃度對丹酚酸B單體的傳質(zhì)系數(shù)影響相對明顯。以丹酚酸B單體溶液為參照,擬合溶質(zhì)質(zhì)量濃度與傳質(zhì)系數(shù)的冪函數(shù)方程,結(jié)果見表11,DI中間體中以分子態(tài)形式存在的丹酚酸B比例為11.44%,此時在75%及85%醇沉的轉(zhuǎn)移率分別為73.5%、52.08%。

    圖11 密度1.05 g?mL?1時丹酚酸B ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表11 密度1.05 g?mL?1時丹酚酸B傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 11 Mass transfer coefficient, power function equation and molecular state percent of salvianolic acid Bat density of 1.05 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 2.0813.6210.41k=12.399 0 C00.142 7k=9.668 0 C00.093 211.44 1.0412.669.63R2=0.992 3R2=0.994 9 0.5211.319.08 0.2610.178.56

    隨著溶液密度由1.05 g/mL增加至1.15 g/mL時,丹酚酸B單體及中間體溶液的線性方程及冪函數(shù)方程均成立(2>0.98),結(jié)果見圖12。且DI中間體中丹酚酸B的分子態(tài)比例為9.43%,相對穩(wěn)定(表12)。相應的醇沉中的轉(zhuǎn)移率變化也均小于6%,說明密度增加,丹酚酸B分子態(tài)比例出現(xiàn)輕微下降,也引起了醇沉過程中轉(zhuǎn)移率的相應下降。

    DI中間體的密度增加至1.25 g/mL時,丹酚酸B單體及中間體溶液的線性方程及冪函數(shù)方程均成立(2>0.98),結(jié)果見圖13。丹酚酸B的分子態(tài)比例由11.44%下降至1.13%(表13)。在DI中間體的濃縮過程中,因熱處理丹酚酸B轉(zhuǎn)化為原兒茶醛、丹參素鈉等成分[13],引起質(zhì)量濃度下降,會減輕醇沉帶來的成分損失。但相較于密度1.05 g/mL,丹酚酸B在75%和85%醇沉的轉(zhuǎn)移率下降均超過17%,說明丹酚酸B分子態(tài)比例下降改變其在醇沉液中的溶解度,促進了不溶物對丹酚酸B的吸附和包裹,導致轉(zhuǎn)移率明顯下降。

    圖12 密度1.15 g?mL?1時丹酚酸B ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表12 密度1.15 g?mL?1時丹酚酸B傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 12 Mass transfer coefficient, power function equation and molecular state percent of salvianolic acid B at density of 1.15 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 1.369.758.46k=9.580 0 C00.063 3k=8.219 0 C00.090 69.43 0.689.377.92R2=0.998 7R2=0.999 7 0.348.967.46 0.178.557.00

    圖13 密度1.25 g?mL?1時丹酚酸B ln[(1-Ro)Jv/Ro]對Jv的相關(guān)性

    表13 密度1.25 g?mL?1時丹酚酸B傳質(zhì)系數(shù)、冪函數(shù)方程及分子態(tài)比例

    Table 13 Mass transfer coefficient, power function equation and molecular state percent of salvianolic acid B at density of 1.25 g?mL?1

    質(zhì)量濃度/(mg·mL?1)傳質(zhì)系數(shù)/(×10?6m·s?1)冪函數(shù)方程分子態(tài)比例/% 單一成分中間體單一成分中間體 0.745.433.44k=5.594 0 C00.098 3k=3.544 0 C00.118 31.13 0.375.093.11R2=0.997 5R2=0.991 3 0.194.712.93 0.094.442.67

    根據(jù)DI中間體中4種指標成分“密度-分子態(tài)比例-醇沉轉(zhuǎn)移率”的相關(guān)性,醇沉濃度由75%增加至85%,按照轉(zhuǎn)移率下降比例,丹參素鈉>丹酚酸B>迷迭香酸>原兒茶醛,與分子態(tài)比例直接相關(guān),呈現(xiàn)出分子態(tài)比例下降,醇沉轉(zhuǎn)移率下降的規(guī)律,說明存在狀態(tài)為影響生產(chǎn)工序間成分傳遞規(guī)律的主導因素。其次,在分子態(tài)比例相對接近的迷迭香酸和丹酚酸B中,醇沉轉(zhuǎn)移率迷迭香酸>丹酚酸B,中間體中迷迭香酸增加,丹酚酸B下降,但丹酚酸B醇沉損失更為明顯,說明在相對分子質(zhì)量偏低的迷迭香酸在醇沉液中溶解度更高,可以減輕醇沉損失。在密度為1.25 g/mL時,DI中間體溶液中丹參素鈉和丹酚酸B分子態(tài)比例接近,隨著醇質(zhì)量濃度增加至85%,丹參素鈉轉(zhuǎn)移率同比下降36.03%,丹酚酸B下降43.80%,因此在醇沉過程中,成分存在狀態(tài)主導醇沉轉(zhuǎn)移率。

    3 討論

    臨床上DI常用于治療冠心病、心絞痛,其中丹參中的酚酸類成分為主要藥效成分[14],在丹參提取、濃縮過程中,丹酚酸類成分因熱處理分解成丹參素鈉、原兒茶醛、迷迭香酸,中間體中成分含量也隨之升高。在納濾分離不同密度下的中間體溶液時,丹參素鈉、原兒茶醛、迷迭香酸的傳質(zhì)系數(shù)并未隨著質(zhì)量濃度升高而增加,表現(xiàn)出與溶解-擴散分離模型不相符的分離現(xiàn)象。而隨著DI中間體的密度增加,膜通量明顯下降直接改變了溶質(zhì)的傳質(zhì)行為,引起傳質(zhì)系數(shù)下降。同時,上述3種成分的截留率也呈下降趨勢,溶質(zhì)的傳質(zhì)行為符合納濾溶解-擴散分離模型。

    醇沉是中藥制藥精制的常用技術(shù)手段[15],在生產(chǎn)過程中通過調(diào)節(jié)制劑中間體的密度、溫度、pH等參數(shù)控制指標性成分轉(zhuǎn)移率的差異[16-18],而與轉(zhuǎn)移率息息相關(guān)的成分存在狀態(tài),因缺少檢測分析手段,無法真正意義上實現(xiàn)醇沉過程的標準化控制,引起生產(chǎn)批次間質(zhì)量的差異性。本實驗基于納濾傳質(zhì)數(shù)學模型,圍繞“密度-成分狀態(tài)-成分轉(zhuǎn)移率”的相關(guān)性,為中藥酚酸類成分醇沉前溶液中存在狀態(tài)定量計算方法提供了研究策略,為中藥制藥醇沉參數(shù)的精準化控制提供技術(shù)支撐,從而提升制劑質(zhì)量的均一性。

    利益沖突 所有作者均聲明不存在利益沖突

    [1] 熊詩慧, 邱婷, 王學成, 等. 基于質(zhì)量標志物(Q-Marker)的中藥制劑質(zhì)量過程控制方法與策略 [J]. 中草藥, 2021, 52(9): 2534-2540.

    [2] 伍振峰, 林瑞華, 王學成, 等. 基于中藥制藥工程質(zhì)量觀的質(zhì)量控制模式研究 [J]. 中國中藥雜志, 2022, 47(1): 271-278.

    [3] Fernando I P S, Dias M K H M, Madusanka D M D,. Step gradient alcohol precipitation for the purification of low molecular weight fucoidan fromand its UVB protective effects [J]., 2020, 163: 26-35.

    [4] 丁鴻, 徐芳芳, 杜慧, 等. 基于決策樹算法的熱毒寧注射液金銀花青蒿醇沉過程質(zhì)量控制研究 [J]. 中草藥, 2021, 52(19): 5836-5844.

    [5] 蔣美林, 張學瑜, 邵峰, 等. 醇料比對雙黃連制劑醇沉效果及沉淀物形態(tài)影響 [J]. 中草藥, 2020, 51(19): 4954-4959.

    [6] Ishida T, Kitagaki R, Watanabe R,. A study of molecular architectural dynamics of crosslinked urethane during photo-aging by two-dimensional infrared correlation spectroscopy [J]., 2020, 179: 109242.

    [7] Mai N N S, Otsuka Y, Goto S,. Effects of polymer molecular weight on curcumin amorphous solid dispersion; at-line monitoring system based on attenuated total reflectance mid-infrared and near-infrared spectroscopy [J]., 2021, 61: 102278.

    [8] Kwon B, Mandal T, Elkins M R,. Cholesterol interaction with the trimeric HIV fusion protein gp41 in lipid bilayers investigated by solid-state NMR spectroscopy and molecular dynamics simulations [J]., 2020, 432(16): 4705-4721.

    [9] 李存玉, 蔣佳麗, 鄒雨岑, 等. 基于傳質(zhì)數(shù)學模型研究丹參-枳實復方中成分存在狀態(tài)與納濾分離機制 [J]. 中草藥, 2021, 52(15): 4544-4551.

    [10] 李存玉, 吳鑫, 馬赟, 等. 夏枯草中迷迭香酸的納濾分離行為及富集工藝對比 [J]. 中藥材, 2018, 41(8): 1947-1949.

    [11] 李存玉, 支興蕾, 李賀敏, 等. 生物堿類成分的納濾傳質(zhì)行為與分子結(jié)構(gòu)的相關(guān)性 [J]. 藥學學報, 2018, 53(12): 1963-1967.

    [12] Cabrera M I, Luna J A, Grau R J A. Modeling of dissolution-diffusion controlled drug release from planar polymeric systems with finite dissolution rate and arbitrary drug loading [J]., 2006, 280(1/2): 693-704.

    [13] 彭九嫚, 柴堯, 潘旻, 等. 丹參中酚酸類成分在不同工藝條件下轉(zhuǎn)化關(guān)系 [J]. 中國實驗方劑學雜志, 2016, 22(5): 67-70.

    [14] 郝晨偉, 李正翔, 張銘慧, 等. 丹參及其配伍制劑治療冠心病的研究進展 [J]. 中草藥, 2021, 52(13): 4096- 4106.

    [15] Tai Y N, Shen J C, Luo Y,. Research progress on the ethanol precipitation process of traditional Chinese medicine [J]., 2020, 15: 84(2020).

    [16] 徐芳芳, 畢宇安, 王振中, 等. 過程分析技術(shù)在中藥注射劑生產(chǎn)過程中的應用研究進展 [J]. 中草藥, 2016, 47(14): 2563-2567.

    [17] Gong X C, Wang S S, Li Y,. Separation characteristics of ethanol precipitation for the purification of the water extract of medicinal plants [J]., 2013, 107: 273-280.

    [18] Gong X C, Li Y, Qu H B. Removing tannins from medicinal plant extracts using an alkaline ethanol precipitation process: A case study of Danshen Injection [J]., 2014, 19(11):18705-18720.

    Alcohol-precipitation mechanism of Danshen Injection based on molecular state quantitative calculation model

    LI Cun-yu1, 2, 3, LIN Ya-juan4, LI Ming-ming5, ZHANG Lian1, LI Shuo1, PENG Guo-ping1, 2, 3, ZHI Xing-lei1

    1. School of Pharmacy, Nanjing University of Chinese Medicine, Nanjing 210023, China 2. Jiangsu Province Engineering Research Center of Classical Prescription, Nanjing 210023, China 3. Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization, Nanjing 210023, China 4. Affiliated Hospital of Nanjing University of Chinese Medicine, Nanjing 210023, China 5. Jiangsu Shenlong Pharmaceutical Co., Ltd., Dongtai 224200, China

    To explore the transfer rules of phenolic acids in the alcohol precipitation process of Danshen Injection (DI) based on molecular state quantitative calculation model.In the experiment, sodium danshensu, protocatechuic aldehyde, rosmarinic acid and salvianolic acid B were selected as indexes for exploring the alcohol-precipitation mechanism. Based on the mathematical model of nanofiltration mass transfer, membrane flux and rejection were collected to calculate the solute mass transfer coefficients, and the quantitative calculation model of molecular state was established with molecular monomer components as reference of DI. Furthermore, the correlation among the transfer rate, molecular state and molecular weight was analyzed with the densities of 1.05, 1.15, 1.25 g/mL of concentrated solution under the alcohol precipitation of 75% and 85%.The power function equation of molecular state quantitative calculation was set up successfully with the regression coefficients (2> 0.99), as the density of intermediate in DI was 1.05—1.25 g/mL, the molecular proportions of sodium danshensu, protocatechuic aldehyde, rosmarinic acid and salvianolic acid B were 1.14%—5.25%, 66.70%—80.25%, 10.51%—19.11% and 1.13%—11.44%, respectively. The correlation analysis of “density- molecular state ratio- alcohol precipitation transfer rate” showed that alcohol precipitation concentration increased from 75% to 85%, and the decreasing rate of transfer rate was in the order of sodium danshensu > salvianolic acid B > rosmarinic acid > protocatechualdehyde. Therefore, the molecular state of phenolic acids was the main factor which affected the transfer rate of alcohol precipitation, and the relative molecular weight was correlated with the transfer rate.The quantitative calculation method of molecular state was established for phenolic acids in the intermediate of DI, and the refining mechanism of alcohol precipitation was preliminarily expounded which was helpful for standardized control of preparation alcohol precipitation process.

    molecular state;et; Danshen Injection; alcohol-precipitation; sodium danshensu; protocatechuic aldehyde; rosmarinic acid; salvianolic acid B; nanofiltration; power function equation

    R283.6

    A

    0253 - 2670(2022)20 - 6431 - 12

    10.7501/j.issn.0253-2670.2022.20.015

    2022-05-09

    國家自然科學基金資助項目(82274106);江蘇省自然科學基金面上項目(BK20211303);江蘇省333高層次人才培養(yǎng)工程;2021年度康緣中藥學院創(chuàng)新創(chuàng)業(yè)項目(kyxysc06);2022年國家級大學生創(chuàng)新創(chuàng)業(yè)訓練計劃項目(103152022105)

    李存玉(1985—),男,博士,副教授,碩士生導師,研究方向為膜分離原理及產(chǎn)業(yè)應用。Tel: (025)86798186 E-mail: 300632@njucm.edu.cn

    支興蕾(1981—),女,博士,副教授,碩士生導師,主要從事中藥皂苷類超濾分離機制研究。Tel: (025)86798186 E-mail: 300040@njucm.edu.cn

    [責任編輯 鄭禮勝]

    猜你喜歡
    兒茶素鈉冪函數(shù)
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(2)
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(1)
    冪函數(shù)、指數(shù)函數(shù)、對數(shù)函數(shù)(1)
    挨打的蜜峰
    熒光分光光度法測定水中示蹤劑熒光素鈉示蹤劑
    兒茶治好了我的鼻炎
    看圖說話,揭開冪函數(shù)的廬山真面目
    收濕斂瘡的兒茶
    復方兒茶酊的生藥學鑒定及質(zhì)量控制
    中成藥(2016年8期)2016-05-17 06:08:39
    貝前列素鈉對全腦缺血/再灌注大鼠腦損傷的作用及機制
    窝窝影院91人妻| 国产精品久久久人人做人人爽| 日本精品一区二区三区蜜桃| 香蕉国产在线看| 午夜成年电影在线免费观看| 欧美性长视频在线观看| 色在线成人网| 成人亚洲精品一区在线观看| 亚洲熟女精品中文字幕| 高清在线国产一区| 欧美精品亚洲一区二区| av欧美777| 欧美日韩乱码在线| 久久精品国产99精品国产亚洲性色 | 亚洲欧美日韩高清在线视频| 亚洲熟女毛片儿| 啦啦啦在线免费观看视频4| 超碰成人久久| 午夜福利视频在线观看免费| 亚洲第一青青草原| 一级,二级,三级黄色视频| 91麻豆av在线| 无限看片的www在线观看| 人人妻,人人澡人人爽秒播| 国产精品九九99| 一区在线观看完整版| 精品第一国产精品| 成人特级黄色片久久久久久久| 午夜福利一区二区在线看| 精品国产一区二区三区久久久樱花| 一级作爱视频免费观看| 中文字幕av电影在线播放| 老司机靠b影院| 亚洲男人天堂网一区| 在线天堂中文资源库| а√天堂www在线а√下载 | 在线永久观看黄色视频| 国产精品 国内视频| 成熟少妇高潮喷水视频| 亚洲专区国产一区二区| 一区二区日韩欧美中文字幕| 高清毛片免费观看视频网站 | 亚洲 欧美一区二区三区| 一二三四在线观看免费中文在| 久久狼人影院| 男女下面插进去视频免费观看| 自拍欧美九色日韩亚洲蝌蚪91| netflix在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| xxxhd国产人妻xxx| 精品电影一区二区在线| 极品教师在线免费播放| 色婷婷av一区二区三区视频| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久人人人人人| 国产xxxxx性猛交| 国产av精品麻豆| 99久久精品国产亚洲精品| 亚洲中文日韩欧美视频| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| 91麻豆av在线| 午夜两性在线视频| 黄色视频不卡| 日本一区二区免费在线视频| 欧美乱色亚洲激情| 久久狼人影院| 日韩大码丰满熟妇| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲高清精品| 大陆偷拍与自拍| av片东京热男人的天堂| 精品卡一卡二卡四卡免费| av天堂在线播放| 国产男女超爽视频在线观看| 一级片'在线观看视频| 美女国产高潮福利片在线看| 黄片小视频在线播放| 国产精品偷伦视频观看了| 人人澡人人妻人| 欧美日韩黄片免| 天堂√8在线中文| 午夜精品在线福利| 久久精品国产综合久久久| 免费看a级黄色片| 日韩欧美免费精品| svipshipincom国产片| 亚洲专区国产一区二区| 成年人免费黄色播放视频| 欧美日韩一级在线毛片| 国产亚洲精品一区二区www | 黄色毛片三级朝国网站| 欧美精品高潮呻吟av久久| 亚洲av第一区精品v没综合| 色综合欧美亚洲国产小说| 成年人午夜在线观看视频| av欧美777| 又大又爽又粗| 美女福利国产在线| 亚洲欧美色中文字幕在线| 亚洲精品中文字幕一二三四区| 亚洲国产中文字幕在线视频| 黄色 视频免费看| av有码第一页| 国产精品亚洲av一区麻豆| 色婷婷久久久亚洲欧美| 国产精品免费一区二区三区在线 | 又黄又爽又免费观看的视频| 99精品欧美一区二区三区四区| 欧美人与性动交α欧美精品济南到| 免费一级毛片在线播放高清视频 | 国产成人啪精品午夜网站| 久久久国产一区二区| 精品久久久久久久毛片微露脸| bbb黄色大片| 91老司机精品| 中文亚洲av片在线观看爽 | 久久久国产成人免费| 大片电影免费在线观看免费| 热re99久久国产66热| 一进一出抽搐gif免费好疼 | 日本a在线网址| 后天国语完整版免费观看| 欧美黄色淫秽网站| 日韩免费高清中文字幕av| 亚洲成人免费av在线播放| 国产一区二区激情短视频| 久久精品国产亚洲av高清一级| 少妇的丰满在线观看| 欧美性长视频在线观看| 欧美av亚洲av综合av国产av| 亚洲精品美女久久av网站| 麻豆av在线久日| 亚洲avbb在线观看| 欧美 亚洲 国产 日韩一| 亚洲av成人不卡在线观看播放网| 日韩欧美三级三区| 日本a在线网址| 99精品久久久久人妻精品| svipshipincom国产片| 九色亚洲精品在线播放| 国产淫语在线视频| 亚洲av电影在线进入| 久久精品亚洲av国产电影网| 精品午夜福利视频在线观看一区| 午夜久久久在线观看| 国产成+人综合+亚洲专区| 亚洲全国av大片| 欧美老熟妇乱子伦牲交| 狠狠婷婷综合久久久久久88av| 后天国语完整版免费观看| 国产激情久久老熟女| 好看av亚洲va欧美ⅴa在| 亚洲专区中文字幕在线| 一进一出好大好爽视频| 天天操日日干夜夜撸| 免费不卡黄色视频| 亚洲色图综合在线观看| 建设人人有责人人尽责人人享有的| 亚洲精品久久成人aⅴ小说| 青草久久国产| 高清欧美精品videossex| 午夜福利一区二区在线看| 如日韩欧美国产精品一区二区三区| 亚洲av成人一区二区三| 另类亚洲欧美激情| 波多野结衣一区麻豆| 男男h啪啪无遮挡| 国产成人系列免费观看| a级毛片在线看网站| 亚洲在线自拍视频| 男女午夜视频在线观看| 视频在线观看一区二区三区| 美国免费a级毛片| av视频免费观看在线观看| 欧美人与性动交α欧美精品济南到| av电影中文网址| 少妇被粗大的猛进出69影院| 超碰97精品在线观看| 亚洲人成77777在线视频| 亚洲avbb在线观看| av天堂在线播放| 色精品久久人妻99蜜桃| 搡老熟女国产l中国老女人| 精品人妻在线不人妻| 久9热在线精品视频| 一本综合久久免费| 日本黄色日本黄色录像| 国产精品欧美亚洲77777| 欧美 日韩 精品 国产| 少妇裸体淫交视频免费看高清 | 久久国产精品人妻蜜桃| av中文乱码字幕在线| 欧美日韩中文字幕国产精品一区二区三区 | 欧美人与性动交α欧美精品济南到| √禁漫天堂资源中文www| 久久久久国产精品人妻aⅴ院 | 亚洲欧美激情综合另类| 高潮久久久久久久久久久不卡| 久久久精品区二区三区| 另类亚洲欧美激情| 色精品久久人妻99蜜桃| 日韩欧美在线二视频 | 亚洲av熟女| 国产精品秋霞免费鲁丝片| 免费观看精品视频网站| 老司机靠b影院| 99精国产麻豆久久婷婷| 50天的宝宝边吃奶边哭怎么回事| 亚洲片人在线观看| 操美女的视频在线观看| 少妇粗大呻吟视频| 亚洲国产欧美一区二区综合| 韩国av一区二区三区四区| a级片在线免费高清观看视频| 精品久久久久久电影网| 国产一区二区三区视频了| 国产亚洲精品久久久久5区| 淫妇啪啪啪对白视频| 欧美av亚洲av综合av国产av| 热99re8久久精品国产| 水蜜桃什么品种好| 天堂中文最新版在线下载| 黄色a级毛片大全视频| 一级a爱片免费观看的视频| 免费观看a级毛片全部| 欧美精品啪啪一区二区三区| 露出奶头的视频| 欧美精品av麻豆av| 午夜福利欧美成人| 久久久久久亚洲精品国产蜜桃av| 成年版毛片免费区| 18禁国产床啪视频网站| 97人妻天天添夜夜摸| av天堂在线播放| 天天添夜夜摸| 国产亚洲欧美精品永久| 一边摸一边抽搐一进一出视频| 国产真人三级小视频在线观看| 国产成人精品久久二区二区91| www.精华液| 精品少妇一区二区三区视频日本电影| 1024香蕉在线观看| tocl精华| 美女福利国产在线| 国产日韩欧美亚洲二区| 中文字幕人妻丝袜制服| 国产男女内射视频| 亚洲专区国产一区二区| 亚洲成人国产一区在线观看| 国产欧美日韩一区二区精品| 免费一级毛片在线播放高清视频 | 日本wwww免费看| 欧美日韩精品网址| 日韩欧美免费精品| 精品国产国语对白av| 精品乱码久久久久久99久播| 精品欧美一区二区三区在线| 午夜免费成人在线视频| 无人区码免费观看不卡| 超碰97精品在线观看| 69av精品久久久久久| 热re99久久国产66热| 99热只有精品国产| 欧美精品av麻豆av| 黄色女人牲交| www.精华液| 涩涩av久久男人的天堂| 无人区码免费观看不卡| ponron亚洲| 高清毛片免费观看视频网站 | 久久久精品区二区三区| 在线观看午夜福利视频| 97人妻天天添夜夜摸| 亚洲九九香蕉| 午夜日韩欧美国产| 视频区欧美日本亚洲| xxx96com| 国产主播在线观看一区二区| 熟女少妇亚洲综合色aaa.| 久久人人97超碰香蕉20202| 黄片大片在线免费观看| 精品少妇一区二区三区视频日本电影| 亚洲国产欧美一区二区综合| 9热在线视频观看99| 极品教师在线免费播放| 欧美日韩国产mv在线观看视频| 一二三四社区在线视频社区8| 欧美中文综合在线视频| 精品国产一区二区三区四区第35| 国产单亲对白刺激| 精品亚洲成a人片在线观看| 精品熟女少妇八av免费久了| 国产成人欧美| 制服诱惑二区| 国产乱人伦免费视频| 天堂动漫精品| tocl精华| 亚洲成人手机| 丁香欧美五月| 久久精品91无色码中文字幕| xxxhd国产人妻xxx| 精品亚洲成a人片在线观看| 亚洲黑人精品在线| xxxhd国产人妻xxx| 精品亚洲成国产av| 最近最新免费中文字幕在线| 窝窝影院91人妻| 国产成人系列免费观看| 又大又爽又粗| 精品卡一卡二卡四卡免费| 精品一区二区三卡| 精品福利永久在线观看| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区蜜桃| 高清黄色对白视频在线免费看| 欧美久久黑人一区二区| 国产亚洲精品第一综合不卡| 99久久国产精品久久久| 美女视频免费永久观看网站| 欧美成狂野欧美在线观看| 久热爱精品视频在线9| 欧美在线一区亚洲| 激情在线观看视频在线高清 | 黄网站色视频无遮挡免费观看| 久久久久久久国产电影| 国产欧美日韩综合在线一区二区| ponron亚洲| 亚洲色图 男人天堂 中文字幕| 女同久久另类99精品国产91| 这个男人来自地球电影免费观看| 一级黄色大片毛片| 十八禁人妻一区二区| 久久精品亚洲熟妇少妇任你| 最近最新中文字幕大全免费视频| 亚洲国产精品合色在线| 亚洲av欧美aⅴ国产| 99国产精品一区二区三区| 久久国产精品大桥未久av| 中文字幕av电影在线播放| 亚洲精品美女久久久久99蜜臀| 久久国产精品人妻蜜桃| 99精品欧美一区二区三区四区| 亚洲欧洲精品一区二区精品久久久| 久久精品熟女亚洲av麻豆精品| 老司机午夜福利在线观看视频| 色在线成人网| 国产成人免费观看mmmm| 在线视频色国产色| 大片电影免费在线观看免费| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| 国产成人系列免费观看| 欧美另类亚洲清纯唯美| 久久久精品免费免费高清| 欧美中文综合在线视频| 久久天堂一区二区三区四区| 久久人人爽av亚洲精品天堂| 欧美黑人精品巨大| 国产国语露脸激情在线看| 久热爱精品视频在线9| 亚洲精品中文字幕一二三四区| 亚洲av第一区精品v没综合| 国产亚洲精品久久久久5区| 一级毛片女人18水好多| 精品国产一区二区三区四区第35| 国产精品美女特级片免费视频播放器 | 欧美日韩国产mv在线观看视频| 免费在线观看完整版高清| 日本wwww免费看| 国产精品一区二区精品视频观看| 久久中文看片网| 午夜老司机福利片| 丰满迷人的少妇在线观看| 动漫黄色视频在线观看| 大陆偷拍与自拍| 久99久视频精品免费| 国产精品一区二区在线不卡| www.自偷自拍.com| 久久精品国产亚洲av高清一级| 亚洲国产中文字幕在线视频| 精品亚洲成国产av| 日韩一卡2卡3卡4卡2021年| 777久久人妻少妇嫩草av网站| 人人妻人人澡人人爽人人夜夜| 亚洲情色 制服丝袜| 十八禁高潮呻吟视频| 国产色视频综合| 国产精品一区二区精品视频观看| 国产欧美亚洲国产| 欧美日韩亚洲国产一区二区在线观看 | 每晚都被弄得嗷嗷叫到高潮| 免费人成视频x8x8入口观看| 日日爽夜夜爽网站| 亚洲av美国av| 亚洲欧洲精品一区二区精品久久久| 啦啦啦免费观看视频1| 在线av久久热| 精品亚洲成a人片在线观看| 99精品久久久久人妻精品| 亚洲自偷自拍图片 自拍| 国产av一区二区精品久久| 在线av久久热| 国产亚洲精品第一综合不卡| 99香蕉大伊视频| 一二三四在线观看免费中文在| 51午夜福利影视在线观看| 一级毛片精品| avwww免费| 人妻久久中文字幕网| 亚洲第一av免费看| tocl精华| 美女国产高潮福利片在线看| 日本a在线网址| 日日夜夜操网爽| 男女午夜视频在线观看| 又大又爽又粗| 高清欧美精品videossex| 国产一区二区三区综合在线观看| 老司机福利观看| 999久久久国产精品视频| 丝袜人妻中文字幕| 亚洲自偷自拍图片 自拍| 女性生殖器流出的白浆| 久久人人爽av亚洲精品天堂| 人妻 亚洲 视频| 色婷婷av一区二区三区视频| 免费在线观看影片大全网站| 视频在线观看一区二区三区| av不卡在线播放| 国产亚洲欧美在线一区二区| 久久午夜综合久久蜜桃| 久久狼人影院| 女人被狂操c到高潮| 亚洲全国av大片| 欧美精品一区二区免费开放| 又紧又爽又黄一区二区| 久久人妻熟女aⅴ| 午夜日韩欧美国产| 19禁男女啪啪无遮挡网站| 淫妇啪啪啪对白视频| 少妇被粗大的猛进出69影院| 成人影院久久| 美女高潮到喷水免费观看| 亚洲一区二区三区不卡视频| 久久亚洲真实| 婷婷成人精品国产| 夜夜夜夜夜久久久久| 啦啦啦免费观看视频1| 亚洲精品国产色婷婷电影| 亚洲精品国产精品久久久不卡| 91大片在线观看| 亚洲专区中文字幕在线| 麻豆国产av国片精品| 身体一侧抽搐| 在线观看免费日韩欧美大片| netflix在线观看网站| 一个人免费在线观看的高清视频| 一区二区三区国产精品乱码| 日韩三级视频一区二区三区| 亚洲欧美日韩高清在线视频| 可以免费在线观看a视频的电影网站| 国产精品综合久久久久久久免费 | 欧美中文综合在线视频| av线在线观看网站| 亚洲一区二区三区欧美精品| 久久精品亚洲熟妇少妇任你| 国产一区二区三区在线臀色熟女 | 精品卡一卡二卡四卡免费| 捣出白浆h1v1| 精品第一国产精品| av视频免费观看在线观看| 欧美一级毛片孕妇| 国产亚洲精品第一综合不卡| 人人妻人人澡人人看| 飞空精品影院首页| 在线观看免费日韩欧美大片| 一二三四在线观看免费中文在| 国产成人精品无人区| 99国产精品一区二区三区| 大陆偷拍与自拍| 看黄色毛片网站| 久久香蕉国产精品| 精品国内亚洲2022精品成人 | 别揉我奶头~嗯~啊~动态视频| 亚洲精品在线观看二区| 欧美国产精品一级二级三级| www.999成人在线观看| 黑人操中国人逼视频| 夜夜爽天天搞| 国产欧美日韩一区二区三区在线| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 成人永久免费在线观看视频| 超色免费av| 黑人巨大精品欧美一区二区蜜桃| 国产精品久久电影中文字幕 | 日韩欧美在线二视频 | 热99国产精品久久久久久7| 国产成人影院久久av| 91字幕亚洲| 高清av免费在线| 国产精品免费一区二区三区在线 | 亚洲av成人av| av不卡在线播放| 久久精品亚洲av国产电影网| 91成人精品电影| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区激情短视频| 岛国在线观看网站| 欧美日韩视频精品一区| tube8黄色片| 欧美久久黑人一区二区| 天天添夜夜摸| 国产精品亚洲一级av第二区| 很黄的视频免费| av视频免费观看在线观看| 真人做人爱边吃奶动态| 老鸭窝网址在线观看| 每晚都被弄得嗷嗷叫到高潮| 成人三级做爰电影| 久久久久精品人妻al黑| 不卡av一区二区三区| 久久性视频一级片| 天天躁狠狠躁夜夜躁狠狠躁| 精品一区二区三区av网在线观看| 精品国产乱码久久久久久男人| 国产亚洲精品久久久久久毛片 | 黄色女人牲交| 在线国产一区二区在线| 中亚洲国语对白在线视频| 中文字幕制服av| 成人国语在线视频| 在线观看一区二区三区激情| 人妻久久中文字幕网| 制服人妻中文乱码| 极品少妇高潮喷水抽搐| 操美女的视频在线观看| 波多野结衣av一区二区av| 精品国产一区二区久久| 丝袜在线中文字幕| 如日韩欧美国产精品一区二区三区| 999久久久精品免费观看国产| 亚洲人成电影观看| 丰满饥渴人妻一区二区三| 国产av又大| 美女视频免费永久观看网站| 丰满人妻熟妇乱又伦精品不卡| 一区二区日韩欧美中文字幕| 午夜成年电影在线免费观看| 他把我摸到了高潮在线观看| 欧美大码av| 午夜福利,免费看| 久久青草综合色| 久久久久久久午夜电影 | 国产精品电影一区二区三区 | 精品卡一卡二卡四卡免费| 中国美女看黄片| 国产日韩欧美亚洲二区| 91老司机精品| 久久久久久久精品吃奶| 嫩草影视91久久| 欧美精品av麻豆av| 国产亚洲精品久久久久5区| 身体一侧抽搐| 亚洲国产中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 又紧又爽又黄一区二区| 欧美黄色淫秽网站| 天堂√8在线中文| 久久精品国产综合久久久| 天堂中文最新版在线下载| 亚洲国产中文字幕在线视频| 在线观看免费视频日本深夜| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品古装| 波多野结衣一区麻豆| 999久久久国产精品视频| 成人黄色视频免费在线看| 九色亚洲精品在线播放| 亚洲av日韩精品久久久久久密| 精品久久久精品久久久| 69av精品久久久久久| 亚洲欧美激情在线| 操出白浆在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 成人18禁在线播放| 99re6热这里在线精品视频| 男男h啪啪无遮挡| 精品免费久久久久久久清纯 | av福利片在线| 极品少妇高潮喷水抽搐| 一级片'在线观看视频| 宅男免费午夜| 免费一级毛片在线播放高清视频 | 亚洲国产欧美日韩在线播放| x7x7x7水蜜桃| 国产又色又爽无遮挡免费看| xxxhd国产人妻xxx| 亚洲人成77777在线视频| 欧美日韩瑟瑟在线播放| 麻豆av在线久日| 欧美亚洲日本最大视频资源| 美女 人体艺术 gogo| 无遮挡黄片免费观看| 午夜久久久在线观看| 国产精品免费一区二区三区在线 | 身体一侧抽搐| 精品久久久久久电影网| 免费少妇av软件| 国产片内射在线| 99精品欧美一区二区三区四区| 国产成人精品久久二区二区免费| 亚洲av熟女| 999久久久精品免费观看国产|