亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        巨大芽孢桿菌ω-轉(zhuǎn)氨酶催化合成手性β-氨基醇

        2022-01-17 08:27:12楊曉曉張朝峰常宏宏張建棟
        高?;瘜W工程學報 2021年6期

        楊曉曉, 李 靜, 張朝峰, 常宏宏, 張建棟

        巨大芽孢桿菌-轉(zhuǎn)氨酶催化合成手性-氨基醇

        楊曉曉, 李 靜, 張朝峰, 常宏宏, 張建棟

        (太原理工大學 生物醫(yī)學工程學院, 山西 太原 030024)

        針對目前用于合成手性-氨基醇的()-轉(zhuǎn)氨酶稀缺的問題,對來自巨大芽孢桿菌(SC6394)的 ()--轉(zhuǎn)氨酶(BMTA)基因密碼子進行優(yōu)化,在大腸桿菌中成功表達,并對BMTA進行純化表征。結(jié)果表明:BMTA對()--氨基醇具有較高活力,對()--氨基醇沒有活力;BMTA在pH為7.5、溫度為55 ℃下酶活最高,在pH為6.5~8.0,溫度為4~30 ℃內(nèi)具有較好的穩(wěn)定性。應用BMTA對不同外消旋-氨基醇進行動力學拆分,底物轉(zhuǎn)化率高達50%,對映體過量值(ee)為50%~99%。應用BMTA不對稱還原胺化不同-羥基酮,產(chǎn)物()--氨基醇轉(zhuǎn)化率高達58%~96%,ee>99%。

        -轉(zhuǎn)氨酶;手性-氨基醇;羥酮;動力學拆分;不對稱還原胺化

        1 前言

        手性-氨基醇是一類非常重要的化合物,可以作為手性助劑或手性配體用于不對稱催化反應[1-2],也可作為手性中間體或砌塊用于合成多種手性藥物[3-4]?;瘜W法合成手性-氨基醇[5-7]普遍存在合成工藝復雜、反應條件苛刻、選擇性不高、催化劑昂貴等缺點。生物法作為一種綠色可持續(xù)方法[8],已成為國際研究熱點。目前,生物法合成手性-氨基醇主要包括轉(zhuǎn)氨酶動力學拆分和不對稱還原胺化[9]、脂肪酶動力學拆分[10]、胺脫氫酶不對稱還原[11]、酮還原酶不對稱還原[12],以及級聯(lián)生物催化環(huán)氧化物不對稱開環(huán)等[13]。轉(zhuǎn)氨酶作為一種高效催化劑在手性胺合成方面已經(jīng)得到廣泛應用[14],但關于手性-氨基醇尤其是()--氨基醇合成的報道較少。目前僅有來自紫色色桿菌()的CVTA被報道可以還原胺化羥酮合成()--氨基醇,但產(chǎn)物轉(zhuǎn)化率不到40%[15]。近年來,對來自河弧菌()[16]和魯杰氏菌(TM1040)[17]的轉(zhuǎn)氨酶(vfTA和3FCR)進行了定向進化,得到的突變酶被用于合成()--氨基醇,但反應效率仍有待提高。因此,篩選新的高活力高選擇性轉(zhuǎn)氨酶用于手性-氨基醇的合成具有重要的價值。

        本研究對來自巨大芽孢桿菌(SC6394)[18]的一個()-轉(zhuǎn)氨酶(BMTA)基因序列進行了密碼子優(yōu)化,在大腸桿菌(BL21)中進行了表達純化及酶學性質(zhì)表征;應用BMTA對外消旋-氨基醇進行了動力學拆分,得到()--氨基醇,同時考察了BMTA不對稱還原胺化-羥基酮合成()--氨基醇的效率。

        2 實驗部分

        2.1 材料與試劑

        菌種與質(zhì)粒:BL21感受態(tài)細胞購于天根生化公司,質(zhì)粒pET-28a(+)購于安諾倫生物科技有限公司。

        培養(yǎng)基(LB培養(yǎng)基):質(zhì)量濃度為10.0 g×L-1胰蛋白胨、5.0 g×L-1酵母提取物、10.0 g×L-1氯化鈉(NaCl),固體培養(yǎng)基添加質(zhì)量濃度為15.0 g×L-1瓊脂粉。

        Terrific Broth培養(yǎng)基(TB培養(yǎng)基):質(zhì)量濃度為12.0 g×L-1胰蛋白胨、24.0 g×L-1酵母提取物、2.31 g×L-1磷酸二氫鉀、12.54 g×L-1磷酸氫二鉀、體積分數(shù)為4 mL×L-1甘油。其中涉及的試劑均為生化試劑級別,購于生工生物工程(上海)有限公司。

        試劑:外消旋-氨基醇((±)-1)、()-1和()-1購于安耐吉化學有限公司和蘇州愛瑪特科技有限公司,-羥基酮(2)購于安耐吉化學有限公司和百靈威科技有限公司。本研究中所涉及的化學藥品均為分析純,其他化學試劑均可從市面上購買獲得。

        2.2 實驗方法

        2.2.1-轉(zhuǎn)氨酶重組菌的構(gòu)建

        -轉(zhuǎn)氨酶BMTA基因來自巨大芽孢桿菌(SC6394),對BMTA基因密碼子進行優(yōu)化并交公司合成(通用生物公司)。將目的基因與表達載體pET28a連接,獲得重組質(zhì)粒pET28a-BMTA。將重組質(zhì)粒pET28a-BMTA導入BL21感受態(tài)細胞中,獲得重組(BMTA)。

        2.2.2 BMTA在大腸桿菌中的表達及純化

        將(BMTA)接種于質(zhì)量濃度為50 μg×mL-1卡那霉素的LB液體培養(yǎng)基中,溫度為37 ℃、轉(zhuǎn)速為200 r×min-1培養(yǎng)7 h,取1 mL的培養(yǎng)液于50 mL的TB培養(yǎng)基(質(zhì)量濃度為50 μg×mL-1的卡那霉素)中進行擴大培養(yǎng)。當培養(yǎng)液在波長為600 nm的吸光值(OD600)達到0.5~0.6時,向培養(yǎng)基中加入濃度為0.5 mmol×L-1的異丙基--D-硫代半乳糖苷(isopropy--D-thiogalactoside,IPTG),20 ℃、200 r×min-1的條件下繼續(xù)培養(yǎng)12~14 h,離心收集菌體,將菌體用磷酸緩沖液(濃度為100 mmol×L-1,pH=8.0)洗滌2次,最后用磷酸緩沖液將細胞懸浮,置于冰上,400 W下進行超聲破碎10 min(工作時間4 s,間歇時間4 s),將破碎液在溫度為4 ℃、轉(zhuǎn)速為3 000 r×min-1的條件下離心30 min,獲得BMTA的粗酶液,通過鎳柱進行純化,SDS-PAGE分析純化后的蛋白。

        2.2.3 BMTA的酶活力檢測及酶活單位定義

        1 mL反應液中有濃度為100 mmol×L-1磷酸鈉緩沖液(pH=8.0)、10 mmol×L-1底物、10 mmol×L-1丙酮酸鈉、0.2 mmol×L-15'-磷酸吡哆醛(pyridoxal 5'-phosphate,PLP)和0.1 mL酶液,30 ℃反應10 min后,取樣300 μL,NaCl飽和,加入10 μL氫氧化鈉(濃度為10 mol×L-1)調(diào)節(jié)溶液pH>10.0,加入300 μL乙酸乙酯(含內(nèi)標物20 mmol×L-1十二烷)萃取,有機相加無水硫酸鈉干燥后,氣相檢測生成-羥基酮的量。1 min內(nèi)催化1 μmol產(chǎn)物生成所需要的酶量定義為1個酶活單位(U)。以牛血清白蛋白(BSA)為標準,采用Bradford法測定酶液中蛋白質(zhì)的濃度[19]。

        2.2.4 pH和溫度對BMTA酶活性的影響

        在不同pH值的緩沖溶液中(磷酸鈉緩沖液(100 mmol×L-1,pH為6.0~8.0)、Tris×HCl緩沖液(100 mmol×L-1,pH為8.0~9.0)和甘氨酸-氫氧化鈉緩沖液(50 mmol×L-1,pH為9.0~10.0))測定BMTA對()-1c酶活。在最佳pH條件下,在不同溫度(25~65 ℃)下檢測BMTA對()-1c的酶活(參照2.2.3節(jié))。

        2.2.5 pH和溫度對BMTA穩(wěn)定性的影響

        將BMTA置于2.2.4節(jié)不同pH的緩沖液中,于4 ℃下保存,不同時間取樣檢測該酶的殘余酶活。將BMTA置于磷酸鈉緩沖液(100 mmol×L-1,pH=7.5)中,分別在4~50 ℃的水浴中孵育,不同時間取樣檢測該酶的殘余酶活(參照2.2.3節(jié))。

        2.2.6 BMTA動力學拆分外消旋-氨基醇

        5 mL反應體系中有濃度為100 mmol×L-1磷酸鈉緩沖液(pH=7.5)、5~20 mmol×L-1底物((±)-1a-l)、0.2 mmol×L-1PLP、5~20 mmol×L-1丙酮酸鈉、質(zhì)量濃度為10 mg×mL-1BMTA,在30 ℃、200 r×min-1的條件下反應12~24 h,取樣300 μL,NaCl飽和,加入10 μL NaOH (10 mol×L-1)調(diào)節(jié)溶液pH>10.0,加入300 μL乙酸乙酯(含內(nèi)標物20 mmol×L-1十二烷)萃取,有機相加無水硫酸鈉干燥后,氣相檢測生成-羥基酮的量以及底物的對映體過量百分數(shù)(electrical engineering, ee),反應路線如圖1所示。

        圖1 ω-轉(zhuǎn)氨酶BMTA拆分外消旋β-氨基醇

        2.2.7 BMTA不對稱還原胺化-羥基酮

        5 mL反應體系中有濃度為100 mmol×L-1磷酸鈉緩沖液(pH=7.5)、10~20 mmol×L-1-羥基酮底物(2a-f)、0.2 mmol×L-1PLP、200~400 mmol×L-1L-丙氨酸、質(zhì)量濃度為10 mg×mL-1BMTA,在30 ℃、200 r×min-1的條件下反應10~20 h,取樣300 μL,NaCl飽和,加入10 μL NaOH(10 mol×L-1)調(diào)節(jié)溶液pH>10.0,加入300 μL乙酸乙酯(含內(nèi)標物20 mmol×L-1十二烷)萃取,有機相加無水硫酸鈉干燥后,氣相檢測生成-氨基醇的量以及產(chǎn)物的ee值,反應路線如圖2所示。

        圖2 ω-轉(zhuǎn)氨酶BMTA催化α-羥基酮不對稱還原胺化

        2.2.8 BMTA動力學拆分(±)-1c制備()-1c

        100 mL反應體系中有濃度為100 mmol×L-1磷酸鈉緩沖液(pH 7.5)、20 mmol×L-1底物((±)-1c)、0.2 mmol×L-1PLP、20 mmol×L-1丙酮酸鈉、質(zhì)量濃度為10 mg×mL-1BMTA,在30 ℃、200 r×min-1的條件下反應24 h,反應液中加NaCl飽和,加入鹽酸使溶液pH<2,100 mL乙酸乙酯連續(xù)萃取3次去除所提取的-羥基酮。加入NaOH溶液(10 mol×L-1)使溶液pH>10,100 mL乙酸乙酯重復萃取3次,有機相加無水硫酸鈉干燥,過濾后旋轉(zhuǎn)蒸發(fā)除去乙酸乙酯,得到的產(chǎn)物真空干燥過夜,計算產(chǎn)物得率。

        2.2.9 分析方法

        使用氣相色譜分析轉(zhuǎn)化率和ee值,具體方法按照Zhang等的方法進行[9]。

        3 結(jié)果與討論

        3.1 BMTA的表達、純化和酶活檢測

        如圖3(a)所示,在蛋白質(zhì)相對分子量r=53 000處出現(xiàn)與BMTA理論大小一致的目的條帶,證明BMTA在大腸桿菌中成功表達,并且大多為可溶性表達(如圖3(a)所示,lane 2),僅產(chǎn)生少量包涵體蛋白(圖3(a),lane 3)。粗酶液經(jīng)鎳柱純化后,得到單一的目的蛋白條帶(圖3(b),lanes 7~8)。對純酶酶活分析發(fā)現(xiàn),BMTA對()-1c有較高酶活(1.1 U×mg-1),對()-1c無活性。Nobili等[16]通過對來自的轉(zhuǎn)氨酶vfTA進行蛋白質(zhì)改造來提高其催化活力,對()-苯甘氨醇的催化活力僅為0.53 U×mg-1。BMTA作為一種()-選擇性轉(zhuǎn)氨酶,可立體選擇性轉(zhuǎn)化()--氨基醇,這一現(xiàn)象可用Cahn-Ingold-Prelog規(guī)則來解釋[17]。

        圖3 ω-轉(zhuǎn)氨酶BMTA表達及純化的SDS-PAGE電泳圖

        3.2 pH對BMTA活性及穩(wěn)定性的影響

        如圖4(a)所示,BMTA的最合適pH=7.5,當pH=6.5~8.0時,BMTA相對酶活穩(wěn)定在80% 以上,當pH>8.5或者pH<6.5時,酶活迅速下降。如圖4(b)所示,BMTA在pH=7.0的條件下較穩(wěn)定,該條件下孵育24 h,殘余酶活仍保持在74% 以上;在pH=7.5的條件下孵育24 h,殘余酶活保持在65% 以上。當pH>8.5或pH<6.5時,BMTA酶活損失較大。BMTA在pH=7.0~8.0較穩(wěn)定,而多數(shù)已報道的轉(zhuǎn)氨酶的最適pH偏高,如Shin等[20]報道的來自JS17的轉(zhuǎn)氨酶vfTA的最適pH=9.2,當pH>9.0時,酶穩(wěn)定性迅速下降。

        圖4 pH對ω-轉(zhuǎn)氨酶BMTA活性及穩(wěn)定性的影響

        3.3 溫度對BMTA活性及穩(wěn)定性的影響

        如圖5(a)所示,BMTA的最合適反應溫度為55 ℃,高于60 ℃或低于50 ℃條件下,酶活迅速下降。BMTA溫度穩(wěn)定性結(jié)果如圖5(b)所示,BMTA在溫度為30 ℃下孵育24 h,殘余酶活仍保持在50% 以上。在50 ℃下孵育9 h后酶活基本喪失,顯示該酶在低于40 ℃溫度下具有較好的穩(wěn)定性。該結(jié)果與已報道的大多數(shù)轉(zhuǎn)氨酶相似,如來自結(jié)核分枝桿菌()的()-轉(zhuǎn)氨酶MVTA[9],當溫度高于40 ℃后,該酶酶活迅速下降。

        圖5 溫度對ω-轉(zhuǎn)氨酶BMTA活性及穩(wěn)定性的影響

        3.4 BMTA的底物特異性

        BMTA對不同-氨基醇的酶活如表1所示,對鏈狀-氨基醇()-1a的活性較低(0.152 U×mg-1),而對含支鏈的鏈狀-氨基醇()-1b活性較高(1.083 U×mg-1);BMTA對含苯環(huán)-氨基醇()-1c的活性最高(1.107 U×mg-1),對苯環(huán)對位含有鹵素取代基的底物(()-1d-f)的活性呈下降趨勢,為0.658~0.353 U×mg-1;BMTA對苯環(huán)對位含─CF3取代基的底物(()-1h)有較高活性(0.624 U×mg-1);當?shù)孜锏谋江h(huán)間位被鹵素取代后(()-1j-k),BMTA的活性也呈現(xiàn)下降趨勢,為0.758~0.346 U×mg-1。BMTA對所有()-底物(1a-l)沒有活性,可見BMTA對-氨基醇具有極高的立體選擇性。

        表1 BMTA對不同底物的酶活性

        3.5 BMTA動力學拆分外消旋β-氨基醇

        如表2所示,BMTA對底物(±)-1a活性較低,反應24 h后,底物轉(zhuǎn)化率僅有31.5%,ee值為50.1%;對其他底物(±)-1b-l,BMTA的活性較高,反應12~24 h,底物完全被拆分,底物轉(zhuǎn)化率為50%左右,ee>98%。

        表2 BMTA拆分外消旋β-氨基醇

        3.6 BMTA不對稱還原胺化α-羥基酮

        如表3所示,BMTA對6種-羥基酮進行不對稱還原胺化反應。反應10 h,底物2b-2d的轉(zhuǎn)化率可達93% 以上,而底物2e轉(zhuǎn)化率可達85.3%,產(chǎn)物ee>99%。進一步提高底物2c濃度至20 mmol×L-1,反應20 h,底物轉(zhuǎn)化率仍可達95.2%。BMTA對2a和2f活力較低,反應10 h,轉(zhuǎn)化率分別為58% 和60%。與BMTA相比,目前已報道轉(zhuǎn)氨酶催化-羥基酮還原胺化活力普遍較低,如Nobili等[16]通過對來自的轉(zhuǎn)氨酶vfTA進行蛋白質(zhì)改造來提高其催化活力,但產(chǎn)物()-苯甘氨醇得率僅有60%。

        表3 BMTA不對稱還原胺化α-羥基酮

        3.7 BMTA動力學拆分(±)-1c制備(S)-1c

        在100 mL的反應體系中對20 mmol×L-1(±)-1c (275 mg)進行動力學拆分,反應24 h,轉(zhuǎn)化率可達50%。經(jīng)純化后獲得110.4 mg ()-1c,得率為40.2%,純度>98%,ee>99%。

        4 結(jié)論

        對來自巨大芽孢桿菌(SC6394)的一個()-選擇性-轉(zhuǎn)氨酶BMTA基因密碼子進行優(yōu)化,并成功在大腸桿菌中進行高效表達,對酶進行了純化和表征。研究發(fā)現(xiàn),在溫度為30 ℃、pH=7.5的條件下,且在PLP和丙酮酸鈉存在時,BMTA可以拆分外消旋-氨基醇合成()--氨基醇;在30 ℃、pH=7.5的條件下,且在PLP和L-丙氨酸存在時,BMTA可不對稱還原胺化-羥基酮合成()--氨基醇。此外,在100 mL的反應體系中對外消旋苯甘氨醇進行了動力學拆分,成功制備了()-1c,得率高達40.2%,ee>99%。本研究證明了BMTA對-氨基醇類化合物具有較高的活力和極好的立體選擇性,對手性β-氨基醇類藥物中間體工業(yè)化生產(chǎn)具有潛在的應用價值。

        [1] AGER D J, SCHAAD D R, PRAKASH I.1,2-Amino alcohols and their heterocyclic derivatives as chiral auxiliaries in asymmetric synthesis [J]. Chemical Reviews, 1996, 27(2): 835-876.

        [2] MAO P, YANG L R, XIAO Y M,. Synthesis of chiral imidazole amino alcohols and their application in the asymmetric copper-catalyzed henry reaction [J]. Chinese Journal of Organic Chemistry, 2019, 39(2): 443-448.

        [3] RADL S, STACH J, PISA O,. An improved synthesis of elvitegravir [J]. Journal of Heterocyclic Chemistry, 2016, 53(6): 1738-1749.

        [4] GUO J, ZHU M Y, WU T X,. Discovery of indolin-2-one derivatives as potent PAK4 inhibitors: Structure-activity relationship analysis, biological evaluation and molecular docking study [J]. Bioorganic & Medicinal Chemistry, 2017, 25(13): 3500-3511.

        [5] 李鑫磊, 許孝良, 李小年.-氨基酸還原制備-氨基醇研究進展[J]. 現(xiàn)代化工, 2017, 37(4): 28-32.

        LI X L, XU X L, LI X N. Progress of the preparation of-amino alcohols by reduction of-amino acid [J]. Modern Chemical Industry,2017, 37(4): 28-32.

        [6] BODKIN J A, MCLEOD M D. The sharpless asymmetric aminohydroxylation [J]. Journal of the Chemical Society Perkin Transactions, 2002, 34(24): 2733-2746.

        [7] WANG C, LUO L, YAMAMOTO H. Metal-catalyzed directed regio- and enantioselective ring-opening of epoxides [J]. Accounts of Chemical Research,2016, 49(2): 193-204.

        [8] 許娜, 關今韜, 張頌紅, 等. 晶膠固載副干酪乳桿菌及合成苯乳酸的實驗研究[J]. 高?;瘜W工程學報, 2020, 34(1): 149-154.

        XU N, GUAN J T, ZHANG S H,. Preparation ofcell-immobilized cryogel biocatalysts for biosynthesis of phenyllactic acid [J]. Journal of Chemical Engineering of Chinese Universities, 2020, 34(1): 149-154.

        [9] ZHANG J D, ZHAO J W, GAO L L,. Enantioselective synthesis of enantiopure-amino alcohols via kinetic resolution and asymmetric reductive amination by a robust transaminase from[J]. Journal of Biotechnology, 2019, 290: 24-32.

        [10] ROUF A, GUPTA P, AGA M A,. Cyclic--amino alcohols: Preparation and enzymatic kinetic resolution [J]. Tetrahedron Asymmetry, 2011, 22(24): 2134-2143.

        [11] WANG H, QU G, LI J K,. Data mining of amine dehydrogenases for the synthesis of enantiopure amino alcohols [J]. Catalysis Science & Technology, 2020, 10(17): 5945-5952.

        [12] PATEL R N, BANERJEE A, HOWELL J M,. Microbial synthesis of (2,3)-(-)-N-benzoyl-3-phenyl isoserine ethyl ester-a taxol side-chain synthon [J].Tetrahedron: Asymmetry, 1993, 4(9): 2069-2084.

        [13] ZHANG J D, YANG X X, JIA Q,. Asymmetric ring opening of racemic epoxides for enantioselective synthesis of ()--amino alcohols by a cofactor self-sufficient cascade biocatalysis system [J]. Catalysis Science & Technology, 2019, 9(1): 70-74.

        [14] FUCHS M, FARNBERGER J E, KROUTIL W. The industrial age of biocatalytic transamination [J]. European Journal of Organic Chemistry, 2015, 2015(32): 6965-6982.

        [15] SEHL T, SIMON R C, HAILES H C,. TTC-based screening assay for-transaminases: A rapid method to detect reduction of 2-hydroxy ketones [J]. Journal of Biotechnology, 2012, 159(3):188-194.

        [16] NOBILI A, STEFFEN-MUNSBERG F, KOHLS H,. Engineering the active site of the amine transaminase fromfor the asymmetric synthesis of aryl-alkyl amines and amino alcohols [J]. ChemCatChem, 2015, 7(5): 757-760.

        [17] PAVLIDIS I V, WEISS M S, GENZ M,. Identification of ()-selective transaminases for the asymmetric synthesis of bulky chiral amines [J]. Nature Chemistry, 2016, 8(11): 1076-1082.

        [18] HANSON R L, DAVIS B L, CHEN Y,. Preparation of ()-amines from racemic amines with an ()-amine transaminase from[J]. Advanced Synthesis & Catalysis, 2008, 350(9): 1367-1375.

        [19] BRADFORD M M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding [J]. Analytical Biochemistry, 1976, 72(1/2): 248-254.

        [20] SHIN J S, YUN H, JANG J W,. Purification, characterization, and molecular cloning of a novel amine: Pyruvate transaminase fromJS17 [J]. Applied Microbiology and Biotechnology, 2003, 61(5): 463-471.

        Synthesis of chiral-amino alcohols by-transaminases from

        YANG Xiao-xiao, LI Jing, ZHANG Chao-feng, CHANG Hong-hong, ZHANG Jian-dong

        (College of Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030024, China)

        ()-transaminases available for chiral-amino alcohol synthesis are limited. A ()--transaminase (BMTA) gene (fromSC6394) codons was optimized and overexpressed in, and it was purified and characterized. The results show that the BMTA had high activity toward ()--amino alcohols, but no activity toward ()--amino alcohols. The maximum activity of BMTA was detected at pH 7.5 and 55 ℃, and it had good stability at pHs from 6.5 to 8.0 and temperatures from 4 ℃ to 30 ℃. Kinetic resolution of a set of racemic-amino alcohols by BMTA resulted in 50% conversion of substrates with 50%-99% ee of ()--amino alcohols. Asymmetric reductive amination of several-hydroxy ketones by BMTA resulted in 58%-96% conversions and >99% ee of ()--amino alcohols.

        -transaminase; chiral-amino alcohols;-hydroxy ketones; kinetic resolution; asymmetric reduction amination

        1003-9015(2021)06-1035-06

        Q555.2

        A

        10.3969/j.issn.1003-9015.2021.06.011

        2020-10-16;

        2020-12-29。

        國家自然科學基金(21772141);山西省青年基金 (201701D221042)。

        楊曉曉(1996-),女,山西臨汾人,太原理工大學碩士生。

        張建棟,E-mail:zhangjiandong@tyut.edu.cn

        国产v综合v亚洲欧美大天堂| 中国人妻与老外黑人| 婷婷五月六月综合缴情| 欧美一欧美一区二三区性| 国产亚洲一区二区三区夜夜骚| 亚洲日本一区二区在线| 久久久亚洲欧洲日产国码aⅴ | 伊人一道本| 久久国产劲爆内射日本| 亚洲av综合av一区| 国产无遮挡又黄又爽免费网站 | 国内精品久久久久国产盗摄 | 免费国产自拍视频在线观看| 男人的天堂一区二av| 欧美人妻日韩精品| 亚洲中文无码久久精品1| 亚洲黄色官网在线观看| 日韩av在线播放人妻| 最近在线更新8中文字幕免费| 久久久国产精品樱花网站| 日本岛国一区二区三区| 亚洲综合一区中文字幕| 99精品欧美一区二区三区| 加勒比精品久久一区二区三区| 91老司机精品视频| 亚洲最黄视频一区二区| 丰满少妇按摩被扣逼高潮| 亚洲精品一区二区| 国产av无码专区亚洲av| 国产福利小视频91| 一区二区三区在线少妇| 我爱我色成人网| 亚洲日韩一区二区一无码| 日本高清无卡一区二区三区| 国产大屁股视频免费区| 69久久夜色精品国产69| 成人精品免费av不卡在线观看 | 成人欧美一区二区三区黑人| 亚洲av无码精品色午夜蛋壳| 嗯啊 不要 啊啊在线日韩a| av成人一区二区三区|