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

    Electrochemical and in silico approaches for liver metabolic oxidation of antitumor-active triazoloacridinone C-1305☆

    2020-09-04 09:33:24AgnieszkPoteDorotZelszczykZofiMzersk
    Journal of Pharmaceutical Analysis 2020年4期

    Agnieszk Pote?g,Dorot˙Zelszczyk,Zo fi Mzersk

    aDepartment of Pharmaceutical Technology and Biochemistry,Faculty of Chemistry,Gdańsk University of Technology,Gabriela Narutowicza St.11/12,Gdańsk,80-233,Poland

    bDepartment of Organic Chemistry,Faculty of Pharmacy,Jagiellonian University,Medyczna St.9,Krakw,30-688,Poland

    Keywords:

    Antitumor triazoloacridinone

    P450-catalyzed reactions

    Electrochemistry/mass spectrometry

    In silico metabolism Prediction

    Liver microsomal assay

    ABSTRACT

    5-Dimethylaminopropylamino-8-hydroxytriazoloacridinone(C-1305)is a promising antitumor compound developed in our laboratory.A better understanding of its metabolic transformations is still needed to explain the multidirectional mechanism of pharmacological action of triazoloacridinone derivatives at all.Thus,the aim of the current work was to predict oxidative pathways of C-1305 that would reflect its phase I metabolism.The multi-tool analysis of C-1305 metabolism included electrochemical conversion and in silico sites of metabolism predictions in relation to liver microsomal model.In the framework of the first approach,an electrochemical cell was coupled on-line to an electrospray ionization mass spectrometer.The effluent of the electrochemical cell was also injected onto a liquid chromatography column for the separation of different products formed prior to mass spectrometry analysis.In silico studies were performed using MetaSite software.Standard microsomal incubation was employed as a reference procedure.We found that C-1305 underwent electrochemical oxidation primarily on the dialkylaminoalkylamino moiety.An unknown N-dealkylated and hydroxylated C-1305 products have been identified.The electrochemical system was also able to simulate oxygenation reactions.Similar pattern of C-1305 metabolism has been predicted using in silico approach.Both proposed strategies showed high agreement in relation to the generated metabolic products of C-1305.Thus,we conclude that they can be considered as simple alternatives to enzymatic assays,affording time and cost efficiency.

    1.Introduction

    Tumor diseases are common and usually rapidly progressive disorders known worldwide[1].Although current medicine is developing very dynamically,the possibilities of treating tumors are still very limited.The success of antitumor therapy with the chemotherapeutic agent depends on the selection of the medication characterized by well-balanced profile between efficacy and safety[2].In this light,the studies on the possible metabolic pathways of a potent drug are extremely important for assessing its bioavailability,activity and toxicity in the human body.Insufficient knowledge about metabolic transformations of the drug is one of the most common causes of failure and the lack of complete results during clinical trials[3].

    In recent years the tools for the prediction of metabolic status of drugs have been developed and improved.The liver microsomal system is one of the best characterized models for drug biotransformation research which is simple to use[4,5].This subcellular fraction contains major drug-metabolizing enzymes,including the cytochrome P450(P450)family capable of catalysing the oxidative biotransformation of most drugs,including antitumor agents,and other xenobiotics[6].Thus,the compound incubation with liver microsomes generates oxidative phase I metabolites.However,the efficiency of the microsomal approach strongly depends on the specific expression of different P450 isoforms in each organism.Consequently,some differences between the metabolites predicted from liver microsome experiments and the metabolites found in the patient's body fluids after the intake of a drug are usually observed[6,7].Moreover,a targeted drug may be metabolized in the human body through different pathways and processes.

    In our studies we considered that in enzymatic metabolism systems,the processes of purification and separation prior to identification and quantification of metabolic products are expensive,laborious,and time-consuming[4].To replace these procedures,otherapproaches,which are able tosupport and improve drug metabolism studies,were developed and they can be classified in two categories:(i)electrochemical simulation of P450-mediated reactions and(ii)in silico prediction of drug metabolism.Electrochemistry(EC)is one of the classical methods to induce oxidationreduction reactionscatalyzed by liverP450 isoforms[8,9].Compared to the enzymatic assays,there are many advantages of this approach,such as low cost,fast response,simplicity,clean system,and ease of automation.These are universally valid and make electrochemical method extremely powerful in drug metabolism studies[10,11].Nowadays,electrochemical cells coupled online to analytical techniques such as mass spectrometry(MS)with or without previous liquid chromatographic(LC)separation(Fig.1)hold a great potential for simulation experiments and are intensively developed.EC(/LC)/MS techniques may provide exhaustive information about the properties of the electrochemically generated transformation products,even though the enzymatic mechanism of transformation differs from the electrochemical oxidative pathway[9,12,13].In turn,in silico screening strategy aims at predicting the most probable sites of metabolism for the investigated compound,without any need for experimental data.Therefore,this approach provides the main metabolite structures which can be formed in various human tissues(liver,skin,brain,and lungs).Then,the metabolites predicted are listed with a likelihood ranking[14,15].

    The aim of the present research was to investigate C-1305(5-dimethylaminopropylamino-8-hydroxytriazoloacridinone)oxidative transformations by the application of three different strategies for studying drug metabolism:experimental EC,in silico prediction,and incubation with rat liver microsomes(RLMs).The unknown products of phase I metabolism were subsequently identified by using electrospray ionization(ESI)MS.

    C-1305(a structure presented in the frame in Fig.1)is the most potent derivative among the antitumor-active triazoloacridinones developed in our laboratory[16].It exhibited significant and clearly differentiated cytotoxicity against a number of tumor cell lines of human and animal origin in tissue cultures[16,17],as was also shown in vitro screening scheme(panel of 60 different human tumor cell lines)in the National Cancer Institute(Bethesda,MD,USA).Moreover,C-1305 displayed high antitumor potency against P388 leukemia in mice and towards a wide spectrum of experimental transplantable solid tumors in animals such as colon and breast cancers[17,18].C-1305 is the close structural analogue of the clinically tested imidazoacridinone antitumor agent,C-1311[19-21];however,some differences in their chemical structure may determine the more advantageous pharmacological properties of the first one.Both compounds were synthesized as a result of rational drug design strategy based on several circumstances connected with combined characteristics of anthracyclines and acridines essential for biological activity[22,23].Overall,the specific structure of C-1305 determines its mode of action that involves several routes as physicochemical interaction with DNA[24],influence on cell cycle progression[25],induction of cell death[25],and others,in particular topoisomerase II inhibition[26,27]and covalent DNA crosslinking[24].The leading concept of the antitumor effect of C-1305,similar to the case of C-1311,is considered as the relationship between its biological activity and the susceptibility to enzymatic oxidation.

    The results regarding the capability of the EC to simulate the oxidative metabolism of C-1305 are discussed here for the first time.Then,electrochemical findings related to phase I metabolites of C-1305 were compared with those predicted from in silico manner and a classical enzymatic approach employing RLMs.Elucidation of the mechanism of the electrochemical oxidation would provide a deeper insight into the chemistry and EC of triazoloacridinones.What is more,the identification of the obtained products would enable their comparison with the C-1305 metabolites eliminated from the body.The results obtained from the proposed multi-tool analysis of C-1305 metabolism would be of great help in designing new compounds with a better pharmacokinetic profile,or in designing prodrugs where the compound needs to be metabolized in order to become active.

    2.Materials and methods

    2.1.Chemicals,reagents,and enzymes

    C-1305[16]was synthesized as dihydrochloride in our laboratory.The compound was of more than 98% pure as determined by liquid chromatography(LC)and nuclear magnetic resonance analyses.The following chemicals were purchased from Sigma-Aldrich(St.Louis,MO,USA):dipotassium phosphate(K2HPO4),and formic acid(HCOOH),magnesium chloride anhydrous(MgCl2),monopotassium phosphate(KH2PO4)and potassium hydroxide(KOH).Methanol(CH3OH;gradient grade for LC)andβ-nicotinamide adenine dinucleotide 2′-phosphate tetrasodium salt(NADPH)were obtained from Merck KGaA(Darmstadt,Germany).Ammonium formate(HCOONH4;reagent grade)was ordered from Fisher Scientific(Loughborough,UK).All other commercially available chemicals and reagents were of the highest possible grade available.Ultrapure water of 18 MΩ?cm of resistivity,used in all the experiments,was passed through a Milli-Q water purification system from Merck KGaA(Darmstadt,Germany).

    Fig.1.Schematic set-up of the system used for the electrochemical generation and identification of C-1305 oxidation products.Molecular structure and atomic numbering of antitumor-active triazoloacridinone C-1305 is shown in the frame.LC,liquid chromatography.

    Pooled RLMs from untreated,male Sprague-Dawley rats(protein concentration,20 mg/mL;P450 content,680 pmol/mg protein)were purchased from Tebu-bio(Le Perray-En-Yvelines,France).

    2.2.Instrumentation

    The simulation of the oxidative metabolism of C-1305 was accomplished in an amperometric electrochemical thin-layer cell equipped with a disc glassy carbon(GC)working electrode(φ=8 mm;A=0.502 cm2)and a Pd/H2reference electrode(reactor cell;Antec Leyden,Zoeterwoude,the Netherlands).Carbon-loaded polytetra fluoroethylene(PTFE)served as auxiliary electrode.The cell potentials were applied using a ROXY EC System(Antec Leyden).All potentials mentioned in this work were based on the reference electrode.The software used for controlling EC was Dialogue(Antec Leyden).

    The outlet of the electrochemical cell was interfaced into an ESI source of a quadrupole-time of flight(Q-TOF)mass spectrometer(Agilent Technologies,Santa Clara,CA,USA)inlet for on-line EC/MS analysis using PEEK tubing.For controlling the MS,MassHunter software(Agilent Technology)was used.The electrochemically generated oxidation products were also off-line injected onto the LC column,separated,and detected by ESI-MS.LC separations were performed with Waters Associates HPLC system(Waters Co.,Milford,MA,USA).It was equipped with a model 600 E system controller,a model 7725i Rheodyne injector,and a model 2996 photodiode array detector(DAD)controlled with Millennium software(Waters Co.).

    2.3.Electrochemical simulation of C-1305 oxidative metabolism

    The setup used for generation of C-1305 products by electrochemical oxidation and subsequent direct(on-line)MS or off-line LC/MS analysis is shown in Figs.1,2,3,4,5 and 6.For electrochemical oxidation,10μM C-1305 in 0.1% HCOOH in water/CH3OH electrolyte(50/50,V/V)was conducted through the EC cell using a dual piston syringe pump model SP2-ROXY(Antec Leyden)at a flow rate of 20μL/min.The effluent from the EC was injected directly into the ESI-MS interface or was collected for further LC/ESI-MS analysis.

    2.4.Rat liver microsomal incubations

    A mixture of microsomal protein(RLMs)and C-1305,dissolved in 0.1 M potassium phosphate buffer solution(adjusted to pH 7.4 with 1 M KOH),was pre-warmed for 5 min at 37°C in a shaking water bath.MgCl2and NADPH were added to the incubation mixture(to a total volume of 100μL),which was then further incubated for 60 min at 37°C.The final concentrations were as follows:2 mg/mL RLMs,0.2 mM C-1305(added from 2 mM stock solution in 0.1 M potassium phosphate buffer solution,pH 7.4),0.5 mM MgCl2,and 2 mM NADPH(added from 10 mM stock solution in 0.1 M potassium phosphate buffer solution,pH 7.4).The incubation was terminated by adding ice-cold CH3OH(1:1,V/V)to the incubation mixture for the precipitation of proteins.The sample was then vortexed and placed in ice for 10 min.After centrifugation at 10,000×g for 15 min,an aliquot of the supernatant(150μL)was then analysed directly by reversed-phase high-performance LC(RP-HPLC)with UV—Vis detection at 420 nm and/ordiode arrayand multiple wavelength detection,and monitored by MS.For each incubation,a negative control was carried out without the cofactor NADPH to prevent P450 catalysis.All assays were conducted in at least triplicate.

    Table 1 EC and ESI-MS parameters as applied in direct EC/ESI-MS experiments for determination of accurate masses of product ions.

    2.5.In silico prediction of C-1305 metabolism in liver

    The MetaSite software tool(version 5.1.1;Molecular Discovery Ltd.,Hertfordshire,UK)is a computational algorithm that predicts the most likely metabolic transformations of the compound related to phase I reactions in the liver.Moreover,it provides the structure of the metabolites,formed with a ranking derived from the site of metabolism predictions,with exact molecular weight(MW)and calculated logarithm of partition coefficient(cLogP)to help and complement the experimental elucidation of metabolite structures.

    A 2-D structure of the C-1305 was imported into the interface of MetaSite to predict metabolic soft spots for liver metabolism and structures of metabolites in the liver related to P450-mediated reactions.Only metabolites with a molecular mass higher than 150 Da and with a likelihood ranking≥50% were considered.

    2.6.EC and ESI-MS conditions

    ESI-MS detection of C-1305 products generated by EC was performed in the positive ion mode in the full scan mode(mass-tocharge ratio,m/z 100-1000).The detailed EC and ESI-MS parameters are presented in Tables 1 and 2.

    2.7.LC conditions

    Chromatographic separations of the parent drug and its transformation products were carried out using a reversed-phase 5-μm Suplex pKb-100 analytical column(C18)(Supelco Inc.,Bellefonte,PA,USA)with the following dimensions:250 mm length×4.6 mm i.d.,5-μm particle size.The column was operated at a room temperature of 25°C.The flow rate utilized was 1 mL/min.For all chromatographic separations,eluent A of the mobile phase was a solution of 0.05 M HCOONH4in water(adjusted to pH 3.4 with HCOOH).Eluent B was CH3OH.A linear gradient from 15% to 85% B in A was kept for 25 min,then followed by a linear gradient from 80% to 100% B in A for 3 min.A rate of 100% of solvent B was kept for 1.5 min before returning to initial conditions within 0.5 min.These were then kept for 10 min to allow for column equilibration.

    3.Results and discussion

    Fig.2.The representative MS spectrum of C-1305(m/z 338)before(reactor cell off;the inset at the top-left corner)and after(reactor cell on)reaction in the electrochemical cell(positive-ion mode).

    Fig.3.The representative ion mass intensity-time curves of 10μM C-1305 oxidation at a GC working electrode(positive-ion mode).The m/z ratios shown correspond to the protonated[M+H]+C-1305 and its oxidation products(see legend).The signal is dependent on the voltage used in the EC cell(two and half cell cycles are shown).Experimental conditions:potential range 0—2.5 V;scan rate 10 mV/s,continuous;T=21 °C;φ GC working electrode 8 mm.

    Any drug development process must proceed through several stages in order to produce a product that is efficacious and safe.Hence,one of the most important topics is identifying metabolic pathways of a new drug.It requires specific methods involved in the generation and characterization of various types of drug metabolites that determine its bioavailability,activity and toxicity profile in humans.In this work,three different strategies for the investigation of C-1305 oxidative metabolism were used.In the following part,the results of electrochemical simulation are discussed towards their correlation with the results from a conventional microsomal experiment and in silico model for the prediction of potential liver metabolism reactions.The C-1305 products of electrochemical and enzymatic transformations were separated and analysed by RP-HPLC with UV—Vis detection and/or diode array and multiple wavelength detection.For structure elucidation of the detected compounds,ESI-MS experiments were performed.

    3.1.Electrochemical generation of C-1305 oxidation products

    Fig.4.The representative EICs of C-1305 and its oxidation products found in an electrochemical simulation(EC)(A)and in an incubation mixture with rat liver microsomes(RLMs)(B).The m/z ratios shown correspond to the protonated[M+H]+C-1305 and its oxidation products.Peak names correspond to compounds presented in Table 2.

    The focus of this study lied in the investigation of the metabolic pathway of C-1305 using the purely instrumental EC method.No metabolism data for this compound is yet available from such a matrix-free environment.In the beginning,in order to obtain a comprehensive overview of C-1305 oxidation products,the electrochemical method was optimized by adjusting several parameters.The obtained results(Tables S1 and S2)indicated that type of electrolyte(composition,pH),working electrode material,flow rate,and potential range influenced the number and/or amounts of products formed during the electrochemical treatment.Finally,the highest electrochemical conversion of C-1305 into its expected products was achieved for 0.1% HCOOH in water/CH3OH electrolyte(50:50,V/V)with a flow rate of 20μL/min and with GC as working electrode material operated in the range of potential 0—2.5 V.GC electrode,widely used in EC,can be applied over an extended potential range when oxidation reactions occur.Thus,due to its high overpotential for oxygen,it cannot be excluded that partial electrolysis of the aqueous solvent also took place.The working electrode surface was cleaned between each oxidation step by flushing with a solution of 50% CH3OH in water.When disassembling the working electrode after an oxidation,no adsorption residue was visible.

    In order to obtain products generated electrochemically from C-1305,the parent compound was introduced into the electrochemical cell at a flow rate of 20μL/min for a designated time period.The scan mode was used for continuous synthesis of C-1305oxidation products.Fig.2 illustrates the MS spectra of C-1305 before(cell off,control measurement;Fig.2 inset)and after(cell on)reaction in the electrochemical cell.When the potentialwas not applied to the EC system,the MS spectrum shows the peak at m/z 338,that is attributed to the C-1305 ion[M+H]+.Additionally,in a solution of C-1305 without electrochemical oxidation,we also observed an intense ion at m/z 310 which might correspond to C-1305 derivative.In turn,the intensityof these peaks decreased with the increasing potential and depended on the potential value the peaks of the appropriate oxidation products appeared.The obtained results were displayed as ion mass intensity-time curves,which made it possible to perform fast screening for electrochemically generated products by showing the ion signal intensity of the relevant product(y-axis)depending on the time of the electrochemical conversion(x-axis)(Fig.3),which reflected the changes in electrochemical potential.After electrochemical reaction seven possible products of C-1305 were detected.A growth in the intensity of signals coming from oxidation products and a drop in the intensity of signal from parent ion were observed.

    Table 2 Molecular formulas of C-1305 and its products found after electrochemical oxidation on a GC working electrode in a potential range of 0—2.5 V versus Pd/H2and predicted in incubations with RLMs or by MetaSite.The structural calculations are based on exact masses,with the relative mass deviations between measured and theoretical m/z values being provided in ppm.For peak identification refer to Fig.4A.

    3.2.Separation and identification of the electrochemical products of C-1305

    For characterization of C-1305 electrochemical products in more detail regarding their polarity,a separation method based on RPHPLC has been developed.The extracted ion chromatograms(EICs)of C-1305 and its generated electrochemically potential products are shown in Fig.4A.These datawere then compared with LC/ESI-MS data obtained from liver microsomal incubation.With the application of off-line EC RP-HPLC/ESI-MS method,C-1305 and its six electrochemical products were identified.Table 2 provides a list of the products of electrochemical conversion of C-1305 characterized by different m/z values,retention time and the proposed transformations of C-1305 molecule.The deviation of measured m/z to the theoretical m/z,given as relative mass deviation(ppm)=((m/zexperimental— m/zcalculated)/m/zcalculated)× 106,was for all products≤4 ppm.

    P450 is the most important enzyme system that catalyses a variety of oxidative transformations of xenobiotics and endogenous compounds(phase I reactions).P450s activate molecular oxygen and transfer a single oxygen atom to a substrate while the other oxygen atom ends up inwater molecule[28].These oxygen transfer reactions include the following:aromatic/aliphatic hydroxylations,epoxidation,and heteroatom oxidations.Additionally,P450s can also promote heteroatom-dealkylations,dehalogenation,and dehydrogenation[29].As can be seen from the literature,the electrochemical system is able to successfully simulate the majority ofP450-catalyzed mixed-function oxidation reactions[9-13].During the direct electrochemical oxidations,the electrode behaves as an oxidant which can be adjusted by the applied potential in order to perform charge transfer reactions[8].The proposed electrochemical conversion of the investigated triazoloacridinone,summarised in the reaction scheme shown in Fig.5,includes mainly N-dealkylation(P5 at m/z 310 and P7 at m/z 324)and single O gain via hydroxylation(P6 at m/z 354).

    N-Dealkylation is a typical biotransformation pathway in P450-catalyzed oxidation reactions that usually inactivates xenobiotics,including secondary and tertiary amines found in many drugs,and facilitates their elimination[30].It can be readily simulated by the EC system.The electrochemical mechanism resembles to a large extent a single-electron transfer mechanism mediated by P450.It consists of an initial electron abstraction,followed by deprotonation and the subsequent abstraction of another electron that results in an iminium intermediate.Following hydrolysis,the intermediate decomposes into the final N-dealkylated product[31].In the case of C-1305,N-demethylated and N-bis-demethylated products were obtained in the electrochemical cell.The chromatographic peak of the latter(P5,m/z 310)was eluted with the peak corresponding to that of substrate(C-1305,m/z 338),which indicates that these compounds may have similar polarity.In turn,product P7,which had been previously found in EC analysis,was not detected after chromatographic separation,likely due to the slightly inferior limit of detection and/or its limited chemical stability.It should be noted that reactive iminium intermediate formed during N-dealkylation process may be associated with the mechanism-based inactivation of P450 enzymes reported for C-1305[32].

    Another remarkable oxidative capability of P450 is hydroxylation[30].This process can be easily induced electrochemically at positive potentials[12].The presence of an ion at m/z 354 corresponding to an increase in mass by 16 Da confirmed that successful electrochemical hydroxylation of C-1305 molecule took place.Further,two different isomers of this product(peaks P6a and P6b in Fig.4A)were found after chromatographic separation of the postelectrochemical mixture.Both are characterized by retention time slightly higher than that of C-1305.

    Fig.5.Overview about all relevant products generated after electrochemical oxidation of C-1305 with proposed structure.Observed m/z value is shown.Numbers correspond to compounds presented in Table 2.

    Due to the presence of an additional hydroxyl group there may be more polar substances than substrate.Therefore,it is suspected that O gain might have occurred at various positions of the C-1305 molecule by aliphatic and aromatic hydroxylation.The electrochemical oxidation potentials for hydroxylation of aliphatic C-H bonds are generally very high[12],but the application of a potential of above 2 V might have enabled oxidation reaction in this case.We suppose that the product of such hydroxylation was obtained with lower yield.For product P6,a putative structure of diol(quinol)is also considered.Electrochemical oxidation of aromatic ring is initiated by electron abstraction.This electron transfer requires a lower positive potential than for aliphatic hydrocarbons due to resonance stabilization of the aromatic radical cation.In order to be oxidized electrochemically within the potential limits of water,the aromatic ring must be activated by an electron-donating group.This requirement was made by 8-hydroxyl group present in the triazoloacridinone moiety.It can be expected that analogously as for compound C-1311[33],the activation of the triazoloacridinone ring system at position ortho to this functional group will be the most preferential.Furthermore,product P3 observed at m/z 352,indicating the mass gain of two hydrogens,was very likely generated from product P6 by dehydrogenation.Thus,the above hydroxylations might have proceed immediately to respective less polar carbonyl or quinone derivatives.

    In our study we also observed signals at m/z 366(P1),368(P4),and 382(P2)which could be supposed to be the potential products resulting from the combination of different oxidative pathways(e.g.,hydroxylation,dehydrogenation)of the investigated antitumor agent.In the elution profile their peaks appeared before the C-1305 peak.When considering the structure of the product P4,it can be assumed that the addition of a single methoxy group achieved by EC could have taken place here.We found that its formation was highly dependent on the electrolyte used.So,it is suggested that this product may result from the interaction between the substrate to be oxidized and CH3OH being a component of the electrolyte solution.

    3.3.Metabolism study in RLMs

    In order to obtain the metabolite profiling of drug compounds,the incubationwith livermicrosomes is a conventional,well known method[4,5].Therefore,this approach was performed on C-1305 to compare the metabolic products with those detected after electrochemical conversion of the studied compound.The representativeEICsofC-1305 productsobtained from microsomal incubations are shown in Fig.4B.In the sample taken from rat liver microsomal incubation with C-1305,compared to the purely instrumental non-enzymatic approach,the presence of only three metabolites of the parent drug:C-1305+O-2H(peak P*3a at m/z=352),C-1305-2CH2(peak P*5 at m/z=310),and C-1305+O(peak P*6a at m/z=354)was observed.No additional metabolites were detected in the LC/ESI-MS analysis of the microsomal sample.Although liver microsomes contain the physiologically relevant combinations of drug-metabolizing enzymes,the expression of the different P450 isoforms may be very diverse in each organism.This may explain why some products of C-1305 generated electrochemically were not produced in liver microsomes.

    The peaks of products P*3b,P*5,and P*6b were characterized byretention time identical tothat of products P3a,P5,and P6b from EC analysis.The first two were of low intensities,whereas the amount of this latter increased with longer incubation time.Based on the experimental mass ion data,the same observations as for the previously identified products of electrochemical conversion were made for the P3a and P5 structures.Thus,N-dealkylation of alkylamine,hydroxylation and dehydrogenation seem to be the main pathways in both EC and enzymatic oxidations of C-1305.The generation of metabolite P*b presumably involves the insertion of a single oxygen atom that,as it was proposed above,may take place under the formation of a hydroxyl group.However,the antitumor agent selected for this study was a good substrate for human recombinant flavin-containing monooxygenases(FMOs),FMO1 and FMO3[34].FMOs can oxidize a wide array of heteroatoms,particularly soft nucleophiles,such as amines,in the presence of an oxygen,an NADPH cofactor,and an FAD prosthetic group[35].Assuming that the active enzymes of RLMs are FMOs,it was suggested that peak P*6b at m/z 354 may be an amine-N-oxide.Nitrogen oxidation is also possible in the EC system,but unlike formation of S-or P-oxides,N-oxide formation can only be detected in low yields[12,36].

    3.4.In silico prediction of C-1305 liver metabolism

    MetaSite software allowed to predict the main sites of liver metabolism for C-1305 and its main metabolic reactions(Fig.6).The most probable sites of P450 attack in the C-1305 molecule proved the dialkylaminoalkylamino moiety and the carbon atom at position 9 of the triazoloacridinone ring.In view of this,five top ranked phase I metabolite structures were reported.

    One of the major metabolites predicted was that derived by the cleavage of the N-C bond leading to N-demethylation(M1;100% of likelihood).The identical product ion P7 at m/z 324 was also identified by EC/ESI-MS analysis.In turn,further degradation of the side chain may result in a complete N-dealkylated metabolite M4(50% of likelihood)which,in contrast to M1,was not generated from electrochemical oxidation. Another metabolic reaction predicted for C-1305,common to the P450 catalytic activity,was also dehydrogenation in the propylamino chain(M5;50% of likelihood).The N-dealkylation and dehydrogenation reactions in the dialkylaminoalkylamino side chain,analogous to those described above,were also predicted by in silico analysis for imidazoacridinone derivative,C-1311[37].Such substituent is flexible and,after metabolic activation,it may be able to induce interstrand covalent crosslinks in DNA of tumor cells playing a key role in biological(antitumor)activity of these compounds[24,38,39].Such mechanism mediating the cytotoxicity was also demonstrated in the case of the clinically used antitumor drugs,mitoxantrone and ametantrone[22].Further,MetaSite indicated the possibility of the formation of two isomeric hydroxylated metabolites(MW=353.15)with various ranking that can be linked with product ion P6 at m/z 354 found after electrochemical simulation.The most probable hydroxylation site in the C-1305 molecule appears to be the ortho position to the hydroxyl group in the aromatic ring,giving metabolite M2(66% of likelihood).The resulting diol(quinol)and its putative quinone derivative,which may correspond to product ion P3 obtained at m/z 352 in the EC analysis,are the preffered structures for the substitution by nucleophiles existing in a living organism(i.e.,reduced glutathione,thiols of proteins or purine and pyrimidine bases of DNA).The formation of stable adducts may alter biological functions of cellular biomolecules,which ultimately leads to a toxic response[5].Thus,both diol and quinone metabolites may be responsible for cytotoxic and antitumor actions of triazoloacridinones.Moreover,they may be related to the observed inactivation of P450 1A2 and 3A4 by C-1305 by the covalent apoprotein modification of P450 isoenzymes[32].In addition to aromatic hydroxylation,aliphatic hydroxylation is also possible,as in the case of metabolite M3(50% of likelihood).

    The ability of the drugs to penetrate various biological membranes,tissues and barriers is determined by the hydrophilic or hydrophobic properties of molecules,for which the partition coefficient,P,of a molecule is used[40,41].In the case of all MetaSite predicted metabolites of C-1305 formed by liver P450s,the calculated cLogP values were in the range of 0.90—1.88.It means that they are rather more hydrophilic substances than the parent compound(cLogP=2.06),so they may be preferentially localized in hydrophilic compartments,such as blood serum,and easily excreted from the body.

    To sum up,in silico study has demonstrated the possibility to theoretically predict the main metabolites of antitumor active triazoloacridinone C-1305.In silico approach proposed here appears very promising,allowing the identification of main metabolites of C-1305,which are extensively metabolized as demonstrated by electrochemical(non-enzymatic)and microsomal(enzymatic)experiments performed.

    4.Conclusions

    Fig.6.Liver metabolism prediction by MetaSite for C-1305.Ranking of the metabolites(A)and plot of the most probable sites of C-1305 metabolism and metabolites predicted(≥50% of likelihood)(B).The functional groups in compound that most likely will be metabolized by liver P450s are marked:the darker the color and the greater the circle of marked functional group—the higher the probability of metabolism to occur.MW and cLogP are predicted values according to MetaSite software tool(version 5.1.1;Molecular Discovery Ltd.,Hertfordshire,UK).MW,molecular weight;cLogP,the logarithm of partition coefficient.

    This work demonstrated three different strategies for the investigation of oxidative metabolism of antitumor triazoloacridinone C-1305.In particular,we presented here the first application of electrochemical approach for generation of potential oxidation products of the studied compound within a matrix-free environment.We found that a dialkylaminoalkylamino moiety was most susceptible to electrochemical oxidation site in the C-1305 molecule.The proposed N-dealkylated,hydroxylated,and dehydrogenated products turned out to be typical for the processes catalyzed by P450 is oenzymes.Next,the results achieved from EC/ESI-MS were correlated with those from liver microsomal experiments.On the one hand,each method was found to provide unique features not identified with those of others.However,despite some differences in the mechanisms of electrochemical and enzymatic oxidation,the products detected after electrochemical and enzymatic conversion of C-1305 were generally well comparable.This underlines the high potential of EC approach as a fast screening tool in the prediction of metabolic transformations in drug development.In the final step,MetaSite in silico analysis reported five top ranked phase I metabolite structures for C-1305 liver metabolism that helped in the interpretation of the results provided by electrochemical and enzymatic approaches.The value of this strategy consists in the possibility of gaining very useful information on the most probable metabolites to be screened in all those cases where the exact metabolism pathways of a compound are unknown.According to the above findings,the electrochemical and in silico methods,which were shown to predict metabolism patterns,can be used as simple and fast alternatives for metabolic(enzymatic)assays,affording time and cost efficiency.

    Conflicts of interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    This work was supported by the National Science Center(Poland)(2012/07/D/NZ7/03395).The authors wish to thank Professor Agata Kot-Wasik(Department of Analytical Chemistry,Faculty of Chemistry,Gdańsk University of Technology,Poland),for help in the field of MS analyses.

    Appendix A.Supplementary data

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

    成年女人永久免费观看视频| 村上凉子中文字幕在线| 日韩人妻高清精品专区| 亚洲av美国av| 欧美成人a在线观看| x7x7x7水蜜桃| 波多野结衣巨乳人妻| 999久久久精品免费观看国产| 久久精品国产自在天天线| 亚洲性夜色夜夜综合| 久久久国产成人免费| 国产av一区在线观看免费| 欧美精品啪啪一区二区三区| 色噜噜av男人的天堂激情| 免费大片18禁| 欧美黄色淫秽网站| 91在线观看av| 好男人电影高清在线观看| 99热只有精品国产| 男人和女人高潮做爰伦理| 亚洲第一电影网av| av天堂中文字幕网| 国产蜜桃级精品一区二区三区| 日日干狠狠操夜夜爽| 桃色一区二区三区在线观看| 久久精品国产亚洲av涩爱 | 亚洲精品成人久久久久久| 亚洲成人久久性| 精品一区二区三区人妻视频| 三级国产精品欧美在线观看| 国产野战对白在线观看| 国产主播在线观看一区二区| 亚洲成人久久性| 三级男女做爰猛烈吃奶摸视频| 99热精品在线国产| 国产成人a区在线观看| 99久久成人亚洲精品观看| 麻豆成人午夜福利视频| www.999成人在线观看| 变态另类成人亚洲欧美熟女| 老熟妇乱子伦视频在线观看| 国产精品99久久久久久久久| 国产精品乱码一区二三区的特点| 国产成+人综合+亚洲专区| 国产伦一二天堂av在线观看| 国产在线精品亚洲第一网站| 少妇的丰满在线观看| 亚洲欧美日韩高清专用| 日本a在线网址| 丰满乱子伦码专区| 少妇的逼好多水| 久久中文看片网| 精品乱码久久久久久99久播| 精品国产三级普通话版| 成人三级黄色视频| 男女午夜视频在线观看| 精华霜和精华液先用哪个| 97人妻精品一区二区三区麻豆| 久久久久亚洲av毛片大全| 国产三级中文精品| 91久久精品电影网| 精品一区二区三区视频在线观看免费| 99视频精品全部免费 在线| 国产精品一区二区三区四区免费观看 | 国产高清有码在线观看视频| 精品久久久久久久末码| 啪啪无遮挡十八禁网站| av女优亚洲男人天堂| 亚洲最大成人手机在线| 十八禁人妻一区二区| 午夜福利在线观看吧| 国内毛片毛片毛片毛片毛片| 长腿黑丝高跟| 国产精品av视频在线免费观看| 欧美日本亚洲视频在线播放| 日韩欧美国产一区二区入口| 亚洲电影在线观看av| 亚洲国产精品合色在线| 欧美高清成人免费视频www| 午夜精品一区二区三区免费看| 亚洲第一欧美日韩一区二区三区| 亚洲国产色片| 亚洲 国产 在线| 精品不卡国产一区二区三区| 成人高潮视频无遮挡免费网站| 91九色精品人成在线观看| 久久久精品大字幕| 精品久久久久久久久久久久久| 最新在线观看一区二区三区| 久久午夜亚洲精品久久| 国产欧美日韩一区二区精品| 亚洲国产色片| 欧美成狂野欧美在线观看| 中国美女看黄片| 嫁个100分男人电影在线观看| 男女视频在线观看网站免费| 最近最新中文字幕大全免费视频| 99国产综合亚洲精品| 淫秽高清视频在线观看| 免费高清视频大片| 九色成人免费人妻av| 国产伦精品一区二区三区视频9 | 在线观看一区二区三区| 日韩成人在线观看一区二区三区| 婷婷精品国产亚洲av在线| 搡老岳熟女国产| 村上凉子中文字幕在线| 亚洲真实伦在线观看| 日韩人妻高清精品专区| 欧美+亚洲+日韩+国产| 日本免费一区二区三区高清不卡| 国产三级中文精品| 亚洲国产高清在线一区二区三| 波野结衣二区三区在线 | 搞女人的毛片| 欧美绝顶高潮抽搐喷水| 欧美乱色亚洲激情| 久久久精品大字幕| 两个人看的免费小视频| 国产 一区 欧美 日韩| 尤物成人国产欧美一区二区三区| 法律面前人人平等表现在哪些方面| 欧美一级a爱片免费观看看| 色老头精品视频在线观看| 国产精品99久久久久久久久| 亚洲国产日韩欧美精品在线观看 | 亚洲无线在线观看| 免费观看人在逋| 制服丝袜大香蕉在线| 国产精品自产拍在线观看55亚洲| ponron亚洲| 国产伦一二天堂av在线观看| 国产成人啪精品午夜网站| 欧美极品一区二区三区四区| 亚洲av免费高清在线观看| 久久久久久国产a免费观看| av国产免费在线观看| 国产精品av视频在线免费观看| 日韩欧美三级三区| 亚洲内射少妇av| 一个人免费在线观看电影| 久久精品91无色码中文字幕| 亚洲最大成人手机在线| 在线看三级毛片| 九色国产91popny在线| 两个人看的免费小视频| 国产精品一区二区免费欧美| 黑人欧美特级aaaaaa片| xxx96com| 99精品在免费线老司机午夜| 99热精品在线国产| 一级毛片高清免费大全| 国产亚洲精品综合一区在线观看| 成人高潮视频无遮挡免费网站| 国语自产精品视频在线第100页| 欧美日韩综合久久久久久 | 国产成人系列免费观看| 午夜福利欧美成人| 免费看a级黄色片| 少妇人妻一区二区三区视频| 国产成人系列免费观看| 草草在线视频免费看| 九色成人免费人妻av| 久久久久久久午夜电影| 99在线人妻在线中文字幕| 90打野战视频偷拍视频| 午夜福利免费观看在线| 亚洲七黄色美女视频| 亚洲精品影视一区二区三区av| 18禁裸乳无遮挡免费网站照片| 日韩亚洲欧美综合| 亚洲人成电影免费在线| 在线看三级毛片| 国产精品永久免费网站| 精品久久久久久久人妻蜜臀av| 无限看片的www在线观看| 夜夜看夜夜爽夜夜摸| 又粗又爽又猛毛片免费看| 少妇人妻一区二区三区视频| 久久香蕉国产精品| 性色avwww在线观看| 91麻豆av在线| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影 | 精品熟女少妇八av免费久了| av黄色大香蕉| 亚洲精品一卡2卡三卡4卡5卡| 国产av在哪里看| 琪琪午夜伦伦电影理论片6080| 国产单亲对白刺激| 悠悠久久av| 日本免费a在线| 九色国产91popny在线| 国产精品综合久久久久久久免费| 亚洲18禁久久av| 国产伦人伦偷精品视频| 男女视频在线观看网站免费| 两人在一起打扑克的视频| 亚洲人成网站高清观看| 国产免费男女视频| 身体一侧抽搐| 他把我摸到了高潮在线观看| 亚洲精品国产精品久久久不卡| 国产一区二区三区在线臀色熟女| 久久天躁狠狠躁夜夜2o2o| 中亚洲国语对白在线视频| 国产亚洲av嫩草精品影院| 精品无人区乱码1区二区| 男女那种视频在线观看| 亚洲人成网站高清观看| av专区在线播放| 中文资源天堂在线| 欧美日韩瑟瑟在线播放| 国内精品美女久久久久久| 午夜a级毛片| 美女高潮喷水抽搐中文字幕| 午夜老司机福利剧场| 两人在一起打扑克的视频| 亚洲第一欧美日韩一区二区三区| 两个人看的免费小视频| 特级一级黄色大片| 精品电影一区二区在线| 国产不卡一卡二| 长腿黑丝高跟| 一个人免费在线观看电影| 国产欧美日韩一区二区精品| 亚洲人成网站在线播| 亚洲国产日韩欧美精品在线观看 | 亚洲18禁久久av| 69av精品久久久久久| 亚洲人与动物交配视频| 成人一区二区视频在线观看| 亚洲专区国产一区二区| 精品电影一区二区在线| 亚洲美女视频黄频| 99在线人妻在线中文字幕| 啪啪无遮挡十八禁网站| 别揉我奶头~嗯~啊~动态视频| 夜夜夜夜夜久久久久| 日本三级黄在线观看| 淫妇啪啪啪对白视频| 久久精品国产综合久久久| 久久久久久人人人人人| 人人妻,人人澡人人爽秒播| 久久中文看片网| 免费观看人在逋| 美女 人体艺术 gogo| 最近最新中文字幕大全电影3| 无限看片的www在线观看| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久亚洲av鲁大| 免费av毛片视频| 亚洲狠狠婷婷综合久久图片| 日韩高清综合在线| 女人被狂操c到高潮| 日韩大尺度精品在线看网址| 麻豆成人av在线观看| 久久久久亚洲av毛片大全| 男女床上黄色一级片免费看| 超碰av人人做人人爽久久 | 成年人黄色毛片网站| 国产久久久一区二区三区| 人人妻人人看人人澡| 国产免费一级a男人的天堂| 久久久成人免费电影| 亚洲欧美一区二区三区黑人| 97超级碰碰碰精品色视频在线观看| 最近最新中文字幕大全电影3| www日本在线高清视频| 亚洲精品在线观看二区| 中文字幕人妻丝袜一区二区| 日韩高清综合在线| 久久精品国产99精品国产亚洲性色| 国产激情欧美一区二区| 99久久九九国产精品国产免费| 此物有八面人人有两片| 免费av毛片视频| 制服丝袜大香蕉在线| 高清毛片免费观看视频网站| 欧美色欧美亚洲另类二区| 观看免费一级毛片| 在线十欧美十亚洲十日本专区| 精品国产超薄肉色丝袜足j| 国产一区在线观看成人免费| 亚洲av成人av| 男女视频在线观看网站免费| 黄色片一级片一级黄色片| 日日摸夜夜添夜夜添小说| 成年女人看的毛片在线观看| 国产精品一区二区三区四区久久| 国产伦人伦偷精品视频| 最近在线观看免费完整版| 欧美另类亚洲清纯唯美| 欧美色视频一区免费| 欧美极品一区二区三区四区| 成人精品一区二区免费| 亚洲精品在线美女| 99久久精品热视频| а√天堂www在线а√下载| 久久久久性生活片| 男人舔女人下体高潮全视频| 国产一区二区在线av高清观看| 中亚洲国语对白在线视频| 欧美大码av| 内射极品少妇av片p| 亚洲人成网站在线播放欧美日韩| 久久国产乱子伦精品免费另类| 女生性感内裤真人,穿戴方法视频| 国模一区二区三区四区视频| 色尼玛亚洲综合影院| 婷婷亚洲欧美| 一级a爱片免费观看的视频| 免费看美女性在线毛片视频| 久久草成人影院| 欧美日本视频| 男人舔奶头视频| 哪里可以看免费的av片| 99riav亚洲国产免费| 最新美女视频免费是黄的| 老熟妇仑乱视频hdxx| 宅男免费午夜| 欧美xxxx黑人xx丫x性爽| 桃色一区二区三区在线观看| 免费一级毛片在线播放高清视频| 精品国产超薄肉色丝袜足j| 九九久久精品国产亚洲av麻豆| 丝袜美腿在线中文| 久久久国产成人精品二区| 小蜜桃在线观看免费完整版高清| 亚洲精品在线美女| 免费av毛片视频| 免费搜索国产男女视频| 欧美成人免费av一区二区三区| 欧美日韩国产亚洲二区| 日韩精品青青久久久久久| 亚洲熟妇熟女久久| 国产欧美日韩一区二区精品| 小蜜桃在线观看免费完整版高清| 99久久久亚洲精品蜜臀av| 亚洲av日韩精品久久久久久密| 成人无遮挡网站| 色综合欧美亚洲国产小说| 搞女人的毛片| 成人三级黄色视频| 日韩欧美免费精品| 美女 人体艺术 gogo| 窝窝影院91人妻| 每晚都被弄得嗷嗷叫到高潮| 亚洲av美国av| 男插女下体视频免费在线播放| 午夜久久久久精精品| 欧美最新免费一区二区三区 | 国产在视频线在精品| 国产亚洲精品久久久久久毛片| 日韩成人在线观看一区二区三区| 99国产综合亚洲精品| 久久九九热精品免费| 亚洲男人的天堂狠狠| 成年版毛片免费区| 欧美av亚洲av综合av国产av| 麻豆国产97在线/欧美| 精品久久久久久久末码| 国产在视频线在精品| 网址你懂的国产日韩在线| 一a级毛片在线观看| 天天添夜夜摸| 欧美黄色淫秽网站| 久久久久国产精品人妻aⅴ院| 99久久精品热视频| 精品久久久久久久毛片微露脸| 中国美女看黄片| 99久久综合精品五月天人人| 国产激情欧美一区二区| 99国产极品粉嫩在线观看| 久久6这里有精品| 精品午夜福利视频在线观看一区| 欧美不卡视频在线免费观看| avwww免费| 国产精品久久久久久人妻精品电影| 校园春色视频在线观看| 欧美日韩精品网址| 国产精品久久久久久精品电影| 欧美成人一区二区免费高清观看| www.熟女人妻精品国产| 精品不卡国产一区二区三区| 宅男免费午夜| 男人和女人高潮做爰伦理| 精品国内亚洲2022精品成人| 老熟妇仑乱视频hdxx| 欧美一区二区国产精品久久精品| 国产成人aa在线观看| 午夜视频国产福利| 亚洲真实伦在线观看| 18禁在线播放成人免费| 国产69精品久久久久777片| 成年版毛片免费区| 亚洲av免费在线观看| 久9热在线精品视频| 国产亚洲精品一区二区www| 国产高清激情床上av| 久久久久久久精品吃奶| 在线播放国产精品三级| 手机成人av网站| 3wmmmm亚洲av在线观看| 天天躁日日操中文字幕| 人妻丰满熟妇av一区二区三区| 一区二区三区高清视频在线| 噜噜噜噜噜久久久久久91| 99久久综合精品五月天人人| 国产极品精品免费视频能看的| 九色成人免费人妻av| 黄色片一级片一级黄色片| 成人永久免费在线观看视频| 久久精品国产清高在天天线| 国产淫片久久久久久久久 | 女人十人毛片免费观看3o分钟| 人妻丰满熟妇av一区二区三区| 人人妻人人澡欧美一区二区| 在线观看美女被高潮喷水网站 | 级片在线观看| 久久精品亚洲精品国产色婷小说| 草草在线视频免费看| 日韩有码中文字幕| 熟妇人妻久久中文字幕3abv| 色综合亚洲欧美另类图片| 桃红色精品国产亚洲av| 久久久久免费精品人妻一区二区| 黄色女人牲交| 51国产日韩欧美| 国内精品一区二区在线观看| 午夜两性在线视频| 男女下面进入的视频免费午夜| 午夜老司机福利剧场| 真人做人爱边吃奶动态| 亚洲国产精品999在线| 久久天躁狠狠躁夜夜2o2o| 国内精品久久久久精免费| av天堂在线播放| 最近在线观看免费完整版| 夜夜看夜夜爽夜夜摸| 久久精品91无色码中文字幕| 色老头精品视频在线观看| 成人特级av手机在线观看| 国内久久婷婷六月综合欲色啪| 国产精品影院久久| 午夜福利视频1000在线观看| 9191精品国产免费久久| 国产主播在线观看一区二区| 欧美日韩福利视频一区二区| 日本熟妇午夜| 亚洲精品色激情综合| 国内精品久久久久久久电影| 少妇高潮的动态图| 日本精品一区二区三区蜜桃| e午夜精品久久久久久久| 熟女电影av网| 免费大片18禁| 日日夜夜操网爽| 午夜a级毛片| 无人区码免费观看不卡| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| 日韩大尺度精品在线看网址| 人人妻人人澡欧美一区二区| 久久久色成人| 国产黄片美女视频| 亚洲av一区综合| 黑人欧美特级aaaaaa片| 国产不卡一卡二| 成年人黄色毛片网站| 一个人看视频在线观看www免费 | 禁无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 90打野战视频偷拍视频| 亚洲天堂国产精品一区在线| 亚洲自拍偷在线| 国产欧美日韩一区二区精品| 又黄又爽又免费观看的视频| 国产三级黄色录像| 在线播放国产精品三级| 国产野战对白在线观看| 亚洲最大成人手机在线| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 日本三级黄在线观看| 日本与韩国留学比较| 两人在一起打扑克的视频| 国产三级中文精品| 桃红色精品国产亚洲av| 日韩国内少妇激情av| 哪里可以看免费的av片| 欧美成人免费av一区二区三区| 麻豆久久精品国产亚洲av| 丝袜美腿在线中文| a级一级毛片免费在线观看| 久久精品国产亚洲av香蕉五月| 午夜激情欧美在线| 波野结衣二区三区在线 | 午夜福利在线观看吧| 国产伦在线观看视频一区| 国产成年人精品一区二区| 婷婷丁香在线五月| 成人国产一区最新在线观看| 一本久久中文字幕| 日韩精品中文字幕看吧| 听说在线观看完整版免费高清| 国产成人av激情在线播放| 国产一区在线观看成人免费| 国产中年淑女户外野战色| 亚洲国产高清在线一区二区三| 一夜夜www| 99久久99久久久精品蜜桃| 亚洲精品国产精品久久久不卡| 最近最新中文字幕大全电影3| 男人舔女人下体高潮全视频| 欧美黄色淫秽网站| 99久久久亚洲精品蜜臀av| 母亲3免费完整高清在线观看| 18+在线观看网站| 男女床上黄色一级片免费看| 国产aⅴ精品一区二区三区波| 国产真人三级小视频在线观看| 99热6这里只有精品| 悠悠久久av| 午夜福利在线观看吧| 日韩精品中文字幕看吧| 制服丝袜大香蕉在线| 欧美高清成人免费视频www| 免费高清视频大片| 国产aⅴ精品一区二区三区波| 中文字幕人成人乱码亚洲影| e午夜精品久久久久久久| 悠悠久久av| 美女cb高潮喷水在线观看| 黄色日韩在线| 久久6这里有精品| 国产精品99久久久久久久久| 嫩草影院精品99| 淫妇啪啪啪对白视频| 国产成年人精品一区二区| xxx96com| 免费电影在线观看免费观看| 国产成人啪精品午夜网站| 欧美在线黄色| 免费在线观看成人毛片| 日本免费一区二区三区高清不卡| 精品一区二区三区视频在线 | 国产一区二区激情短视频| 欧美一级a爱片免费观看看| 宅男免费午夜| 国产久久久一区二区三区| 3wmmmm亚洲av在线观看| 少妇的丰满在线观看| 此物有八面人人有两片| 成人性生交大片免费视频hd| 国产免费一级a男人的天堂| 久久性视频一级片| 在线天堂最新版资源| 久久天躁狠狠躁夜夜2o2o| 99久久九九国产精品国产免费| 一夜夜www| 婷婷亚洲欧美| 亚洲精华国产精华精| 国产精品国产高清国产av| 搡老妇女老女人老熟妇| 欧美乱码精品一区二区三区| 麻豆国产97在线/欧美| 高潮久久久久久久久久久不卡| 国产精品久久久人人做人人爽| 深夜精品福利| 国产在线精品亚洲第一网站| 黄片小视频在线播放| 亚洲av电影不卡..在线观看| 美女大奶头视频| 欧美黄色淫秽网站| 免费看a级黄色片| 国产视频内射| 午夜福利在线在线| 亚洲人成伊人成综合网2020| 老熟妇仑乱视频hdxx| 最近在线观看免费完整版| 人妻夜夜爽99麻豆av| 禁无遮挡网站| 国产精品久久久久久久电影 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲片人在线观看| av黄色大香蕉| 美女高潮喷水抽搐中文字幕| 国产又黄又爽又无遮挡在线| 国产精品,欧美在线| 国产成人a区在线观看| 亚洲精品粉嫩美女一区| 人人妻人人看人人澡| 国产一区二区三区在线臀色熟女| 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 黄色日韩在线| 俺也久久电影网| 成人特级黄色片久久久久久久| 亚洲最大成人手机在线| 悠悠久久av| 亚洲五月婷婷丁香| 女同久久另类99精品国产91| 一边摸一边抽搐一进一小说| 夜夜看夜夜爽夜夜摸| av女优亚洲男人天堂| 中文字幕av成人在线电影| 在线看三级毛片| 日韩欧美精品v在线| 老司机午夜福利在线观看视频| 国产不卡一卡二| 高清在线国产一区| 国产av不卡久久| 免费av观看视频| 在线观看美女被高潮喷水网站 |