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

    Curcumin encapsulated dual cross linked sodium alginate/montmorillonite polymeric composite beads for controlled drug delivery

    2021-05-12 06:30:16SreeknthReddySuhJithendrMdhviChowdojiRo
    Journal of Pharmaceutical Analysis 2021年2期

    O.Sreeknth Reddy,M.C.S.Suh,*,T.Jithendr,C.Mdhvi,K.Chowdoji Ro

    aDepartment of Chemistry,Sri Krishnadevaraya University,Ananthapuramu,515003,India

    bDepartment of Polymer Science and Technology,Sri Krishnadevaraya University,Ananthapuramu,515003,India

    Keywords:

    Curcumin

    Sodium alginate

    Montmorillonite

    Microbeads

    Drug delivery

    A B S T R A C T

    The aim of the present work is fabrication of dual cross linked sodium alginate(SA)/montmorillonite(MMT)microbeads as a potential drug vehicle for extended release of curcumin(CUR).The microbeads were prepared using in situ ion-exchange followed by simple ionotropic gelation technique.The developed beads were characterized by Fourier transform infrared spectroscopy(FTIR),differential scanning calorimetry(DSC),thermogravimetric analysis(TGA),X-ray diffraction(X-RD)and scanning electron microscopy(SEM).The effect of MMT on encapsulation efficiency of CUR and intercalation kinetics was investigated.Dynamic swelling study and in vitro release study were investigated in simulated intestinal fluid(pH 7.4)and simulated gastric fluid(pH 1.2)at 37°C.Results suggested that both the swelling and in vitro release studies were influenced by the pH of test media,which might be suitable for intestinal drug delivery.The release mechanism was analyzed by fitting the release data into Korsmeyer-Peppas equation.

    1.Introduction

    In today's pharmaceutical formulation,polymer based hydrogel plays a vital role in several biomedical applications.Hydrogels have many advantages over other drug carriers due to their water intake capacity,biocompatibility and biodegradability[1,2].However,hydrogels have a few drawbacks like unrestrained release rate and non-controlled swelling properties,which lead to several side effects.To control the release rate and swelling properties of hydrogels,certain substances have been incorporated in the hydrogels such as clay minerals,surface coating with other polymers such as chitosan[3]and poly-L-lysine[4].Controlled release drugs have many advantages over conventional forms such as reducing release rate,minimizing side effects and keeping drug concentration at effective levels in plasma.

    For the last few decades,clay minerals have been extensively used in pharmaceutical formulations as lubricants,disintegrants,diluents,pigments and binders.Since ancient times,clay materials have been used to treat food poisoning,aches and pains,infections,and mineral deficiencies[5].In recent years,researchers show much interest in montmorillonite(MMT)in pharmaceutical formulations,because of good adsorbance ability,large Specific area,high cation-exchange capacity,standout adhesive ability,intercalate large molecules into the space and drug carrying capability[6].MMT is a naturally occurring inorganic material,composed of silica tetrahedral sheets layered between alumina octahedral sheets at a ratio of 2:1,respectively[7].MMT intercalates with drug molecules by adsorption process;during intercalation drug molecules replace the hydroxyl groups present on the outer surface and inter layer space of clay material[8].Many researchers have reported that the mechanisms involved in the interaction between drug molecules and MMT are hydrogen bonding,ion exchange,van der Waals interaction,hydrophilic/hydrophobic interaction and so on[9-11].Among all the mechanisms,ion exchange interactions have received considerable attention between the two components of the hybrid system[12].Ion-exchange interaction might take place by mixing solid substrates(namely,ion exchangers)with ionic drugs in solution.In biological fluids,“counterions”can displace the drug from the substrate and deliver it into the body.The exchanger may be then eliminated or biodegraded[13].

    Sodium alginate(SA)is an anionic polysaccharide comprising 1→4 linkedβ-(D)-guluronic(G)andα-(L)-mannuronic(M)acid residues arranged as-M-G-sequences randomly[14].SA is one of the most popular anionic polymers used in biomedical applications due to its versatile properties such as biocompatibility,hydrophilicity,biodegradability,non-toxicity and good potentiality in drug delivery applications[15,16].SA forms three-dimensional network hydrogel by electrostatic attraction between the guluronic acid residues through the exchange of sodium ions of guluronic acid with the multivalent cations such as Ca2+,Mg2+,Ba2+,Sr2+and makes egg-box structure[17].However,the calcium alginate hydrogels have few limitations such as non-controlled swelling properties and drug release profile due to the inherent stiffness of the Ca2+alginate[18].In order to control the release rate and swelling properties,calcium alginates beads are blended with other polymers like chitosan,locust bean gum and gelatin to form a complex structure through hydrogen bonding,ion-ion,and dipoleion interactions.A few articles have reported that MMT containing alginate polymer blend microbeads acts as good carriers for drug delivery applications[19-21].

    Curcumin(CUR)(Fig.1),an active bio-compound obtained from the yellow spice of Curcuma longa(Indian spice),is a natural pigment with antibacterial,antifungal antiviral,anti-inflammatory,anti-malarial,antioxidant,anti-mutagenic effects and wound healing properties,and it enhances anti-tumour activity against different types of cancer cells[22-25].However,therapeutic use of CUR is limited due to its poor water solubility,short life and low bioavailability[26,27].Generally simple ionotropic gelation techniques are used to encapsulate the hydrophilic drugs in hydrogel beads,but this technique gives low encapsulation efficiency of CUR due to its low water solubility.Therefore,in order to increase the drug encapsulation efficiency of CUR,MMT clay material was used in the present study;the muco-adhesiveness of MMT clay material facilitates the intercalation of drug molecules,which inturn increases the bioavailability of CUR.

    Fig.1.Chemical structure of curcumin.

    In the present study,we focused on fabricating heteroionic SA/MMT beads by doping calcium-alginate matrix with different multivalent cations(Mg2+,Ba2+and Al3+).The incorporation of multivalent ions into calcium-alginate matrix would modify swelling property and release the rate of bio-active molecules.The developed microbeads were characterized by different techniques such as Fourier transform infrared spectroscopy(FTIR),X-ray diffraction(X-RD),differential scanning calorimetry(DSC),thermo gravimetric analysis(TGA)and scanning electron microscopy(SEM).The swelling study and in vitro drug release kinetics were performed in both simulated intestinal fluid and simulated gastricfluid at 37°C and the results are presented here.

    2.Experimental

    2.1.Materials

    MMT was purchased from Sigma-Aldrich(USA).SA and calcium chloride,magnesium chloride,barium chloride and aluminium chloride were purchased from SD Fine Chemicals,Mumbai,India.CUR was purchased from Loba Chemicals,Mumbai,India.Water used was of high purity grade after double distillation.

    2.2.Preparation of SA/MMT microbeads

    200mg of SA was weighed and dissolved in water under constant stirringovernight.To this solution,200mg of MMT was added and stirred well.A required amount of CUR was added and stirred to obtain a homogeneous solution.Afterwards the suspension was placed in sonicator for 5min to get homogenous suspension.The resulting suspension was slowly dropped into different ratios of CaCl2,MgCl2,AlCl3and BaCl2solution(as given in Table 1),where the spherical beads formed instantly were kept for 40 min.The obtained wet beads were collected by decantation,washed three times with double distilled water to remove the drug attached on the bead surface,and finally were dried in air overnight at room temperature.

    Table 1 Formulation and composition of all samples used for the studies.

    2.3.Characterization methods

    2.3.1.Intercalation kinetics

    To estimate the maximum time requiredfor intercalation of CUR with MMT,50 mg of CUR and 100mg of MMT were weighed and dissolved in 20 mL of double distilled water with continuous stirring at 37°C.At regular intervals of time(15,30,60,90,120,240 and 360min),the drug solution was filtered and the concentration of CUR was assayed using UV spectrophotometerat fixedλmaxvalue of 470.00 nm.

    2.3.2.Fourier transform infrared spectroscopy(FTIR)

    FTIR spectra of SA,CUR,MMT,placebo microbeads,and drug loaded microbeads were measured with an FTIR spectrophotometer(model Bomem MB-3000,with Horizon MB?FTIR software).The samples were mixed with KBr and compressed into pellets and scanned from 400 to 4000cm-1to find out the possible chemical interactions between polymers and drug.

    2.3.3.Differential scanning calorimetry(DSC)

    DSC curves of CUR,MMT,placebo microbeads,and drug loaded microbeadswererecordedusingthermogravimetry analyzer(Rheometric Scientific,Model DSC-SP,UK).The analysis was performed by heating the sample from 40 to 600°C at the heating rate of 10°C/min under nitrogen atmosphere.

    2.3.4.Thermogravimetric analysis(TGA)

    TGA of CUR,MMT,placebo microbeads,and drug loaded microbeads was carried out using thermogravimetry analyzer(Rheometric Scientific,Model DSC-SP,UK).About 5-7mg of sample was placed into alumina crucible and the thermo grams were recorded between 40°C to 600°C at a heating rate of 10°C/min under nitrogen atmosphere.

    2.3.5.X-ray diffraction(X-RD)

    The X-ray diffraction of CUR,placebo microbeads,and drug loaded microbeads was performed by a wide angle X-ray scattering diffractometer(Panalytical X-ray Diffractometer,model-X'pert Pro)with CuKα radiation(λ =1.54060)at a scanning rate of 10°/min to determine the crystallinity.

    2.3.6.Scanning electron microscopy(SEM)

    The morphologicalcharacterization ofmicrobeads was observed using SEM(JOEL MODEL JSM 840A)with an accelerated voltage of 20 kV equipped with an EDAX detector.

    2.3.7.Swelling measurements

    The swelling behaviour of different formulations was determined gravimetrically in simulated intestinal fluid(pH 7.4)and simulated gastric fluid(pH 1.2)at 37°C.The percentage of equilibrium swelling degree was calculated using the following equation:

    Where Wsis the weight of swollen beads and Wdis the weight of dry beads.

    2.4.Encapsulation efficiency(EE)

    A known mass of drug loaded microbeads(40mg)was immersed into 100 mL of phosphate buffer solution(pH 7.4,containing 5% absolute ethyl alcohol)for 24h and then the suspension was agitated with an agate mortar to ensure the complete extraction of CUR from the microbeads and filtered through filter paper.The drug solution was analyzed by ultraviolet(UV)spectrophotometer(LabIndia,Mumbai,India)at theλmaxof 470.00nm with placebo microbeads used as a blank correction.Concentration of drug was determined using calibration curve constructed by a series of CUR standard solutions.Percentage of EE was determined using the following formula.

    2.5.In vitro drug release studies

    In vitro drug release studies of different formulations were carried out by using a dissolution tester(Lab India,Mumbai,India)containing eight baskets each with 900 mL of phosphate buffer solution(PBS)maintained at 37°C,at a rotation speed of 50 rpm to replicate intestinal fluid (pH 7.4)and gastric fluid (pH 1.2)atmosphere,respectively.100mg of the CUR loaded microbeads were taken in a dialysis bags for the drug release studies.At regular intervals of time,5 mL aliquot samples were withdrawn,and analyzed using UV spectrophotometer at fixedλmaxvalue of 470.00 nm,and the released drug amount was obtained by using concentration versus absorbance calibration curves.The withdrawn aliquot samples werereplenished with equal volumes of PBS to simulate physiological conditions.The sink conditions were maintained throughout the release study.

    3.Results and discussion

    3.1.Intercalation kinetics

    From Fig.2,it was observed that 15.11% of CUR was intercalated with MMT electrostatically within 90 min and remained constant up to 360 min.Therefore,we should keep the time of 90 min for interaction between CUR and MMT to avoid partial interaction in the following experiments.

    Fig.2.Effect of time for intercalation of CUR with MMT.

    3.2.FTIR spectral analysis

    The FTIR spectra of CUR,MMT,placebo microbeads,and drug loaded SA/MMT microbeads are presented in Fig.3.The FTIR spectra of CUR show a characteristic broad peak at 3496 cm-1,which corresponds to phenolic O-H stretching vibrations,a peak at 2923 cm-1is assigned to aromatic C-H stretching vibrations,a peak at 1596 cm-1corresponds to stretching vibration of benzene ring skeleton,a peak at 1513 cm-1corresponds to mixed(C=O)and(C=C)vibration,a peak at 1272cm-1is assigned to Ar-O stretching vibrations.The FTIR spectra of MMTshow distinct peaks at 3615 and 3545cm-1that are assigned to O-H stretching vibrations of Si-OH and Al-Al-OH,a peak at 3371cm-1corresponds to H-O-H stretching vibrations of interlayer water,a peak at 1627cm-1corresponds to H-O-H bending mode of adsorbed water,peaks at 1381,1016 and 931cm-1correspond to(Si-O-Si)stretching frequency and peaks at 794 and 509 cm-1are assigned to(O-Si-O)and(Al-Si-O)bending vibrations[28-30].In case of placebo microbeads a peak at 3371cm-1is responsible for O-H stretching vibrations,a peak at 1616 cm-1is assigned to C=O stretching vibrations,and a peak at 1393 cm-1corresponds to COO-symmetric vibration.On comparing the drug loaded SA/MMT microbeads and pristine MMT,the Si-OH stretching of MMT at 3615 cm-1disappeared in the spectra of drug SA/MMT microbeads,which confirmed that active sites of polymer matrix interacted with MMT.On comparing the placebo microbeads and drug loaded SA/MMT microbeads,the O-H stretching vibration of SA at 3371 cm-1was shifted to lower side in dug loaded SA/MMT microbeads,which confirms that intermolecular hydrogen bonding takes place between polymer matrix and MMT clay[19].And also a peak at 678cm-1indicates the presence of MMTcomposition in the matrix.The peak at 1616 cm-1in placebo microbeads was shifted to 1606 cm-1in drug SA/MMT loaded microbeads,which indicates that drug should have been loaded in the microbeads.

    Fig.3.FTIR spectrum of(A)CUR,MMT,placebo microbeads,and(B)drug SA/MMT loaded microbeads(Ca,CaMg,CaBa and CaAl).

    3.3.DSC analysis

    To investigate the interaction of polymer matrix with MMT as well as crystalline nature of encapsulated CUR in the polymer matrix,DSC analysis was performed and the DSC curves are displayed in Fig.4.The thermogram of placebo microbeads(Fig.4b)shows a peak at 208°C,whereas drug loaded SA/MMT microbeads(Fig.4c)also showed the same peak as that of placebo microbeads but a slight variation was observed.And also another peak was observed in drug loaded SA/MMT microbeads at 164°C.This confirmed that the interaction took place between active sites of polymer matrix and MMT,because MMT(Fig.4d)also showed the same peak but a slight variationwasobserved.Fig.4a shows a sharp peak at 189°C for pure CUR,whereassuch peak was notobservedin drug loaded SA/MMT microbeads(Fig.4c),which confirmed that the drug was molecularly dispersed in the microbeads.

    Fig.4.DSC curves of(a)CUR,(b)placebo microbeads,(c)drug loaded SA/MMT microbeads(CaMg formulation)and(d)pristine MMT.

    3.4.TGA analysis

    The TGA thermograms of CUR(a),placebo microbeads(b),drug loaded SA/MMT microbeads(c)and MMT(d)are displayed in Fig.5.As shown in Fig.5a,CUR should remain stable up to 169°C;after that it followed mass loss and became maximum at 393°C due to total degradation of the compound.In case of MMT(Fig.5d),weight loss of 11% at 94 and 154°C corresponds to dehydration of water molecules absorbed in pores and between the silicate layers[31],followed by weight loss of 4% in the region of 160-700°C,corresponds to loss of structural water.The thermal decomposition of placebo microbeads(Fig.5b)occurs in three consecutive steps.The first weight loss of 17% was found between 41 and 176°C,which is ascribed to the evaporation of adsorbed water.The second weight loss step was observed between 180 and 299°C with a loss of 23%,indicating the decomposition of the salt from the polymer network.The last step was observed in the region of 304-600°C with a loss of 11% due to decomposition of polymer network which resulted in the formation of sodium carbonate[21].In the case of drug loaded SA/MMT microbeads(Fig.5c)four weight loss steps were observed.The first weight loss of 11% was observed in the region of 42-181°C,followed by weight loss of 11% between the region of 185-254°C due to degradation of CUR.The next two steps with weight loss of 10% and 17% were found in the region of 259-304°C and 309-600°C,which corresponds to decomposition of polymer network.The TGA results suggest that CUR loaded SA/MMT matrix showed an overall improvement in the thermal stability of microbeads.

    Fig.5.TGA curves of(a)CUR,(b)placebo microbeads,(c)drug loaded SA/MMT microbeads(CaMg formulation)and(d)MMT.

    3.5.X-RD analysis

    To investigate the molecular dispersion of CUR in polymer matrix,X-RD analysis was performed and the diffractograms of CUR(a),placebo microbeads(b)and CaMg microbeads(c)are presented in Fig.6.The diffractogram of CUR(Fig.6a)shows characteristic peaks in the 2θ region of 12-28°because of its crystallinity.However,these characteristic peaks disappeared in drug loaded SA/MMT microbeads(Fig.6c),which indicates that the CUR changes its state from crystalline to amorphous nature.These results suggest that CUR was molecularly dispersed in the polymeric network.

    Fig.6.X-RD patterns of(a)CUR,(b)placebo microbeads,and(c)drug loaded SA/MMT microbeads(CaMg formulation).

    3.6.SEM analysis

    To study the morphology of microbeads,SEM analysis was performed and the topographical images are displayed in Fig.7.From Fig.7,it can be clearly observed that the microbeads have rough surface with visible wrinkles.In all formulations it was observed that high roughness is present,which indicates the presence of MMT platelets on outer surface of microbeads.On comparing the Ca(Fig.7A)and CaMg(Fig.7B)formulations,more porous nature was observed in CaMg than in Ca.The porous nature is due to the presence of Mg2+cations in the network,which forms less rigid structure.The reason is that Mg2+ions are less strongly coordinated with the alginate than Ca2+ions and therefore more porous nature was observed.This was further correlated with the swelling and drug release studies.A similar observation was reported by Sanchez-Ballester et al.[18].In contrast,the formulations CaBa(Fig.7C)and CaAl(Fig.7D)show less porous nature,which is due to formation of more rigid structure.From the results of SEM images,theaveragesizeofmicrobeadswasfound tobe 800-1100μm.

    Fig.7.SEM images show effect of MMT and divalent ions on polymer network:(A)Ca,(B)CaMg,(C)CaBa and(D)CaAl.

    3.7.Energy-dispersive X-ray spectra(EDS)analysis

    To evaluate the elemental composition of Ca,CaMg,CaBa and CaAl microbeads,EDS analysis was performed and their spectra are presented in Fig.8.From Fig.8,it can be clearly observed that in all formulations Al and Si peaks were observed,which confirmed that MMT presented in all formulations.EDS spectra of CaMg,CaBa and CaAl show Mg,Ba and Al peaks,respectively.This confirmed that multivalent ions were incorporated in the calcium alginate network.

    Fig.8.EDS spectra of(A)Ca,(B)CaMg,(C)CaBa and(D)CaAl.

    3.8.Swelling measurements

    Swelling degree is one of the important factors of SA beads for control release in drug delivery systems[32].In order to determine the swelling degree of all formulations,swelling experiments were performed in simulated intestinal fluid(pH 7.4)and simulated gastric fluid(pH 1.2)at 37°C and the results are presented in Fig.9.From Fig.9,it is clearly observed that the swelling degree is more at pH 7.4 than at pH 1.2.Therefore,SA beads are good promising carriers to deliver drug molecules in the intestines and to avoid gastric release of drugs.

    From Fig.9 it is observed that the swelling degree of all formulations differed based on the amount of multivalent ions present in the beads.The maximum swelling degree of all profiles was observed at60 min and thereafterslightly decreased.On comparing all the formulations,CaMg beads showed high percentage of water uptake because the Mg2+ions formed less rigid network than the all other formulations;therefore,the water molecules entered the beads through the pores;thereby the percentage of swelling increased.On the other hand,the CaBa and CaAl showed less percentage of water uptake because the Ba2+and Al3+formed a tight rigid network,such that the polymer network was unable to allow a large number of water molecules,so the swelling decreased.

    Fig.9.Swelling studies of Ca,CaMg,CaBa and CaAl at(A)pH 7.4 and(B)pH 1.2 at 37°C.

    3.9.Encapsulation efficiency(EE)

    The EE of CUR varied for all formulations depending on the presence of multivalent ions in beads and the results are listed in Table 2.The EE decreased in CaBa and CaAl formulations,because the network structure of these formulations was more rigid,thereby swelling to less extent,which resulted in low EE values.In Ca and CaMg formulations,the EE increased because it formed less rigid network,thereby swelling to more extent,which resulted in high EE values.

    Table 2 Release kinetics parameters at pH 7.4 and encapsulation efficiency(EE)of all samples.

    3.10.In vitro drug release studies

    The in vitro release studies were carried out in simulated intestinal fluid(pH 7.4)and simulated gastric fluid(pH 1.2)at 37°C and the results are displayed in Fig.10.From Fig.10,it is observed that the release rate is higherat pH 7.4 than at pH 1.2,because at pH 7.4 the carboxylic groups show less interactions with buffer media,allowing the network to be more loose;hence the entrapped drug molecules easily leaches out from the network.From Fig.10,it is observed that similar to swelling behaviour,in vitro release rate also shows the same type of results.In the case of Ca and CaMg profiles,CaMg shows high release than Ca,which is due to its more porous structure and also formation of weaker interaction between Mg2+ions and alginates.This would contribute to CaMg beads’disintegrating more readily in the dissolution medium and drug molecules easily leach out from the polymer matrix,whereas in the case of CaBa and CaAl the drug release rate decreases because it forms more rigid structure;hence CUR molecules have less chance to escape out from the matrix.

    Fig.10.CUR release of Ca,CaMg,CaBa and CaAl at(A)pH 7.4 and(B)pH 1.2 at 37°C.

    The obtained in vitro drug release data of CUR loaded microbeads at pH 7.4 were fitted into different mathematical models such as zero order, first order and Higuchi.The release rate constant and correlation coefficient of all formulations are shown in Table 2.The correlation coefficient values of CUR loaded microbeads followed neither zero order nor first order,but the correlation coefficient values were close to the Higuchi model.Therefore,the drug release kinetics of CUR loaded microbeads followed Higuchi model.According to Higuchi model,the release of drug from the microbeads involves the penetration of liquid into the matrix and dissolving of the drug,which then diffuses the drug into the exterior liquid through pores or intestinal channels.However,this type of phenomenon is observed in hydrophilic matrix system.Therefore,the drug release rate of CUR loaded microbeads shows the phenomenon of swelling and erosion of the polymer simultaneously.

    To understand the drug release mechanism,the data obtained from in vitro drug release studies in PBS(7.4)were fitted into the following Korsmeyer-Peppas equation[33].

    Where,Mtis the cumulative release of CUR at time t,Mαis the total amount of CUR in the matrix,k is a characteristic release constant of the drug-polymer system and n is the release exponent indicating the type of drug release mechanism.The results of n and k are listed in Table 2.For spherical drug carriers,if n<0.43,the drug diffuses from the polymer matrix according to Fickian diffusion;if 0.430.85,the mode of drug release follows the super Case-II diffusion.In the present data n values were obtained in the range of 0.550-0.572,which indicates non-Fickian type of diffusion process.

    4.Conclusion

    In the present work,dual cross-linked SA/MMT microbeads were fabricated by ionotropic gelation method for the extended release of CUR.FTIR confirmed the interaction between drug molecules and polymer network.DSC,TGA and X-RD studies confirmed the chemical stability and molecular level dispersion of CUR in microbeads,respectively.SEM studies confirmed the presence of MMT platelets on the outer surface of microbeads.The effect of the incorporation of multivalent ions on crosslinking of SA chains,swelling studies and drug release studies were thoroughly examined.The CaMg beads showed highest percentage of water uptake at pH 7.4 compared to other formulations due to the formation of less rigid network,whereas the CaAl beads showed less percentage of water uptake due to the formation of tight rigid network.The drug release profiles at pH 7.4 were correlated with the swelling results,and the CaMg and CaAl beads showed the highest and lowest drug release rates,respectively.The in vitro results were fitted into peppas equation and the results showed that the drug mechanism followed non-Fickian type of diffusion process.Based on all these results,the developed microbeads are considered to be potentially good drug carriers for extended release of CUR.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    The authors C.Madhavi and K.Chowdoji Rao thank UGC-BSR,New Delhi,India,for the financial support provided.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jpha.2020.07.002.

    欧美激情久久久久久爽电影| 只有这里有精品99| 亚洲国产欧美在线一区| 精品久久久精品久久久| 国产在线一区二区三区精| 日韩国内少妇激情av| 欧美日韩一区二区视频在线观看视频在线 | 成年女人看的毛片在线观看| 成人国产麻豆网| 在线观看美女被高潮喷水网站| 国产 精品1| 国产一区有黄有色的免费视频| 国产精品一区二区三区四区免费观看| 欧美精品国产亚洲| 国产免费一区二区三区四区乱码| 精品亚洲乱码少妇综合久久| 欧美日韩综合久久久久久| eeuss影院久久| 九色成人免费人妻av| 亚洲色图综合在线观看| 天堂俺去俺来也www色官网| 久热这里只有精品99| 欧美日韩一区二区视频在线观看视频在线 | 久久久精品欧美日韩精品| 在线免费十八禁| 亚洲人成网站高清观看| 国产女主播在线喷水免费视频网站| av免费在线看不卡| 一级毛片我不卡| 免费黄网站久久成人精品| 免费看光身美女| 男人舔奶头视频| 2018国产大陆天天弄谢| 久久久久国产精品人妻一区二区| 亚洲精品一二三| 交换朋友夫妻互换小说| 国国产精品蜜臀av免费| 秋霞在线观看毛片| kizo精华| 少妇人妻一区二区三区视频| 午夜福利视频精品| 蜜桃久久精品国产亚洲av| 成人高潮视频无遮挡免费网站| 欧美日韩在线观看h| 听说在线观看完整版免费高清| 亚洲人与动物交配视频| 欧美+日韩+精品| 欧美亚洲 丝袜 人妻 在线| 欧美另类一区| 欧美3d第一页| 久久亚洲国产成人精品v| 男人舔奶头视频| 亚洲最大成人av| 亚洲成人久久爱视频| 久久精品熟女亚洲av麻豆精品| 国产成人精品久久久久久| 亚洲一级一片aⅴ在线观看| 能在线免费看毛片的网站| 亚洲精品乱码久久久v下载方式| 国产成人精品福利久久| 久久久久久久久大av| 大片免费播放器 马上看| 蜜桃久久精品国产亚洲av| 国产黄片美女视频| 国产欧美另类精品又又久久亚洲欧美| 下体分泌物呈黄色| 免费av观看视频| 国产黄频视频在线观看| 亚洲国产精品成人综合色| 免费看日本二区| 不卡视频在线观看欧美| 日本wwww免费看| 欧美成人a在线观看| 九九爱精品视频在线观看| 99久久九九国产精品国产免费| 少妇丰满av| 国产成人免费无遮挡视频| 亚洲美女搞黄在线观看| 国产有黄有色有爽视频| 激情五月婷婷亚洲| 婷婷色av中文字幕| 欧美最新免费一区二区三区| 亚洲国产精品专区欧美| 在线 av 中文字幕| 国产在线男女| 国产精品偷伦视频观看了| 18禁在线播放成人免费| 亚洲久久久久久中文字幕| 成人免费观看视频高清| 国产成人91sexporn| 禁无遮挡网站| 18+在线观看网站| 蜜臀久久99精品久久宅男| 国产精品爽爽va在线观看网站| 99热全是精品| 国产精品.久久久| 久久99热这里只有精品18| 欧美成人一区二区免费高清观看| 日本一本二区三区精品| 2022亚洲国产成人精品| 免费观看无遮挡的男女| 一级爰片在线观看| 亚洲精品久久久久久婷婷小说| 99九九线精品视频在线观看视频| 狠狠精品人妻久久久久久综合| h日本视频在线播放| 亚洲国产精品专区欧美| 少妇人妻精品综合一区二区| 久久久久国产网址| 亚洲av免费在线观看| 在线 av 中文字幕| 免费av观看视频| 97人妻精品一区二区三区麻豆| 女人十人毛片免费观看3o分钟| 免费黄网站久久成人精品| 国产熟女欧美一区二区| 91久久精品国产一区二区三区| 国产视频内射| 三级男女做爰猛烈吃奶摸视频| 国产精品一及| 久久久久国产精品人妻一区二区| 久久久久久久久久久免费av| 国国产精品蜜臀av免费| 又爽又黄无遮挡网站| 精品酒店卫生间| 国产成人91sexporn| 在线 av 中文字幕| 看免费成人av毛片| 亚洲自偷自拍三级| 麻豆乱淫一区二区| 免费大片黄手机在线观看| 草草在线视频免费看| 亚洲欧美清纯卡通| 日韩成人伦理影院| 亚洲一级一片aⅴ在线观看| 国产 精品1| 国产伦在线观看视频一区| 又大又黄又爽视频免费| 亚洲国产精品成人久久小说| 免费看光身美女| 亚洲av男天堂| 久久久久精品性色| 在线精品无人区一区二区三 | 大片免费播放器 马上看| 亚洲av成人精品一区久久| 神马国产精品三级电影在线观看| 国模一区二区三区四区视频| 中文资源天堂在线| 免费电影在线观看免费观看| 99热这里只有是精品在线观看| 国产精品.久久久| 亚洲aⅴ乱码一区二区在线播放| 欧美成人午夜免费资源| 久久久久久九九精品二区国产| 成人漫画全彩无遮挡| 三级国产精品片| 黄色一级大片看看| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩另类电影网站 | 国产成人freesex在线| 少妇人妻 视频| 国产老妇伦熟女老妇高清| 国产成人午夜福利电影在线观看| 国产精品av视频在线免费观看| 色婷婷久久久亚洲欧美| 99久久九九国产精品国产免费| 亚洲精品色激情综合| 51国产日韩欧美| 国产成人freesex在线| 国产精品国产三级国产av玫瑰| 国产精品一区www在线观看| 国产伦理片在线播放av一区| 各种免费的搞黄视频| 午夜激情福利司机影院| 欧美区成人在线视频| av一本久久久久| 中文字幕亚洲精品专区| 亚洲av中文av极速乱| 亚洲欧美日韩卡通动漫| 日韩一区二区视频免费看| 国产精品无大码| 最近手机中文字幕大全| 日韩国内少妇激情av| 爱豆传媒免费全集在线观看| 国产精品久久久久久久久免| 肉色欧美久久久久久久蜜桃 | 国产黄片美女视频| 我要看日韩黄色一级片| 亚洲av成人精品一区久久| av播播在线观看一区| 午夜老司机福利剧场| 久热这里只有精品99| 久久精品久久久久久噜噜老黄| 久久久久网色| 亚洲经典国产精华液单| 亚洲av电影在线观看一区二区三区 | 91午夜精品亚洲一区二区三区| 国产成人免费无遮挡视频| 免费人成在线观看视频色| 3wmmmm亚洲av在线观看| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站| 91在线精品国自产拍蜜月| 男女那种视频在线观看| 日韩中字成人| 久久精品夜色国产| 亚洲国产高清在线一区二区三| 噜噜噜噜噜久久久久久91| 九九爱精品视频在线观看| 男女啪啪激烈高潮av片| 亚洲四区av| 天天一区二区日本电影三级| 免费不卡的大黄色大毛片视频在线观看| 少妇熟女欧美另类| 日本黄大片高清| 91在线精品国自产拍蜜月| 青春草视频在线免费观看| 哪个播放器可以免费观看大片| 亚洲欧美清纯卡通| 伦精品一区二区三区| 在线观看三级黄色| 亚洲精品中文字幕在线视频 | 免费人成在线观看视频色| eeuss影院久久| 熟女人妻精品中文字幕| 亚洲欧美精品自产自拍| 综合色丁香网| 欧美高清成人免费视频www| 精品久久久久久电影网| 美女高潮的动态| 日本黄大片高清| 丝袜美腿在线中文| 日日摸夜夜添夜夜添av毛片| 亚洲欧洲日产国产| 天堂网av新在线| av免费观看日本| 国产人妻一区二区三区在| 久久女婷五月综合色啪小说 | 中文字幕久久专区| www.色视频.com| 韩国av在线不卡| 九九在线视频观看精品| 久久精品久久精品一区二区三区| 久久亚洲国产成人精品v| 爱豆传媒免费全集在线观看| 中文字幕亚洲精品专区| 欧美日韩在线观看h| 久久久久网色| 亚洲电影在线观看av| 国产精品女同一区二区软件| 国产又色又爽无遮挡免| 午夜激情久久久久久久| 国产色爽女视频免费观看| 亚洲成人精品中文字幕电影| 亚洲欧美一区二区三区国产| 国内精品美女久久久久久| 亚洲色图综合在线观看| 欧美激情国产日韩精品一区| 男女国产视频网站| 成人综合一区亚洲| 国产一区亚洲一区在线观看| 亚洲精品一区蜜桃| 欧美日韩视频精品一区| 天堂俺去俺来也www色官网| 三级国产精品欧美在线观看| 久久午夜福利片| 尾随美女入室| 嫩草影院精品99| 男女那种视频在线观看| 99久久中文字幕三级久久日本| 干丝袜人妻中文字幕| 永久免费av网站大全| 欧美日韩国产mv在线观看视频 | 97人妻精品一区二区三区麻豆| 国产 精品1| 你懂的网址亚洲精品在线观看| 天美传媒精品一区二区| 欧美日本视频| 一级毛片电影观看| 免费人成在线观看视频色| 亚洲av.av天堂| 中文字幕制服av| 国产精品国产三级专区第一集| 毛片女人毛片| 美女脱内裤让男人舔精品视频| av在线播放精品| 国产在线男女| 久久久精品94久久精品| 男人添女人高潮全过程视频| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 中文天堂在线官网| 最近手机中文字幕大全| 美女脱内裤让男人舔精品视频| 久久6这里有精品| 欧美高清成人免费视频www| 国产国拍精品亚洲av在线观看| av线在线观看网站| 国产精品秋霞免费鲁丝片| 久久精品久久久久久久性| 两个人的视频大全免费| 国内精品宾馆在线| 国产免费福利视频在线观看| 亚洲av免费高清在线观看| 成人免费观看视频高清| 男男h啪啪无遮挡| 国产精品不卡视频一区二区| 中国国产av一级| 精品视频人人做人人爽| 国产成人精品婷婷| 青春草视频在线免费观看| 中国三级夫妇交换| 又黄又爽又刺激的免费视频.| 最新中文字幕久久久久| 欧美日韩亚洲高清精品| 亚洲精品视频女| 一本久久精品| av.在线天堂| 毛片一级片免费看久久久久| 国产精品一区二区三区四区免费观看| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放| 综合色丁香网| 天天一区二区日本电影三级| 一区二区三区四区激情视频| 国产成人freesex在线| av在线天堂中文字幕| 久久这里有精品视频免费| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品国产av蜜桃| 街头女战士在线观看网站| 51国产日韩欧美| 男的添女的下面高潮视频| 日韩亚洲欧美综合| 在线看a的网站| 久久国内精品自在自线图片| 夫妻性生交免费视频一级片| 女人十人毛片免费观看3o分钟| 亚洲在久久综合| 舔av片在线| 秋霞伦理黄片| 国产一区二区亚洲精品在线观看| 欧美变态另类bdsm刘玥| 中国国产av一级| kizo精华| 午夜亚洲福利在线播放| 精品少妇黑人巨大在线播放| 久久久久久久精品精品| 国产一区二区三区综合在线观看 | 亚洲精品456在线播放app| 蜜臀久久99精品久久宅男| 在线免费观看不下载黄p国产| 日韩欧美 国产精品| 久久久久网色| 久久精品人妻少妇| 建设人人有责人人尽责人人享有的 | 一级毛片黄色毛片免费观看视频| 亚洲最大成人手机在线| 国内精品美女久久久久久| 亚洲最大成人手机在线| 91久久精品电影网| 国产精品99久久99久久久不卡 | 蜜桃久久精品国产亚洲av| 在线观看一区二区三区| 高清在线视频一区二区三区| 亚洲av免费在线观看| 又黄又爽又刺激的免费视频.| 丰满人妻一区二区三区视频av| 伊人久久国产一区二区| 亚洲精品色激情综合| av在线老鸭窝| 亚洲av免费在线观看| 亚洲人成网站在线播| 婷婷色麻豆天堂久久| 一区二区三区乱码不卡18| 夜夜爽夜夜爽视频| 国产午夜精品久久久久久一区二区三区| 春色校园在线视频观看| 亚洲无线观看免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品一区蜜桃| 亚洲精品亚洲一区二区| 女人十人毛片免费观看3o分钟| 欧美潮喷喷水| 久久久久久九九精品二区国产| 观看美女的网站| 亚洲精品一区蜜桃| 亚洲激情五月婷婷啪啪| 久久久精品欧美日韩精品| 中文字幕人妻熟人妻熟丝袜美| 男人舔奶头视频| 丰满少妇做爰视频| av在线老鸭窝| 日本色播在线视频| 特大巨黑吊av在线直播| 久久久成人免费电影| 99久久人妻综合| 国产精品一及| 亚洲丝袜综合中文字幕| 综合色av麻豆| 天天一区二区日本电影三级| 一区二区三区免费毛片| 我的老师免费观看完整版| 国产 精品1| 午夜福利在线在线| 女的被弄到高潮叫床怎么办| 亚洲美女视频黄频| 亚洲人成网站高清观看| 久久精品综合一区二区三区| 国产精品99久久久久久久久| 免费观看无遮挡的男女| 日本免费在线观看一区| 亚洲国产精品成人久久小说| 自拍偷自拍亚洲精品老妇| 91精品国产九色| 少妇被粗大猛烈的视频| 日日啪夜夜爽| 日韩av免费高清视频| 久热这里只有精品99| 听说在线观看完整版免费高清| 伦精品一区二区三区| 免费观看无遮挡的男女| 肉色欧美久久久久久久蜜桃 | av国产免费在线观看| 高清毛片免费看| 精品国产乱码久久久久久小说| 欧美丝袜亚洲另类| 夫妻性生交免费视频一级片| 自拍欧美九色日韩亚洲蝌蚪91 | 看非洲黑人一级黄片| 亚洲人与动物交配视频| 五月开心婷婷网| 免费不卡的大黄色大毛片视频在线观看| 亚洲人与动物交配视频| 内射极品少妇av片p| 精品午夜福利在线看| 韩国av在线不卡| 搞女人的毛片| 国产一区亚洲一区在线观看| 欧美3d第一页| 欧美三级亚洲精品| 中国国产av一级| 亚洲av在线观看美女高潮| 一区二区三区精品91| 色婷婷久久久亚洲欧美| av.在线天堂| 乱系列少妇在线播放| 汤姆久久久久久久影院中文字幕| 晚上一个人看的免费电影| 视频区图区小说| 乱码一卡2卡4卡精品| 麻豆国产97在线/欧美| av在线亚洲专区| 国产成人免费无遮挡视频| 高清午夜精品一区二区三区| 午夜激情福利司机影院| 丰满少妇做爰视频| 国产女主播在线喷水免费视频网站| 久久人人爽人人片av| 日日啪夜夜撸| 秋霞在线观看毛片| 亚洲精品日韩在线中文字幕| 亚洲最大成人av| 欧美成人a在线观看| 免费看不卡的av| xxx大片免费视频| 又爽又黄a免费视频| 国产精品偷伦视频观看了| 最近的中文字幕免费完整| 两个人的视频大全免费| 最后的刺客免费高清国语| 亚洲色图av天堂| 少妇裸体淫交视频免费看高清| 人人妻人人看人人澡| 国产亚洲午夜精品一区二区久久 | 国产精品99久久久久久久久| 噜噜噜噜噜久久久久久91| 久久这里有精品视频免费| 中文精品一卡2卡3卡4更新| av黄色大香蕉| 少妇丰满av| 如何舔出高潮| 国产一区二区三区av在线| 天堂俺去俺来也www色官网| 国产极品天堂在线| 男人爽女人下面视频在线观看| 亚洲第一区二区三区不卡| 性色avwww在线观看| 国产精品一区二区性色av| 久久精品久久久久久久性| 久久久成人免费电影| 久久这里有精品视频免费| 免费观看性生交大片5| 精品国产一区二区三区久久久樱花 | 97在线视频观看| 99re6热这里在线精品视频| 一区二区三区乱码不卡18| 亚洲精品自拍成人| 亚洲精品国产av蜜桃| 国产亚洲91精品色在线| 久久久久久久久久成人| 免费观看的影片在线观看| 只有这里有精品99| 中文在线观看免费www的网站| a级毛片免费高清观看在线播放| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 亚洲精品日韩av片在线观看| 成人美女网站在线观看视频| 日本av手机在线免费观看| av.在线天堂| 天堂俺去俺来也www色官网| 99热网站在线观看| 亚洲人成网站在线观看播放| 免费观看无遮挡的男女| 久久久久久九九精品二区国产| 国产精品蜜桃在线观看| 2018国产大陆天天弄谢| 精品酒店卫生间| 国产精品久久久久久精品电影小说 | 美女国产视频在线观看| 国产又色又爽无遮挡免| 丰满少妇做爰视频| 亚洲精品,欧美精品| 乱码一卡2卡4卡精品| 青春草亚洲视频在线观看| .国产精品久久| 国产精品一区www在线观看| av国产免费在线观看| 男人添女人高潮全过程视频| 亚洲,一卡二卡三卡| 永久网站在线| 亚洲成人一二三区av| 亚洲自偷自拍三级| 日韩,欧美,国产一区二区三区| 好男人在线观看高清免费视频| 97人妻精品一区二区三区麻豆| 精品久久久噜噜| 国产黄片视频在线免费观看| 亚洲精品aⅴ在线观看| 三级国产精品欧美在线观看| 欧美xxxx黑人xx丫x性爽| 日日啪夜夜撸| 老女人水多毛片| av在线天堂中文字幕| 亚洲国产高清在线一区二区三| 热re99久久精品国产66热6| av免费在线看不卡| 在线精品无人区一区二区三 | 精品视频人人做人人爽| 国产成人免费无遮挡视频| 国产乱人视频| 高清毛片免费看| 国产精品爽爽va在线观看网站| 国产老妇女一区| 韩国av在线不卡| 丰满人妻一区二区三区视频av| 伦精品一区二区三区| 老司机影院成人| 少妇丰满av| 波多野结衣巨乳人妻| 欧美xxxx黑人xx丫x性爽| 只有这里有精品99| 欧美一级a爱片免费观看看| 亚洲色图综合在线观看| 国产乱人视频| videos熟女内射| 国产高清国产精品国产三级 | 国产熟女欧美一区二区| 在线 av 中文字幕| 亚洲精品久久久久久婷婷小说| 久久久久久久大尺度免费视频| 99热网站在线观看| 国产日韩欧美在线精品| 麻豆国产97在线/欧美| kizo精华| 日本爱情动作片www.在线观看| 内射极品少妇av片p| 激情 狠狠 欧美| 国产精品一区二区三区四区免费观看| 亚洲国产最新在线播放| 亚洲第一区二区三区不卡| 亚洲欧美中文字幕日韩二区| 久久热精品热| 久久99热6这里只有精品| 在线观看一区二区三区| videos熟女内射| 精品久久国产蜜桃| 欧美激情国产日韩精品一区| 天天躁日日操中文字幕| 国产老妇女一区| 久久久色成人| 国产精品福利在线免费观看| 国产久久久一区二区三区| 亚洲av福利一区| 日韩成人伦理影院| 毛片一级片免费看久久久久| 在线看a的网站| 日本黄色片子视频| 国产成人91sexporn| 欧美潮喷喷水| 日本av手机在线免费观看| 精品国产一区二区三区久久久樱花 | 精华霜和精华液先用哪个| 在线观看人妻少妇| av一本久久久久| 亚洲色图av天堂| 天堂俺去俺来也www色官网| 高清av免费在线| 最近中文字幕2019免费版| 欧美 日韩 精品 国产| 日韩强制内射视频| 大话2 男鬼变身卡| 免费黄网站久久成人精品| 一级毛片久久久久久久久女| 国产永久视频网站|