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

    Salvianolic acid B dry powder inhaler for the treatment of idiopathic pulmonary fibrosis

    2022-09-06 13:50:22

    aState Key Laboratory of Component-Based Chinese Medicine,Tianjin University of Traditional Chinese Medicine,Tianjin 301617,China

    bEngineering Research Center of Modern Chinese Medicine Discovery and Preparation Technique,Ministry of Education,Tianjin 301617,China

    cHaihe Laboratory of Modern Chinese Medicine,Tianjin 301617,China

    dCollege of Chinese Medicine,Tianjin University of Traditional Chinese Medicine,Tianjin 301617,China

    eDepartment of Pharmacy,Suzhou TCM Hospital Affiliated to Nanjing University of Chinese Medicine,Suzhou,Jiangsu 215009,China

    fEndocrine and Metabolic Disease Center,The First Affiliated Hospital,and College of Clinical Medicine of Henan University of Science and Technology,Medical Key Laboratory of Hereditary Rare Diseases of Henan,Luoyang Sub-center of National Clinical Research Center for Metabolic Diseases,Luoyang,Henan 471003,China.

    Keywords:Salvianolic acid B Dry powder inhaler Idiopathic pulmonary fibrosis Pulmonary administration

    ABSTRACT Idiopathic pulmonary fibrosis (IPF) is a serious and fatal pulmonary inflammatory disease with an increasing incidence worldwide.The drugs nintedanib and pirfenidone,are listed as conditionally recommended drugs in the“Evidence-Based Guidelines for the Diagnosis and Treatment of Idiopathic Pulmonary Fibrosis”.However,these two drugs have many adverse reactions in clinical application.Salvianolic acid B (Sal B),a water-soluble component of Salvia miltiorrhiza,could alleviate bleomycin-induced peroxidative stress damage,and prevent or delay the onset of IPF by regulating inflammatory factors and fibrotic cytokines during the disease’s progression.However,Sal B is poorly absorbed orally,and patient compliance is poor when administered intravenously.Therefore,there is an urgent need to find a new non-injection route of drug delivery.In this study,Sal B was used as model drug and L-leucine (LL) as excipient to prepare Sal B dry powder inhaler (Sal B-DPI) by spray drying method.Modern preparation evaluation methods were used to assess the quality of Sal B-DPI.Sal B-DPI is promising for the treatment of IPF,according to studies on pulmonary irritation evaluation,in vivo and in vitro pharmacodynamics,metabolomics,pharmacokinetics,and lung tissue distribution.

    1.Introduction

    Idiopathic Pulmonary Fibrosis (IPF) is a common alveolar disease in which normal lung tissue is replaced by excessive interstitial cells and extracellular matrix,resulting in the destruction of alveolar structures,the gradual densification of the alveolar region,and ultimately respiratory failure and even death [1-6].The age of onset of IPF is mostly between 50 and 70 years old.With the aging of the population,the incidence of IPF worldwide is increasing year by year.According to the statistics of the World IPF Joint Association,there are about 3.2 million people suffering from IPF and 1.22 million new cases every year [6].Studies have shown that environmental factors (inhalation of dust,metal powder and wood chips in the air),viral infection(influenza virus,hepatitis C virus,etc.),gastroesophageal reflux and smoking are the main potential pathogenic factors for IPF.Since December 2019,there has been an outbreak of Coronavirus disease 2019(COVID-19) in Wuhan,Hubei Province.Patients with severe COVID-19 may have left a certain degree of changes in lung damage repair for a period of time,such as the IPF [7-9].The autopsy results of the first case of COVID-19 showed severe pulmonary fibrosis,although the pulmonary fibrosis and parenchymal lesions were not as severe as those caused by the SARS virus [9].Current studies have shown that the pathological changes were due to the combined effects of oxidative stress,inflammatory response and fibrosis[2-5].

    The commonly used drugs for clinical treatment of IPF mainly include tyrosine kinase inhibitors,cytokine inhibitors,glucocorticoids,collagen synthesis inhibitors,immunosuppressants,phosphodiesterase inhibitors,anticoagulants,endothelin antagonists and antioxidants [2].There was no definitive and effective treatment other than lung transplantation,and without lung transplantation,the mortality rate for IPF patients was 50% at 3 years and 80% at 5 years [2-5],respectively.However,lung transplantation had problems such as scarcity of lung donor,high risk and cost.Currently,in the evidence-based guidelines for the diagnosis and treatment of IPF,only nitanib,a tyrosine kinase inhibitor,and pirfenidone,a cytokine inhibitor,are listed as conditional recommended drugs[10].However,the above two drugs have large clinical dosage,high treatment cost,and are prone to drug resistance and more adverse reactions (diarrhea,rash,insomnia and dizziness,etc.) [11-12].Therefore,it is of great economic value and social significance to continue to search for safe and effective drugs for the treatment of IPF.

    Traditional Chinese medicine (TCM) is an integral part of modern medicine.It is also an important resource of natural medicines and plays a broader role in the treatment of human diseases[13].In recent years,there have been many studies on TCM prevention or treatment of IPF[13-14].Salvia miltiorrhiza is one of herbal medicines commonly used in the treatment of IPF [15].Its water-soluble component,salvianolic acid,is the main active component to promote blood circulation and remove blood stasis.It not only has a protective effect on the cardiovascular system,but also plays an important role in the treatment of pulmonary diseases such as IPF and chronic obstructive pulmonary disease[14].Salvianolic acid B(Sal B),the main component of salvianolic acid,could relieve bleomycin-induced peroxide-stress injury and regulate Tumor necrosis factor-α(TNF-α) and other inflammatory cytokines and fibrosis related cytokines to prevent or treat pulmonary fibrosis [16-17].However,Sal B is poorly absorbed orally and rapidly eliminatedin vivo.It has obvious liver firstpass effect and hepatointestinal circulation,and its absolute bioavailability is between 0.02% and 5.56% [18-20].The way of injection administration can make the drugs directly enter the body to play its role,and the therapeutic effect is very good.However,it needs long-term and frequent intravenous administration to maintain the therapeutic effect,which leads to poor patient compliance due to the inconvenience of the way of use.Therefore,there is an urgent need to find a new non-injection route of drug delivery.

    Pulmonary delivery as a non-injectable drug delivery route can target drug delivery to the lungs,avoiding low systemic toxicity due to intestinal absorption and liver metabolism[21].The lungs have a large drug absorption surface area (about 100-140m2),relatively high blood flow,low enzyme activity and thin epithelial cells,which have unique advantages for the treatment of lung diseases[22].Compared with nebulizer(NEB)and metered dose inhaler(MDI),dry powder inhaler(DPI)has the advantages of good drug stability,convenient use and carrying,no need for ejectors,larger drug delivery doses,and good patient compliance[23].At present,the preparation methods of DPI mainly include grinding method,spray freeze drying method,spray drying method,airflow crushing method and supercritical fluid technology [24].DPI mainly includes carrier type and carrier-free type,among which lactose is the only carrier approved by the Food and Drug Administration(FDA)for the preparation of it.At present,lactose with various physicochemical properties in the market can be used for the preparation of DPI with good safety.However,lactose as a carrier has certain limitations (some drugs may have strong absorption of lactose,and lactose is easy to deposit in the oropharynx during administration,resulting in a very small amount of drugs entering the respiratory tract;longterm use of lactose-containing DPI may cause dental caries,etc.) and it is suitable for low doses of drugs.The deposition of drug particles of DPI without carrier type in the lung is not affected by the separation of drug and carrier surface,and can avoid the potential toxic and side effects of carrier on the body,so it is suitable for high dose drugs [25-26].Among them,L-leucine (LL) is one of the most widely used excidients,which is enriched on the surface of drug particles,reduces the interparticle force,improves the fluidity of dry powder,and improves the dispersion and atomization properties of dry powder [26-28].In addition,some studies have shown that it can enhance the anti-hygroscopicity of DPI to improve the stability of the drug [26-28].At present,the preparation of DPI using LL as excidient by spray drying technology has been widely used [26-28].In this study,Sal B was used as the model drug and LL was used as the exciphant to prepare Sal B-DPI by spray drying.The quality of Sal B-DPI was evaluated by modern preparation evaluation methods.The lung irritation evaluation,in vivoandin vitropharmacodynamics,metabolomics,pharmacokinetics and lung tissue distribution were studied to prove the feasibility and rationality of lung administration of Sal B-DPI in the treatment of IPF.

    2.Materials and methods

    2.1.Materials

    Sal B (Lot#:P18J9F65817,purity 98.0%),Lleucine(Lot#:H15J7H17735) were purchased from Shanghai Yuanye Biotechnology Co.,Ltd (Shanghai,China).Transforming growth factor-β1 (TGF-β1) was purchased from R&D System (USA).DMEM medium and fetal bovine serum were purchased from Gibco (New York,USA).ⅠCollagen typeα1 (COL1A1),ⅢCollagen typeα1 (COL3A1),Intercellular adhesive molecule-1,ICAM-1),Arginase type 1 (Arg-1),Fibronectin (FN),Inducible nitric oxide synthase(iNOS) primers were purchased from Shanghai Sangon Bioengineering Co.,Ltd.(Shanghai,China).IL-8 ELISA kit was purchased from Shanghai Sinovac Biotechnology Co.,Ltd.(Shanghai,China).Total protein quantitative,LDH,SOD activity and MDA assay kit were purchased from Nanjing Jiancheng Biotechnology Co.,Ltd.(Nanjing,China).Bleomycin was purchased from Nippon Chemical Co.,Ltd.Vegat(nidanib sulfonic acid soft capsule)(C31H33N5O4,150 mg/capsule,Lot#:900,678) from Catalent Germany Eberbach GmbH (Germany).IL-1β,IL-6,IL-4,IL-18 and TGF-β1 ELISA kits were purchased from Wuhan Bode Bio-engineering Co.,Ltd.(Wuhan,China).IFN-γ,MPO and HYP ELISA kits were purchased from Sinovac Biotech Co.,Ltd.(Shanghai,China).Protocatechuic aldehyde(PA),used as internal standard (IS),(Lot#:w06-2-5,purity 98.0%) were purchased from Zhongxin Pharmaceutical Research Center (Tianjin,China).Heparin sodium (batch no.325D024,titer ≥140 units/mg)was purchased from Suolai Bao Technology Co.Ltd (Beijing,China).Formic acid (batch no.K3JHG.LS,purity 98%)was purchased from Chemical Industry Co.Ltd(Tokyo,Japan).HPLC-grade acetonitrile,methanol and phosphoric acid(Lot#:130,271) were purchased from Fisher Scientific Inc.(Fairlawn,USA).Analytical-grade hydrochloric acid was purchased from Jiangtian Chemical Technology Co.Ltd (Tianjin,China).Ascorbic acid (batch no.102045-201,504,purity 98%) was purchased from the National Institutes for Food and Drug Control (Beijing,China).HPLC-grade water was purified using a Milli-Q Reagent Water system(Millipore,Milford,MA,USA).

    2.2.Cell lines and animals

    NIH-3T3 mouse embryonic fibroblast was obtained from American Type Culture Collection (Maryland,USA).The cells were cultured in high glucose DMEM medium containing 10% fetal bovine serum,100 mg/l streptomycin and 100 kU/l penicillin,pH 7.2-7.4.The cells were incubated at 37°C with 5.0% CO2atmosphere.Male Sprague-Dawley rats(weighing 170-200 g;Beijing Weitong Lihua Experimental Animal Technology Co.,Ltd,Beijing,China;license SCXK 2016-0006) were kept in an environmentally controlled breeding room for 1 week before the experiments.Rats were fasted for 12 h and supplied with unlimited water prior to the experiment.The protocol and any amendments or procedures involving the care or use of animals in this study were in accordance with the regulations for animal experimentation issued by the State Committee of Science and Technology of China and approved by TJUTCM’s Institutional Animal Care and Use Committee(document number TCM-LAEC 2,020,055).

    2.3.Production of powder formulations by spray dryer

    Feed solutions were prepared by dissolving Sal B and LL in water with known mass ratios of Sal B to LL.The feed concentration(total solids content)was maintained at 1.50%.A B-90 nano spray-dryer (Büchi Labortechnik AG,Falwil,Switzerland) was operated under the following parameters[29]: inlet temperature of 80 °C,air flow of 100 l/min,pump speed of 23%,spray rate of 60%,and internal pressure of 34 Mbar.After the above parameters were stabilized,the spray would begin.The spray-dried powders were stored in a desiccator containing silica gel at room temperature without light until used.

    2.4.Quality evaluation of spray-dried powders

    2.4.1.Particle sizing of spray-dried powders

    The volumetric diameter and span of the sample powders were analyzed by a laser particle size analyzer(HELOS(Hi214)&OASISDRY,R1;Sympatec GmbH,Goslar,Germany) [30].The parameters of D10,D50,D90and volume surface mean diameter were separately determined,representing the sizes at which 10%,50%,and 90%of the particles were smaller than the remaining particles and the mean volume diameter.

    2.4.2.Scanning electron microscopy(SEM)

    SEM (JSM-7500F;JEOL Int.,Tokyo,Japan) was employed to examine particle morphology of the powder formulations at 10 kV.The samples were coated with a layer of gold and dried in a vacuum[26].

    2.4.3.X-ray powder diffraction(XRD)

    Powder crystallinity was evaluated by XRD (D/MAX-2550 V;Rigaku Int.,Tokyo,Japan),using Cu Ka radiation at a voltage of 40 kV and a current of 30 mA.Data were recorded from 5°-50°(2θ)using a scan speed of 2 per min and a step size of 0.02[26].

    2.4.4.Dynamic vapor sorption(DVS)

    A DVS instrument (Surface Management Systems,London,United Kingdom) was used to derive the moisture sorption isotherms of sample powders.The chamber was kept at 25°C under continuous nitrogen flow.30-50 mg of each sample was dried at 0 RH for 1 h before being exposed to one cycle of 0-90% RH,with 10% RH increments.After the sample gradually reached the moisture absorption balance at each relative humidity,the next level of relative humidity test would be automatically performed[26].

    2.4.5.Content determination

    The content of Sal B in spray-dried powders of various mass ratios were determined by HPLC according to the chromatographic conditions in the Supplementary material.

    2.4.6.In vitro aerosol performances

    A next-generation pharmaceutical impactor (NGI-1302;Copley Scientific,Nottingham,UK) was used to evaluate the aerosolization behavior of spray-dried powders with a United States Pharmacopoeia (USP) induction port.An ethanol solution containing a certain amount of benzyl and glycerin was added to the surface of each collection cup and then evaporated at ambient temperature for 2 h to avoid particle bounce or re-entrainment.Spray-dried powders were loaded into an Easyhaler multiple dose powder inhaler (Orion Int.,Espoo,Finland).The test was carried out (4 s,60 l/min flow rate).Five shots were fired into Next Generation Impactor(NGI) for each formulation.Thoroughly rinse all stages with a predetermined volume of water after aerosolization.The collected solutions were analyzed by HPLC as content determination.

    The parameters of aerosol properties were calculated by CITDAS software (version 3.00;Copley Scientific): fine particle dose (FPD),fine particle fraction (FPF),mass median aerodynamic diameter (MMAD) and geometric standard deviation (GSD).FPD refers to the dose of drug particles with aerodynamic diameter of<5 μm,whereas FPF is defined as the ratio of FPD to the total recovery dose [31].In addition,delivery efficiency refers to the ratio of the amount of powder collected in each part of the NGI to the total dose.

    2.5.Preliminary stability investigation

    The Sal B-DPI(80SalB20LL) was stored in a desiccator containing silica gel at room temperature without light.The temperature and humidity corresponding to this condition are 20±2 °C and 20% ± 2%,respectively.Samples were taken on Day 1 and 35 for particle size distribution,morphology,XRD analysis,content determination andin vitroaerodynamics analysis.

    2.6.Anti-pulmonary fibrosis in vitro

    2.6.1.MTT assay

    MTT assay was used to verify whether these isolates had toxic side effects on normal fibroblasts and the effect on vitality of fibroblasts[32].NIH-3T3 mouse embryonic fibroblasts(Normal fibroblasts)were plated in 96-well plates,and each compound was set to a concentration gradient of nidanib (0.25,0.5,1,2.5,5 and 10 μM) and Sal B-DPI(80SalB20LL) (25,50,100,125,150,175 and 200 μg/ml).After incubation for 48 h,20 μl MTT solution(5 mg/ml)was added to each well,and the culture was continued for 4 h.After that,110 μl DMSO was added to each well,and the absorbance value of each well was measured at 490 nm with a microplate analyzer.

    2.6.2.Sirius red method

    The effect of Sal B-DPI on collagen deposition in NIH-3T3 cells induced by TGF-β1 was determined by Sirius red method [33].The optimal stimulating dose of TGF-β1 was 5 ng/ml.SB431542,a selective inhibitor of TGF-β1 signaling,was used as a positive control at a concentration of 5 μM.Each compound was set to a concentration gradient of nidanib(0.25 μM and 0.5 μM)and Sal B-DPI(25,50 and 125 μg/ml).After incubation for 48 h,Sirius red 100 μl was added to each well and kept away from light at room temperature for 4 h.The detection was performed at 540 nm.

    2.6.3.Scratch test

    NIH-3T3 mouse embryonic fibroblasts were coated with 12-well plates and cultured for 24 h and then scratched [34].After 24 h of synchronization,5 ng/ml TGF-β1 and drugs(SB431542,nidanib,Sal B-DPI) were added and incubated for 48 h.Photographs were taken with an inverted microscope at 0,24 and 48 h after drug addition,respectively.

    2.6.4.qPCR detection[31-34]

    NIH-3T3 mouse embryonic fibroblasts were coated with 6-well plates and cultured for 24 h.5 ng/ml TGF-β1 and Sal BDPI(50 μg/ml) were added and incubated for 48 h.RNA was extracted from cells for qPCR,and the expressions of COL1A1,COL3A1,ICAM-1,FN,iNOS and Arg-1 were detected to analyze the activation degree of NIH-3T3.Primer sequences are shown in Table S1.

    2.7.Pulmonary irritation

    2.7.1.Grouping and administration

    18 experimental rats were fed adaptive for 3 d The rats were divided into 3 groups by random number expression method:normal control group(NC),sham operation group(sham),and pulmonary administration of Sal B-DPI(80SalB20LL)group(Sal B-DPI),with 6 rats in each group.SD rats in Sal B-DPI group were anesthetized and fixed by intraperitoneal injection of 10%chloral hydrate(2 ml/kg).SD rats were administrated with 10 mg/kg of Sal B-DPI by the method described below[35-36].The method was as follows: Sal B-DPI was loaded into a DP-4R dry powder injector and equipped with an SC-X sample chamber dilator (Penn Pancentury,Inc.,Wyndmoor,PA).The delivery tube of the inhaler had been carefully introduced into the trachea before the first trachea bifurcated.1.5 ml air corresponding to the tidal volume of SD rats was pumped to Sal B-DPI at a relatively steady rate for three consecutive times.The sham group was administered in the same manner as Sal B-DPI group,except that no Sal B-DPI was administered.

    2.7.2.Sample collection and preservation

    The sample collection method is optimized based on the previous research,especially the collection of broncho alveolar lavage fluid(BALF) [14,16].After weighing,the rats in the three groups were anesthetized by intraperitoneal injection of 10% chloral hydrate (2 ml/kg).The rats were sacrificed after blood collection through abdominal aorta.After the tracheal branch was found,the left lung of the rats was ligated with a hemostatic clip.A small incision was cut near the tracheal branch,and a plastic tube was inserted into the right bronchus along the incision,which was immediately fixed with a sterile surgical thread,and the lavage operation was ready to begin.2 ml normal saline at 4°C was taken with a sterile syringe and slowly injected into the right lung of rats through the three-way valve of the retained needle.The normal saline remained in the lung tissue for 1 min,and the lavage fluid was slowly withdrawn.After repeated lavage for 3 times,the lavage fluid was collected and mixed in ice bath.After lavage,the lung tissue was separated,the surface of the lung tissue was rinsed with normal saline and dried with filter paper.The left lung of the rats was soaked in 10% neutral formalin solution and fixed at room temperature for 48 h,then used for pathological observation of lung tissue.

    2.7.3.Determination of total protein,IL-8 concentration and lactate dehydrogenase(LDH)activity in BALF

    Total protein,IL-8 and LDH activity in the supernatant of BALF in each group were quantitatively determined by BCA method,enzymatic linked immunosorbent assay(ELISA)and microplate method,respectively,according to the kit instructions.

    2.8.Preliminary pharmacodynamics study in vivo

    2.8.1.Rat model of IPF

    After the experimental rats were anesthetized by intraperitoneal injection of 10% chloral hydrate (2 ml/kg),0.1 ml bleomycin saline solution (5 mg/kg) was injected into the trachea through the space of the tracheal cartilage ring into the rat’s lungs.Immediately place the rat upright and rotate it for 5 min to make the solution evenly distributed in the rat’s lungs[14,16].

    2.8.2.Grouping and administration

    30 experimental rats were fed adaptively for 3 d They were divided into 5 groups by random number expression,normal control group (NC),bleomycin group (bleo),and pulmonary administration of Sal B-DPI (80SalB20LL) group (bleo+Sal BDPI),intravenous injection of Sal B solution group (bleo+Sal B iv) and intragastric administration of positive drug nidanib group(bleo+nidanib).

    The NC group was injected with equal volume of normal saline in the way of tracheal non-invasive pulmonary administration,and the other groups were injected with bleomycin normal saline solution to induce the rat model of IPF.The model was established as the 0 d of the experiment.Each treatment group was given drug intervention on the 1st d.The bleo+Sal B-DPI group was given Sal B-DPI (10 mg/kg,measured by SalB)by the method described in“2.7.1′′for 28 d,once a day.The bleo group was pushed into air as described above.Rats in the bleo+Sal B iv group were injected with Sal B solution at the same dose as those in the bleo+Sal B-DPI group.The positive drug nidanib was administrated intragastrically with a dose of 31.5 mg/kg (measured by nidanib).During administration,activities,breathing,feeding,hair and body weight of rats in each group were observed carefully every day.

    2.8.3.Non-invasive lung function measurement

    On Day 29,the system of drug administration and respiratory physiological detection for awake animals was started.After the equipment was calibrated,rats in each group were placed into the animal plysmography chamber in batch.When the rats were in normal state and breathing was stable,various indicators of lung function including tidal volume(TV),minute volume (MV).After the measurement,further specimen collection was carried out.TV: The volume of air inhaled or exhaled with each breath is an indicator of lung volume of lung function in rats;MV: The total amount of gas inhaled or exhaled per minute,which is also the product of TV and respiratory frequency,is one of the ventilation indicators of lung function in rats.

    2.8.4.Sample collection and preservation

    On Day 29,rats in each group were sacrificed,and serum,BALF and lung tissue were collected.After weighing,the rats in the three groups were anesthetized by intraperitoneal injection of 10%chloral hydrate(2 ml/kg).The rats were sacrificed after blood collection through abdominal aorta.The procedure for collecting BALF was the same as “2.7.2′′.After lavage,the lung tissue was separated,the surface of the lung tissue was washed with ice normal saline,and the water on the surface of the lung tissue was dried with filter paper.Part of the left lung of the rats was soaked in 10%neutral formalin solution,fixed at room temperature for 48 h,and then used for staining of the lung tissue sections.Part of the left lung was stored at-80 °C for the determination of Hyp,SOD activity and MDA content[37].

    2.9.Metabonomics study of lung tissue in rats

    2.9.1.Processing of lung tissue samples

    According to a weight-to-volume ratio of 1:5 (g/ml),physiological saline was added to the lung tissues of each group in an appropriate amount,and placed in an ice-water bath to prepare a homogenate.Take 200 μl lung tissue homogenate,add 4 times the amount of methanol,vortex for 2 min,and centrifuge at 12 000 rpm at 4 °C for 10 min.Take the supernatant,dry it with nitrogen at room temperature,add 100 μl methanol for reconstitution,and centrifuge at 4°C at 12 000 rpm for 10 min.The supernatant was injected and UPLC-Q-TOF/MS liquid mass spectrometry was used to analyze the difference components in each group.The UPLC-Q-TOF/MS analysis conditions were descripted in the Supplementary material

    2.9.2.Metabonomics data processing[38]

    Import the Raw.data after sample injection and analysis based on QI software (Waters,USA).Based on the software’s automatic identification function,complete basic operations such as feature extraction of batch samples,QC feature standardization,chromatographic peak alignment,peak matching,retention time correction,etc.Export a list of all chromatographic peaks in the form of an Excel table,mainly including: RT,m/z,MS/MS,S/N,standardized peak area and other related parameters.

    Then,on the basis of the original table,analyze the number of“0 values”in each group,and further eliminate and improve irrelevant data according to relevant rules.The candidate differential metabolites were screened based on the VIP value of multivariate statistical analysis,and statistical analysis was carried out by parametric test according to the data type and specific distribution.

    Using the precise molecular weightm/zof high resolution mass spectrometry,it can be found in HMDB database(http://www.hmdb.ca/),Chemspider database (http://www.chemspider.com/),Metlin (http://metlin.scripps.edu) to find possible endogenous components,combined with MS/MS and molecular error (ppm) to determine the final differential molecule.

    In order to illustrate the role of differential metabolites in organisms,this article conducts pathway enrichment analysis for all metabolites on the open source website MetPA(https://www.metaboanalyst.ca/).Screen the biologically significant metabolites and analyze the metabolic pathways where the metabolites accumulate,and find out the potential pathways for Sal B-DPI in the treatment of IPF.

    2.10.Preliminary pharmacokinetics and lung tissue distribution study

    2.10.1.Preliminary pharmacokinetic study

    Twelve male Sprague-Dawley rats (200±20 g) were randomly divided into two groups (n=6 per group).For Sal BDPI group,Sal B-DPI (10 mg/kg,measured as Sal B) was administered to the lung by the method under “2.7.1′′.In another group,the same dose of Sal B solution was intravenously administered.Blood samples (0.5 ml) were immediately collected in heparinized 1.5 ml polythene tubes from the suborbital vein at 0.033,0.083,0.25,0.5,1,2,4,6 and 12 h after administration.Blood samples were immediately centrifuged at 3000 rpm(Thermo Fisher,Massachusetts,USA)for 10 min at 4°C to obtain plasma.All of the plasma samples were kept at-80°C until analysis.

    Each plasma sample (100 μl) was spiked with 10 μl ascorbic acid solution (1.4 mg/ml in water),100 μl IS working solution (0.98 μg/ml PA),20 μl hydrochloric acid(2 mol/l in water) and 1 ml ethyl acetate and briefly vortexmixed for 2 min,centrifuged (12 000 rpm) for 10 min and the supernatant separated.The supernatant was carefully transferred into another tube,dried with nitrogen and diluted with 100 μl methanol-water (80:20,v/v),vortex-mixed for 2 min and centrifuged (12 000 rpm) for 10 min.Finally,3 μl the supernatant was injected into the UPLC-MS/MS system for analysis [39].The specific methods were presented in supplementary materials.The relevant pharmacokinetic parameters were determined using the plasma concentrationtime data and calculated by WinNonlin version 6.4(Pharsight Corporation,Mountain View,CA).The method was developed and validated with respect to specificity,linearity,precision,accuracy,recovery,matrix effect,and stability through the processing procedure and sample storage.

    2.10.2.Lung tissue distribution study

    Fifty-four male Sprague-Dawley rats (200±20 g) were randomly divided into two groups(pulmonary administration group and intravenous administration group).Each group has 9 time points of 0.033,0.083,0.25,0.5,1,2,4,6 and 12 h (n=3 per time point).For Sal B-DPI group,Sal B-DPI(10 mg/kg,measured as Sal B) was administered to the lung by the method under “2.7.1′′.In the another,the same dose of Sal B solution was intravenously administered.After the administration of 0.033,0.083,0.25,0.5,1,2,4,6,and 12 h later,the rats were sacrificed via the abdominal aorta and lung tissue was collected.All of the samples were kept at -80 °C until analysis.

    According to the weight to volume ratio of 1:5 (g/ml),an appropriate amount of lung tissue was added to normal saline,and homogenate was prepared under the condition of ice water bath.The lung tissue homogenate processing method was the same as the plasma sample processing method under“2.10.1′′.The relevant pharmacokinetic parameters were determined using the concentration-time data and calculated by WinNonlin version 6.4 (Pharsight Corporation,Mountain View,CA).The method was developed and validated with respect to specificity,linearity,precision,accuracy,recovery,matrix effect,and stability through the processing procedure and sample storage.

    2.11.Statistical analysis

    Statistical software SPSS 23.0 was used for statistical analysis,and measurement data was expressed as mean ± SD.Oneway ANOVA was used for comparison among groups.P <0.05 indicated significant difference,P <0.01 indicated extremely significant difference.P >0.05 indicated no significant difference.All data were mapped using GraphPad Prism 8 software and Adobe Illustrator CS6 software.

    3.Results and discussion

    3.1.Quality evaluation of spray-dried powders

    3.1.1.Particle size analysis

    The particle size distribution data of powder formulations are summarized in Table 1.Most particles of each formulation were smaller than 5 μm,suitable for DPI formulation [40].The D50values of pray-dried (SD) pure Sal B and pure LL were 1.42 μm and 2.35 μm,while the co-spray dried formulations had fairly similar D50values around 1.10 μm (range 1.00-1.17 μm).

    Table 1-Particle size distributions of various formulations(n=3).

    3.1.2.Particle morphology

    SEM images of various formulations were shown in Fig.1A,which confirmed the particle size data presented in Table 1,with the majority of particles being<5 μm diameter.The surface roughness of SD Sal B particles increased with the amount of LL,most obviously in the sample containing 20% (w/w) LL.The formulations were made up of spherical or nearly spherical particles with a corrugated surface and hollow particles with a smooth surface.This could be due to the enrichment of LL on the surface of the drug particles,thereby reducing the force between them,improving the fluidity,dispersion and atomization performance of DPI [41].The crystallinity of LL in various formulations would be discussed below.

    3.1.3.Crystallinity

    The raw material of LL exhibited a distinct diffraction peak in Fig.1B-1,wheras Sal B did not.SD LL was crystalline,with peaks around 6°,12°,20°,25°,31°,34° and 38° In contrast,powders containing both Sal B showed only a broad peak at about 20°,suggesting that the powder was mostly amorphous(Fig.1B-2).During long-term storage,amorphous powders may gradually change to crystalline form .As a result,storage conditions for Sal B-DPI should be monitored(temperature,humidity,etc.).

    Fig.1-Quality evaluation of spray-dried powders.(A)SEM images of SD Sal B-LL particles.(B)X-ray powder diffraction patterns of(B-1)raw powders and(B-2)spray-dried formulations containing Sal B and LL.(C)Hygroscopicity diagram of spray-dried formulations containing Sal B and LL.(D)The effect of amounts of LL on in vitro aerosolization performances of Sal B-LL composite powders(n=3).

    3.1.4.DVS

    At 0-90%RH,SD Sal B alone absorbed a significant amount of water.According to hygroscopicity classification systems,SD Sal B is very hygroscopic which is partly due to the amorphous form of Sal B,as confirmed by XRD.In contrast,SD LL alone was crystalline and non-hygroscopic,with less than 1%water absorbed.Thus,the contribution of water uptake of LL in the composite powders would be negligible.As LL was added to the formulations,the water uptake at elevated RH was decreased with increasing amounts of LL.It was discovered that LL could enhance the anti-hygroscopicity of DPI to a certain extent,thereby improving the stability of the drug.The hygroscopic behavior of 80SalB20LL and 75SalB25LL was similar when RH was 0-80%,and the hygroscopic behavior of 80SalB20LL was slightly stronger when RH was 80%-90%(Fig.1C).

    3.1.5.Content determination

    The content determination results of Sal B-DPI with various fomulations were shown in Table 2.The content of Sal B is 99.32%,94.24%,89.69%,84.76%,79.39% and 74.56%,respectively.

    Table 2-Content of Sal B in various formulations and results of in vitro deposition properties.

    3.1.6.In vitro aerosol performances

    The results ofin vitrodeposition properties of various formulations were shown in Table 2 and Fig.1D.The delivery efficiencies of SalB-DPI with various fomulations were all greater than 70%,and the FPF value ranged from 0.48% to 49.71%,with the FPF value reaching the maximum when the LL content was 20%.With the addition of LL increasing,thein vitrodeposition properties of DPI were improved [42].Among the six formulations mentioned above,80SalB20LL had the most obvious surface roughness,particle size less than 3 μm,relatively strong hygroscopic resistance,the highest delivery efficiency and the largest FPF value,and the best deposition effectin vitro.Therefore,80SalB20LL was selected for the preliminary stability test and subsequent test.

    On the basis of previous studies [41],the diameter of DPI was reduced from 8 μm to 3 μm by replacing the nozzle of B-90 spray dryer,and the particle size distribution was more uniform,so as to further improve the deposition efficiency of drugsin vitro.Similarly,the FPF value of 80SalB20LL increased from 37.55% to 49.71%.Other studies on Salvia miltiorrhiza related dry powder inhalers mainly adopted grinding method,with drug delivery efficiency less than 86.00% and FPF value less than 46.00% [43].In this experiment,the drug delivery efficiency of 80SalB20LL dry powder inhaler was 90.43%,and the FPF value was 49.71%.Moreover,the preparation method and prescription composition were convenient and simple,which provided reference for the preparation of DPI of traditional Chinese medicine.

    3.2.Preliminary stability investigation

    The particle size of Sal B-DPI(80SalB20LL)stored in the brown dryer at room temperature on Day 1 and 35 was less than 3 μm(Fig.S1A).With the increase of time,the surface roughness of particles did not increase significantly (Fig.S1B).The 80SalB20LL stored in the brown dryer at room temperature on Day 1 and 35 was mainly amorphous and did not crystallized(Fig.S1C).The content of Sal B in Sal B-DPI stored at room temperature on Day 1 and 35 was 79.39%±0.70%and 79.08%±0.39%,respectively,without significant change (Fig.S1D).Thein vitrodeposition properties of Sal B-DPI stored at room temperature on Day 1 and 35 were shown in Fig.S1E.The deposition of Sal B at various levels of NGI was shown in Fig.S1F.The delivery efficiency of Sal B-DPI on Day 1 and 35 were 90.13%±0.46%and 89.85%±0.82)%,respectively,and the FPF values were 49.05%±0.30%and 47.96%±0.95%,respectively.After 35 d,the deposition ratein vitrodecreased slightly,but there was no significant difference.The above indicated that Sal B-DPI stored in the brown dryer at room temperature was stable for 35 d

    3.3.Anti-pulmonary fibrosis in vitro

    3.3.1.MTT assay

    NIH-3T3 cells are commonly used cell lines forin vitroexperiments in cell-level experiments to investigate fibrosis.During the development of fibrosis,some fibroblasts differentiate into myofibroblasts and produce and accumulate excessive collagen,and the deposition of collagen will destroy the structure and function of lung tissue,thus triggering fibrosis [44].In this study,SalB-DPI(80SalB20LL) containing 79.39% Sal B was selected as therapeutic agent.MTT assay was used to detect the effects of nidanib and Sal B-DPI on NIH3T3 cell proliferation.The results showed that nidanib(<1 μM) and Sal B-DPI (<125 μg/ml d no significant effect on cell viability.Therefore,nidanib (0.25 and 0.5 μM) and Sal B-DPI (25,50 and 125 μg/ml were selected for subsequent experiments(Fig.S2).

    Fig.2-Effect of TGF-β1 on collagen content of NIH-3T3 cells(A)and effects of TGF-β1 on NIH-3T3 cell migration(B).###P <0.001 vs control group,???P <0.001 vs TGF-β1 group.

    3.3.2.Effect of TGF-β1-induced collagen content in NIH-3T3 cells

    In the injury and repair stage of pulmonary inflammation,macrophages would participate in the regulation of inflammatory response,secreting transforming growth factor TGF-β1,etc.,which further stimulated the proliferation of fibroblasts,resulting in remodeling of extracellular matrix(ECM),excessive deposition of collagen,and the development of irreversible pulmonary fibrosis.In this study,Sirius red staining was used to detect intracellular collagen deposition after TGF-β1-induced activation to confirm the successful induction of the cell model and to use it as a cell model for collagen deposition in IPF [44].As can be seen from Fig.2A,collagen content in cells significantly increased (red area significantly increased) 48 h after TGF-β1 induction,and the occurrence of this process was inhibited by the positive drug group (SB431542).Compared with TGF-β1 group,nidanib significantly inhibited the increase of collagen content in a concentration-dependent manner.When the concentration of Sal B-DPI was greater than 50 μg/ml it had a significant inhibitory effect on collagen content in NIH-3T3 cells in a concentration-dependent manner.The inhibitory effect of 125 μg/ml Sal B-DPI on collagen content was similar to that of 0.5 μM nidanib.

    3.3.3.Effect of TGF-β1-induced NIH-3T3 cell migration

    As shown in Fig.2B,NIH-3T3 cells migrated significantly after 48 h induced by TGF-β1 compared with the blank group,and SB431542 significantly inhibited the process.Compared with the TGF-β1 group,both nidanib and Sal B-DPI significantly inhibited the migration of NIH-3T3 cells in a concentration dependent manner.

    3.3.4.Effect of TGF-β1-induced NIH-3T3 cell activation

    Current studies have shown that IPF was associated with stages of oxidative stress,inflammatory response and fibrosis.In recent years,studies on the pathogenesis of IPF have not been limited to the pathogenesis of alveolar inflammation,which may be caused by endogenous and exogenous cellular stress,resulting in fibroblast activation,alveolar structural damage and ECM protein deposition,etc.,and further causing gas exchange disorders and respiratory failure.In the course of its pathogenesis,macrophages run through the whole pathogenesis,and different phenotypes have different effects on it [45].The polarization of macrophages is regulated by microbial signals,tissue-specific induction and cytokines,and the polarization of different subtypes is regulated by different signaling pathways.In different stages of pulmonary fibrosis,the two types of cells are polarized to different subtypes,and their polarized phenotypes include M1-type macrophages and M2-type macrophages[45].M1-type macrophages secrete iNOS as one of the nitric oxide synthase (NOS),and under normal circumstances,iNOS is almost not expressed in all tissues and cells,but under pathological conditions,it can be induced and activated by some cytokines,leading to the increase of NO content,and then affecting some pathological processes[46].i NOS and its product NO promote the occurrence of pulmonary fibrosis by participating in the oxidative stress injury of lung tissue.As an important enzyme in arginine metabolism,Arg-1 catalyzes the denitrification of arginine into ornithine,which is an essential precursor molecule for collagen synthesis.Under the regulation of type 2 helper T cell(Th2)cytokines,Arg-1 can be highly expressed in type M2 macrophages.It was found that overexpression of Arg-1 significantly increased bleomycin-induced pulmonary fibrosis in mice [47].The TGF-β1/Smad signaling pathway plays an important role in the process of cell polarization,and the accumulation of TGF-β1-mediated ECM precipitation is a Smad3-dependent process.FN,COL1A1 and COL3A1 are the main components of ECM.Excessive deposition of FN,COL1A1 and COL3A1 can lead to scarring,pulmonary dysfunction,respiratory failure,and ultimately death[48-49].ICAM-1 belongs to a family of immunoglobulins that regulate cell-to-cell and cell-ECM interactions.The main receptor of ICAM-1 was Lymphocyte function-associated antigen-1 (LFA-1),which was mainly expressed in the inflammatory cells such as eosinophils,neutrophils and T lymphocytes.The interaction between ICAM-1 and LFA-1 caused the above inflammatory cells to adhere to vascular endothelial cells and infiltrate the diseased bronchus.In addition,during lung injury,ICAM-1 was highly expressed in pulmonary vascular endothelial cells,leading to adhesion of inflammatory cells to the vascular endothelium and mediating its transendothelial metastasis.At the same time,the high expression of ICAM-1 in the alveolar epithelium also led to the connection of ICAM-1 with activated T cells and multinucleated leukocytes in the alveolar lumen,which led to the development of alveolar inflammation,and further increased the expression of ICAM-1 in inflammatory cells and alveolar epithelium[35].

    In order to detect the mRNA expression of TGF-β1 on fibrosis markers of NIH-3T3 cells,fibrosis related factors such as COL1A1,COL3A1,ICAM-1,FN,iNOS and Arg-1 were detected respectively.As shown in Fig.3,TGF-β1 could significantly increase the expression of fibrosis related factors in NIH-3T3 cells,but the intervention of Sal B-DPI (50 μg/ml)significantly reduced the expression of fibrosis related factors.

    Fig.3-Effects of TGF-β1-induced fibrosis related factors in NIH-3T3 cells.(A-F)mRNA expression of COL1A1,COL3A1,ICAM-1,FN,iNOS,and Arg-1.#P <0.05,## P <0.01,### P <0.001,#### P <0.0001 vs control group,?P <0.05,??P <0.01 vs TGF-β1 group.

    3.4.Pulmonary irritation

    The rats in normal control group (NC),sham operation group(sham) and Sal B-DPI group all responded sensitively with good mental status,stable breathing,good diet,fat body,shiny fur and significant weight gain.

    The lung tissues of rats in the NC group,sham group and Sal B-DPI group were intact and clear,with thin alveolar walls,no bleeding points,and no inflammation in the alveolar cavity cells infiltrated and exudated,with less collagen deposition(Fig.4A).

    When lung injury and inflammation occurred,inflammatory cells,cytokines and related enzymes in the BALF were abnormally increased.IL-8 is a chemotactic cytokine and an important mediator of inflammatory diseases.IL-8 can chemotactic and activate neutrophils,and then participate in the pulmonary inflammatory response and local lung injury [50].LDH is widely present in the cytoplasm of all tissues and cells of the body.In normal lung tissues,LDH is mainly derived from the active secretion of lung epithelial cells and the decomposition and release of exfoliated cells.However,in damaged lung tissue,the mucosal integrity of the airway epithelium is affected,resulting in elevated LDH enzyme expression [51].The total protein concentrations in BALF in NC group and Sal B-DPI group were 8.08±0.76 and 8.99±0.33 mg/ml,respectively,and there was a significant difference between them (P <0.05).There was no significant difference in the concentration of IL-8 and LDH activity in BALF between the three groups,indicating that the noninvasive tracheal administration is safe,and the lung injury is relatively small.Sal B-DPI may cause slight irritation to the lungs(Fig.4B).

    Fig.4-Pulmonary irritation.(A)Hematoxylin&eosin(H&E)staining of rat lung tissue(×200).(B)The concentrations of total protein,IL-8 and LDH in BALF.?P <0.05 vs NC group.

    3.5.Preliminary pharmacodynamics study in vivo

    As a glycoside antibiotic,bleomycin is mainly used in the clinical treatment of tumors,but there are obvious adverse reactions of pulmonary interstitial fibrosis,so it is often used to induce pulmonary fibrosis animal models.The histological changes of bleomycin-induced rat model of IPF are very similar to those of human IPF and can well reflect the pathological characteristics of human pulmonary fibrosis[52].In the experiment of establishing animal model,lung injury after bleomycin administration mainly experienced three different stages:(1) The first 7 d were in the stage of acute injury and inflammation.Due to more epithelial cell damage and increased vascular permeability,a large number of inflammatory cytokines were released and the number of inflammatory cells increased;(2) From the 7th d to the 14th d,it was in the transitional stage from inflammation to active pulmonary fibrosis,during which the inflammatory response gradually subsided,accompanied by a gradual increase in myofibroblast generation and collagen deposition;(3) From Day 14 to 28,chronic pulmonary fibrosis occurs,with obvious forms of alveolar and septal fibrosis,characterized by an increase in the number of myofibroblasts and the deposition of extracellular matrix.Generally,about 28 d after bleomycin administration,the degree of collagen deposition in the lungs reaches its peak.Bleomycin was widely used to establish animal models of IPF because it can cause pulmonary inflammation and fibrosis in the short term[52].In this study,the rats model of IPF was induced by one-time infusion of bleomycin normal saline through trachea,and the whole course of administration was used for intervention.

    The rats in the NC group are sensitive,and have good mental state,stable breathing,good diet,plump body,shiny fur,and significant weight gain.In the early stage,the modeled rats in other groups were poor in spirit,short of breath,wheezing,slow movements,sparse fur,dull fur,and severe shedding.The rats in the bleo group are in the worst condition,with a significant decrease in food intake,a decrease in drinking water,and an insignificant weight gain.After 7 d in the other treatment groups,the state of rats gradually improved,the breathing became stable,and the weight gradually increased,which was significantly better than the bleo group,but slightly worse than NC group.

    Low TV value indicates insufficient lung ventilation,while high TV value indicates excessive ventilation.MV value is related to TV and respiratory rate.With a small TV value,a higher respiratory rate is needed to ensure adequate ventilation.Pulmonary diseases such as IPF are often characterized by shallow and fast breathing,decreased TV value and increased respiratory rate.The decrease of TV value and MV value can be used as early diagnostic indicators of IPF.As shown in Fig.5A,TV value of rats in bleo group decreased,which was significantly different from that in NC group (P <0.01).Compared with bleo group,bleo+nidanib group was significantly higher (P <0.01),with a very significant difference (P <0.01),and bleo+Sal B-DPI group also increased,with significant difference (P <0.05).While bleo+Sal B iv group also increased,but there was no significant difference.The MV of bleo group was significantly lower than that of NC group(P <0.01).Compared with bleo group,bleo+nidanib group and bleo+Sal B-DPI group were significantly higher (P <0.01),and bleo+Sal B iv group were significantly higher(P <0.05).

    As shown in Fig.5B,the lung tissue structure of NC group rats is clear and intact,with thin alveolar walls,no bleeding points,and normal morphological structure.There was no infiltration and exudation of inflammatory cells and less collagen deposition in the alveolar cavity.Bleo group showed alveolar wall congestion,inflammatory cell infiltration and emphysema.After nidanib,Sal B-DPI and intravenous injection of Sal B solution,the degree of alveolar wall congestion,inflammatory cell infiltration and emphysema in the lung tissue of rats were lower than that of bleo group,bleo+Sal B-DPI group and bleo+nidanib group had better therapeutic effects.

    IPF is a common alveolar disease in which healthy lung tissue is replaced by excessive interstitial cells and extracellular matrix,resulting in the destruction of alveolar structures and the gradual densification of alveolar areas,eventually leading to respiratory failure and even death[1-6].At present,studies have shown that IPF was associated with stages of oxidative stress,inflammatory response and fibrosis[2-5].Oxidative stress caused by oxidation/antioxidant imbalance is one of the important factors in the pathogenesis of IPF,and MDA content can reflect the degree of lipid peroxidation in the body.MPO is mainly secreted and catalyzed by neutrophils to produce excessive oxidizing substances,leading to oxidative stress and oxidative tissue damage.SOD can scavenge oxygen free radicals and prevent cells from being damaged by free radicals.TGF-β1 is an important pro-fibrosis factor,which can inhibit its expression,thereby inhibiting the proliferation and activation of fibroblasts,inhibiting collagen synthesis,and reducing extracellular matrix deposition.In the process of pulmonary fibrosis,the increased expression of a variety of inflammatory factors will also cause inflammation and fibrosis.The imbalance of Th1/Th2 in the BALF plays an important role in the pathogenesis of pulmonary fibrosis,and various cytokines secreted by Th1 cells and Th2 cells are involved in the formation of pulmonary fibrosis.Th2 cytokines,such as IL-1β,IL-4,IL-6 and IL-18,can all cause alveolar inflammation and damage in the early stage of pulmonary fibrosis,leading to the activation and proliferation of fibroblast cells and collagen synthesis,leading to matrix protein deposition and production of fibrous tissue [53].IFN-γis a Th1 cytokine that inhibits the proliferation of fibroblasts and the formation of fibrous tissue [54].Hyp is one of the twelve nonessential amino acids that is characteristic of collagen and is a major component of the extracellular matrix during the development of fibrosis in various organs[55].

    The MDA content in the lung tissue homogenate of rats in the bleo group was significantly higher than that in NC group(P<0.01),and the serum MPO content was also significantly higher than that in NC group(P <0.01).The content of MDA and MPO in bleo+nidanib group,bleo+Sal B-DPI group,and bleo+Sal B iv group can be significantly reduced.Bleo+Sal B-DPI group and bleo+nidanib group had a more significant reduction in MDA and MPO content.Compared with NC group,the SOD activity in the lung tissue homogenate of the bleo group was significantly reduced(P <0.01).The SOD activity of bleo+nidanib group,bleo+Sal B-DPI group,and bleo+Sal B iv group increased significantly,but the improvement effect was more significant in bleo+Sal B-DPI group and bleo+nidanib group,as shown in Fig.5C.

    Compared with NC group,the contents of TGF-β1 protein,IL-1β,IL-4,IL-6 and IL-18 in bleo group were significantly increased (P <0.01),and the contents of IFN-γwere significantly decreased (P <0.01).Compared with bleo group,the contents of TGF-β1 protein,IL-1β,IL-4,IL-6 and IL-18 were significantly decreased in bleo+nidanib group and bleo+Sal B-DPI group,and the contents of IFN-γwere significantly increased(P <0.01).TGF-β1 protein,IL-6 (P <0.01) and IL-1β,and IL-18 (P <0.05) were significantly decreased in bleo+Sal B iv group,but IL-4 and IFN-γcontents were not significantly different from those in the bleo group,as shown in Fig.5D.

    Compared with NC group,the content of Hyp in lung tissue homogenate of bleo group was significantly increased(P <0.01).Compared with bleo group,the levels of Hyp were significantly lower in bleo+nidanib group,bleo+Sal B-DPI group and bleo+Sal B iv group.However,bleo+Sal B-DPI group and bleo+nidanib group had a more significant reduction effect(P <0.01),as shown in Fig.5E.

    Fig.5-Preliminary pharmacodynamics study in vivo.(A)Pulmonary function indexes of rats(Mean±SD,n=6).(B)Hematoxylin&eosin(H&E)staining of rat lung tissue(×200).(C)Comparison of MPO contents in serum,SOD and MDA contents in lung homogenate(Mean±SD,n=6).(D)The concentrations of inflammatory factors and IFN-γ in BALF(Mean±SD,n=6).(E)Comparison of Hyp content in lung tissue(Mean±SD,n=6).## P <0.01 vs NC group,?P <0.05,??P <0.01 vs bleo group.

    3.6.Metabonomics study of lung tissue in rats

    Multivariate statistical analysis is mostly used to deal with the characteristics of complex sample classification by dimensionality reduction.Its core idea is to condense multiple variables into one or two with the strongest characteristics through variance changes.As shown in Fig.S3,A1 showed the scatter diagram of PCA analysis of multiple groups.Each sample in each group presents a certain natural aggregation state under the projection of principal component analysis,and this characteristic is further strengthened in OPLS-DA analysis (Fig.S3A2).This indicates that there are certain similarities among each sample group,that is,common characteristics,and there are significant differences among each group.

    To further search for and characterize the difference components,OPLS-DA between the two groups was used to screen potential metabolic markers,as shown in Fig.12B-12E.The R2Y and Q2parameters of OPLS-DA were all close to 1.The replacement test analysis showed that the model established by OPLS-DA in each group was good without overfitting (Table S2),which could be used for further screening of candidate metabolic markers based on VIP value.In the multivariate statistical analysis of the negative ion mode,compared with the model group,the Q2parameter of OPLSDA in the pulmonary administration of Sal B-DPI group(F)and intravenous injection of Sal B solution group(W)was less than 0.9 but greater than 0.5,which was still an effective model and could be further analyzed.The specific results were shown in Fig.S4 and Table S3.

    In the OPLS-DA discriminant model based on multivariate statistical analysis,VIP parameters were used as the Cut off for screening differential endogenous small molecules through supervision and identification.On the basis of VIP>1,Metlin and HMDB online databases were used to identify the sources (endogenous,exogenous,drug-derived,etc.) and the distribution of body fluids (serum,urine,bile,tissue,etc.) of each feature,screen the endogenous ion peak of lung tissue,and confirm the final difference components by MS/MS and PPM.The HMDB database was used to compare the primary excimer ion peaks of the differential components according to the precise molecular weight,and the final differential components were confirmed by MS/MS and ppm.T=5.88 min,m/z496.3406) was used as an example to illustrate the identification process of the compounds.Firstly,m/z496.3406 was used to conduct an accurate molecular weight search in the HMDB database,ppm was set to 15,and the ion summation forms wereM+H,M+Na,M+NH4,M+K,M+Li,etc.,for matching.Then,on the basis of matching results,exogenous components were eliminated and only endogenous components were retained.The molecular composition was preliminarily determined to be C24H50NO7P,ppm=1.61<<15.Then,based on MassLYNX and QI workstations,the precise molecular weight of the component was searched,and its MS/MS at RT=5.87 was obtained.As shown in Fig.S5,the mass spectrum fragments ofm/z104.1078,125.0009 and 184.0743 correspond to-C5H13NO,-C9H17,-C5H13NO4P,and basically determine thatm/z496.3406 and RT=5.87 were LYSOPC(16:0/0:0).

    As shown in Fig.S6,a total of 23 differential endogenous metabolites were identified based on UPLC-Q-TOF/MS in the positive and negative ion mode.After administering the positive drug nintedanib to rats,the 23 disease markers can be recalled to 7 [3-Methyl-2-butene-1-thiol,Indole-3-carboxaldehyde,LysoPC (P-18:0/0:0)),SM (d18:1/16:1),Indoleacrylic acid,L-Cystine and Hypoxanthine] (P<0.05);after pulmonary administration of Sal B-DPI to rats,23 disease markers can be recalled 7 [LysoPC (14:0/0:0),LysoPC(P-18:0/0:0),PE (16:0/15:0),PC (18:4/18:1),SM (d18:1/16:1),LCystine and Hypoxanthine] (P <0.05);after administration of Sal B solution to rats via the tail vein route,23 disease markers can be recalled 3[PE(16:0/15:0),PC(22:6/18:2)and L-Cystine](P<0.05).Bleo+Sal B-DPI group is equivalent to bleo+nintedanib group.In addition,through metabolic pathway analysis,it was found that the differential metabolites in bleo+Sal B-DPI group were mainly enriched in four pathways:glycerophospholipid metabolism,linolenic acid metabolism,α-linolenic acid metabolism and glycosylphosphatidylinositol anchor biosynthesis(Fig.S7).

    3.7.Preliminary pharmacokinetics and lung tissue distribution study

    Duration curves of Sal B at different times after pulmonary administration of Sal B-DPI(Sal B-DPI group) and intravenous injection of Sal B solution(Sal B iv group) were shown in Fig.S8A.According to the blood concentration of Sal B,WinNonlin 6.4 software was used to calculate the pharmacokinetic parameters,as shown in Table 3.Tmaxof Sal B-DPI group and Sal B iv were 0.083 and 0.033 h,respectively.TheCmaxvalues were 8012.26±751.84 and 16,404.17±2237.56 ng/ml,respectively,and the differences were extremely significant(P <0.01).AUC0-∞was 3841.29±515.73 and 6554.95±1485.24 h?ng/ml,respectively,with extremely significant difference(P<0.01).The absolute bioavailability of Sal B was 59.67%±6.30%for Sal B-DPI group.

    Table 3-Pharmacokinetic parameters of Sal B in rats plasma after intravenous injection of Sal B solution and pulmonary administration of Sal B-DPI(Mean±SD,n=6).

    Table 4-Pharmacokinetic parameters of Sal B in rat lung tissue after intravenous injection of Sal B solution and pulmonary administration of Sal B-DPI(Mean±SD,n=3).

    The distribution of Sal B in the lung tissue after pulmonary administration of Sal B-DPI(Sal B-DPI group)and intravenous injection of Sal B solution(Sal B iv group)were shown in Figure S8B.According to the concentration of Sal B in the lung tissue,WinNonlin 6.4 software was used to calculate the pharmacokinetic parameters in Table 4.Tmaxof Sal B-DPI group and Sal B iv group were 0.033 and 0.033 h,respectively.Cmaxvalues were 142,594.47±3654.22 and 13,150.38±619.02 ng/ml,respectively,with extremely significant differences(P <0.01).AUC0-∞was 80,583.35±4225.84 and 1898.5±134.03 h?ng/ml,respectively.AUC0-∞of Sal B-DPI group was 42.45 times that of Sal B iv group,with significant difference(P <0.01).

    Some studies reported that the oral bioavailability of Sal B was no more than 5.56% [18-20].In previous studies,the absolute bioavailability of Sal B was 19.15%±7.44% [41],while the absolute bioavailability of Sal B was 59.67%±6.30%after pulmonary administration of Sal B-DPI.Compared with the oral bioavailability reported in the study,the oral bioavailability increased by nearly 11 times,and the oral bioavailability increased by more than 3 times compared with the previous experimental results,which should be attributed to the reduction of Sal B-DPI particle size and the improvement ofin vitrodeposition efficiency.From the perspective of systemic absorption,the development of Sal B pulmonary drug delivery system is very promising.The therapeutic effect of drugs is usually directly related to the drug concentration of the target.The concentration of Sal B in lung tissue was still 2940.62±117.04 ng/g after 1 h of pulmonary administration,indicating that dry powder inhalers have certain advantages in the treatment of pulmonary diseases(such as IPF).

    4.Conclusions

    In this study,L-leucine as excipient was used to prepare Sal BDPI successfully by spray drying method,which improved the powder properties,drug delivery efficiency,in vitrodeposition effect of Sal B-DPI,and enhanced its moisture resistance.In the form of DPI,it was targeted to the lungs through tracheal non-invasive lung drug delivery technology,which had no irritation to the lungs,overcomed the shortcomings of poor oral absorption and low bioavailability of Sal B,increased the drug concentration in the lesion site,and improved the therapeutic effect of IPF.The above work provided new ideas and methods for the clinical treatment of IPF and the clinical application of Sal B.

    Conflicts of interest

    The authors have no conflicts of interest to declare.

    Acknowledgements

    This work was supported by Innovation Team and Talents Cultivation Program of National Administration of Traditional Chinese Medicine (No.ZYYCXTD-D-202002) and Scientific Research Project of Tianjin Municipal Education Commission(No.2019KJ083).

    Supplementary materials

    Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.ajps.2022.04.004.

    99riav亚洲国产免费| 国产单亲对白刺激| 国产黄色视频一区二区在线观看 | 性插视频无遮挡在线免费观看| eeuss影院久久| 久久人人精品亚洲av| 免费高清视频大片| 一本久久中文字幕| 一本精品99久久精品77| av天堂中文字幕网| 少妇丰满av| 观看美女的网站| 国产成人aa在线观看| 国产高清三级在线| 成人亚洲欧美一区二区av| 深夜a级毛片| 日韩强制内射视频| 麻豆av噜噜一区二区三区| 人妻丰满熟妇av一区二区三区| 插逼视频在线观看| 日韩中字成人| 免费不卡的大黄色大毛片视频在线观看 | 此物有八面人人有两片| 小说图片视频综合网站| .国产精品久久| 久久久久国产网址| 少妇的逼水好多| 日韩大尺度精品在线看网址| 一级毛片aaaaaa免费看小| 老女人水多毛片| 成人一区二区视频在线观看| 真人做人爱边吃奶动态| 成人高潮视频无遮挡免费网站| 精品一区二区三区人妻视频| 亚洲精品日韩在线中文字幕 | 日韩欧美精品v在线| 最新中文字幕久久久久| АⅤ资源中文在线天堂| 男女视频在线观看网站免费| 麻豆精品久久久久久蜜桃| 欧美成人免费av一区二区三区| 亚洲四区av| 亚洲精品久久国产高清桃花| 哪里可以看免费的av片| 亚洲av中文av极速乱| 亚洲最大成人中文| 免费av观看视频| 又爽又黄a免费视频| 热99re8久久精品国产| 搡老妇女老女人老熟妇| 亚州av有码| 高清日韩中文字幕在线| 国产高清激情床上av| 亚洲自拍偷在线| 日本色播在线视频| av免费在线看不卡| 黄片wwwwww| 99国产极品粉嫩在线观看| 天天躁夜夜躁狠狠久久av| av福利片在线观看| 欧美日本视频| 又黄又爽又刺激的免费视频.| 国产单亲对白刺激| 非洲黑人性xxxx精品又粗又长| 午夜福利在线观看吧| 国产精品女同一区二区软件| 久久精品国产自在天天线| 国产午夜精品久久久久久一区二区三区 | 赤兔流量卡办理| 免费观看精品视频网站| 成熟少妇高潮喷水视频| 一进一出抽搐动态| 欧美成人一区二区免费高清观看| 日本a在线网址| 国产精品日韩av在线免费观看| 国产精品av视频在线免费观看| 无遮挡黄片免费观看| 日韩欧美 国产精品| 九九在线视频观看精品| 在线a可以看的网站| 人妻制服诱惑在线中文字幕| 国产精品野战在线观看| 国产成人a区在线观看| 天堂影院成人在线观看| 亚洲av电影不卡..在线观看| 夜夜夜夜夜久久久久| 亚洲国产精品成人久久小说 | 亚洲高清免费不卡视频| 少妇猛男粗大的猛烈进出视频 | 国产三级在线视频| av.在线天堂| 免费av不卡在线播放| 大香蕉久久网| 熟妇人妻久久中文字幕3abv| 亚洲人成网站在线观看播放| 国产精品久久电影中文字幕| 国产成人a区在线观看| 欧美潮喷喷水| 日韩在线高清观看一区二区三区| av黄色大香蕉| 成年免费大片在线观看| 亚洲最大成人手机在线| 岛国在线免费视频观看| 日韩亚洲欧美综合| 日本欧美国产在线视频| 欧美中文日本在线观看视频| 久久精品国产亚洲av天美| 男人狂女人下面高潮的视频| 女的被弄到高潮叫床怎么办| 欧洲精品卡2卡3卡4卡5卡区| 大又大粗又爽又黄少妇毛片口| 欧美+日韩+精品| 免费无遮挡裸体视频| 一进一出抽搐gif免费好疼| 午夜a级毛片| 日本精品一区二区三区蜜桃| 国产日本99.免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品国产三级国产av玫瑰| 久久午夜亚洲精品久久| 午夜日韩欧美国产| 草草在线视频免费看| 俺也久久电影网| av黄色大香蕉| 亚洲,欧美,日韩| 欧美绝顶高潮抽搐喷水| 高清午夜精品一区二区三区 | 在线观看av片永久免费下载| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久久亚洲| 国产亚洲91精品色在线| 久久综合国产亚洲精品| 成人无遮挡网站| 少妇人妻精品综合一区二区 | 丰满的人妻完整版| a级毛片免费高清观看在线播放| 中国美女看黄片| 又爽又黄a免费视频| 国产一级毛片七仙女欲春2| 久久精品影院6| 秋霞在线观看毛片| 中文亚洲av片在线观看爽| 国产黄a三级三级三级人| 少妇熟女欧美另类| 少妇高潮的动态图| 99热只有精品国产| 亚洲精品456在线播放app| av福利片在线观看| 久久久久久大精品| 午夜免费激情av| 国产久久久一区二区三区| 中文字幕免费在线视频6| 国产又黄又爽又无遮挡在线| 日韩一区二区视频免费看| 国产精品一区二区三区四区免费观看 | 黑人高潮一二区| 99在线视频只有这里精品首页| 久久久欧美国产精品| 欧美bdsm另类| 乱系列少妇在线播放| 国产精品不卡视频一区二区| 色综合亚洲欧美另类图片| 1024手机看黄色片| 国产精品一区二区三区四区免费观看 | 亚洲国产色片| 97超碰精品成人国产| 精品国产三级普通话版| 亚洲最大成人手机在线| 欧美日韩国产亚洲二区| 国产成年人精品一区二区| 国产成人aa在线观看| 国产免费男女视频| 欧美精品国产亚洲| 亚洲国产色片| 国产精品野战在线观看| 99久久九九国产精品国产免费| 久久久欧美国产精品| 欧美激情久久久久久爽电影| 伦精品一区二区三区| 九九在线视频观看精品| 女同久久另类99精品国产91| 国语自产精品视频在线第100页| 成人永久免费在线观看视频| 赤兔流量卡办理| 日本熟妇午夜| 国产精品不卡视频一区二区| 看片在线看免费视频| 久久久久国产网址| 亚洲一区二区三区色噜噜| 神马国产精品三级电影在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产av麻豆久久久久久久| 久久久a久久爽久久v久久| 国产色爽女视频免费观看| 亚洲av成人精品一区久久| 午夜视频国产福利| 内射极品少妇av片p| 日本成人三级电影网站| 国产国拍精品亚洲av在线观看| 欧美一区二区国产精品久久精品| 麻豆国产av国片精品| 免费看日本二区| 精品久久久噜噜| 大香蕉久久网| 色在线成人网| 嫩草影院精品99| 国产人妻一区二区三区在| 欧美一级a爱片免费观看看| 国产精品久久久久久久电影| 日日撸夜夜添| 99久国产av精品国产电影| 欧美日韩国产亚洲二区| 女生性感内裤真人,穿戴方法视频| 国产视频内射| 伊人久久精品亚洲午夜| 日韩精品有码人妻一区| 中文字幕熟女人妻在线| 亚洲va在线va天堂va国产| 可以在线观看的亚洲视频| 国产 一区 欧美 日韩| 欧美人与善性xxx| 国产成人精品久久久久久| 全区人妻精品视频| 18禁黄网站禁片免费观看直播| 欧美性感艳星| 91精品国产九色| 看十八女毛片水多多多| 久久精品国产自在天天线| 日日干狠狠操夜夜爽| 级片在线观看| 少妇丰满av| 亚洲欧美中文字幕日韩二区| 亚洲国产高清在线一区二区三| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 久久国产乱子免费精品| 日本撒尿小便嘘嘘汇集6| 国产老妇女一区| 99九九线精品视频在线观看视频| 日日摸夜夜添夜夜添小说| 精品国产三级普通话版| 老司机午夜福利在线观看视频| 亚洲中文字幕一区二区三区有码在线看| 国产精品精品国产色婷婷| 国产精品一区二区免费欧美| 久久99热这里只有精品18| 成年免费大片在线观看| 亚洲精品日韩在线中文字幕 | 欧美中文日本在线观看视频| 两个人的视频大全免费| 99久久无色码亚洲精品果冻| 欧美绝顶高潮抽搐喷水| av在线天堂中文字幕| 久久人妻av系列| 亚洲自拍偷在线| 亚洲欧美日韩高清专用| 日本黄色视频三级网站网址| 一边摸一边抽搐一进一小说| 国产毛片a区久久久久| 亚洲最大成人av| 观看免费一级毛片| 精品一区二区免费观看| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 狂野欧美白嫩少妇大欣赏| av.在线天堂| 国产高清有码在线观看视频| 搡老熟女国产l中国老女人| 国产高清激情床上av| 一进一出抽搐动态| 国产精品综合久久久久久久免费| 99久久精品国产国产毛片| 搞女人的毛片| 日日撸夜夜添| 一级黄片播放器| 午夜日韩欧美国产| 亚洲成人av在线免费| eeuss影院久久| 九九热线精品视视频播放| videossex国产| 日韩大尺度精品在线看网址| 91久久精品电影网| 精品久久久久久久久av| 精品人妻熟女av久视频| 在线播放国产精品三级| 大香蕉久久网| 人人妻,人人澡人人爽秒播| 国产一区二区激情短视频| 欧美成人一区二区免费高清观看| 亚洲图色成人| 搞女人的毛片| 神马国产精品三级电影在线观看| 精品乱码久久久久久99久播| 亚洲七黄色美女视频| 午夜精品在线福利| 国产精品永久免费网站| 天天一区二区日本电影三级| 三级国产精品欧美在线观看| 国产蜜桃级精品一区二区三区| 一级毛片久久久久久久久女| 最近手机中文字幕大全| 直男gayav资源| 天天一区二区日本电影三级| 2021天堂中文幕一二区在线观| 亚洲18禁久久av| 亚洲婷婷狠狠爱综合网| 精品无人区乱码1区二区| 日本一本二区三区精品| 免费无遮挡裸体视频| 成人三级黄色视频| 中国国产av一级| 日韩 亚洲 欧美在线| 国产大屁股一区二区在线视频| 欧美激情久久久久久爽电影| 国产精品福利在线免费观看| 成人特级av手机在线观看| 男人和女人高潮做爰伦理| 中文字幕人妻熟人妻熟丝袜美| 欧美一区二区国产精品久久精品| 搡老妇女老女人老熟妇| 啦啦啦啦在线视频资源| 午夜免费男女啪啪视频观看 | 久久午夜福利片| 欧洲精品卡2卡3卡4卡5卡区| 国产精品国产高清国产av| 国产 一区精品| av在线播放精品| 成人二区视频| 内射极品少妇av片p| 精品人妻视频免费看| 18禁在线播放成人免费| a级毛色黄片| 亚洲成av人片在线播放无| 中文字幕av成人在线电影| 长腿黑丝高跟| 麻豆国产av国片精品| 国内精品久久久久精免费| 欧美日韩综合久久久久久| 亚洲精品一区av在线观看| 国产成人福利小说| 亚洲精品久久国产高清桃花| 欧美国产日韩亚洲一区| 亚洲经典国产精华液单| 久久亚洲精品不卡| 毛片一级片免费看久久久久| 内地一区二区视频在线| .国产精品久久| 色综合亚洲欧美另类图片| 午夜福利在线观看吧| 国产乱人偷精品视频| 国产精品亚洲一级av第二区| 我的女老师完整版在线观看| 最后的刺客免费高清国语| 久久精品国产鲁丝片午夜精品| 在线观看66精品国产| 久久这里只有精品中国| 美女内射精品一级片tv| 校园春色视频在线观看| 国产毛片a区久久久久| 日韩精品有码人妻一区| 亚洲五月天丁香| 国产大屁股一区二区在线视频| 国产老妇女一区| 黑人高潮一二区| 亚洲欧美中文字幕日韩二区| 日韩一区二区视频免费看| 神马国产精品三级电影在线观看| 激情 狠狠 欧美| 99热这里只有是精品在线观看| 伦精品一区二区三区| 深夜a级毛片| 一夜夜www| 亚洲国产精品久久男人天堂| 日韩在线高清观看一区二区三区| 国产高清视频在线播放一区| 成人特级av手机在线观看| 成年免费大片在线观看| 欧美激情国产日韩精品一区| 精华霜和精华液先用哪个| 六月丁香七月| 国产在线精品亚洲第一网站| 国产高清视频在线播放一区| 日韩强制内射视频| 中出人妻视频一区二区| 亚洲av五月六月丁香网| 丰满的人妻完整版| 亚洲av免费高清在线观看| 午夜精品一区二区三区免费看| 在线观看av片永久免费下载| 国产精品精品国产色婷婷| 成年女人毛片免费观看观看9| 黄色一级大片看看| 国产精品久久久久久精品电影| 日本欧美国产在线视频| 天堂av国产一区二区熟女人妻| 自拍偷自拍亚洲精品老妇| 最近中文字幕高清免费大全6| 亚洲美女搞黄在线观看 | 听说在线观看完整版免费高清| 91精品国产九色| av福利片在线观看| 国产欧美日韩精品亚洲av| 国产中年淑女户外野战色| 国产成人freesex在线 | 久久精品国产亚洲av天美| 永久网站在线| 狂野欧美激情性xxxx在线观看| 欧美3d第一页| av福利片在线观看| 国产黄a三级三级三级人| 日本三级黄在线观看| 寂寞人妻少妇视频99o| 成人午夜高清在线视频| 日韩在线高清观看一区二区三区| 美女免费视频网站| 国产亚洲精品久久久com| 日日摸夜夜添夜夜添小说| 成人三级黄色视频| 午夜福利视频1000在线观看| 亚洲人成网站在线观看播放| 欧美bdsm另类| 亚洲欧美日韩高清专用| 联通29元200g的流量卡| 久久人人精品亚洲av| 亚洲国产精品成人久久小说 | 天天一区二区日本电影三级| 亚洲无线在线观看| 国产又黄又爽又无遮挡在线| 午夜福利18| 51国产日韩欧美| 我的老师免费观看完整版| 夜夜爽天天搞| 成人精品一区二区免费| 国产av不卡久久| 中文字幕人妻熟人妻熟丝袜美| 亚洲成人中文字幕在线播放| 一级黄片播放器| 全区人妻精品视频| 亚洲婷婷狠狠爱综合网| 又黄又爽又免费观看的视频| 免费观看人在逋| 青春草视频在线免费观看| 丰满人妻一区二区三区视频av| 亚洲国产精品国产精品| 国产蜜桃级精品一区二区三区| 免费黄网站久久成人精品| 国产午夜精品久久久久久一区二区三区 | 日韩av在线大香蕉| 国产成人一区二区在线| 日本精品一区二区三区蜜桃| 欧美另类亚洲清纯唯美| 欧美又色又爽又黄视频| 99久国产av精品国产电影| 国产毛片a区久久久久| 国产精品女同一区二区软件| 波多野结衣巨乳人妻| 长腿黑丝高跟| 国产高清三级在线| 久久久久久久久久黄片| 精品99又大又爽又粗少妇毛片| 我的女老师完整版在线观看| 99久久精品一区二区三区| 欧美日韩乱码在线| 女同久久另类99精品国产91| 可以在线观看的亚洲视频| 亚洲精品一区av在线观看| 国产91av在线免费观看| 波多野结衣高清作品| 亚洲成av人片在线播放无| 欧美zozozo另类| 亚洲电影在线观看av| 最近的中文字幕免费完整| 自拍偷自拍亚洲精品老妇| 高清毛片免费看| 免费高清视频大片| 久久久久久久久久成人| 俺也久久电影网| 又粗又爽又猛毛片免费看| 99久国产av精品| 国产精品免费一区二区三区在线| 亚洲人成网站在线观看播放| 色播亚洲综合网| 寂寞人妻少妇视频99o| a级毛色黄片| 男女做爰动态图高潮gif福利片| 人妻丰满熟妇av一区二区三区| av在线天堂中文字幕| 婷婷精品国产亚洲av| 久久午夜亚洲精品久久| 午夜福利18| 国产在视频线在精品| 亚洲精品粉嫩美女一区| 亚洲国产精品合色在线| 精品欧美国产一区二区三| 蜜臀久久99精品久久宅男| 亚洲精品乱码久久久v下载方式| 久久精品影院6| 熟妇人妻久久中文字幕3abv| 18禁裸乳无遮挡免费网站照片| 91午夜精品亚洲一区二区三区| 色av中文字幕| 99在线人妻在线中文字幕| 欧美绝顶高潮抽搐喷水| 国产精华一区二区三区| 午夜影院日韩av| 亚洲不卡免费看| 九九久久精品国产亚洲av麻豆| 国产aⅴ精品一区二区三区波| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| avwww免费| 中文字幕免费在线视频6| 欧美bdsm另类| 精品午夜福利视频在线观看一区| 国产高清激情床上av| 内地一区二区视频在线| 免费电影在线观看免费观看| 如何舔出高潮| 精品久久久久久久久亚洲| 麻豆av噜噜一区二区三区| 亚洲美女黄片视频| 亚洲精品456在线播放app| 五月伊人婷婷丁香| 一卡2卡三卡四卡精品乱码亚洲| 午夜a级毛片| 久久久久久久久久久丰满| 久久久成人免费电影| 一本精品99久久精品77| 国产精品国产三级国产av玫瑰| 亚洲欧美成人精品一区二区| 最后的刺客免费高清国语| 嫩草影视91久久| 亚洲第一电影网av| 成人特级av手机在线观看| 黄色日韩在线| 热99re8久久精品国产| 男女做爰动态图高潮gif福利片| 天堂av国产一区二区熟女人妻| 夜夜爽天天搞| 亚洲av中文av极速乱| 久久久久精品国产欧美久久久| 精品久久国产蜜桃| 亚洲av二区三区四区| 国产成人aa在线观看| 婷婷六月久久综合丁香| 亚洲精华国产精华液的使用体验 | 麻豆乱淫一区二区| 久久精品国产清高在天天线| 国产午夜精品论理片| 成人无遮挡网站| 亚洲一区二区三区色噜噜| 国产精品美女特级片免费视频播放器| 日韩中字成人| 天堂av国产一区二区熟女人妻| 成人特级av手机在线观看| 97超碰精品成人国产| 一区二区三区四区激情视频 | 亚洲熟妇熟女久久| 国产精品国产高清国产av| 国产美女午夜福利| 免费av毛片视频| 91狼人影院| 男人和女人高潮做爰伦理| 黄片wwwwww| 国产在线男女| 久久久久九九精品影院| 久久人人精品亚洲av| 国内精品久久久久精免费| 狠狠狠狠99中文字幕| 午夜福利高清视频| 99久久成人亚洲精品观看| 国产精品女同一区二区软件| av在线蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 我要看日韩黄色一级片| 国产极品精品免费视频能看的| 欧美成人免费av一区二区三区| 国产一区二区亚洲精品在线观看| 成人亚洲欧美一区二区av| 亚洲一区高清亚洲精品| 少妇的逼好多水| 亚洲经典国产精华液单| 桃色一区二区三区在线观看| 国产男人的电影天堂91| 国语自产精品视频在线第100页| 亚洲av免费高清在线观看| 国产精品福利在线免费观看| 久久精品国产鲁丝片午夜精品| 波多野结衣高清作品| 国产乱人偷精品视频| 成人亚洲欧美一区二区av| 99久国产av精品国产电影| 日韩av不卡免费在线播放| 晚上一个人看的免费电影| 久久精品影院6| 欧美绝顶高潮抽搐喷水| 国产亚洲精品久久久com| 日本与韩国留学比较| 亚洲精品乱码久久久v下载方式| 天堂动漫精品| 欧美高清成人免费视频www| 一个人观看的视频www高清免费观看| 成人av一区二区三区在线看| 欧美成人免费av一区二区三区| 国产高清视频在线观看网站| 中文资源天堂在线| 中国美女看黄片| 日本免费一区二区三区高清不卡| 男人狂女人下面高潮的视频| 我的女老师完整版在线观看| 在线播放无遮挡| 中国美女看黄片| 又粗又爽又猛毛片免费看| 亚洲国产精品久久男人天堂| 精品久久久久久久久av|