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

    Simulation analysis and pharmacokinetic evaluation of atractylide I for the treatment of diabetic nephropathy

    2022-11-03 06:35:54ShuQingCaoYuJingZhangMengChenYuanYuanZhangChenJieGuoYaoYao
    Medical Data Mining 2022年4期

    Shu-Qing Cao,Yu-Jing Zhang,Meng Chen,Yuan-Yuan Zhang,Chen-Jie Guo,Yao Yao

    1School of Forensic Medicine,Shanxi Medical University,Jinzhong 030600,China.2School of Public Health,Shanxi Medical University,Taiyuan 030001,China.3School of Clinical Medicine,Shanxi Medical University,Taiyuan 030000,China.

    Abstract Objective:The study aims to investigate the molecular mechanism and druggability of Atractylide I (AT-I) for the treatment of diabetic nephropathy (DN) using network pharmacology,molecular docking,and pharmacokinetic analysis.Method: Multiple databases and service platforms were used to predict and screen the potential genes of AT-I for the prevention and treatment of DN.Further,gene ontology and kyoto encyclopedia of genes and genomes enrichment analyses were performed on the genes to evaluate the mechanism of AT-I against DN.The protein interaction network was constructed by the STRING database,and the key gene targets of AT-I against DN were screened by Cytoscape software.The binding activity of AT-I with the essential target proteins was analyzed by molecular docking technique.Finally,the pharmacokinetic parameters of AT-I were evaluated to determine its druggability.Result:AT-I has 45 drug action targets intersecting with DN targets,which may be the potential targets of AT-I for the prevention and treatment of DN.The enrichment analysis indicated that gene expression regulation,insulin resistance,hypoxia inducible factor-1 signaling pathway,tumor necrosis factor signaling pathway,and PI3K-Akt signaling pathway determine the biological mechanism of AT-I in the prevention and treatment of DN.Further,network pharmacology and molecular docking results indicated that key genes such as AKT1,MAPK8,EGFR,MAPK1,MMP9,and MTOR bind well with AT-I,while AKT1 is probably the primary target gene of AT-I against DN.Pharmacokinetic analysis showed that AT-I conforms to Lipinski’s rule of five drug-like properties,has good gastrointestinal absorption,can pass the blood-brain barrier,and has superior biological properties.Conclusion:AT-I has desirable druggability and multi-target,multi-pathway,and multi-mechanism activity against DN.This study provides an important basis for future research on the activity of AT-I against DN.

    Keywords: Atractylide I;diabetic nephropathy;inflammatory response;AKT1;pharmacokinetics;network pharmacology

    Introduction

    DN is one of the most common and serious microvascular complications of diabetes mellitus.In recent years,DN incidence has been on the rise,with a gradual increase in the diabetic population in China and aging [1].Microproteinuria and reduced glomerular filtration rate are the major diagnostic indicators of early-stage DN[2].A rapid decline in renal function can lead to a gradual development of DN into end-stage renal disease.Due to the complex pathogenesis of DN and the lack of specific drugs,the current clinical prevention measures of DN include reduction in cardiovascular risk,controlling blood glucose,controlling blood pressure,and inhibiting the renin-angiotensin-aldosterone system [3].Despite current therapeutic approaches,DN onset and progression risk remains high.Therefore,the search for more effective anti-DN drugs is imperative.

    Chinese medicine has a long history of treating DN.Based on its clinical manifestations,Chinese medicine classifies DN under“edema”,“kidney dissipation”,“asthenia”,“urinary turbidity”,etc.,and has significant advantages over common Western medicine treatment [4].Atractylide,the rhizome of Atractylide,species of genus Atractylodes in the family Asteraceae,is bitter,sweet,and spleen and stomach-friendly.It can tonify the qi to strengthen the spleen,dry dampness,and promote water retention [5].Wenshun et al.reviewed the modern DN literature and based on the statistical analysis of the Ancient and Modern Medical Case Cloud platform,concluded that Atractylide ranked ninth in terms of frequency of use (with 42 times)for the treatment of early DN,and was one of the major herbal medicines used in TCM for DN treatment [6].A study pointed out that metformin combined with Atractylide powder from Scutellaria baicalensis could significantly reduce the levels of Interleukin (IL)-6 and tumor necrosis factor (TNF)-α in DN patients and improve the insulin sensitivity of the body through anti-inflammatory effects,thus improving DN[7].It was observed that Atractylide Dihuang decoction could significantly reduce urinary microalbumin in DN patients,and there was no toxic clinical side.The efficacy was shown to be better combined with Western medical foundations than conventional Western medical treatment alone [8].Atractylide is rich in volatile oil and polysaccharides,and volatile oils of Atractylide have been attributed to multiple pharmacological activities such as antioxidant,anti-inflammatory,immune modulation,and glucose metabolism [9].A research study by Qinhua et al.highlighted that AT-I could promote significant changes in the cytokine expression levels of inflammatory macrophages,and its anti-inflammatory activity was stronger than that of Atractylide-III [10].CQ et al.screened the anti-inflammatory component of Atractylide by leukocyte membrane chromatography and concluded that AT-I had significant anti-inflammatory effects in both acute and chronic inflammation models in mice [11].The anti-inflammatory activity of AT-I was similar to that of dexamethasone.However,it is noteworthy that there are only a few studies on the anti-DN effects of single drug Atractylide and its active ingredient AT-I.The mechanism of action of AT-I against DN is not well understood.

    Several factors are involved in the pathogenesis of DN,including genetic factors,environmental factors,and pathophysiological factors such as inflammatory response,oxidative stress,and renal hemodynamic alterations [12].The traditional one-drug,one-target model is no longer suitable for the treatment of complex diseases such as DN.In this context,the multi-component,multi-pathway,and multi-target characteristics of traditional Chinese medicine bring a ray of hope for the prevention and treatment of DN.Therefore,network pharmacology and computer-aided drug design methods are used in this study to understand the mechanism of action of AT-I against DN.The pharmacokinetic characteristics of AT-I can provide a scientific basis for the clinical development of Atractylide for the prevention and treatment of DN (Figure 1).

    Materials and methods

    AT-I target acquisition

    In this study,the drug targets of AT-I were evaluated using two strategies,namely,database search and molecular conformational relationship search.

    The targets and SMILE codes of AT-I were retrieved from TCMSP(https://tcmspw.com/tcmsp.php) and Herb (http://herb.ac.cn/) with“Atractylide” as the keyword.The 3D chemical structure of AT-I was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/).The energy minimized structure of AT-I was obtained using the minimum energy module of the MM2 procedure in ChemBio3D Ultra 14.0 software.The structure was saved in mol2 format and then imported into the Pharm Mapper server(http://www.lilab-ecust.cn/pharmmapper). The pharmacophore matching method was used to predict the potential targets of AT-I.The top 100 human target proteins were selected according to the “Norm Fit” score in descending order.Further,the SMILE code of AT-I was entered into the Swiss Target Prediction web server(http://www.swisstargetprediction.ch/) to predict the potential targets based on the molecular structure of AT-I.The UniProt database(https://www.uniprot.org/) was used to identify the genes for the above target proteins and integrate the above four sources to obtain AT-I targets.

    AT-I target gene acquisition for DN treatment

    For identification of genes,the keyword “diabetic nephropathy” was used,and the species was limited to “Homo sapiens”.The target genes associated with DN were obtained from DisGeNET (https://www.disgenet.org/),GeneCards (https://www.genecards.org/),and CTD (https://ctdbase.com) databases.The top 300 genes with high scores were selected and combined after screening out duplicates to develop a DN gene collection.The DN genes thus obtained and AT-I drug genes were provided as inputs into Venny 2.1.0(https://bioinfogpcnb.csic.es/),and the intersections of the two sets of genes were the targets of Atractylide against DN.

    Bioenrichment analysis

    The gene ontology enrichment analysis and kyoto encyclopedia of genes and genomes pathway enrichment analysis were performed using DAVID (https://david.ncifcrf.gov/) to identify Atractylide targets against DN.The species was limited to “Homo sapiens”,and the biological functions and signaling pathways withP<0.001 were screened out.Thus,the potential biological mechanism of AT-I against DN was evaluated.

    Screening of key targets of AT-I against DN and construction of protein-protein interaction network

    AT-I targets against DN were imported into the STRING database(https://string-db.org/),and protein interaction analysis was conducted by setting the species as “Homo sapiens”.The protein interaction network was constructed by fixing the confidence score as“medium confidence=0.4” and other default settings.The tab-separated values file was saved and imported into Cytoscape 3.9.1 software.The top 10 genes were selected as the key target genes of AT-I for DN control using the maximal clique centrality algorithm in the CytoHubba plug-in.

    Molecular docking validation

    For further verification of AT-I affinity for crucial target genes,the 3D structures of key target genes were identified and downloaded from PDB (http://www.rcsb.org/) using the original mol2 format of AT-I.Then,the key proteins were pre-processed using AutoDockTools for dehydration,hydrogenation,detection of molecular center nodes,etc.After the identification of the key protein as a high molecular weight receptor with AT-I as a low molecular weight ligand for molecular docking,the results were visualized using Pymol software.

    Pharmacokinetic evaluation

    Figure 1 Flowchart of research design.AT-I,Atractylide I;DN,diabetic nephropathy;PPI,protein-protein interaction;GO,gene ontology;KEGG,kyoto encyclopedia of genes and genomes.

    The SMILES code of AT-I was imported into the SwissADME server(http://www.swissadme.ch/) for analyzing and evaluating the pharmacokinetic parameters of AT-I such as druggability,blood-brain barrier permeability,gastrointestinal absorption,and hepatic pharmacase inhibition.

    Result

    Anti-DN gene screening for AT-I

    The AT-I drug targets obtained from TCMSP,Herb,Swiss Target Prediction,and Pharm Mapper databases were combined and de-weighted to obtain 158 AT-I targets.Then based on the relevant score,the top 300 human DN genes were obtained from GeneCards,DisGeNET,and CTD databases.After combining and de-weighting,621 human DN major target genes were obtained.The above two sets of genes were then intersected to obtain 45 anti-DN genes for AT-I(Figure 2).

    Bioenrichment analysis

    Figure 2 Potentially useful AT-I genes for DN prevention and treatment.AT-I,Atractylide I;DN,diabetic nephropathy.

    Figure 3 Bar chart showing gene ontology bioenrichment analysis

    For elucidation of the biological processes and enrichment pathways involved in the action of AT-I genes against DN,a bioenrichment analysis was conducted on 45 potential genes.The top ten smallest biological processes,cellular components,and molecular function withP< 0.001 were visualized for gene ontology functional enrichment analysis (Figure 3).The enrichment results indicated that the major DN genes were involved in the negative regulation of the apoptotic process,positive regulation of transcription from RNA polymerase II promoter,and negative regulation of cellular response to insulin stimulation.The negative regulation of the apoptotic process,positive regulation of transcription from RNA polymerase II promoter,negative regulation of cellular response to insulin stimulus,and positive regulation of gene expression are biological processes.The cellular components and molecular function are primarily related to the cell surface receptor complex,steroid hormone receptor activity,insulin binding,and phosphatidylinositol 3-kinase binding.Thus,it is implied that AT-I can interfere with DN progression by affecting gene expression and hormone binding.

    The top 20 smallest kyoto encyclopedia of genes and genomes enrichment entries withP<0.001 were visualized (Figure 4).The enrichment results indicated that AT-I mainly inhibited DN development through the hypoxia inducible factor-1 (HIF-1) signaling pathway,pathways involved in cancer development,TNF signaling pathway,insulin resistance,PI3K-Akt signaling pathway,vascular endothelial growth factor signaling pathway,etc.(Table 1).

    Anti-DN key target screening and protein-protein interaction construction for AT-I

    For elucidation of the interactions between the major target genes and for extracting the key genes of AT-I for DN prevention and treatment,45 major target genes were imported into the String database for protein-protein interaction analysis.The genes with disconnected nodes were hidden,and 45 nodes with 378 edges were obtained.This indicated that the 45 commonly associated genes were interlinked by 378 edges.After saving the file in the tab-separated values format and mapping it in Cytoscape software,the top 10 key genes(Table 2) were filtered using the maximal clique centrality algorithm in the CytoHubba plugin (Figure 5).

    Molecular docking validation

    The molecular docking validation of the obtained 10 key gene targets of AT-I was performed using AutoDockTools software,and the binding activities of their high molecular weight receptors and low molecular weight ligands are shown in Table 3.It can be observed that the binding energy of all AT-I and 10 key anti-DN genes were less than or equal to -5.0 kcal-mol-1,i.e.,AT-I has a strong binding with the key genes(Figure 6).

    Figure 4 Scatter plot of kyoto encyclopedia of genes and genomes enrichment results.VEGF,vascular endothelial growth factor;HIF-1,hypoxia inducible factor-1;TNF,tumor necrosis factor;AMPK,adenosine 5’-monophosphate-activated protein kinase.

    Figure 5 PPI diagram of key genes.PPI,protein-protein interaction.

    Figure 6 Molecular docking map of AT-I and key genes.AT-I,Atractylide I.

    Table 1 Kyoto encyclopedia of genes and genomes pathway enrichment analysis

    Table 2 Key genes involved in the prevention and treatment of DN by AT-I

    Table 3 Results of AT-I molecular docking

    Pharmacokinetic evaluation

    The SMILE code CC1=C2CC3C(=C)CCCC3(C=C2OC1=O)C of AT-I was imported into the SwissADME platform,and its prototype pharmacokinetic parameters were analyzed.As shown in the bioavailability radar plot (Figure 7B),it can be observed that the six physicochemical properties of Atractylide,namely,lipophilicity,size,polarity,solubility,saturation,and flexibility are ideal.AT-I conforms to the five rules of Lipinski-like drugs and has good druggability(Table 4).

    Discussion

    DN has become the leading cause of chronic kidney disease in China,and its growth has adversely affected the quality of life of patients and imposed a great burden on their families and society [13].As a core drug in the treatment of DN in TCM,the lack of research to understand the mechanism of action and druggability of a single drug and its primary active ingredient against DN has restricted the clinical use of Atractylide and the discovery of effective novel anti-DN drugs to a certain extent.

    Using the data from multiple platforms,a total of 45 major targets of AT-I against DN were obtained in this study.The bioenrichment analysis indicated that AT-I can exert anti-DN activity through multi-target and multi-pathway characteristics.HIF is a major regulator of the body’s adaptation to hypoxia and can interact with signaling pathways such as Notch and NF-kB for directly stimulating the expression of hundreds of target genes to maintain cellular oxygen homeostasis [14].In a study including 1,165 US type 1 diabetic patients with and without DN,researchers determined a protective association between the hypoxia inducible factor-1 alpha (HIF1A)Pro582Ser polymorphism and DN in male subjects by identifying polymorphism of the HIF1A gene through in vivo and in vitro experiments [15].Moreover,the animal experiments indicated that conditional knockout of HIF-1α exacerbates renal tubular injury in DN mice [16].This phenomenon did not occur when there was an overexpression of HIF-1α and in the HO-1 agonist groups,thus implying the alleviation of renal tubular injury in DN mice by HIF-1α via HO-1-mediated mitochondrial dynamics.

    On the other hand,AT-I attenuated tissue injury by inhibiting hypoxia-induced HIF-1α and downregulating the Akt signaling pathway [17].Diabetes significantly increases the risk of cancer in diabetic patients when compared to non-diabetic patients and leads to more than one-fold higher cancer mortality than non-diabetic patients[18,19].A study conducted by Luo et al.showed that the disturbance of calcium metabolism in DN rats was associated with their susceptibility to kidney cancer [20].Dong et al.analyzed the co-expressed differential genes between DN and kidney cancer using bioinformatics and observed a significant association between DN and kidney cancer [21].In our previous study,the differentially expressed genes,namely,VEGFA and MMP9,in DN were also significantly upregulated in renal clear cell carcinoma.Moreover,the study also indicated that high MMP9 expression levels were associated with a poorer prognosis for patients.Interestingly,AT-I has been shown to inhibit cell proliferation and induce cancer cell death.It has also been clinically tested for the treatment of cancer-induced cachexia [9].Insulin resistance (IR) refers to a decrease in the sensitivity of insulin-dependent receptor cells or disruption of insulin-mediated signaling pathways.IR plays an essential role in the development of diabetes and its complications [22].In addition,the PI3K-AKT pathway,an important signaling pathway for IR production,inhibits hepatic gluconeogenesis and promotes glucose transporters 4 to transfer glucose into the cell,thereby causing IR when AKT is inactivated,inhibited,or continuously activated [23].A study by Han et al.indicated that Atractylide extracts promoted glucose transport by increasing glucose transporters 4,PI3K,and IRS-1 levels in addition to organismal lipid differentiation [24].However,the role of AT-I was not determined in that study.A meta-analysis suggested that TNF-α-308G/A polymorphism may enhance DN susceptibility in Asian patients with type 2 diabetes [25].A research study by Xu et al.showed that the expression levels of TNF-α and IL-6 were significantly higher in DN patients when compared to those of normal human subjects [26].When DN-affected rats were treated with TNF-α inhibitors for 12 weeks,it was found that TNF-α inhibition significantly reduced the expression of NLRP3 inflammatory vesicles and IL-6 and L-17A mRNA in the renal tubules,thereby resulting in reduced albuminuria levels and glomerular and tubular damage in DN rats [27].Moreover,it has been observed that AT-I can significantly inhibit the production of TNF-α,IL-6,and IL-1β through the TLR4/NF-kB signaling pathway,thus acting as a protective barrier[28,29].Lastly,due to the lack of relevant literature,it is difficult to directly understand the mechanism of action and effect of AT-I against DN.However,based on the available literature,it can be predicted that multiple pathways are involved in the action of AT-I against DN.

    Figure 7 Structural formula of AT-I (A) and radar diagram of bioavailability (B).AT-I,Atractylide I.

    Table 4 Pharmacokinetic parameters of AT-I

    Next,the algorithm was used to screen the ten key genes of AT-I against DN,and their AT-I binding activity was verified by molecular docking.The results indicated that all the key genes have a strong binding activity with AT-I,and there is some similarity between the active sites of AT-I and the inhibitors of key genes.This further illustrates the reliability and authenticity of the results from this study,with AKT1 found to have the strongest binding activity with AT-I,and is presumed to be the possible major target of AT-I against DN.AKT is a downstream effector of the PI3K pathway that regulates many cellular functions,including glucose metabolism,glycogen synthesis,cell proliferation,death,and hypertrophy.It is also involved in several signaling pathways [30].AKT is an important mediator of insulin action,which regulates glucose uptake in muscles,adipocytes,liver,and other tissues by promoting the translocation of glucose transporter proteins from intracellular stores to the plasma membrane in response to insulin stimulation [31].A study conducted by Xing et al.indicated that repeated low-dose radiation activated AKT phosphorylation in kidney tissues of DN mice,which in turn,significantly ameliorated diabetes-induced renal inflammatory damage [32].Further,a study done by Rane et al.also showed that a prolonged hyperglycemic environment in the body downregulated the level of AKT activation and accelerated the renal tubular cell apoptosis via enhanced p38 MAPK activation,which in turn,contributed to the DN development [33].AKT1,one of the isoforms of AKT,regulates neuronal insulin signaling,glucose uptake,and insulin resistance.Numerous experimental studies have indicated that AKT1 inhibition leads to increased levels of uncoupled mitochondrial respiration and oxidative stress found in the renal tubules,while activation of AKT1 signaling prevents the development of acute kidney damage and subsequent chronic kidney disease [34,35].However,it should be noted that the role of AKT1 in the development of DN was confirmed in our previous study and other network pharmacology studies [36,37].Although previous studies have shown that mixed herbs in combination with Atractylide lead to good DN control,there is still a lack of studies on the activity of Atractylide and AT-I on DN and AKT1.It is important to experimentally verify whether AKT1 is the primary target of AT-I against DN.However,the present study undoubtedly highlights the importance of Atractylide as an active ingredient.In addition,other key genes also have a major impact on DN.In a study that included 650 Chinese Han patients with type 2 DN and 580 patients with only diabetes mellitus,genetic variants at the rs2010963 and rs69947 loci of the VEGFA gene were found to increase the risk of DN possibly [38].Moderate amounts of VEGFA are beneficial for maintaining the glomerular structure;however,excess VEGFA leads to abnormal angiogenesis and thereby promotes the development of DN [39].AT-I downregulates vascular endothelial growth factor expression,and VEGFA,one of the isoforms of vascular endothelial growth factor,may improve DN by inhibiting renal tissue angiogenesis[40].A study by Zhu et al.highlighted that the ALB gene encoding serum albumin levels were associated with chronic vascular complications in Chinese patients with type 2 diabetes [41].A network pharmacology study indicated that AT-I can act on the ALB gene and then participate in the regulation of metabolic syndrome by the Atractylide-Porin drug pair and improving glucolipid metabolism and insulin sensitivity [42].AT-I was shown to inhibit the expression of EGFR,and the activation of EGFR has a protective effect on the kidney [43].However,sustained activation exacerbates DN progression,and the pathophysiology of DN involves several factors such as altered hemodynamics,metabolic disturbances,and inflammatory responses [44].CASP3-dependent scorch death pathway has attracted attention as a newly discovered regulatory mechanism involved in DN progression,and studies have highlighted that AT-I can inhibit CASP3 activity while retaining mitochondrial function[45].The analogs to CASP3,MTOR,and AKT1 are together involved in the PI3K/AKT/mTOR signaling pathway,and it has been observed that the activation of MTOR reduces SOD activity and GSH levels and activates CASP3,which in turn,promotes oxidative stress-induced DN thylakoid cell apoptosis [46-48].MAPK1 and MAPK8 are MAPK family proteins that play a key role in regulating gene expression and cytoplasmic functional activities.AT-I can inhibit apoptosis and oxidative stress by inhibiting the MAPK pathway and inhibiting the activities of JAK2/STAT3,PI3K/Akt,p38 MAPK,and Wnt/β-linked protein to reduce renal fibrosis[49,50].Although there are only a few studies on the mechanism of AT-I against DN,the major key targets of AT-I against DN screened by the data platform and algorithm in this study directly or indirectly corroborated with the previous studies,thereby implying that the action of AT-I against DN involves multiple targets.

    Finally,the pharmacokinetics of AT-I was evaluated by simulation.It was found that AT-I conforms to Lipinski’s rule of five drug-like properties,has good druggability and gastrointestinal absorption,and can cross the blood-brain barrier.Therefore,it is expected to be an excellent drug against DN.However,due to its inhibitory effect on some hepatic enzymes,utmost care should be taken when used in combination with other drugs.

    Lastly,this study utilizes network pharmacology and computer-aided drug design methods to overcome the limitations of traditional experiments,such as time-consuming processes,high cost,and fewer data.This study investigates the pharmacological mechanism of AT-I against DN,which provides a theoretical basis for further research on the active ingredients of Atractylide and the development of novel anti-DN drugs in the future.

    国产真实伦视频高清在线观看| 下体分泌物呈黄色| 久久久久久久久久久丰满| 国产精品99久久久久久久久| 欧美精品一区二区大全| 亚洲精品色激情综合| 亚洲色图综合在线观看| 成人亚洲精品av一区二区| 免费看不卡的av| 精品熟女少妇av免费看| 日韩三级伦理在线观看| 久久久色成人| 一级毛片aaaaaa免费看小| 日韩成人伦理影院| 麻豆精品久久久久久蜜桃| 91久久精品国产一区二区三区| 美女xxoo啪啪120秒动态图| 久久久久九九精品影院| 色婷婷久久久亚洲欧美| 超碰97精品在线观看| 联通29元200g的流量卡| 国产一区二区三区综合在线观看 | 欧美少妇被猛烈插入视频| 国产黄a三级三级三级人| 在线观看国产h片| 2018国产大陆天天弄谢| 丰满人妻一区二区三区视频av| 午夜激情久久久久久久| 听说在线观看完整版免费高清| av国产精品久久久久影院| 一本久久精品| 少妇 在线观看| 欧美日韩亚洲高清精品| videos熟女内射| 亚洲精品一二三| 亚洲经典国产精华液单| 国产成人精品久久久久久| 国产亚洲5aaaaa淫片| 99久久精品热视频| 91精品伊人久久大香线蕉| 日本黄色片子视频| 丝袜喷水一区| 久久亚洲国产成人精品v| 国产毛片a区久久久久| 国产有黄有色有爽视频| 91午夜精品亚洲一区二区三区| 26uuu在线亚洲综合色| 在线观看三级黄色| 综合色丁香网| 最近手机中文字幕大全| 亚洲欧美日韩另类电影网站 | 免费观看的影片在线观看| 99热全是精品| 国产亚洲av片在线观看秒播厂| 在线精品无人区一区二区三 | 在线观看国产h片| 韩国av在线不卡| 欧美日本视频| 成人国产麻豆网| 国产成人a∨麻豆精品| 亚洲不卡免费看| 亚洲色图av天堂| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 国产一区二区亚洲精品在线观看| 看十八女毛片水多多多| 伊人久久精品亚洲午夜| 91精品国产九色| 成年av动漫网址| 国产精品99久久久久久久久| 亚洲激情五月婷婷啪啪| 五月玫瑰六月丁香| 亚洲高清免费不卡视频| 欧美 日韩 精品 国产| 91狼人影院| 欧美精品人与动牲交sv欧美| 色吧在线观看| 久久人人爽av亚洲精品天堂 | 国语对白做爰xxxⅹ性视频网站| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 国产精品一区二区在线观看99| 国产精品国产三级专区第一集| 国产欧美另类精品又又久久亚洲欧美| av在线播放精品| 久久精品久久久久久久性| 中文字幕久久专区| 国产精品一区二区在线观看99| 中文欧美无线码| 夫妻午夜视频| 97热精品久久久久久| 亚洲国产精品成人久久小说| 久久久久久久亚洲中文字幕| 色综合色国产| 精品人妻视频免费看| 色哟哟·www| 毛片女人毛片| 国产欧美日韩一区二区三区在线 | 午夜福利在线观看免费完整高清在| 只有这里有精品99| 777米奇影视久久| 国产免费视频播放在线视频| 成人二区视频| 国产一区有黄有色的免费视频| 伦精品一区二区三区| 亚洲av日韩在线播放| 国产伦精品一区二区三区视频9| 欧美性感艳星| 亚洲av免费高清在线观看| 精品亚洲乱码少妇综合久久| 黄色配什么色好看| 婷婷色麻豆天堂久久| 国产成人精品婷婷| 激情 狠狠 欧美| 日韩免费高清中文字幕av| 免费看日本二区| 99久久精品一区二区三区| 久久久久久久久久人人人人人人| 亚洲自拍偷在线| 2018国产大陆天天弄谢| 王馨瑶露胸无遮挡在线观看| 免费看日本二区| 18+在线观看网站| 大话2 男鬼变身卡| 日韩精品有码人妻一区| 久久人人爽av亚洲精品天堂 | 国产精品99久久99久久久不卡 | 婷婷色综合www| 嘟嘟电影网在线观看| 小蜜桃在线观看免费完整版高清| 人体艺术视频欧美日本| 五月伊人婷婷丁香| 久久ye,这里只有精品| 亚洲四区av| 日本三级黄在线观看| 搞女人的毛片| 亚洲精品久久久久久婷婷小说| 自拍欧美九色日韩亚洲蝌蚪91 | 最近中文字幕高清免费大全6| 久久精品国产亚洲av天美| 成人亚洲欧美一区二区av| 国产成年人精品一区二区| 日韩 亚洲 欧美在线| 亚洲在久久综合| 少妇人妻 视频| 久久国内精品自在自线图片| 日韩国内少妇激情av| 国产成人91sexporn| 在线 av 中文字幕| 国产成人一区二区在线| 免费不卡的大黄色大毛片视频在线观看| 在线免费十八禁| 日韩成人伦理影院| 黑人高潮一二区| 秋霞在线观看毛片| 插逼视频在线观看| 亚洲精品国产色婷婷电影| 亚洲国产成人一精品久久久| 亚洲av欧美aⅴ国产| 青青草视频在线视频观看| 日本三级黄在线观看| 纵有疾风起免费观看全集完整版| 久久人人爽av亚洲精品天堂 | 亚洲av电影在线观看一区二区三区 | 亚洲三级黄色毛片| 99热国产这里只有精品6| 最近手机中文字幕大全| 国产成人freesex在线| 99久久精品热视频| 久久99热6这里只有精品| 三级经典国产精品| 嫩草影院精品99| 日本三级黄在线观看| 中文天堂在线官网| 国产精品久久久久久av不卡| 人人妻人人看人人澡| 精品国产乱码久久久久久小说| 免费观看的影片在线观看| kizo精华| 日日撸夜夜添| 蜜桃亚洲精品一区二区三区| av网站免费在线观看视频| 久久ye,这里只有精品| 99久久人妻综合| 欧美精品人与动牲交sv欧美| 国模一区二区三区四区视频| 人人妻人人爽人人添夜夜欢视频 | 99视频精品全部免费 在线| 男女边吃奶边做爰视频| 九九爱精品视频在线观看| 纵有疾风起免费观看全集完整版| 亚洲高清免费不卡视频| 国产伦在线观看视频一区| 七月丁香在线播放| 极品教师在线视频| 日本午夜av视频| 国产免费一区二区三区四区乱码| 亚洲性久久影院| 中国三级夫妇交换| 中文字幕免费在线视频6| 大话2 男鬼变身卡| 欧美激情在线99| 建设人人有责人人尽责人人享有的 | 在线免费十八禁| 伦精品一区二区三区| 一本色道久久久久久精品综合| 人妻系列 视频| 老司机影院毛片| 亚洲最大成人av| 亚洲国产精品国产精品| 免费看a级黄色片| 最近手机中文字幕大全| 99久久九九国产精品国产免费| 国产成人福利小说| 国产色爽女视频免费观看| 街头女战士在线观看网站| 久久久久九九精品影院| 伦理电影大哥的女人| 麻豆国产97在线/欧美| 国产欧美日韩一区二区三区在线 | 欧美日韩亚洲高清精品| 边亲边吃奶的免费视频| 精品久久久精品久久久| 麻豆成人av视频| 亚洲av免费在线观看| 男人狂女人下面高潮的视频| 国产爱豆传媒在线观看| 99久久精品国产国产毛片| 欧美人与善性xxx| 各种免费的搞黄视频| 干丝袜人妻中文字幕| 精品熟女少妇av免费看| 一区二区三区免费毛片| 久久久久国产网址| 我要看日韩黄色一级片| 精品视频人人做人人爽| 国产精品一及| 蜜桃久久精品国产亚洲av| 国产黄色免费在线视频| 禁无遮挡网站| 免费av不卡在线播放| 青春草亚洲视频在线观看| 中文在线观看免费www的网站| 国产午夜精品久久久久久一区二区三区| 亚洲欧美清纯卡通| 免费av不卡在线播放| 亚洲av电影在线观看一区二区三区 | 婷婷色综合大香蕉| 建设人人有责人人尽责人人享有的 | 伊人久久精品亚洲午夜| 精品久久国产蜜桃| 精品久久久噜噜| 26uuu在线亚洲综合色| 少妇人妻精品综合一区二区| 成人毛片60女人毛片免费| 色播亚洲综合网| 国产午夜精品一二区理论片| 日本午夜av视频| 精品久久久久久久末码| 精品亚洲乱码少妇综合久久| 人人妻人人澡人人爽人人夜夜| av免费在线看不卡| 亚洲精品中文字幕在线视频 | 日韩欧美一区视频在线观看 | 免费观看a级毛片全部| 水蜜桃什么品种好| 两个人的视频大全免费| 国产亚洲精品久久久com| 国产精品99久久久久久久久| 国产成人福利小说| 日韩视频在线欧美| 精品久久久久久电影网| 欧美成人午夜免费资源| 久久精品国产自在天天线| 在线免费观看不下载黄p国产| 国产av码专区亚洲av| 日本爱情动作片www.在线观看| 成人午夜精彩视频在线观看| 少妇的逼好多水| 中文在线观看免费www的网站| 五月伊人婷婷丁香| 蜜桃亚洲精品一区二区三区| 日韩欧美精品v在线| 日韩 亚洲 欧美在线| 免费观看在线日韩| 亚洲高清免费不卡视频| 久久久久久久久久久免费av| 女人十人毛片免费观看3o分钟| 免费看不卡的av| 久久人人爽人人爽人人片va| 日韩视频在线欧美| 久久久久久伊人网av| 成人亚洲精品av一区二区| 亚洲真实伦在线观看| 人妻少妇偷人精品九色| 亚洲精品aⅴ在线观看| 亚洲国产最新在线播放| 波野结衣二区三区在线| 亚洲精品影视一区二区三区av| 欧美+日韩+精品| 丰满乱子伦码专区| 波多野结衣巨乳人妻| 亚洲自偷自拍三级| 免费少妇av软件| 国产精品久久久久久精品电影| 一级av片app| 久久久午夜欧美精品| 国产片特级美女逼逼视频| 免费av不卡在线播放| 国产高清有码在线观看视频| 日韩,欧美,国产一区二区三区| 舔av片在线| 午夜视频国产福利| 一级毛片电影观看| 国产成人午夜福利电影在线观看| 91久久精品国产一区二区三区| 久久久久久久久大av| 一级av片app| 一区二区三区精品91| 乱系列少妇在线播放| 亚洲电影在线观看av| 毛片一级片免费看久久久久| 伦理电影大哥的女人| 日本爱情动作片www.在线观看| 日韩成人伦理影院| 久久人人爽av亚洲精品天堂 | 国产亚洲5aaaaa淫片| 熟女人妻精品中文字幕| 男女边摸边吃奶| 亚洲欧美中文字幕日韩二区| 中文天堂在线官网| 亚洲欧美一区二区三区国产| 日韩av在线免费看完整版不卡| 精品一区二区三区视频在线| 1000部很黄的大片| 乱系列少妇在线播放| 纵有疾风起免费观看全集完整版| 国产 精品1| 国产成人精品久久久久久| 亚洲综合精品二区| 高清午夜精品一区二区三区| 欧美成人一区二区免费高清观看| 成人亚洲欧美一区二区av| av卡一久久| 亚洲色图av天堂| 91久久精品国产一区二区成人| 亚洲在久久综合| 青青草视频在线视频观看| 欧美老熟妇乱子伦牲交| 久久综合国产亚洲精品| 不卡视频在线观看欧美| 69av精品久久久久久| 欧美激情久久久久久爽电影| 毛片一级片免费看久久久久| 亚洲成人av在线免费| 午夜福利视频精品| 欧美日韩综合久久久久久| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 男人和女人高潮做爰伦理| 夜夜看夜夜爽夜夜摸| 国语对白做爰xxxⅹ性视频网站| 在线a可以看的网站| 国产免费一区二区三区四区乱码| 熟女电影av网| 老司机影院成人| 身体一侧抽搐| 最近2019中文字幕mv第一页| 蜜臀久久99精品久久宅男| 久久久久性生活片| 亚洲熟女精品中文字幕| 亚洲久久久久久中文字幕| 少妇被粗大猛烈的视频| 99久久人妻综合| 精品久久久久久电影网| 国产精品麻豆人妻色哟哟久久| 高清日韩中文字幕在线| 亚洲在线观看片| 在线观看三级黄色| 人妻 亚洲 视频| 夫妻午夜视频| 色5月婷婷丁香| 午夜亚洲福利在线播放| 国产成人精品福利久久| 一区二区三区免费毛片| 午夜福利在线观看免费完整高清在| 国产亚洲av嫩草精品影院| 国产精品久久久久久久电影| 久久久久网色| 日韩av在线免费看完整版不卡| 99热6这里只有精品| 久久99精品国语久久久| 女的被弄到高潮叫床怎么办| 国产精品秋霞免费鲁丝片| 日日啪夜夜爽| 亚洲精品色激情综合| 大又大粗又爽又黄少妇毛片口| 高清av免费在线| 精品视频人人做人人爽| 搡老乐熟女国产| 亚洲性久久影院| 国产黄片视频在线免费观看| 国产黄片美女视频| 美女cb高潮喷水在线观看| 欧美xxⅹ黑人| 亚洲人与动物交配视频| 精品人妻视频免费看| 狠狠精品人妻久久久久久综合| 男女下面进入的视频免费午夜| 熟妇人妻不卡中文字幕| 97超视频在线观看视频| 可以在线观看毛片的网站| 亚洲精品国产色婷婷电影| 色5月婷婷丁香| 嫩草影院新地址| 成人免费观看视频高清| 最近手机中文字幕大全| 天天躁夜夜躁狠狠久久av| 久久精品国产亚洲网站| 亚洲精品日韩av片在线观看| 久久精品国产鲁丝片午夜精品| 亚洲一级一片aⅴ在线观看| 免费少妇av软件| 国产一区二区亚洲精品在线观看| 亚洲性久久影院| 欧美zozozo另类| 国产白丝娇喘喷水9色精品| 成人无遮挡网站| 亚洲欧美一区二区三区黑人 | 欧美成人一区二区免费高清观看| 三级国产精品欧美在线观看| 人妻系列 视频| 精华霜和精华液先用哪个| 内地一区二区视频在线| 亚洲自偷自拍三级| 亚洲精华国产精华液的使用体验| 亚洲欧美日韩卡通动漫| 插逼视频在线观看| 亚洲欧美一区二区三区国产| 成人黄色视频免费在线看| 亚洲人成网站在线播| 国产成人免费观看mmmm| 日韩欧美一区视频在线观看 | 国产亚洲午夜精品一区二区久久 | 一区二区三区免费毛片| 九草在线视频观看| 久久综合国产亚洲精品| 亚洲三级黄色毛片| 久久久久久伊人网av| 男人狂女人下面高潮的视频| 久久久成人免费电影| 国产成人精品一,二区| 国产亚洲91精品色在线| 日本一本二区三区精品| 草草在线视频免费看| 在线免费十八禁| 国产精品蜜桃在线观看| 看免费成人av毛片| 色视频在线一区二区三区| 国产探花极品一区二区| 人人妻人人澡人人爽人人夜夜| 一级片'在线观看视频| 丝袜喷水一区| 丰满人妻一区二区三区视频av| 在线亚洲精品国产二区图片欧美 | 国产熟女欧美一区二区| 欧美zozozo另类| 国产一区二区三区综合在线观看 | 看非洲黑人一级黄片| 麻豆精品久久久久久蜜桃| 交换朋友夫妻互换小说| 日韩免费高清中文字幕av| 色综合色国产| 禁无遮挡网站| 亚洲精品,欧美精品| 99热这里只有精品一区| 国产成人精品婷婷| 亚洲在久久综合| 一级毛片我不卡| 男的添女的下面高潮视频| 干丝袜人妻中文字幕| 精品一区二区免费观看| 欧美潮喷喷水| 女人十人毛片免费观看3o分钟| 亚洲电影在线观看av| 99热全是精品| 久久久成人免费电影| 国产伦精品一区二区三区四那| 校园人妻丝袜中文字幕| 你懂的网址亚洲精品在线观看| 国产黄色视频一区二区在线观看| 国产一区二区在线观看日韩| 精品一区二区三卡| 日韩成人av中文字幕在线观看| 九九在线视频观看精品| 大香蕉97超碰在线| 亚洲精品视频女| 中文字幕制服av| 国产精品人妻久久久影院| 中文资源天堂在线| 国语对白做爰xxxⅹ性视频网站| 男女边摸边吃奶| 欧美极品一区二区三区四区| 大片电影免费在线观看免费| 亚洲国产日韩一区二区| 日韩亚洲欧美综合| 国产精品无大码| 观看免费一级毛片| 成人漫画全彩无遮挡| 国产男女超爽视频在线观看| 97精品久久久久久久久久精品| 又大又黄又爽视频免费| 国产欧美日韩精品一区二区| 国产av国产精品国产| 亚洲av日韩在线播放| 成人二区视频| 亚洲国产欧美在线一区| 久久精品国产亚洲av涩爱| 性色avwww在线观看| 亚洲精品亚洲一区二区| 亚洲av在线观看美女高潮| 欧美少妇被猛烈插入视频| 亚洲色图av天堂| 精品人妻熟女av久视频| 日韩不卡一区二区三区视频在线| 美女被艹到高潮喷水动态| 国产av国产精品国产| 午夜福利网站1000一区二区三区| 国产v大片淫在线免费观看| 亚洲av福利一区| 日韩国内少妇激情av| 久久久色成人| 3wmmmm亚洲av在线观看| 国产欧美日韩精品一区二区| 一个人观看的视频www高清免费观看| 男的添女的下面高潮视频| 狠狠精品人妻久久久久久综合| 精品久久久久久久久av| 网址你懂的国产日韩在线| 亚洲欧美一区二区三区国产| 国产伦精品一区二区三区视频9| 波多野结衣巨乳人妻| 亚洲精品影视一区二区三区av| 男女边摸边吃奶| 伊人久久国产一区二区| 听说在线观看完整版免费高清| 国产av码专区亚洲av| 啦啦啦啦在线视频资源| 亚洲精品色激情综合| 一个人看视频在线观看www免费| 男人和女人高潮做爰伦理| av国产免费在线观看| 国产探花在线观看一区二区| 少妇猛男粗大的猛烈进出视频 | 国产精品蜜桃在线观看| 精品久久国产蜜桃| 三级男女做爰猛烈吃奶摸视频| 日韩一区二区三区影片| 黑人高潮一二区| 天天一区二区日本电影三级| 国产 精品1| 久久久成人免费电影| 欧美老熟妇乱子伦牲交| 亚洲最大成人手机在线| 黄片wwwwww| 日韩欧美一区视频在线观看 | 一区二区三区四区激情视频| 欧美精品人与动牲交sv欧美| av线在线观看网站| 深爱激情五月婷婷| 成年女人看的毛片在线观看| 极品教师在线视频| 国产一区二区亚洲精品在线观看| 在线a可以看的网站| 国产亚洲一区二区精品| 欧美精品国产亚洲| 久久久久久久久久久免费av| 国产精品国产三级国产专区5o| 爱豆传媒免费全集在线观看| 免费大片黄手机在线观看| 国产黄色视频一区二区在线观看| 国产成人免费观看mmmm| 不卡视频在线观看欧美| 水蜜桃什么品种好| 我的女老师完整版在线观看| 国产成人freesex在线| 国产 精品1| 国产精品一及| 亚洲国产日韩一区二区| 狂野欧美白嫩少妇大欣赏| 国产成人免费观看mmmm| 亚洲国产日韩一区二区| 少妇人妻精品综合一区二区| 亚洲一级一片aⅴ在线观看| 国产欧美亚洲国产| 国产中年淑女户外野战色| 色吧在线观看| 另类亚洲欧美激情| 日本爱情动作片www.在线观看| 免费观看av网站的网址| 99热全是精品| 搡女人真爽免费视频火全软件| 日韩视频在线欧美| 欧美性猛交╳xxx乱大交人| 免费大片黄手机在线观看| 男插女下体视频免费在线播放| 欧美精品一区二区大全| 女人被狂操c到高潮| 九九久久精品国产亚洲av麻豆| 免费少妇av软件| 一级二级三级毛片免费看| 黄色日韩在线| 黄色视频在线播放观看不卡| 亚洲熟女精品中文字幕| 观看免费一级毛片|