• <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.

    欧美人与善性xxx| 欧美精品一区二区大全| 免费av毛片视频| 18禁在线播放成人免费| 久久99热6这里只有精品| 成人美女网站在线观看视频| 国产老妇女一区| 日本与韩国留学比较| 午夜激情久久久久久久| 亚洲av免费高清在线观看| 日日啪夜夜爽| 伊人久久精品亚洲午夜| 久久久久精品久久久久真实原创| 一夜夜www| 国产麻豆成人av免费视频| 亚洲人成网站高清观看| 欧美三级亚洲精品| 美女国产视频在线观看| 99久久精品国产国产毛片| av线在线观看网站| 欧美最新免费一区二区三区| 国产精品一及| av国产久精品久网站免费入址| 亚洲av不卡在线观看| 亚洲,欧美,日韩| 亚洲av中文av极速乱| 亚州av有码| 国产精品久久久久久av不卡| 亚洲人成网站在线播| 亚洲在久久综合| 自拍偷自拍亚洲精品老妇| 有码 亚洲区| 2021少妇久久久久久久久久久| 中文字幕人妻熟人妻熟丝袜美| 波野结衣二区三区在线| 免费电影在线观看免费观看| 中文乱码字字幕精品一区二区三区 | 蜜桃久久精品国产亚洲av| 熟妇人妻久久中文字幕3abv| 在线观看一区二区三区| 亚洲国产最新在线播放| 亚洲人成网站高清观看| 亚洲精品自拍成人| 一级毛片久久久久久久久女| 国产91av在线免费观看| 欧美区成人在线视频| 日本wwww免费看| 亚洲国产精品成人综合色| 亚洲国产成人一精品久久久| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 熟女人妻精品中文字幕| 啦啦啦韩国在线观看视频| 国产一级毛片七仙女欲春2| 丰满人妻一区二区三区视频av| 能在线免费看毛片的网站| 色尼玛亚洲综合影院| 亚洲经典国产精华液单| 久久久久久伊人网av| 欧美日韩国产mv在线观看视频 | 韩国高清视频一区二区三区| 日韩国内少妇激情av| 久久国产乱子免费精品| 日本av手机在线免费观看| 女人被狂操c到高潮| 人妻一区二区av| 国产精品无大码| 日本免费a在线| 麻豆精品久久久久久蜜桃| 欧美日本视频| 又大又黄又爽视频免费| 美女cb高潮喷水在线观看| 在线观看人妻少妇| 亚洲激情五月婷婷啪啪| 在线观看美女被高潮喷水网站| 麻豆av噜噜一区二区三区| 色5月婷婷丁香| 日本与韩国留学比较| 国产黄色小视频在线观看| 亚洲精品中文字幕在线视频 | 色哟哟·www| 日本av手机在线免费观看| 麻豆精品久久久久久蜜桃| videos熟女内射| 内射极品少妇av片p| 欧美日韩在线观看h| 免费观看性生交大片5| 18禁动态无遮挡网站| 国产白丝娇喘喷水9色精品| 亚洲欧洲国产日韩| 国产精品熟女久久久久浪| 天堂网av新在线| 国产高潮美女av| 激情五月婷婷亚洲| 日本黄色片子视频| 国内精品一区二区在线观看| 97超视频在线观看视频| 久久久久精品性色| 国产亚洲精品av在线| 精品人妻视频免费看| 69av精品久久久久久| 欧美成人午夜免费资源| 国产精品一区二区性色av| 久久久国产一区二区| 国产精品美女特级片免费视频播放器| 国产探花极品一区二区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久人人爽人人片av| 亚洲无线观看免费| www.色视频.com| 日本一二三区视频观看| 日韩在线高清观看一区二区三区| 中文字幕av成人在线电影| 搡老乐熟女国产| 亚洲一区高清亚洲精品| 大陆偷拍与自拍| 亚洲18禁久久av| 人妻少妇偷人精品九色| 久久精品国产亚洲av涩爱| 最新中文字幕久久久久| 中文乱码字字幕精品一区二区三区 | 国产成人91sexporn| 18禁在线播放成人免费| 中文字幕制服av| 国产高清不卡午夜福利| 国产精品人妻久久久影院| 国产人妻一区二区三区在| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 国产av国产精品国产| 国产在线一区二区三区精| 在线天堂最新版资源| 一边亲一边摸免费视频| 亚洲国产成人一精品久久久| 91久久精品电影网| 91久久精品国产一区二区成人| 高清在线视频一区二区三区| 亚洲av中文av极速乱| 久久97久久精品| 精品亚洲乱码少妇综合久久| 边亲边吃奶的免费视频| 日本爱情动作片www.在线观看| 国产91av在线免费观看| 日韩av不卡免费在线播放| 在线观看美女被高潮喷水网站| 欧美成人午夜免费资源| 日日啪夜夜撸| 亚洲国产最新在线播放| 亚洲无线观看免费| 五月玫瑰六月丁香| 久热久热在线精品观看| 亚洲电影在线观看av| 综合色av麻豆| 极品少妇高潮喷水抽搐| 亚洲精品日韩在线中文字幕| 最后的刺客免费高清国语| 精品国内亚洲2022精品成人| 日本与韩国留学比较| 亚洲精品影视一区二区三区av| 免费少妇av软件| 国产黄色免费在线视频| 精品久久久久久久末码| 亚洲精华国产精华液的使用体验| 久久久久久久久大av| 人妻少妇偷人精品九色| 三级男女做爰猛烈吃奶摸视频| 亚洲av福利一区| 99热这里只有是精品在线观看| 看十八女毛片水多多多| 久久久精品94久久精品| 午夜久久久久精精品| 中文天堂在线官网| 六月丁香七月| 亚洲无线观看免费| 亚洲自偷自拍三级| 中文字幕亚洲精品专区| av网站免费在线观看视频 | 80岁老熟妇乱子伦牲交| 婷婷六月久久综合丁香| 国产精品国产三级国产专区5o| 亚洲色图av天堂| 一级片'在线观看视频| 干丝袜人妻中文字幕| 亚洲成人av在线免费| 国产精品三级大全| 久久精品国产鲁丝片午夜精品| 两个人的视频大全免费| 一二三四中文在线观看免费高清| 一级毛片aaaaaa免费看小| 97人妻精品一区二区三区麻豆| 日韩欧美精品免费久久| 日韩国内少妇激情av| 国产欧美另类精品又又久久亚洲欧美| 男人狂女人下面高潮的视频| 国内少妇人妻偷人精品xxx网站| 最近中文字幕高清免费大全6| 高清在线视频一区二区三区| 亚洲熟女精品中文字幕| 能在线免费观看的黄片| 亚洲最大成人中文| 亚洲国产欧美在线一区| 99热这里只有精品一区| 日本与韩国留学比较| 亚洲va在线va天堂va国产| 看黄色毛片网站| 草草在线视频免费看| 成年人午夜在线观看视频 | a级毛片免费高清观看在线播放| 99久国产av精品| 国产午夜福利久久久久久| 视频中文字幕在线观看| 免费黄频网站在线观看国产| 最新中文字幕久久久久| 97精品久久久久久久久久精品| 色综合色国产| 国产一区二区亚洲精品在线观看| 最近最新中文字幕大全电影3| 国产三级在线视频| 日本一二三区视频观看| 久久久久久久久久久免费av| 久久久精品94久久精品| 又黄又爽又刺激的免费视频.| av播播在线观看一区| 能在线免费看毛片的网站| 中国美白少妇内射xxxbb| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成人亚洲精品av一区二区| 三级国产精品欧美在线观看| 18禁在线无遮挡免费观看视频| 国产伦精品一区二区三区四那| 欧美bdsm另类| 欧美3d第一页| 久久精品人妻少妇| 日韩亚洲欧美综合| 国产精品一区二区在线观看99 | 一级爰片在线观看| 国产伦精品一区二区三区视频9| 欧美xxxx性猛交bbbb| 三级经典国产精品| 免费观看性生交大片5| 可以在线观看毛片的网站| 哪个播放器可以免费观看大片| 亚洲国产精品成人久久小说| 精品不卡国产一区二区三区| 一夜夜www| 成人一区二区视频在线观看| 亚洲av成人精品一二三区| 最近2019中文字幕mv第一页| 国内精品美女久久久久久| 久久99热这里只频精品6学生| 黄色欧美视频在线观看| 99热网站在线观看| 国产91av在线免费观看| 日本与韩国留学比较| 在线 av 中文字幕| 久久久久久久久久人人人人人人| 欧美性感艳星| 免费人成在线观看视频色| 精品久久久久久久久亚洲| 国产片特级美女逼逼视频| 大香蕉久久网| 青青草视频在线视频观看| 国产日韩欧美在线精品| 国产成人91sexporn| 建设人人有责人人尽责人人享有的 | 一本久久精品| 少妇熟女aⅴ在线视频| 国产 亚洲一区二区三区 | 三级毛片av免费| 你懂的网址亚洲精品在线观看| 久久99热这里只频精品6学生| 国产伦精品一区二区三区四那| 美女国产视频在线观看| 亚洲图色成人| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情福利司机影院| 久久久久久九九精品二区国产| 国产成人a∨麻豆精品| 久久久欧美国产精品| 色综合亚洲欧美另类图片| 99热这里只有是精品在线观看| 永久网站在线| 国产亚洲av嫩草精品影院| 国产一区二区三区综合在线观看 | 一级毛片 在线播放| 建设人人有责人人尽责人人享有的 | 久久久久性生活片| av网站免费在线观看视频 | 国产综合懂色| 人人妻人人澡欧美一区二区| 好男人在线观看高清免费视频| 欧美97在线视频| 国产av在哪里看| 2021少妇久久久久久久久久久| 大香蕉久久网| 欧美zozozo另类| 精品久久久久久久人妻蜜臀av| 国产精品一及| 精品99又大又爽又粗少妇毛片| 只有这里有精品99| 国产精品美女特级片免费视频播放器| 麻豆久久精品国产亚洲av| 国产精品久久久久久av不卡| 亚洲色图av天堂| 乱人视频在线观看| a级一级毛片免费在线观看| 国产 亚洲一区二区三区 | 国产单亲对白刺激| 汤姆久久久久久久影院中文字幕 | 欧美不卡视频在线免费观看| 欧美性猛交╳xxx乱大交人| 免费观看精品视频网站| 美女黄网站色视频| 久久久久久久国产电影| 中文欧美无线码| 国产精品美女特级片免费视频播放器| 成人毛片60女人毛片免费| 男女视频在线观看网站免费| 成年版毛片免费区| 一夜夜www| 国产午夜精品一二区理论片| 淫秽高清视频在线观看| 青春草视频在线免费观看| 国产伦精品一区二区三区视频9| 天堂俺去俺来也www色官网 | 人妻少妇偷人精品九色| 欧美三级亚洲精品| 禁无遮挡网站| 精品熟女少妇av免费看| 欧美xxxx黑人xx丫x性爽| 熟妇人妻久久中文字幕3abv| 日本爱情动作片www.在线观看| 波多野结衣巨乳人妻| 国产精品国产三级国产av玫瑰| 亚洲国产成人一精品久久久| 美女大奶头视频| 亚洲欧美一区二区三区国产| 日韩欧美 国产精品| 欧美xxxx黑人xx丫x性爽| 国产白丝娇喘喷水9色精品| 免费看美女性在线毛片视频| 国产成人午夜福利电影在线观看| 亚洲精品国产成人久久av| 久久久精品免费免费高清| 免费不卡的大黄色大毛片视频在线观看 | av网站免费在线观看视频 | 老司机影院毛片| 在线观看美女被高潮喷水网站| 80岁老熟妇乱子伦牲交| 色综合站精品国产| 亚洲精华国产精华液的使用体验| 搡老乐熟女国产| 99久久精品一区二区三区| 国产男人的电影天堂91| 国产 一区 欧美 日韩| 精品久久久久久久久亚洲| 男人狂女人下面高潮的视频| 亚洲欧美精品专区久久| 成年免费大片在线观看| 97人妻精品一区二区三区麻豆| 亚洲精品成人av观看孕妇| 久久热精品热| 久久人人爽人人片av| 91午夜精品亚洲一区二区三区| 午夜精品在线福利| av在线蜜桃| 一本一本综合久久| 肉色欧美久久久久久久蜜桃 | 真实男女啪啪啪动态图| 国产精品人妻久久久影院| 国产亚洲av片在线观看秒播厂 | 成年免费大片在线观看| 亚洲欧美日韩无卡精品| 久久久久久久大尺度免费视频| 人人妻人人澡欧美一区二区| 国产精品一二三区在线看| 人人妻人人澡欧美一区二区| 欧美丝袜亚洲另类| 国产精品国产三级专区第一集| 51国产日韩欧美| 亚洲精品视频女| 久久久久精品性色| 亚洲精品视频女| 久久久久性生活片| 自拍偷自拍亚洲精品老妇| 国产乱来视频区| 只有这里有精品99| 亚洲精品456在线播放app| 人人妻人人澡欧美一区二区| 欧美3d第一页| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人精品欧美一级黄| 国产精品无大码| 国产精品久久久久久久久免| 3wmmmm亚洲av在线观看| 久久久久久久大尺度免费视频| 日韩精品有码人妻一区| 成人亚洲精品av一区二区| 日本色播在线视频| 欧美高清性xxxxhd video| 国内揄拍国产精品人妻在线| av线在线观看网站| 免费电影在线观看免费观看| 免费人成在线观看视频色| 中文精品一卡2卡3卡4更新| 欧美一级a爱片免费观看看| 国产国拍精品亚洲av在线观看| 国产一区有黄有色的免费视频 | 亚洲美女搞黄在线观看| 精品国产三级普通话版| 欧美成人午夜免费资源| 久久99热这里只有精品18| 高清毛片免费看| 亚洲精品成人久久久久久| 一级av片app| 欧美日韩精品成人综合77777| 欧美成人精品欧美一级黄| 欧美日韩在线观看h| 成人亚洲精品av一区二区| 黄片wwwwww| 三级国产精品片| 免费观看a级毛片全部| 亚洲国产成人一精品久久久| 亚洲伊人久久精品综合| 午夜精品一区二区三区免费看| 国产亚洲91精品色在线| 久久亚洲国产成人精品v| 偷拍熟女少妇极品色| 久久久久国产网址| 能在线免费看毛片的网站| av.在线天堂| 久久久色成人| 午夜日本视频在线| 国产v大片淫在线免费观看| 26uuu在线亚洲综合色| av线在线观看网站| 五月玫瑰六月丁香| 91精品一卡2卡3卡4卡| 免费看美女性在线毛片视频| 在线播放无遮挡| 久久久久久国产a免费观看| 亚洲无线观看免费| 免费观看无遮挡的男女| 国产男女超爽视频在线观看| 亚洲国产色片| 自拍偷自拍亚洲精品老妇| 麻豆av噜噜一区二区三区| av国产免费在线观看| 别揉我奶头 嗯啊视频| 日韩大片免费观看网站| 成人欧美大片| 青青草视频在线视频观看| 欧美激情久久久久久爽电影| 禁无遮挡网站| 亚洲激情五月婷婷啪啪| 成人亚洲精品av一区二区| 亚洲国产精品专区欧美| 免费人成在线观看视频色| 91久久精品国产一区二区三区| 高清视频免费观看一区二区 | 免费无遮挡裸体视频| 国语对白做爰xxxⅹ性视频网站| 日本黄大片高清| 最近的中文字幕免费完整| 免费看a级黄色片| 18禁动态无遮挡网站| 亚洲国产精品专区欧美| 亚洲av中文字字幕乱码综合| 一个人免费在线观看电影| 国产高清不卡午夜福利| 亚洲成人久久爱视频| 日本与韩国留学比较| 日本三级黄在线观看| 国产欧美另类精品又又久久亚洲欧美| 在线播放无遮挡| 成人午夜高清在线视频| 国产免费福利视频在线观看| 久久久久性生活片| 99久久精品热视频| 精品一区二区三卡| 一个人免费在线观看电影| h日本视频在线播放| 五月天丁香电影| 日产精品乱码卡一卡2卡三| 蜜桃久久精品国产亚洲av| 97超视频在线观看视频| 亚洲精华国产精华液的使用体验| 女人十人毛片免费观看3o分钟| 在线天堂最新版资源| 熟女人妻精品中文字幕| 少妇丰满av| 老师上课跳d突然被开到最大视频| 天堂√8在线中文| av在线亚洲专区| 国产av国产精品国产| 国产熟女欧美一区二区| 人妻系列 视频| 日本色播在线视频| 成人美女网站在线观看视频| or卡值多少钱| 亚洲经典国产精华液单| 亚洲国产色片| 久久99蜜桃精品久久| 精品人妻一区二区三区麻豆| 亚洲av成人av| 国产伦在线观看视频一区| 国精品久久久久久国模美| 国产男女超爽视频在线观看| 国产美女午夜福利| 亚洲美女视频黄频| 又黄又爽又刺激的免费视频.| 我要看日韩黄色一级片| 看十八女毛片水多多多| 男人爽女人下面视频在线观看| 亚洲欧美精品专区久久| 亚洲美女搞黄在线观看| 97在线视频观看| 亚洲在久久综合| 国模一区二区三区四区视频| 午夜日本视频在线| 中文字幕免费在线视频6| 久久久国产一区二区| 亚洲,欧美,日韩| 最后的刺客免费高清国语| 日日啪夜夜撸| 久久久久久久久久久丰满| 精品久久国产蜜桃| 黄色一级大片看看| 精品人妻熟女av久视频| 国产精品av视频在线免费观看| 免费观看精品视频网站| 久久6这里有精品| av.在线天堂| 亚洲一区高清亚洲精品| 国产精品国产三级国产av玫瑰| 亚洲无线观看免费| 插阴视频在线观看视频| av播播在线观看一区| 亚洲天堂国产精品一区在线| 欧美一区二区亚洲| 亚洲电影在线观看av| a级毛色黄片| 丰满人妻一区二区三区视频av| 日韩在线高清观看一区二区三区| 欧美日韩国产mv在线观看视频 | 日本黄大片高清| 成人二区视频| av又黄又爽大尺度在线免费看| 搞女人的毛片| 免费观看的影片在线观看| 亚洲乱码一区二区免费版| 国产综合懂色| 日韩成人av中文字幕在线观看| 亚洲第一区二区三区不卡| a级毛色黄片| 亚洲av免费在线观看| 高清视频免费观看一区二区 | 黄片无遮挡物在线观看| 国产高清三级在线| 一夜夜www| 亚洲av日韩在线播放| 中文字幕亚洲精品专区| 精品久久久久久久久av| 成年版毛片免费区| 国产黄色视频一区二区在线观看| 搡女人真爽免费视频火全软件| 国产伦精品一区二区三区四那| 亚洲高清免费不卡视频| 久久久亚洲精品成人影院| 99热这里只有精品一区| 搡老乐熟女国产| 国产男女超爽视频在线观看| 成年免费大片在线观看| 最近最新中文字幕免费大全7| 日韩在线高清观看一区二区三区| 美女cb高潮喷水在线观看| 2021天堂中文幕一二区在线观| 亚洲人成网站高清观看| 嫩草影院新地址| 国产高潮美女av| 久久精品国产鲁丝片午夜精品| or卡值多少钱| 成年女人看的毛片在线观看| 国产69精品久久久久777片| 午夜激情久久久久久久| 欧美不卡视频在线免费观看| 最近2019中文字幕mv第一页| 国产精品一及| av在线亚洲专区| 中文乱码字字幕精品一区二区三区 | 在线 av 中文字幕| 26uuu在线亚洲综合色| 久久久欧美国产精品| 18禁裸乳无遮挡免费网站照片| 亚洲熟女精品中文字幕| 国产免费又黄又爽又色| 在线观看免费高清a一片| 成人毛片a级毛片在线播放| 伊人久久国产一区二区| 深爱激情五月婷婷| 免费观看在线日韩| 亚洲久久久久久中文字幕| 国产一区二区亚洲精品在线观看| 欧美区成人在线视频| 建设人人有责人人尽责人人享有的 | 国产又色又爽无遮挡免| 丰满少妇做爰视频| 汤姆久久久久久久影院中文字幕 | 国产不卡一卡二| 99久久精品一区二区三区| 久久久久九九精品影院| 免费观看a级毛片全部| 久久久精品免费免费高清| 最近手机中文字幕大全| 国产一区二区三区av在线|