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

    Licorice-saponin A3 is a broad-spectrum inhibitor for COVID-19 by targeting viral spike and anti-inflammation

    2024-03-21 05:51:14YngYiWenzheLiKengLiuHengXueRongYuMengZhngYngOujieBoXinyunLiJingjingFnYuxiHungJingWngXiomengShiJunhuLiHongpingWeiKunhuiXingLinjieLiRongZhngXinZhoXueQioHngYngMinYe
    Journal of Pharmaceutical Analysis 2024年1期

    Yng Yi , Wenzhe Li , Keng Liu , Heng Xue , Rong Yu , Meng Zhng ,Yng-Oujie Bo ,Xinyun Li ,Jingjing Fn ,Yuxi Hung ,Jing Wng ,Xiomeng Shi ,Junhu Li ,Hongping Wei ,Kunhui Xing ,Linjie Li ,Rong Zhng ,Xin Zho ,Xue Qio , Hng Yng ,g,**, Min Ye ,b,*

    a State Key Laboratory of Natural and Biomimetic Drugs, School of Pharmaceutical Sciences, Peking University, Beijing,100191, China

    b Yunnan Baiyao International Medical Research Center, Peking University, Beijing,100191, China

    c CAS Key Laboratory of Pathogen Microbiology and Immunology, Institute of Microbiology, Chinese Academy of Sciences, Beijing,100101, China

    d CAS Key Laboratory of Special Pathogens and Biosafety, Center for Biosafety Mega-Science, Wuhan Institute of Virology, Chinese Academy of Sciences,Wuhan, 430071, China

    e University of Chinese Academy of Sciences, Beijing,100049, China

    f Department of Microbiology and Infectious Disease Center, School of Basic Medical Sciences, Peking University Health Science Center, Peking University,Beijing,100191, China

    g Hubei Jiangxia Laboratory, Wuhan, 430000, China

    h Shenzhen Children's Hospital, Shenzhen, Guangdong, 518036, China

    Keywords:COVID-19 Inflammation Licorice-saponin A3 Receptor-binding domain (RBD)SARS-CoV-2

    A B S T R A C T Currently, human health due to corona virus disease 2019 (COVID-19) pandemic has been seriously threatened.The coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike (S)protein plays a crucial role in virus transmission and several S-based therapeutic approaches have been approved for the treatment of COVID-19.However, the efficacy is compromised by the SARS-CoV-2 evolvement and mutation.Here we report the SARS-CoV-2 S protein receptor-binding domain (RBD)inhibitor licorice-saponin A3(A3)could widely inhibit RBD of SARS-CoV-2 variants,including Beta,Delta,and Omicron BA.1, XBB and BQ1.1.Furthermore, A3 could potently inhibit SARS-CoV-2 Omicron virus in Vero E6 cells, with EC50 of 1.016 μM.The mechanism was related to binding with Y453 of RBD determined by hydrogen-deuterium exchange mass spectrometry (HDX-MS) analysis combined with quantum mechanics/molecular mechanics (QM/MM) simulations.Interestingly, phosphoproteomics analysis and multi fluorescent immunohistochemistry (mIHC) respectively indicated that A3 also inhibits host inflammation by directly modulating the JNK and p38 mitogen-activated protein kinase (MAPK) pathways and rebalancing the corresponding immune dysregulation.This work supports A3 as a promising broad-spectrum small molecule drug candidate for COVID-19.

    1.Introduction

    Currently, 6.87 million deaths of corona virus disease 2019(COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been reported till March 2023 [1,2].While many vaccines have been developed for clinical use [3,4],new waves of infections due to SARS-CoV-2 variants have been challenged.These variants may ablate antibody binding and neutralization, resulting in diminished efficacy of antibodies derived from natural infections or vaccinations [5].

    The spike (S) protein receptor-binding domain (RBD) of SARSCoV-2 is a significant target to develop antiviral therapies [6,7].Given the crucial role of RBD(aa 319-541)in recognizing the host cell receptor angiotensin converting enzyme 2 (ACE2), mutations may remarkably affect its binding to ACE2[8], thus decreasing the protective ability of current vaccines and create a huge challenge in inhibiting diverse SARS-CoV-2 variants[9].The first variant with an S mutation of D614G appeared in Germany and showed enhanced infectivity [10].A number of SARS-CoV-2 variants, identified as variants of concern (VOCs), have been reported afterwards.First,Alpha (B.1.1.7) reported in September 2020 in Southeastern England spread quickly in 160 countries as of June 2021 [11].In late 2020, Beta (1.351) and Gamma (P.1) both containing three mutations at RBD emerged respectively.Beta(1.351)has since spread to 113 countries with great concern about its immune escape[12,13].After Delta (B.1.617.2) variant [14], a newly recognized variant Omicron (B.1.1.529), containing 26-32 aa mutations in the S protein,is the predominant variant globally with several sub-variants merging[15].Compared to Delta,Omicron is more destructive due to the fast spread rate, immune escape potential, and high transmissibility.

    Although antiviral therapies targeting S protein have been developed,various studies showed that neutralization efficiency of antibodies against SARS-CoV-2 VOCs was less effective [16-18].Beta and Gamma variants with mutations K417 N/T, E484K, and N501Y at RBD decreased the neutralization capacity of antibodies activated by infection with vaccination, and showed resistance to neutralization by N-terminal domain-specific monoclonal antibodies [17,19].The Delta variant with L452R mutation causes structural changes of RBD and promotes the interaction with ACE2,and has over 60% more transmissibility than Alpha.Recently, Omicron and its sub-variants carrying over 15 mutations in RBD are spreading rapidly in countries with high levels of population immunity, which demonstrates the limited antiviral potency of vaccination[20].Thus,it will be of huge benefit to develop effective and broad-spectrum RBD inhibitors to suppress SARS-CoV-2 variants infection.

    Thus far,several attractive therapies have been reported to block SARS-CoV-2 viral entry.For instance, the peptidomimetic compound N-0385, which abrogates the proteolytic activity of TMPRSS2,could inhibit the entry of SARS-CoV-2 variants into Calu-3 cells [21].A novel fusion inhibitor of 5-Helix could target the S2 subunit HR2 domain of S protein to inhibit infection by pseudotyped SARS-CoV-2 and variants Delta and Omicron [22,23].However,no potent broad-spectrum small molecule inhibitors for SARSCoV-2 variants have been reported.

    In addition to the direct coronavirus-induced injury to the respiratory system, more evidence reported the hyperinflammatory reaction and immune response triggered by SARS-CoV-2 infection.Several studies have shown that inflammation-related signaling pathways, as well as cytokines can be affected by SARS-CoV-2 and thus induce hyper inflammatory reactions[24-26].The number and function of various subtypes of immune cells are severely affected in patients with COVID-19,which are associated with disease severity and can lead to severe complications[27].These evidences suggest that regulation of host inflammatory levels and immune response may be effective strategies for COVID-19 treatment.

    Herein, we report the natural product licorice-saponin A3 (A3)and glycyrrhetinic acid (GA) could widely inhibit RBDs of SARSCoV-2 Beta, Delta, and Omicron variants.GA showed inhibitory activities against Omicron on Vero E6 cells with an inhibition rate of 26.5% at 1 μM, and A3 showed an EC50of 1.016 μM.The antiviral mechanisms were elucidated by theoretical calculations and experiments.The anti-inflammatory activities of A3 were also evaluated.

    2.Materials and methods

    2.1.Cells and animals

    Vero E6 cells (Vero C1008; CRL 1586, American Type Culture Collection, Manassas, VA, USA) were stored in high glucose Dulbecco's modified Eagle medium (DMEM) (GIBCO, Thermo Fisher,Waltham, MA, USA) with penicillin (100 IU/mL), 10% fetal bovine serum (FBS, GIBCO, Thermo Fisher) and streptomycin (100 μg/mL)at 37°C in a 5% CO2environment.SARS-CoV-2 RBDs including prototype,Beta,Delta and Omicron BA.1/XBB/BQ1.1 were provided by Prof.George Fu Gao(Institute of Microbiology,Chinese Academy of Sciences,Beijing,China).ACE2 was from Novoprotein Technology Co., Ltd (Shanghai, China).SARS-CoV-2 Omicron BA.2(CSTR:16,533.06IVCAS6.7600) was propagated in Vero E6 and titrated by standard plaque assay following the standard procedure.All experiments related to authentic SARS-CoV-2 were conducted at the Zhengdian Biosafety Level 3(BSL3)facility of Wuhan Institute of Virology,China.GA(CAS No.471-53-4,HPLC≥98%)and A3(CAS No.118325-22-7, HPLC≥98%) were obtained from a compound library of the authors’ laboratory [28].PJ-34 (CAS No.344458-19-1,HPLC≥98%) was purchased from Aladdin (Shanghai, China).The protocol for animal experiments was approved by the Animal Care and Use Committee of Peking University Health Science Center(Approval number: LA2021537).

    2.2.Surface plasmon resonance (SPR) assay

    Biacore 8 K was performed for SPR study.RBDs were immobilized on the CM5 sensor chip,and the running buffer was PBS with 0.05% P20 and 5% DMSO (D8371, Solarbio, Beijing, China).RBD proteins were then diluted with sodium acetate solution (10 mM,pH 5.5)to a final concentration of 20 μg/mL.The final immobilized levels for RBDs were ~10,000 RU (response units).In binding studies, compounds were used at calculated concentrations, and flowed with buffer (30 μL/min).Both the contact time and dissociation time were set to 60 s.Biacore evaluation software was used to analyze the data.The binding affinity of ligand-protein was evaluated with the equilibrium dissociation constants(KD).Kinetics analysis and steady state affinity fitting method were used.

    2.3.Enzyme linked immunosorbent assay (ELISA)

    The inhibition of drug candidates against SARS-CoV-2 RBDs was performed as reported in the literature [29].Briefly, the S RBD proteins were diluted with coating liquid (NaHCO3, 2.93 mg/mL;Na2CO3,1.59 mg/mL;pH 9.6)and incubated at 4°C overnight.After washing with a wash solution(PBS with 0.05%Tween-20(V/V)),the microplates were blocked with 2%bull serum albumin(BSA)in the wash solution at 37°C for 40 min.After three times washing, the small molecule compounds were added at a final concentration of 0.625,1.25, 2.5, 5,10, 20, 40, and 80 μM (1 h before ACE2).Then,ACE2 was added and incubated for 30 min before being washed off.Enzyme-linked antibody SA-HRP (1:10,000) was subsequently conducted and tetramethyl benzidine was chosen as the chromogenic agent.Finally, the reaction was stopped with 1 M HCl.At OD450 nm, the enzyme activity was measured.

    2.4.Pseudovirus-based neutralization assay

    In pseudovirus-based neutralization assay [29], Beta, Delta and Omicron BA.1 pseudoviruses were provided by Prof.George Fu Gao(Institute of Microbiology, Chinese Academy of Sciences, China).Vero E6 cells were seeded at a density of 20,000 cells/well in 96-well plates.A3 and GA at different concentrations were respectively added with diluted VSV-SARS-CoV-2-EGFP virus(MOI=0.5 PFU/mL).The mixture was cultivated at 37°C for 1 h and then added into seeded Vero E6 cells.After 36 h of incubation,fluorescence images were measured by a CQ1 confocal image cytometer (Yokogawa, Tokyo, Japan) with a × 10 objective, and analyzed by a Harmony high-content analysis software (PerkinElmer, Akoya, MA, USA).

    2.5.Antiviral assay

    After being cultured with gradient concentrations of drugs overnight, Vero E6 cells in 48-well plates were infected by SARSCoV-2 Omicron BA.2 (MOI = 0.01).After 24 h of infection, RNA was extracted from the cell supernatant using QIAamp viral RNA mini kit(52,906,Qiagen,Duesseldorf,Germany).The viral RNA was quantified by quantitative real time polymerase chain reaction(qRT-PCR) with a reported primer pair [29].

    2.6.Hydrogen-deuterium exchange mass spectrometry (HDX-MS)analysis

    RBD proteins(1 mg/mL),and RBD proteins(1 mg/mL)with GA/A3 (1 mM) were respectively deuterium labeled with 100 mM phosphate buffer(D2O,pD 7.0,20-fold dilution).After 0.083,0.25,1,10,30,60,and 240 min,the quenching buffer(100 mM phosphate,0.5 M tris(2-carboxyethyl)phosphine (TCEP), 4 M GdHCl, pH 2.0)was added to stop the labeling reaction.Samples were then injected and digested online on Waters ENZYMATE BEH pepsin columns(2.1 mm × 30 mm, 5 μm; Milford, MA, USA).The peptides were trapped and desalted on a VanGuard Pre-Column trap (ACQUITY UPLC BEH C18,1.7 μm) for 3 min, eluted from the trap performing 15%acetonitrile at 100 μL/min,then separated by an ACQUITY UPLC BEH C18column (1.0 mm × 100 mm, 1.7 μm, Waters).All mass spectra were acquired on a mass spectrometer (Waters Xevo G2).Peptides from an unlabeled protein were identified by applying ProteinLynx Global Server (PLGS) searches of a protein database including RBD sequences only.Peptides in relative deuterium levels were measured by subtracting the mass of the undeuterated control sample from that of the deuterium-labeled sample [30].

    2.7.Proteomics and phosphoproteomics analysis

    Vero E6 cells were incubated with or without 10 μM A3 for 8 h,and subsequently infected with SARS-CoV-2 for 24 h (MOI = 1).A high intensity ultrasonic processor (Scientz, Ningbo, China) was used to sonicate the samples in lysis buffer(8 mM urea,1%protease inhibitor cocktail).The first digestion was performed overnight with a 1:50 trypsin-to-protein mass ratio.Tandem mass tags(TMT)reagents(Thermo Fisher Scientific)were used to label each channel of peptide, and each channel was incubated for 2 h at room temperature.We desalted the pooled samples using a Strata X C18SPE column(Phenomenex,Torrance,CA,USA)and dried the samples by vacuum centrifugation.Solvent A (0.1% formic acid and 2% acetonitrile in water) was used to dissolve tryptic peptides, which were directly loaded onto a reversed-phase analytical column (25 cm length, 75 μm i.d.).Using solvent B (0.1% formic acid in 90%acetonitrile),peptides were separated from 5%to 25%over 60 min,25%-35% in 22 min, up to 80% in 4 min, and then held at 80% for 4 min.The UPLC system (ThermoFisher Scientific) operated at a constant flow rate of 450 nL/min.Then the separated peptides were characterized by a Q Exactive?HF-X mass spectrometer (ThermoFisher Scientific) with a nano-electrospray ion source[30].

    2.8.Molecular docking and quantum mechanics/molecular mechanics (QM/MM) simulations

    Omicron S protein (PDB ID: 7QTK) and SARS-CoV-2 S protein(PDB ID:6M0J)were used as templates[31,32].The structures of A3 and GA were optimized using density functional theory (DFT) at B3LYP/6-311G(d)[33-35].Molecular docking of molecules and RBD was carried out by AutoDock Vina and AutoDocktools v1.5.6 software [36-38].Furthermore, QM/MM simulations were applied by Gaussian16 to optimize the ligand and residues (400-422 and 453-470) from the docking model.The B3LYP/6-311G in DFT calculations was used to describe the QM atoms from A3 and Y453/Y453F, and the UFF force field was calculated to simulate the MM atoms [39].

    2.9.Western blotting

    Proteomics samples were determined using BCA Protein Content Assay Kit (AKPR017 Boxbio, Beijing, China).After boiling for 5 min, 20 μg of each sample was added to 10% sodium dodecyl sulfate (SDS)-polyacrylamide minigels (EpiZyme, Shanghai, China)and electrophoresed at 80 V (stacking gel) and 120 V (separating gel).A polyvinylidene fluoride (PVDF) membrane (Bio-Rad, Hercules,CA,USA)was then used to transfer proteins out of the gel,and was blocked for 1 h at room temperature in the blocking buffer(Beyotime, Shanghai, China).Primary and secondary antibodies were then incubated overnight at 4°C on the membrane (nuclear factor-kappa-B (NF-κB) (AF5006, Affinity, Liyang, Jiangsu, China),1:1000;PAK1/2/3(2604 T,Cell Signaling,Boston,MA,USA),1:1000;p38α (9218 T, Cell Signaling), 1:1000; c-Jun (BF0245, Affinity),1:1000; JNK2 (DF2357, Affinity),1:1000; STAT1 (AF6300, Affinity),1:1000; MEK4 (4514 T, Cell Signaling), 1:1000).Signals were imaged by Bio-Rad ChemiDoc XRS.GAPDH was the loading control.

    2.10.Lipopolysaccharide(LPS)-induced mice acute lung injury(ALI)model, pathological analysis, and pharmacokinetic study

    Male BALB/c mice (20 g) were obtained from the Experimental Animal Center of Peking University Health Science Center(Beijing,China).Mice in the control,model,and A3 groups were treated with saline solution,LPS(2 mg/kg,i.n.),and LPS(2 mg/kg,i.n.)with A3(10 and 20 mg/kg,i.g.),respectively.After 8 h,blood and lung tissue samples were collected.Small pieces of lung tissues of ALI model were fixed in 4% paraformaldehyde.After dehydration of the specimens in graded ethanol and embedding in paraplast, the sections were cut at 5-μm thickness.Then, hematoxylin and eosin(H&E) staining was performed to demonstrate the histological structures.Images were taken by WISLEAP (WS-10, Changzhou,Jiangsu, China).

    For pharmacokinetic study, A3 (20 mg/kg) was given to BALB/c mice (20 g) by intragastric administration.Blood samples were collected at 1,4,8,10,12,24,and 48 h,respectively,and centrifuged at 15,000 rpm(4°C)for 30 min to obtain the plasma.Then 100 μL of plasma sample, 200 μL of methanol, and 100 μL internal standard solution (5 mg/mL astragaloside IV in methanol) were mixed.The mixture was vortexed and then centrifuged at 15,000 rpm for 30 min.The supernatant was separated, dried, and then dissolved in 100 μL of 50% methanol.The sample was centrifuged(15,000 rpm for 30 min), and then analyzed by Q-Exactive hybrid quadrupole Orbitrap mass spectrometer (Thermo Scientific, San Jose, CA, USA) equipped with a heated electrospray ionization source(HESI).The conditions were as follows: negative ion mode;spray voltage, -3.5 kV; sheath gas flow rate, 45 arb; aux gas flow rate, 10 arb; capillary temperature, 350°C; aux gas heater temperature, 4000°C; and S-lens RF level, 60 V.Samples were separated on a Waters Acquity UPLC BEH C18column(1.7 μm×150 mm)(Waters).The mobile phase consisted of water containing 0.1%formic acid(V/V,A)and acetonitrile containing 0.1%formic acid(V/V,B).The gradient elution program was as follows:0-5 min,10%-30%; 5-8 min, 30%-50% B; 8-15 min, 50%-100% B; 15-18 min,100%-10% B.The flow rate was 0.3 mL/min.The injection volume was 2 μL.A3 and astragaloside IV were determined by multiple reaction monitoring (MRM) in the negative ion mode, with MRM transitions of m/z 983.5 →821.3 and 783.4 →651.4, respectively.

    2.11.Immunohistochemistry (IHC)

    We applied the same lung tissues of ALI model at 5 μm for IHC.Following 10 min of incubation with 3% H2O2, the samples were incubated in citrate buffer (pH 8.0) for 3 min at high heat and 15 min at low heat for antigen retrieval.Sections were incubated with primary antibodies for 2 h after blocking with goat serum(tumour necrosis factor alpha (TNF-α) (Abcam, Cambridge, United Kingdom),1:20;interleukin-6 (IL-6)(Abcam),1:50;interleukin-1β(IL-1β)(Abcam),1:250,interleukin-7(IL-7)(Thermo Fisher),1:100).Phosphate-buffered saline was used to wash the sections three times and anti-horseradish peroxidase (HRP) goat anti-rabbit IgG was used for incubation (PV-6001, ZSGB-BIO, Beijing, China).3,3'-diaminobenzidine (DAB) was incubated for 2 min.The sections were then stained with hematoxylin after the reaction was stopped.Images were taken using WISLEAP(WS-10).

    2.12.Multi fluorescent IHC (mIHC)

    The 5-μm sections obtained from the lung tissues of ALI model were incubated in citrate buffer (pH 6.0, CD3 (Abcam), Foxp3(Abcam) and iNOS (Affinity)) or ethylenediamine tetraacetic acid(EDTA) (pH 9.0, CD4 (Abcam) and CD8 (Abcam)) for 3 min high heat, and 15 min low heat (microwave oven).There were five sequential rounds of staining for each section, including a protein block with goat serum,a primary antibody(CD3,1:200;CD4,1:125;CD8,1:75;iNOS,1:75 and Foxp3,1:375),and a secondary goat antirabbit IgG HRP stain.By using tyramide signal amplification (CD3,fluorescein isothiocyanate(FITC);CD4,Texas;CD8,CY5;iNOS,CY3 and Foxp3, Texas), each horseradish peroxidase-conjugated polymer covalently bound a different fluorophore.To remove bound antibodies before the next step in the sequence,additional antigen retrieval using heated citrate buffer(pH 6.0,CD3,Foxp3 and iNOS)or EDTA (pH 9.0, CD4 and CD8) was performed after this covalent reaction.The sections were counterstained with 4',6-diamidino-2-phenylindole (DAPI) (AlphaTSA, Beijing, China) and mounted with fluorescence mounting medium after completing all six sequential reactions.Vectra Polaris Imaging System version (PerkinElmer,Waltham, MA, USA) was used to image multiplex stained slides.Tissue sections stained with seven-plex fluorescence were scanned using the ×40 objective, with every 10 nm of the emission light spectrum across the range of each emission filter cube.Filter cubes used for multispectral imaging were DAPI (440-680 nm), FITC(520-680 nm), Cy3 (570-690 nm), Texas red (580-700 nm) and Cy5 (670-720 nm).A spectral library containing the emitting spectral peaks of all fluorophores was created using multispectral images obtained from single stained slides for each marker and associated fluorophore.This spectral library was then used to separate each multispectral image cube into its individual components (spectral unmixing) allowing for the color-based identification of all markers of interest in a single image using the inForm 2.6 image analysis software (Akoya, MA, USA).All spectrally unmixed and segmented images were subsequently subjected to a proprietary inForm active learning phenotyping algorithm, allowing for the phenotyping and subsequent quantifications of different cell phenotypes.

    2.13.Limitations of the study

    The complex crystal structure of SARS-CoV-2 RBD with A3 needs to be solved to further interpret the inhibition mechanisms.Moreover,SARS-CoV-2 infected animals treated with A3 should be investigated in future studies to confirm the anti-virus activity of A3 and further explain the anti-inflammatory capability of A3 in treating SARS-CoV-2.The anti-inflammation and immune regulation efficacy of A3 administered at different time points or with different doses in ALI and acute respiratory distress syndrome(ARDS) models are of vital importance in evaluating its druggability, which should also be studied in the future.

    Moreover,it is of vital importance to construct different types of inflammation animal models to confirm the anti-inflammatory activity of A3.Analysis of more detailed subtypes of immune cells will also be helpful in explaining the immune regulation capability of A3,as well as interpreting the correlation between the immune regulation and inflammation inhibition.

    3.Results and discussion

    3.1.Inhibitory activities of A3 and GA against RBDs of SARS-CoV-2

    Potent inhibitory activities of A3 and GA against the prototype SARS-CoV-2 RBD have been reported[29].Herein,we showed that A3 and GA could inhibit RBDs of SARS-CoV-2 variants Beta, Delta,and Omicron BA.1/XBB/BQ1.1 evaluated by measuring the binding capacity between A3/GA and RBDs of the three SARS-CoV-2 variants by ELISA, respectively.Both A3 and GA exhibited strong inhibitory activities against the RBDs of Beta (A3,IC50= 4.99 ± 2.86 μM; GA, IC50= 9.12 ± 2.52 μM), Delta (A3,IC50=8.25±1.06 μM;GA,IC50=12.4±0.06 μM),Omicron BA.1(A3,IC50=6.26±0.33 μM;GA,IC50=2.37±0.16 μM),Omicron XBB(A3,IC50= 4.36 ± 3.82 μM; GA, IC50= 8.47 ± 1.32 μM), and Omicron BQ1.1 (A3, IC50= 22.98 ± 1.91 μM; GA, IC50= 17.93 ± 2.9 μM)(Fig.1A).Subsequently,the binding of A3 and GA with RBDs of Beta,Delta and Omicron variants were determined by SPR,respectively.Not surprisingly,both GA and A3 showed binding capabilities with RBDs of SARS-CoV-2 and the five variants.Particularly, A3 presented increased binding affinity to Omicron variants, with KDvalues of 9.1, 4.66, and 8.28 μM, respectively (Fig.1B).

    3.2.Weak inhibitory activities of A3 and GA against SARS-CoV RBD

    SARS-CoV is another human coronavirus,and possesses a highly similar gene sequence to SARS-CoV-2 [40].Most proteins of these two viruses share high homology(95%-100%)[41].However,there is only 60% similarity between the RBD sequences of SARS-CoV-2 and SARS-CoV (Fig.S1).Considering the abundant licorice compound glycyrrhizic acid (GA-gg) had been reported as an active compound to inhibit replication of SARS-CoV in 2003 [42], we tested the inhibitory activities of A3,GA,and GA-gg,against SARSCoV RBD.Consistently,GA-gg could inhibit SARS-CoV RBD with IC50of 62.95±7.02 μM(Fig.S2).However,both A3 and GA showed weak inhibitory activities against SARS-CoV RBD with IC50>80 μM(Fig.2A), though SPR indicated A3 and GA could target SARS-CoV RBD with KDvalues of 61.9 and 10.9 μM, respectively(Fig.2B).

    Fig.1.Enzyme linked immunosorbent assay(ELISA)and surface plasmon resonance(SPR)assays of licorice-saponin A3(A3)and glycyrrhetinic acid(GA)binding with S receptorbinding domains(RBDs)of different severe acute respiratory syndrome coronavirus 2(SARS-CoV-2)variants.(A)Inhibitory activities of A3(upper)or GA(lower)against the S RBD of SARS-CoV-2 and its variants Beta, Delta,and Omicron BA.1/XBB/BQ1.1 by ELISA.RBD proteins of SARS-CoV-2(0.2 μg/mL),Beta(0.3 μg/mL),Delta(0.2 μg/mL)and Omicron BA.1/XBB/BQ1.1 (0.5 μg/mL) were respectively coated onto 96-well microplates.After incubation with A3 or GA for 1 h, ACE2 at the same concentration as RBD protein was added and incubated with the S protein for 30 min.(B)SPR analysis of A3(upper) or GA(lower)binding to RBDs of SARS-CoV-2,or SARS-CoV-2 variants(Beta,Delta,and Omicron BA.1/XBB/BQ1.1).

    To elucidate the binding mechanisms between A3/GA and SARSCoV RBD, we employed HDX-MS to determine protein conformation changes.Deuterium uptakes of A3 and GA with SARS-CoV RBD increased for majority of the peptides from 330 to 503 (Fig.2C),demonstrating that uptake of the ligand-protein structures was slower than the individual SARS-CoV protein.It is noteworthy that deuterium uptakes of peptides 362-388 and 413-427 increased remarkably, indicating that binding residues of A3 and GA with SARS-CoV RBD may locate in these two peptides (Figs.2C and D).Y453 had been considered as a crucial residue for inhibitory activity of A3 and GA against SARS-CoV-2 RBD(Fig.2E)[29].As Y453 is not located within these two peptides,this may be the reason why A3 and GA could bind to SARS-CoV RBD but showed weak inhibitory activities.

    3.3.A3 inhibits SARS-CoV-2 Omicron variants through binding site Y453 of RBD

    We had reported A3 and GA could potently inhibit SARS-CoV-2 in Vero E6 cells with EC50of 75 nM and 3.17 μM,respectively.One inhibitory mechanism was hydrogen bonding with Y453 of RBD[29].To confirm the significance of Y453, we conducted HDX-MS analysis.The receptor-binding motif (RBM) of S protein to bind with ACE2 consists of residues 437 to 508.According to our experimental results, the uptake of peptides 448-512 which mainly located in the RBM showed obvious changes upon A3 or GA treatment (Figs.2F and G).These changes may be caused by hydrogen binding of A3 or GA with Y453 (Fig.2H) [29].

    The SARS-CoV-2 Omicron BA.1 carries 15 mutations (Fig.3A),but not including Y453.Both A3 and GA showed noticeable inhibitory activities against the RBD of Omicron (Fig.1A).HDX-MS analysis indicated the uptake of peptide 453-470 with A3 and GA was remarkably increased(Figs.3B-D).Subsequently,we used the structure of Omicron RBD(PDB ID:7QTK)to simulate binding with A3.Two hydrogen bonds were respectively formed between A3 and Y453/R408, and between GA and Y453/N417 (Figs.3E and F).A3 showed higher inhibitions against SARS-CoV-2 Beta, Delta, and Omicron BA.1 pseudoviruses than GA (Fig.S3).Particularly, A3 showed an EC50of 1.016 μM against SARS-CoV-2 Omicron in Vero E6 cells (Fig.S4).

    To confirm the significance of Y453 in the binding between A3 and Omicron RBD, we constructed the Y453F mutant.QM/MM calculations demonstrate the binding energy of A3 with peptide 453-470 of Omicron RBD was -26.7 kcal/mol, which was much lower than that with peptide 453-470 of the Omicron Y453F mutant (4.37 kcal/mol) (Fig.3G), indicating a stable confirmation with Omicron.Moreover, the phenolic hydroxyl group of Y453 could form a hydrogen bond with-OH of C(6’)-OOH of A3(2.8 ?),while no hydrogen bond could form for the Y453F mutant(Figs.3H and I).These results confirmed that Y453 was a crucial site in the binding between A3 and Omicron RBD.

    Fig.2.Binding mechanisms of licorice-saponin A3(A3)targeting S receptor-binding domains(RBDs)of severe acute respiratory syndrome coronavirus(SARS-CoV)and SARS-CoV-2.(A)Inhibitory activities of A3(upper)and glycyrrhetinic acid(GA)(lower)against the S RBD of SARS-CoV,performed by enzyme linked immunosorbent assay(ELISA).RBD protein of SARS-CoV (0.3 μg/mL) was coated onto 96-well microplates.After incubation with test drugs (1 h), angiotensin converting enzyme 2 (ACE2) at the same concentration as RBD protein was added and incubated with the S protein for 30 min.(B)Surface plasmon resonance(SPR)analysis of A3(upper)and GA(lower)binding to SARS-CoV RBD.(C)Deuterium uptake differences of SARS-CoV RBD with A3 (upper) and GA (lower), respectively.Sum is total uptakes from 0 s to 4 h.(D) Deuterium uptake plots of two peptides 413-427 and 362-388 of SARS-CoV RBD with A3 and GA.(E)Location of peptides 413-427 and 362-388 in SARS-CoV S RBD.(F)Location of peptides 448-465,471-487,and 487-512 in SARSCoV-2 S RBD.(G)Deuterium uptake differences of SARS-CoV-2 RBD with A3(upper)and GA(lower),respectively.Sum is total uptakes from 0 s to 4 h.(H)Deuterium uptake plots of two peptides 448-465 and 487-512 of SARS-CoV-2 RBD with A3 and GA.

    3.4.A3 regulates the mitogen-activated protein kinase (MAPK)signaling pathway of SARS-CoV-2 infected host cells

    Fig.3.Binding mechanisms of licorice-saponin A3 (A3) and glycyrrhetinic acid (GA) targeting the S receptor-binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron BA.1.(A) 15 mutations of Omicron S RBD compared to wild type SARS-CoV-2.(B) Deuterium uptake differences between Omicron RBD,and Omicron RBD with A3 or GA.Sum is total uptakes from 0 s to 4 h.(C)Deuterium uptake plots of peptides 400-422 and 453-470 of Omicron RBD,Omicron RBD with A3,and Omicron RBD with GA.(D) Mass spectra of peptides 400-422 and 453-470 of Omicron RBD and Omicron RBD with A3/GA.(E) Location of peptides 400-422 and 453-470 in Omicron RBD.(F) Molecular docking of A3 and GA with Omicron RBD, respectively.Hydrogen bonds (red dashes) are shown.(G) Relative energies of A3 with RBD (peptides 400-422 and 453-470), and A3 with RBDY453F (peptides 400-422 and 453-470)computed by quantum mechanics/molecular mechanics (QM/MM).(H,I) Binding modes of A3 with active residues of RBD (H) and RBDY453F (I), respectively.Hydrogen bonds (red dashes) are shown.

    Fig.4.Phosphoproteomics profiling of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-infected cells treated with licorice-saponin A3 (A3).(A) Experimental scheme.Vero E6 cells were infected with SARS-CoV-2(MOI=1)in the presence or absence of 10 μM A3 for 24 h.Proteins were prepared using lysis buffer(8 M urea,1%protease inhibitor cocktail)and trypsin.(B)Down-regulated viral proteins after A3 treatment(fold change >1.3),popularly reported proteins presented in pink color.(C)Inhibitory activity of A3 against nucleoprotein at 20 μM determined by enzyme linked immunosorbent assay(ELISA)kit(Abclonal,https://abclonal.com.cn/),and the positive control was PJ-34[50].(D)Four categories Q1-Q4 were divided by A3/mock ratio (<0.667, 0.667-0.769,1.3-1.5 and >1.5), and the numbers of their phosphoproteomic sites were presented.Mock group:Vero E6 cells infected by SARS-CoV-2 for 24 h; A3 group: Vero E6 cells infected by SARS-CoV-2 in the presence of A3 (10 μM) for 24 h, n = 3.(E) Enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for categories Q1-Q4.(F) Mitogen-activated protein kinase (MAPK) signaling pathway analysis to explain the antiinflammatory mechanism of A3.(G) Expression of proteins participated in inflammation-related MAPK signaling pathway in SARS-CoV-2 infected Vero E6 cells before and after the treatment of A3.Mock group:Vero E6 cells infected by SARS-CoV-2 for 24 h;A3 group:Vero E6 cells infected by SARS-CoV-2 in the presence of A3(10 μM)for 24 h,n=3.(H)Interleukin (IL)-1, tumor necrosis factor (TNF)-α, IL-6, IL-1β, IL-8, IL-7, function of interferons (IFN)-β and IFN-γ of proteomic samples were measured by monkey ELISA kit(MEIMIAN, www.mmbio.cn).Vero E6 cells respectively treated with SARS-CoV-2 in the absence (control), and presence of A3 for 24 h.*P <0.05, **P <0.01, compared with the control group, n = 3.LC-MS/MS: liquid chromatography coupled to tandem mass spectrometry; TMT: tandem mass tags.

    Fig.5.Anti-inflammatory effects of licorice-saponin A3 (A3) in lipopolysaccharide (LPS)-induced mice acute lung injury (ALI) mice model.(A) Experimental scheme to construct LPS-induced ALI mice model.(B) Representative histopathology of lung sections of ALI mice with or without A3 administration.Control group: healthy mice treated with saline solution; Model group: LPS-induced ALI mice without A3 administration; A3 groups: LPS-induced ALI mice with A3 administration in a dose of 10 or 20 mg/kg, i.g., respectively,n=4.The scale bars represent 500 μm for ×5 magnification images and 100 μm for×20 magnification images.(C)Pathology evaluation for the therapeutic effects of A3 in LPSinduced ALI as shown in Fig.5B, n = 4, **P <0.01, ***P <0.001.(D) Tumor necrosis factor (TNF)-α, interleukin (IL)-6, IL-1β and IL-7 of blood samples were measured by mouse enzyme linked immunosorbent assay (ELISA) kit (MEIMIAN, www.mmbio.cn).*P <0.05, **P <0.01 compared with the control group, n = 4.(E) Representative immunohistochemistry (IHC) of lung sections of LPS-induced ALI mice with or without A3 administration for TNF-α, IL-6, IL-1β and IL-7 analysis.

    To explore the effects of A3 on SARS-CoV-2 infected host cells,we conducted proteomics analysis of Vero E6 cells infected by SARS-CoV-2 with the treatment of A3(10 μM,pretreated for 8 h and infected for 24 h) (Fig.4A).The expression levels of 1 protein in virus and 20 proteins in host cells were significantly altered after A3 treatment (Fig.S5A).We further conducted phosphoproteomic analysis,according to the significant cellular regulatory mechanism of protein phosphorylation [43].The results showed 17 downregulated phosphoproteins in virus, together with 1070 downregulated and 31 up-regulated phosphoproteins in host cells (fold change >1.3) (Fig.S5B).Among the 17 down-regulated proteins,aside from the S protein,we identified two other structural proteins including the membrane protein (M) and the nucleoprotein (N),which are both targets for SARS-CoV-2 viral inhibition(Fig.4B)[7].Moreover, one popular functional protein of replicase polyprotein 1 ab(rep)was determined[44].Subsequently,enrichment analysis of biological processes performed by gene ontology annotations(P value <0.05)was conducted,and several biological activities of the proteins were listed(Figs.S6 and S7).In order to evaluate whether A3 could inhibit the down-regulated proteins, we measured the inhibitory activity of A3 against nucleoprotein whose quantity was decreased with an A3/mock ratio of 0.593 by ELISA.As a result,A3 could inhibit nucleoprotein by 41.1%at 20 μM(Fig.4C).These data demonstrated that A3 could not only down-regulate the protein expression, but also inhibit their activities.

    Next, phosphoproteins in the host cells regulated by A3 in phosphoproteomic analysis were studied.GO annotations indicated that 128 proteins were related to immune system (Fig.S8),and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis presented that A3 remarkably regulated inflammation through MAPK signaling pathway in four categories Q1-Q4 divided by A3/mock ratio(Figs.4D and E).

    Among the proteins related to anti-inflammation(Table S1),PAK,p38/p38γ, NF-κB, MEK4, JNK, STAT1 and Jun respectively participated in the c-Jun amino-terminal kinases 1/2/3 and p38(α,β,γ,and δ)MAPK signaling pathways were obviously down-regulated.Thus,we deduced the anti-inflammatory activity of A3 was related to the JNK and p38 MAPK pathways (Fig.4F).Consistently, Western blotting analysis confirmed that PAK1/2,p38(α,β)and NF-κB in the p38 MAPK pathway, and MEK4, JNK, STAT1 and Jun in JNK1/2/3 were remarkably regulated(Figs.4G and S9-S16).Because the expression of IL-6,TNF-α,and IL-1β could be suppressed by inhibiting NF-κB and JNK activation[26],the protein samples of Vero E6 cells treated with SARS-CoV-2 (control), and SARS-CoV-2 combined with A3 were respectively measured by ELISA[30].As a result,cytokines including IL-1β,IL-6,TNF-α,IL-1,IL-7,IL-8,IFN-β,and IFN-γ were all decreased significantly upon A3 treatment (Fig.4H).The above results indicated that A3 possessed anti-inflammatory activities of SARS-CoV-2 infected host cells through regulating key pro-inflammation proteins involved in the MAPK signaling pathway.

    Fig.6.Multi-spectral imaging and quantification analysis of the immune phenotype in lung tissues of lipopolysaccharide(LPS)-induced mice acute lung injury(ALI)mice treated by licorice-saponin A3 (A3).(A, B)Representative images displaying the lung tissues after multispectral imaging(raw image)and after spectral unmixing (composite image).Mice in the model(A)and A3 groups(B)were respectively administered LPS(2 mg/kg,i.t.),and LPS(2 mg/kg,i.t.)coupled with A3(20 mg/kg,i.g.).The markers included CD3(membrane,fluorescein isothiocyanate(FITC),red),CD4(membrane,Texas,green),CD8(membrane,CY5,blue),iNOS(membrane,CY3,cyan)and Foxp3(nuclear,Texas,orange).The scale bars represent 100 μm for the single channel images.(C) Positive rate of the five markers in the control, model, and A3 groups.*P <0.05, **P <0.01, ***P <0.001 compared with the control group, n = 3.

    3.5.Host-directed anti-inflammatory activity of A3 in vivo

    To further evaluate the anti-inflammatory potential and to explore the anti-inflammatory mechanisms of A3 in vivo, we constructed an ALI mice model as demonstrated in Fig.5A [45].After treatment with 2 mg/kg LPS by nasal intubation drip(i.t.),the two groups of mice were respectively treated with A3 at a dose of 10 mg/kg (low dose) or 20 mg/kg (high dose) by intragastric administration (i.g.).Mice treated with saline without LPS administration were used as the control.The mice were sacrificed 8 h after LPS treatment,and histopathology changes of the lung tissues were examined by H&E staining (Figs.5B and S17).After analyzing and scoring the histopathology injuries of the excised lung tissue sections(Fig.5C),we found the injury levels in lung section of the ALI mice treated with A3 were obviously alleviated.Moreover, the high-dose group showed more potent anti-inflammatory effects than the low-dose group.

    Subsequently,the levels of proinflammatory cytokines IL-1β,IL-6, TNF-α and IL-7 of blood samples from ALI mice treated with A3 were determined (Fig.5D).IL-7 was decreased significantly upon A3 treatment.Moreover, high-dose A3 (20 mg/kg) was more effective than the low-dose (10 mg/kg).Then, IHC staining of ALI mice lung section was conducted(Fig.5E),IL-1β,IL-6,TNF-α and IL-7 were remarkably down-regulated after A3 treatment.The results provided solid evidence that A3 could remarkably inhibit inflammation by directly down-regulating the pro-inflammatory cytokines in host cells.

    3.6.Host-directed immunoregulatory activity of A3

    As p38 and JNK play critical roles in the innate and adaptive immune systems [46,47], and immune-modulating approaches are commonly used in inflammation treatment,we studied the infiltrates of differentsubtypesofimmune cells inthe lungtissues ofALI mice.Tcell subpopulations including CD3+total T cells,CD4+T helper cells(Th cells),CD8+cytotoxic T lymphocytes(CTL),and Foxp3+ regulatory T cells (Tregs), as well as iNOS + M1 macrophages of the lung tissue samples from treated or untreated ALI mice as well as control micewere analyzed by multiplex fluorescent immunohistochemistry(mIHC) staining(Fig.6).Multispectral imaging followed by spectral unmixing of five biomarkers in lung tissue samples of the three groups could provide the spatial distribution information and positive rates of each immune cell subtype (Figs.6A, 6B and S18).The results showed the CD4+Th cells,CD8+CTL and iNOS+M1 macrophages,as well as the immune-suppressive Foxp3+Tregs in lungs of ALI mice were generally improved compared with that in the control group without LPS, and were remarkably decreased after A3 treatment (n = 3) (Fig.6C).This finding suggested LPS-induced inflammation caused a broad immune dysregulation and that A3 could rebalance the immune system.

    4.Discussion

    In this work,we report the natural product A3 from the traditional herbal medicine licorice(Gan-Cao)could inhibit RBDs of SARS-CoV-2 variants Alpha, Delta, and Omicron BA.1/XBB/BQ1.1.Particularly, A3 could potently inhibit SARS-CoV-2 Omicron in Vero E6 cells with an EC50of 1.016 μM(Fig.S4).While several peptide or protein therapies could inhibit the entry of SARS-CoV-2 variants [21,22,48], A3 is the first broad-spectrum small molecule inhibitor of RBDs.

    Considering that Y453 of the S protein RBD had been reported as a key amino acid in the hydrogen bonding between A3 and SARSCoV-2 RBD [29], we explored the significance of Y453 in the inhibition of SARS-CoV-2 Omicron variant by A3.HDX-MS analysis and molecular docking revealed that A3 inhibited the Omicron RBD also through binding to Y453.Moreover,QM/MM calculations indicated the binding energy of A3 with peptide 453-470 of Omicron RBD(-26.7 kcal/mol) was much lower than that with Omicron Y453F mutant (4.37 kcal/mol) (Fig.3G), suggesting a stable binding between A3 and Omicron RBD mediated by Y453.

    For severe COVID-19 patients,overproduction of early response proinflammatory cytokines IL-1β, IL-6, and TNF-α can result in cytokine storm with symptoms like fever,cough,and diarrhea[26].In the SARS-CoV-2-induced hyperinflammatory reaction, several reports have shown that the IL-1-IL-6 axis is the most biologically relevant signal[49].Since protein phosphorylation is an important cellular regulatory mechanism, including inflammation [42],phosphoproteomic analysis of Vero E6 cells treated with A3 and then infected by SARS-CoV-2 was conducted to explore the antiinflammatory mechanisms of A3 for COVID-19.Surprisingly,inflammation related proteins in JNK and p38 MAPK signaling pathways, including PAK, p38/p38γ, NF-κB, MEK4, JNK, STAT1 and Jun(Table S1),were all down-regulated by A3.Moreover,A3 could significantly and dose-dependently decrease the levels of IL-1β,IL-6 and TNF-α in both blood samples and lung tissues of mice in the ALI mice model.These data collectively suggested that A3 could directly inhibit inflammation by regulating key proinflammation proteins involved in the JNK and p38 MAPK signaling pathways.

    Dysregulation of host immune responses is also a feature prominently in COVID-19 and immunomodulatory is another COVID-19 treatment strategy besides antiviral agents,mainly used for inhibiting the SARS-CoV-2-induced hyperinflammatory reaction[48].Inspired by this,we explored the potential host-directed therapeutic capabilities of A3 by studying the immune phenotype in lung tissues of LPS-induced ALI mice treated with A3.Spatial visualization and quantification of several subtypes of immune cells infiltrating in lungs of ALI mice, including the immune activable CD4+Th cells,CD8+CTL and iNOS+M1 macrophages,as well as the immune-suppressive Foxp3+ Tregs.Results suggested that the four subtypes of immune cells were generally decreased after A3 treatment(n=3)(Fig.6C),indicating that A3 could restore the immune dysregulation of host induced by inflammation, further demonstrating the anti-inflammatory activities of A3 [48].

    Furthermore,we monitored the pharmacokinetics of A3 in mice after 20 mg/kg administration (i.g.) (Fig.S19), and found that A3 was metabolically stable with a half-life of 7.13 h.Majority of the drugs were present in the circulation in its parent form.Furthermore, A3 showed little toxicity in both the cell line and the mice model, suggesting a fine biocompatibility.

    5.Conclusion

    In this work,we report that the natural product A3 could widely inhibit RBDs of SARS-CoV-2 variants, including Beta, Delta, and Omicron BA.1/XBB/BQ1.1.Furthermore, we find that A3 could directly suppress host inflammation by modulating the JNK and p38 MAPK pathways, thus restoring the corresponding immune dysregulation.Given its good safety and druggability,A3 could be a promising broad-spectrum small molecule drug candidate for the prevention and treatment of COVID-19.

    CRediT author statement

    Yang Yi:Conceptualization, Methodology, Validation, Formal analysis, Investigation, Data curation, Writing - Original draft preparation, Reviewing and Editing, Visualization, Funding acquisition;Wenzhe Li:Methodology, Validation, Formal analysis, Data curation, Writing - Reviewing and Editing, Funding acquisition;Kefang Liu:Conceptualization, Methodology, Validation, Formal analysis,Data curation;Heng Xue:Validation,Formal analysis,Data curation;Rong Yu:Validation;Meng Zhang:Validation;Yang-Oujie Bao:Validation;Xinyuan Lai:Validation;Jingjing Fan:Validation;Yuxi Huang:Validation;Jing Wang:Validation;Xiaomeng Shi:Validation;Junhua Li:Validation;Hongping Wei:Formal analysis, Data curation;Kuanhui Xiang:Formal analysis,Data curation;Linjie Li:Validation;Rong Zhang:Validation;Xin Zhao:Formal analysis, Data curation;Xue Qiao:Formal analysis,Funding acquisition;Hang Yang:Formal analysis, Data curation,Writing-Reviewing and Editing,Supervision,Funding acquisition.Min Ye:Conceptualization,Methodology,Formal analysis,Writing-Reviewing and Editing, Supervision, Funding acquisition.

    Declaration of competing interest

    The authors declare that there are no conflicts of interest.

    Acknowledgments

    This work was supported by National Natural Science Foundation of China (Grant Nos.: 81891010/81891011, 81725023,82003614, 82173950, 31770192, 32070187, 32161133003 and 82003681) and China Postdoctoral Science Foundation (Grant No:2022T150029).We thank Professor George Fu Gao from CAS Institute of Microbiology for helpful and constructive suggestions.We also thank Tao Du,Lun Wang,Jin Xiong,and the entire running team from Zhengdian Biosafety Level 3 Laboratory of Wuhan Institute of Virology for technical support.We thank the Wroclaw Center for Networking and Super Computing for providing generous computer time, and Jingjie PTM BioLab Co., Ltd.(Hangzhou, China) for proteomics mass spectrometry analysis.Appendix A.Supplementary data

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

    99热精品在线国产| 国产精品野战在线观看| 在线观看午夜福利视频| 亚洲无线在线观看| 久久亚洲精品不卡| 日韩人妻高清精品专区| 在线国产一区二区在线| 免费黄网站久久成人精品| 亚洲欧美日韩高清专用| 青青草视频在线视频观看| 最新中文字幕久久久久| 成人特级黄色片久久久久久久| 欧美潮喷喷水| 国产一级毛片在线| 午夜精品在线福利| 深爱激情五月婷婷| 国产高清视频在线观看网站| 亚洲精品亚洲一区二区| 欧美变态另类bdsm刘玥| 麻豆av噜噜一区二区三区| 欧美zozozo另类| 一边亲一边摸免费视频| 国产爱豆传媒在线观看| 搞女人的毛片| 国产午夜精品论理片| 日韩欧美在线乱码| 亚洲精品456在线播放app| 亚洲高清免费不卡视频| 亚洲一级一片aⅴ在线观看| 91久久精品电影网| kizo精华| 午夜精品一区二区三区免费看| 亚洲国产色片| 日日摸夜夜添夜夜爱| 我要看日韩黄色一级片| 99热只有精品国产| 亚洲精品自拍成人| 成人高潮视频无遮挡免费网站| 十八禁国产超污无遮挡网站| 给我免费播放毛片高清在线观看| a级毛片免费高清观看在线播放| 只有这里有精品99| 婷婷精品国产亚洲av| 日本一本二区三区精品| 日韩一区二区三区影片| 成人鲁丝片一二三区免费| 亚洲精品456在线播放app| 日韩 亚洲 欧美在线| 久久6这里有精品| 97超视频在线观看视频| 男人舔奶头视频| 婷婷亚洲欧美| 免费一级毛片在线播放高清视频| 天堂网av新在线| 联通29元200g的流量卡| 午夜精品一区二区三区免费看| 日韩在线高清观看一区二区三区| 成人综合一区亚洲| 天堂网av新在线| 真实男女啪啪啪动态图| 国产乱人偷精品视频| 99久久人妻综合| 国产成人aa在线观看| 99久久无色码亚洲精品果冻| 亚洲最大成人av| av又黄又爽大尺度在线免费看 | 日韩高清综合在线| 国产精品女同一区二区软件| 长腿黑丝高跟| 国产精品三级大全| 国产亚洲91精品色在线| 久久久久久伊人网av| 26uuu在线亚洲综合色| 精品人妻视频免费看| 欧美一区二区精品小视频在线| 一级黄色大片毛片| 网址你懂的国产日韩在线| 国产免费男女视频| 国产在线男女| av在线亚洲专区| 在线国产一区二区在线| 亚洲,欧美,日韩| 欧美激情国产日韩精品一区| 欧美人与善性xxx| 中文资源天堂在线| 一边摸一边抽搐一进一小说| 国产亚洲精品久久久com| 精品免费久久久久久久清纯| 噜噜噜噜噜久久久久久91| 九色成人免费人妻av| 国产成人福利小说| 美女黄网站色视频| 久久精品国产清高在天天线| 日韩中字成人| 欧美日韩在线观看h| 久久欧美精品欧美久久欧美| 日本黄大片高清| 在线a可以看的网站| 亚洲最大成人中文| 中文资源天堂在线| 亚洲三级黄色毛片| 丝袜喷水一区| 国产精品一二三区在线看| 18禁在线无遮挡免费观看视频| 日日摸夜夜添夜夜爱| av在线播放精品| 麻豆国产97在线/欧美| 免费电影在线观看免费观看| 亚洲欧美精品专区久久| 3wmmmm亚洲av在线观看| www日本黄色视频网| 久久精品久久久久久噜噜老黄 | 亚洲在久久综合| 久久久久国产网址| 夫妻性生交免费视频一级片| 国产成人午夜福利电影在线观看| 在线国产一区二区在线| 亚洲av中文字字幕乱码综合| 又黄又爽又刺激的免费视频.| 久久精品国产亚洲av天美| 成熟少妇高潮喷水视频| 一级av片app| 国产v大片淫在线免费观看| 中文字幕制服av| 国产精品久久久久久久久免| 亚洲欧美精品综合久久99| 桃色一区二区三区在线观看| 午夜激情欧美在线| 在线免费十八禁| 国产片特级美女逼逼视频| 男人和女人高潮做爰伦理| 青春草国产在线视频 | 非洲黑人性xxxx精品又粗又长| 插阴视频在线观看视频| 超碰av人人做人人爽久久| 亚洲精品粉嫩美女一区| 亚洲精品粉嫩美女一区| 欧美精品国产亚洲| 狂野欧美激情性xxxx在线观看| 女人十人毛片免费观看3o分钟| 日日摸夜夜添夜夜爱| 欧美色视频一区免费| 日韩 亚洲 欧美在线| 人妻系列 视频| 级片在线观看| 97超碰精品成人国产| 日韩国内少妇激情av| 中文字幕久久专区| 国产成人午夜福利电影在线观看| 亚洲无线观看免费| 婷婷亚洲欧美| 国产亚洲精品久久久久久毛片| 国产一级毛片七仙女欲春2| 国产精品乱码一区二三区的特点| 午夜视频国产福利| 欧美变态另类bdsm刘玥| 国产探花在线观看一区二区| 亚洲在线自拍视频| 观看免费一级毛片| 啦啦啦观看免费观看视频高清| 欧美日韩国产亚洲二区| 国产高清不卡午夜福利| 成人毛片60女人毛片免费| 欧美日本亚洲视频在线播放| 在线观看一区二区三区| 亚洲自偷自拍三级| 嫩草影院入口| 免费大片18禁| 欧美成人a在线观看| 久久久久网色| 久久草成人影院| 国产视频内射| 一本精品99久久精品77| 欧洲精品卡2卡3卡4卡5卡区| 男女那种视频在线观看| 又黄又爽又刺激的免费视频.| 精品久久久久久久久久免费视频| 久久久精品大字幕| 免费观看精品视频网站| 中文字幕久久专区| 欧美日本视频| 高清午夜精品一区二区三区 | 91午夜精品亚洲一区二区三区| 夜夜看夜夜爽夜夜摸| 91午夜精品亚洲一区二区三区| 26uuu在线亚洲综合色| 久久久精品欧美日韩精品| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 简卡轻食公司| 欧美成人免费av一区二区三区| 99久久人妻综合| 小蜜桃在线观看免费完整版高清| 天天一区二区日本电影三级| 在线观看免费视频日本深夜| 男的添女的下面高潮视频| 国产精品福利在线免费观看| 国产不卡一卡二| 国产高清视频在线观看网站| 国产精品精品国产色婷婷| 亚洲va在线va天堂va国产| 韩国av在线不卡| 欧美高清性xxxxhd video| 国内精品久久久久精免费| 国产亚洲精品av在线| 三级经典国产精品| 在线观看免费视频日本深夜| 日本免费一区二区三区高清不卡| 亚洲精品日韩在线中文字幕 | 久久久色成人| 日韩一区二区三区影片| 久久精品国产99精品国产亚洲性色| 一卡2卡三卡四卡精品乱码亚洲| 日韩大尺度精品在线看网址| 中国美白少妇内射xxxbb| 寂寞人妻少妇视频99o| 在线免费观看的www视频| 噜噜噜噜噜久久久久久91| 午夜福利在线观看吧| 欧美色视频一区免费| 一卡2卡三卡四卡精品乱码亚洲| 免费人成在线观看视频色| 国产午夜精品论理片| 小说图片视频综合网站| 国产熟女欧美一区二区| 国产乱人视频| 久久人人精品亚洲av| 搡女人真爽免费视频火全软件| 国产精品国产三级国产av玫瑰| 欧美性猛交╳xxx乱大交人| 亚洲七黄色美女视频| 国产精品美女特级片免费视频播放器| 高清毛片免费看| av专区在线播放| 日本成人三级电影网站| 亚洲自拍偷在线| 男女边吃奶边做爰视频| 18禁在线播放成人免费| 久久国内精品自在自线图片| 久久精品久久久久久久性| 亚洲av不卡在线观看| 黄色欧美视频在线观看| 久久久久久久久久成人| 国产极品精品免费视频能看的| 蜜臀久久99精品久久宅男| 国产一区二区激情短视频| 亚洲最大成人手机在线| 精品人妻一区二区三区麻豆| 久久久久国产网址| 亚洲成人av在线免费| 1000部很黄的大片| 91av网一区二区| 国产一区亚洲一区在线观看| 欧美性猛交黑人性爽| 精品免费久久久久久久清纯| 日韩 亚洲 欧美在线| 大又大粗又爽又黄少妇毛片口| 一个人看的www免费观看视频| 亚洲真实伦在线观看| 国产精品,欧美在线| 亚洲aⅴ乱码一区二区在线播放| 日本黄色片子视频| 亚洲,欧美,日韩| 欧美不卡视频在线免费观看| 国产精品一区二区性色av| 国产高清激情床上av| 深爱激情五月婷婷| 国产精华一区二区三区| 国产精品一二三区在线看| 老女人水多毛片| 内射极品少妇av片p| 精品熟女少妇av免费看| 国产一区二区三区av在线 | 免费大片18禁| 亚洲欧美日韩高清专用| 晚上一个人看的免费电影| 日本一本二区三区精品| 在线观看66精品国产| 蜜臀久久99精品久久宅男| 波多野结衣高清作品| a级毛片免费高清观看在线播放| 一本一本综合久久| 国产v大片淫在线免费观看| 99热精品在线国产| 久久久久久九九精品二区国产| 青春草亚洲视频在线观看| 51国产日韩欧美| 中文在线观看免费www的网站| 午夜免费男女啪啪视频观看| av天堂中文字幕网| 日产精品乱码卡一卡2卡三| 亚洲18禁久久av| av在线天堂中文字幕| 最近手机中文字幕大全| 免费看av在线观看网站| 日本av手机在线免费观看| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式| 亚洲在久久综合| 久久99精品国语久久久| 一级毛片久久久久久久久女| 老女人水多毛片| 人妻夜夜爽99麻豆av| av又黄又爽大尺度在线免费看 | 亚洲经典国产精华液单| 国产精品久久久久久亚洲av鲁大| 国产私拍福利视频在线观看| 国产黄片视频在线免费观看| 国产淫片久久久久久久久| 九九爱精品视频在线观看| 亚洲成av人片在线播放无| 亚洲电影在线观看av| 亚洲成人中文字幕在线播放| 精品久久久久久久久av| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 久久鲁丝午夜福利片| 五月玫瑰六月丁香| 男人和女人高潮做爰伦理| 亚洲国产精品久久男人天堂| 日韩欧美在线乱码| 不卡一级毛片| 久久久久久久久大av| 黄色视频,在线免费观看| 给我免费播放毛片高清在线观看| 欧美极品一区二区三区四区| av天堂中文字幕网| 爱豆传媒免费全集在线观看| 简卡轻食公司| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 嫩草影院精品99| 久久99热6这里只有精品| 国产亚洲精品av在线| 少妇被粗大猛烈的视频| 国产黄片视频在线免费观看| 久久精品久久久久久久性| 亚洲av中文字字幕乱码综合| 色哟哟哟哟哟哟| 日韩欧美三级三区| 高清毛片免费观看视频网站| 久久精品国产鲁丝片午夜精品| 熟女电影av网| 男女下面进入的视频免费午夜| 成人午夜高清在线视频| 国产精品久久久久久久久免| 国产精品人妻久久久久久| 在线观看免费视频日本深夜| 久久亚洲精品不卡| АⅤ资源中文在线天堂| 18禁在线无遮挡免费观看视频| 校园春色视频在线观看| 久久久欧美国产精品| 国产极品精品免费视频能看的| 亚洲精品粉嫩美女一区| 高清毛片免费看| 国产乱人视频| 美女xxoo啪啪120秒动态图| 免费黄网站久久成人精品| 色哟哟哟哟哟哟| 99九九线精品视频在线观看视频| 天堂影院成人在线观看| 成人毛片a级毛片在线播放| av视频在线观看入口| 美女被艹到高潮喷水动态| www.色视频.com| 色视频www国产| 国产亚洲精品久久久久久毛片| 久久久精品欧美日韩精品| 麻豆精品久久久久久蜜桃| 成人高潮视频无遮挡免费网站| 午夜a级毛片| 午夜福利视频1000在线观看| 亚洲不卡免费看| 精品久久久久久久久久久久久| 特级一级黄色大片| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 国产精品女同一区二区软件| 国产av一区在线观看免费| 免费观看在线日韩| 亚洲内射少妇av| 久久精品久久久久久久性| 日韩,欧美,国产一区二区三区 | 两个人的视频大全免费| 嫩草影院新地址| 成人综合一区亚洲| 12—13女人毛片做爰片一| 黄色视频,在线免费观看| 国产亚洲欧美98| 国产黄片美女视频| 99热只有精品国产| 人妻少妇偷人精品九色| 中国美白少妇内射xxxbb| 久久亚洲精品不卡| 亚洲av免费高清在线观看| 天堂av国产一区二区熟女人妻| 欧美日本亚洲视频在线播放| 色播亚洲综合网| 亚洲天堂国产精品一区在线| 国产爱豆传媒在线观看| 黄色视频,在线免费观看| 亚洲av不卡在线观看| av.在线天堂| 噜噜噜噜噜久久久久久91| 国产91av在线免费观看| 日韩av不卡免费在线播放| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美精品综合久久99| 中国美白少妇内射xxxbb| 欧美+亚洲+日韩+国产| 亚洲av第一区精品v没综合| 蜜桃久久精品国产亚洲av| 亚洲成人av在线免费| 亚洲18禁久久av| 啦啦啦啦在线视频资源| 国产精品国产高清国产av| 99热网站在线观看| 插阴视频在线观看视频| 成年av动漫网址| 亚洲欧美精品自产自拍| 国产一区二区亚洲精品在线观看| 亚洲欧美日韩高清专用| av在线亚洲专区| 婷婷色综合大香蕉| 特大巨黑吊av在线直播| 一级黄色大片毛片| 有码 亚洲区| 欧美丝袜亚洲另类| 特级一级黄色大片| 国产成年人精品一区二区| 51国产日韩欧美| 97超碰精品成人国产| 日韩欧美精品v在线| 亚洲精品久久国产高清桃花| 精品久久久久久久人妻蜜臀av| 91aial.com中文字幕在线观看| 夜夜爽天天搞| 欧美变态另类bdsm刘玥| 久久亚洲精品不卡| 欧美xxxx性猛交bbbb| 成人永久免费在线观看视频| 草草在线视频免费看| 亚洲欧美清纯卡通| 国产色婷婷99| 欧美日本亚洲视频在线播放| 国产精品久久久久久av不卡| 亚洲天堂国产精品一区在线| 99久久久亚洲精品蜜臀av| 欧美一区二区国产精品久久精品| 日本黄色片子视频| 国产精品乱码一区二三区的特点| 男女下面进入的视频免费午夜| 禁无遮挡网站| 色5月婷婷丁香| 久久精品国产亚洲av香蕉五月| 麻豆成人午夜福利视频| 日日摸夜夜添夜夜添av毛片| 欧美bdsm另类| 国内少妇人妻偷人精品xxx网站| 日本黄色视频三级网站网址| 精品久久国产蜜桃| 91精品国产九色| av.在线天堂| 舔av片在线| 别揉我奶头 嗯啊视频| 91aial.com中文字幕在线观看| 久久午夜亚洲精品久久| 国产精品1区2区在线观看.| 舔av片在线| 日韩视频在线欧美| 国产伦精品一区二区三区视频9| 国产美女午夜福利| 亚洲一区二区三区色噜噜| 日韩成人av中文字幕在线观看| 久久久久久久亚洲中文字幕| 色噜噜av男人的天堂激情| 熟女人妻精品中文字幕| 一个人看的www免费观看视频| 青春草亚洲视频在线观看| 联通29元200g的流量卡| 长腿黑丝高跟| 国产色爽女视频免费观看| 久久精品久久久久久噜噜老黄 | 国产高清有码在线观看视频| 小蜜桃在线观看免费完整版高清| 国产av不卡久久| 亚洲在线观看片| 热99在线观看视频| 青春草国产在线视频 | 国内揄拍国产精品人妻在线| 国产精品美女特级片免费视频播放器| 三级经典国产精品| 91久久精品国产一区二区三区| 日韩视频在线欧美| 国产私拍福利视频在线观看| 欧美高清性xxxxhd video| 久久精品人妻少妇| 午夜福利在线观看吧| 欧美日本视频| av免费在线看不卡| 晚上一个人看的免费电影| 在线观看av片永久免费下载| 国产成人a∨麻豆精品| 日日摸夜夜添夜夜添av毛片| 亚洲国产欧洲综合997久久,| 九九爱精品视频在线观看| 日产精品乱码卡一卡2卡三| 最后的刺客免费高清国语| 日韩,欧美,国产一区二区三区 | 中文精品一卡2卡3卡4更新| 亚洲国产精品国产精品| av国产免费在线观看| 最近2019中文字幕mv第一页| 99视频精品全部免费 在线| 免费观看人在逋| 国产成人一区二区在线| 九九热线精品视视频播放| 99久久精品一区二区三区| 亚洲精品色激情综合| 成年av动漫网址| 亚洲精品日韩av片在线观看| 久久久国产成人免费| 国产精品.久久久| 观看美女的网站| 日韩成人伦理影院| 成年女人永久免费观看视频| 99久久成人亚洲精品观看| 亚洲欧美精品专区久久| 精品无人区乱码1区二区| 亚洲中文字幕日韩| 一区二区三区四区激情视频 | 26uuu在线亚洲综合色| 韩国av在线不卡| 久久精品91蜜桃| а√天堂www在线а√下载| 99久久精品一区二区三区| 爱豆传媒免费全集在线观看| 精品无人区乱码1区二区| 悠悠久久av| 毛片女人毛片| 成人亚洲欧美一区二区av| av视频在线观看入口| 内射极品少妇av片p| 黄色视频,在线免费观看| 成人综合一区亚洲| 在线免费观看的www视频| 欧美一级a爱片免费观看看| 亚洲一级一片aⅴ在线观看| 久久久精品大字幕| 哪里可以看免费的av片| 成人午夜精彩视频在线观看| 变态另类丝袜制服| 啦啦啦韩国在线观看视频| 成年女人永久免费观看视频| 免费一级毛片在线播放高清视频| 亚洲国产欧美人成| 欧美最黄视频在线播放免费| 国产高清视频在线观看网站| 亚洲人成网站高清观看| 丰满人妻一区二区三区视频av| 亚洲美女搞黄在线观看| 一区福利在线观看| 色尼玛亚洲综合影院| 亚洲综合色惰| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 最近最新中文字幕大全电影3| 色综合亚洲欧美另类图片| 亚洲成av人片在线播放无| 爱豆传媒免费全集在线观看| 中文在线观看免费www的网站| 精品一区二区免费观看| 日本三级黄在线观看| 精品国产三级普通话版| 中出人妻视频一区二区| 亚洲国产精品sss在线观看| 亚洲成a人片在线一区二区| 男女啪啪激烈高潮av片| 国产精品一区www在线观看| 国产成人精品久久久久久| 亚洲四区av| 精品一区二区三区人妻视频| 久久久久网色| 久久热精品热| 99热这里只有是精品在线观看| 亚洲av中文av极速乱| 国产视频内射| 国产 一区精品| 国产精品伦人一区二区| 深夜a级毛片| 内地一区二区视频在线| 国内揄拍国产精品人妻在线| 日本在线视频免费播放| 国产私拍福利视频在线观看| 中文在线观看免费www的网站| 草草在线视频免费看| 日本三级黄在线观看| 久久久久久伊人网av| 亚洲国产精品国产精品| 99久久中文字幕三级久久日本| 国产在线精品亚洲第一网站| 变态另类丝袜制服| 看黄色毛片网站| 日韩成人av中文字幕在线观看| 色综合色国产| 男插女下体视频免费在线播放| 免费观看在线日韩| 成人av在线播放网站| 自拍偷自拍亚洲精品老妇| 国产精品蜜桃在线观看 | 成年女人永久免费观看视频| 少妇的逼水好多| 免费看美女性在线毛片视频| 成人永久免费在线观看视频|