A Complex of Silver(Ⅰ)with N-O Chelating Agent and Phosphine Ligand:Design,Structures,and Biological Activity

2022-08-09 03:50WANGRenShuFENGJingWANGYuanLanOUYANGChenFeng
无机化学学报 2022年8期

WANG Ren-ShuFENG JingWANG Yuan-LanOUYANG Chen-Feng

(1School of Chemistry and Materials Engineering,Liupanshui Normal University,Liupanshui,Guizhou 553004,China)(2Guizhou Key Laboratory of Coal Clean Utilization,Liupanshu,Guizhou 553004,China)(3E.N.T.department,General Hospital of Zhongshan,Liupanshu,Guizhou 553000,China)

Abstract:A new Ag(Ⅰ) complex,[Ag(Qina)(Tpp)2]·1.5H2O(Qina=2-quinolinecarboxylate,Tpp=triphenylphosphine),was synthesized and successfully obtained as a single crystal.The structure of the complex was characterized by single-crystal X-ray diffraction,IR,NMR,and powder X-ray diffraction.The results show that the complex belongs to monoclinic crystal system,C2/c space group,and the unit cell parameters are a=3.195 90(13)nm,b=1.210 96(4)nm,c=2.319 74(7)nm,β=102.166(4)°,V=8.776 0(5)nm3,and Z=8.The supramolecular structures and weak intermolecular force of the complex were deduced by Hirshfeld surface analysis.Meanwhile,the biological activities of the complex were tested by evaluating DNA-binding efficacy,antibacterial activity,and in vitro toxic activity of cancer cells.CCDC:2142361.

Keywords:Ag(Ⅰ)complex;Hirshfeld surface analysis;crystal structures;biological activity

0 Introduction

Since ancient times,ancestors have used silverware for food preservation,silver nitrate for disinfection,and then silver compounds went on the market in the last century as an anti-inflammatory drug for burns represented by silver sulfadiazine[1-2].Silver,as a metal with biological activity and medicinal value,has been flowing throughout the long history of mankind.Because microorganisms have almost no drug resistance to silver compounds,silver compounds can coexist with a variety of antibiotics and occupy a place[3].Firstly,silver ions can penetrate the surface of the cell wall and strongly bind to sulfur-containing enzymes,or they can remove electrons from cell components,destroy the integrity and permeability of the cell wall,and finally interfere with cell metabolism and kill cells[4-6].Therefore,silver can play a role in the pharmaceutical neighborhood.Secondly,due to the existence of a large number of chloride ions in the human body environment,after the activation of silver ions,these chloride ions will combine with them and then convert into AgCl without toxicity[7].As a traditional safe metal,silver is undoubtedly favored by developers of nonplatinum cancer chemotherapeutic drugs.Platinumbased chemotherapy drugs have been criticized in clinical applications because of their high toxicity and side effects[8].

Quinoline compounds have excellent biological activities.As a result,many compounds in the quinoline ring family,such as plasmocin,bulaquine,and tafenoquine,have been used as medicines because of their antimalarial properties[9].From the perspective of structural chemistry,it is found that the 2-quinolinecarboxylate(Qina)ligand molecule has an electron-rich anionic COO-group and a neutral N atom that can provide a lone electron pair.It is a classical nitrogenoxygen chelating ligand,so it is often selected as a metal chelating agent for metal wet selection and analytical chemistry[10-11].More importantly,the Qina ligand also has a wide range of biochemical activity significance[12].The Qina ligand in the quinoline ring family has been reported to be used as a nucleic acid matrix to assist in the DNA laser desorption process[13].In addition,Qina is a metabolite that eventually exists in urine after humans receive tryptophan through the kynurenine metabolic pathway[14-16].Studies have pointed out that the tryptophan kynurenine metabolic pathway can produce Qina because it is closely related to the role of intestinal flora[17-18].At the same time,there is also research evidence that the stem extract ofEphedra Far Easthas a strong inhibitory effect onClostridium difficileandClostridium perfringens[19].These two bacteria are considered to be related to intestinal diseases in animals[20-22].The reason whyEphedrastem extract expresses antibacterial activity is that the extract contains the key substance Qina,which plays a vital role in the antibacterial process[19].More specifically,studies have shown that Qina can also affect the expression of oncogene p53 at the gene and protein molecule level,and then express its inhibitory effect on LS180 human colon adenocarcinoma cells[23].More pragmatically,thiostrepton,an antibiotic originally used as a feed additive,also has Qina functional group in its side chains that enables its molecules to express anticancer and antibacterial activities[24-25].In recent years,the improved synthesis of this antibiotic has often been carried out around the Qina group contained in its side chain[26-27].Correspondingly,there are also reports on the anticancer activity of metal complexes containing Qina ligand[28-29].Besides,related complexes are also used in luminescence[30],catalysis,functional group protection[31],and so on.

In the synthesis of small molecular complexes with mixed ligands,synthesizers usually select Lewis base bidentate chelating agents with nitrogen-containing heterocycles as the auxiliary ligand.However,choosing bidentate or multidentate ligands as the main ligand may cause excessive steric hindrance to the secondary coordination of the auxiliary ligand,resulting in difficulties in the synthesis of target compounds or ligand dissociation[32-33].Recent reports have shown that the spatial obstacles and restrictions between primary and secondary ligands could be avoided by the use of phosphorus-containing monodentate auxiliary ligands in synthesis[34].Triphenylphosphine(TPP)is often used as a reducing agent,and it participates in various organic reactions after being converted to triphenylphosphine oxide in the synthesis process[35].TPP is also a good candidate for monodentate ligands with multiple-aromatic groups,which have good lipophilicity and can effectively penetrate the mitochondrial membrane[36].It has been reported that the use of TPP in the synthesis of silver complexes can contribute to the improvement of the stability and light resistance of the silver complex[37].In addition,the fluorescence activity of Ag-TPP type complexes has also been studied[38].

PI3K kinase is a key signal molecule for cancer cell growth and one of the most active signal pathways in cancer cells[39-40].In addition,PI3K is closely related to the Warburg effect of glycolysis in cancer cells,that is,cancer cells use Warburg metabolism to maintain the activity of the PI3K signaling pathway,to ensure continuous growth and division of cancer cells[41-42].

We previously reported a case of Ag-Qina type linear coordination polymer with good DNA binding efficiency and moderate cancer cytotoxic activity[43].Based on the excellent biological characteristics of Ag,Qina,and TPP,we are interested in the preparation and activity evaluation of small molecule complex with these three components and hope to use the combination of the three components to make them synergistically express more biological activities in the complementary process.

1 Experimental

1.1 Material and instruments

All reagents and solvents were commercially purchased and used without any further treatment.The elemental composition(including C,H,and N)of these samples was determined by a model Elementar Vario Macro cube.The powder X-ray diffraction(PXRD)experiment was completed by using a Shimadzu 6100 diffractometer,which employed CuKα(0.154 056 nm)radiation and was operated at 40 kV and 40 mA in a 2θrange of 5°-80°.The1H NMR spectra of the complex in DMSO-d6solutions were recorded on a Bruker AVANCE Ⅲ 600 spectrometer.The FT-IR spectra were recorded on a Bruker Equinox 55 FT-IR spectrophotometer.In addition,the emission spectra were obtained on a Perkin-Elmer LS55 fluorescence Spectrofluorometer.

1.2 Synthesis of the complex

0.3 mmol of silver acetate and 0.6 mmol of TPP were mixed in a small volume of toluene.The reaction system was heated,and the reaction was stopped when the system solution suddenly became clear.Toluene was removed by vacuum drying with an oil pump,and white bulk Ag-TPP intermediate was collected.The intermediate was stirred with the mixed solvent of ethanol and distilled water(1∶1,V/V)containing Qina at room temperature.After carefully adjusting the pH value to 5.0,the reaction continued for at least 3 h.A large amount of white precipitate that appeared during the reaction was collected by filtration and washed with ethanol.The dried precipitate was recrystallized with methanol and acetonitrile(1∶1,V/V).The product solution was kept in a ventilated darkroom for about a week,and the complex was precipitated out of the solution in the form of colorless massive crystals.The yield wasca.65.5%based on Ag.1H NMR(600 MHz,DMSO-d6):δ8.51(d,J=8.4 Hz,1H),8.34(d,J=8.4 Hz,1H),8.00(d,J=7.72 Hz,1H),7.63(d,J=8.4 Hz,1H),7.55(t,J=7.08 Hz,1H),7.45-7.36(m,7H),7.33-7.14(m,24H)(Fig.S1a,Supporting information).13C NMR(101 MHz,DMSO-d6):δ166.84,155.82,144.87,138.18,133.25(d,J=16.9 Hz),132.49(d,J=23.3 Hz),130.51,130.05,129.10(d,J=9.4 Hz),128.64,128.48,128.19,127.63,122.21(Fig.S1b).Elemental analysis Calcd.for C46H39AgNO3.5P2(%):C 66.43;H 4.72;N 1.68;Found(%):C 65.68;H 4.50;N 1.67.Selected FT-IR bands(KBr,cm-1):3 051ν(Ar—H),1 606ν(C=C),1 585ν(C=C),1 502ν(C=C),1 477ν(C=C),1 606ν(O—C—O)asym,1 351ν(O—C—O)sym,1 428ν(C=N),781ν(phenyl rings),722ν(phenyl rings),688ν(phenyl rings),490ν(C—P)(Fig.S2).

1.3 Single-crystal X-ray crystallography

An Agilent Gemini E single-crystal X-ray diffractometer(MoKα,λ=0.071 073 nm)was applied for polymorphism analysis.Data indexing,integration,and absorption correction were carried out using the CrysAlisPro software package.Besides,the structure solution and refinement were done with the support of full-matrix least-squares procedures using the SHELXL-2014 software package through the OLEX2 suite[44-46].All non-hydrogen atoms were located in successive difference Fourier syntheses and refined by full-matrix least-squares onF2.The locations of all hydrogen atoms are theoretically correct.

CCDC:2142361.

1.4 Hirshfeld surface calculations

The CIF file of the complex molecule analysis was studied by the Crystal Explorer software[47].All the obtained Hirshfeld surfaces were output with high resolution.Moreover,the complex has been analyzed by the coordinate range of 2D fingerprint plots.However,the calculated coordinate range of complex molecules was 0.1 to 0.28 nm.

1.5 Fluorescence spectrum

As a well-known fluorescent molecule,ethidium bromide(EtBr)can embed into the DNA base,which shows fluorescence enhancement after fixation.The effects of complex differ in its concentration because of the competitive binding with calf thymus DNA(ct-DNA)-EtBr(cDNA=5 μmol·L-1,cEtBr=1 μmol·L-1).In this work,a buffer solution with a pH value of 7.4 was prepared,containing 10 nmol·L-1NaCl and 50 mmol·L-1Tris-HCl.Before the fluorescence test,the samples were incubated in the buffer solution for 4 h at 20℃.The excitation wavelength was set as 602.5 nm,while the emission scale was set in a range of 550-750 nm.

1.6 Antibacterial efficacy tests

The cultivation of the bacterial liquid and plate smearing inoculation were completed according to the conventional procedures.Then,sterilized filter paper with a diameter of about 1 cm was used and immersed in a sterile aqueous solution of complexes with concentration gradients of 0.5,1.0,1.5,and 2.0 mmol·L-1.After 20 min,the filter papers were placed on central media plates containingEscherichia coli(E.coli),Staphylococcus aureus(S.aureus),andParatyphoid salmonella(P.salmonella),and the concentration of each triplicate set.Finally,all plates were cultured in biochemical incubators at 37℃for 24 h to complete the record of bacterial inhibition rings.In addition,the inhibitory effects of the complex with different concentrations on bacterial growth could be compared by drawing the bacterial growth inhibition curves.

A bacterial suspension(1 mL)and 10 mL of liquid culture medium were added in parallel to five sterile test tubes.Then,each test tube was sequentially added with 1 mL of the acetonitrile solution of the complex with the concentration of 0,0.5,1.0,1.5,and 2.0 mmol·L-1,and placed in a biochemical incubator at 37℃for constant temperature incubation for 24 h.The culture medium was transferred to the centrifuge tube and centrifuged at 1 000 r·min-1for 10 min.After centrifugation,the supernatant of the bacterial solution was collected.The absorbance of the supernatant at a wavelength of 600 nm was measured using an ultraviolet spectrophotometer.The penicillin sodium solution with the same concentration gradient as the tested complex was set as the control group to compare the antibacterial efficacy of the Ag complex more objectively.

1.7 Cytotoxicity assay

The growth inhibitory effect of the metal complex on the HeLa and MCF-7 cells was evaluated by the MTT assay.Briefly,the cells firstly were cultured in RPMI1640 containing 10%fetal bovine serum in an incubator at 37 ℃ under a CO2atmosphere(φ=5%).Then each well of the 96-well culture plate was precisely inoculated with 100 μL suspension containing 3×104cells and then incubated for 24 h.The solution of the complex in each well was diluted by the addition of 150 μL RPMI 1640 solution.Besides,dimethyl sulfoxide(DMSO)with a final concentration of not higher than 0.1%was added with three replication settings.After incubation for 72 h,the old solution in each well was washed using 200 μL of saline.Then,100 μL of MTT solution(0.5 mg·mL-1)diluted 10 times with RPIM1640 was introduced into each well and then the cells were cultivated for another 3 h.Then,the media containing MTT was removed.Finally,the obtained formazan crystals were dissolved in 150 μL of DMSO at room temperature for 0.5 h.The absorbance performance was evaluated at 492 and 630 nm,respectively.The detailed effects of the cisplatin and carboplatin on the cell growth were investigated as a control team in another experiment.

1.8 Molecular docking

The PI3K enzyme protein crystal file number 3csf was downloaded from the PDB database.The opensource PyRx 0.8 nested AutoDock Vina program was utilized to simulate the molecular docking of the protein and the target compound[48].The visual projection of the 2D effect of the docking result was computed and output by Discovery Studio.3D results of the docking were displayed and output using the Pymol program.

2 Results and discussion

2.1 Crystal structure determination

Single-crystal X-ray analysis shows that the complex crystallizes successfully in the monoclinic system with theC2/cspace group.The data of a colorless single crystal of the complex was collected in aθrange of 2.58°-29.324°with a total of 20250 reflections.There might be a bug with CheckCIF because there is a hydrogen bond O2W—H2W…O3W between the water molecules O2W and O3W.To make the ciffile no warning,the influence of water molecules was avoided in the process of subsequent crystal structure analysis and calculation.As shown in Fig.S3,the peak positions of the simulated and experimental PXRD patterns matched each other.As shown in the thermogravimetric analysis curve in Fig.S4,the complex exhibited three steps of weight loss.The first weight loss corresponding to the removal of lattice water was 3.25%(Calcd.3.46%)from 30 to 120℃,while the second weight loss corresponding to the removal of ligands was 79.88%(Calcd.83.79%)from 120 to 300℃.The final residue is the metal oxide(Calcd.27.87%).The crystal data and structural refinement of the complex are summarized in Table 1.The selected bond lengths and bond angles are listed in Table 2.

Table 1 Crystal data and structure refinements for the complex

Table 2 Selected bond lengths(nm)and angles(°)for the complex

The complex presents a tetracoordinate mode.Each silver(Ⅰ)ion is coordinated by two phosphorus atoms from two TPP ligands and chelated by the oxygen and nitrogen atoms of the Qina ligand.The basic structure of the complex is exhibited in Fig.1.

Fig.1 ORTEP view with 30%probability ellipsoids of the complex

2.2 Hirshfeld surface analysis and structural details of the complex

From thednormanalysis results in Fig.2a,it can be inferred that the red dot region indicates the existence of someπinteractions and hydrogen bond interaction regions between the complex molecules.Fig.2b shows that the results of the shape index can be mutually confirmed withdnormanalysis.In the area where red dots appear indnorm,the shape index will produce surface depressions at the corresponding positions,and it will still be red[49].The shape index also shows that there are strong O…H andπinteractions between molecules.Fig.2c shows the smooth curvedness surface of molecules,which indicates that the intermolecular interactions tightly fit together.

Fig.2 Hirshfeld surface of the complex:dnorm(a),shape index(b),and curvedness(c)

π-stacking is widely present in biological macromolecular systems and often plays a key role in the stable stacking of DNA and RNA bases,the maintenance of polysaccharides and protein-specific structures,and the recognition of protein-ligand interactions[50-51].In addition,supramolecular interactions such as theπ-stacking network and hydrogen bond network between complex molecules often become potential action points for drug molecules to exert force on biological targets[52].

The dissection ofdnorm(Fig.3)shows that there is an edge-to-face type weak C—H…πinteraction between the complex molecules along the crystal axisc.The H15 subordinate to C15 on the aromatic ring of the TPP ligand interacts with the aromatic ring Cg(C40-C45)of another adjacent complex TPP ligand at a distance(H…Cg)of 0.321 3 nm.The closest distance between the H15 in the C—H…πinteraction and theatom of the aromatic ring is 0.279 6 nm,the sec-ondary distance to theatom of the aromatic ring is 0.302 5 nm,and the dihedral angle between the plane C42-H15-C43 and the plane C42-centroid-C43 is 97.646°.The details of the geometric parameters of the weak interaction are shown in Fig.4.This C—H…πinteraction can meet the parameter range defined and required in published papers[53-54].

Fig.3 Details of C—H…π interaction in dnormanalysis

It is through the alternation of intermolecular C—H…πinteractions that the 1D supramolecular configuration of the complex can be formed.The green dotted line in Fig.5 shows this C—H…πinteraction,and the 1D supramolecular chain structure of the complex is also shown in Fig.5.

Fig.5 One-dimensional supramolecular chain structure of the complex

Based on the 1D chain of supramolecular structure,there is a pair of C—H…O type hydrogen bonds between the O1 of the carboxyl functional group of the Qina ligand and the CH group of the side chain.Each supramolecular chain can extend the 2D supramolecular surfaces on thebcplane through this effect.The bond parameters of the hydrogen bond pair are the same,the D…A distance is 0.261 3 nm,and∠DHA=139.078°,which conforms to the range of hydrogen bond parameters of this type[55].The green dotted line in Fig.6 shows the paired hydrogen bonding between the supramolecular chains and the 2D supramolecular surface structure of the complex.

Fig.6 Two-dimensional supramolecular planar structure of the complex

The shape index analysis in Fig.7 shows that a large red concave surface appears above the quinoline fused ring of the Qina ligand,indicating that the complex molecules will haveπ…πstacking interactions here.The curvedness analysis also shows that the molecule surface above the quinoline fused ring is smooth and fitted to the adjacent molecule surface,which further illustrates that theπ…πstacking interactions inferred from the results of the shape index analysis are established[56].

Fig.7 π…π stacking interactions between two Qina molecules:shape index(a)and curvedness(b)

Along the crystal axisa,the 2D supramolecular planes depend on the equidistant doubleπ…πinteractions to maintain each other.Theseπ-action pairs occur between the fused rings of quinoline in the Qina ligand's structure,and finally construct a complete crystal packing structure of the complex molecular system(Fig.8).This can be confirmed by the conclusion of the Hirshfeld surface analysis.The intermolecularπ-interaction surfaces match each other well,and the center-to-center distance is 0.377 6 nm,which is within the regular numerical range[57].

Fig.8 Three-dimensional supramolecular structure of the complex

Two-dimensional fingerprint plots can show in detail the quantitative relationship between the atomic charges between molecules[58].As shown in Fig.9,the interaction ratio of H…C/C…H is as high as 24.1%.The projected shape is a symmetrical sharp flying wing with a bright area in the center,which confirms the existence of C—H…πinteractions between molecules.In addition,the H…O/O…H interactions account for 10.1%.The two wings of the projection map are symmetrical and very sharp,indicating that there is a strong O…H hydrogen bond between the molecules.The C…C interactions in the central area of the shadow map a regular and symmetrical pattern with a light spot in the center,indicating that there areπ…πinteractions between the molecules.

Fig.9 Two-dimensional fingerprint plots of the complex:(a)H…C/C…H interactions,(b)H…O/O…H interactions,(c)C…C interactions

The proportion of charge interaction between all atoms is illustrated in Table 3.The sum of the proportion of charge interaction between all atoms is 100%,indicating that the Hirshfeld surface analysis of the crystal is complete and strict.The analysis of the 2D fingerprint plots can be consistent with the analysis result of the crystal structure.

Table 3 Hirshfeld surface calculations for the complex

2.3 Fluorescence properties

According to the fluorescence spectra of the complex-EB-DNA system(Fig.10),in the EB-DNA system,the sensitivity reached its maximum when the fluorescence wavelength was 602.5 nm.After the introduction of the complex,the fluorescence of the EB-DNA system was quenched.With the gradual increase of the complex concentration,the fluorescence quenching became more obvious.The intensity was quenched from 77.63 to 34.23,with a difference of 43.40.

Fig.10 Fluorescence spectra for the binding of the complex to EB-DNA

The Stern-Volmer equation[59]I0/I=1+rKscan be used to quantitatively evaluate the competitive binding ability of the complex with the EB-DNA system,whereI0andIrepresent the fluorescence intensities with and without the complex,respectively;Ksrefers to the quenching constants andrrepresents the ratio of the gradient concentration of the complex to that of DNA(ccomplex/cDNA).In addition,regardingI0/Iandrscatter plots,the slope of the linear fitting can be obtained by theKsvalue.Additionally,the linear fitting diagram is also reflected in Fig.11.TheKsvalue of the complex was 0.326.In this case,the complex has displaced the EtBr molecules in the EtBr-DNA system to a certain extent.In addition,fluorescence studies show that the main mode of the complex interacting with DNA is intercalation[60].

Fig.11 Stern-Volmer quenching plots of the complex in EB-DNA system

2.4 Antibacterial activity

The complex showed certain inhibitory activity againstS.aureus,and the area of the inhibitory ring increased with the increase of solution concentration.The comparative test shows that the antibacterial effect of the complex was not as good as that of penicillin sodium,as shown in Table 4.Fig.12 and 13 show the bacteriostatic zones of the complex and penicillin sodium againstS.aureusunder different concentration gradients,respectively.In addition,the ligands and complex cannot inhibitP.salmonellaandE.coli.

Table 4 Diameter of the antibacterial ring of tested compounds

Fig.13 Antibacterial activity of penicillin sodium on S.aureus

The inhibitory effects of the complex with different concentrations on bacterial growth can be compared by bacterial growth inhibition curves.According to the bacterial liquid system,the antibacterial action line of the test substance obtained by measuring the OD600(the optical density at 600 nm)values is shown in Fig.14.When the concentration of the complex reached 1.0 mmol·L-1,the inhibitory effect onS.aureuswas obvious at first,and then gradually leveled off with the increase of the concentration of the complex.Consistent with the results of the bacteriostatic ring experiment,the complex showed a mild inhibitory effect onS.aureus,but its activity was inferior to penicillin sodium.

Fig.14 OD600values of the complex and penicillin sodium

2.5 Cytotoxicity in vitro study

The complex was examined regarding thein vitrocytotoxic activity using an MTT assay.According to Table 5,it can be found that IC50 values were evaluated for two human cancer cells,and the lines of HepG2 and SKOV-3 are listed.The results clearly show that the complex exhibited a high level of anti-tumor activity as well as a level of toxicity.In particular,the cytotoxicity of the complex was much better than that of cisplatin and carboplatin.

Table 5 Cytotoxicity of tested compounds against tumor cells after 72 h incubation

2.6 Molecular docking studies

To further investigate the biological activity of the complex,we simulated the possible binding modes of the complex with PI3K kinase.Fig.15 exhibits the possible combination of thermodynamically stable complexes with PI3K,where the complex to PI3K with computed binding affinity was-40.184 kJ·mol-1,respectively.The interaction model between the complex with PI3K can be displayed in 2D mapping pictures(Fig.16).The aromatic ring groups of the auxiliary ligand can contribute to the multipleπ-binding sites for the complex and the enzyme.

Fig.15 Three-dimensional model of the docking between the complex and PI3K

Fig.16 Bonding interactions of the complex with PI3K

3 Conclusions

A complex[Ag(Qina)(Tpp)2]·1.5H2O was synthesized and characterized mainly by X-ray single-crystal diffraction.The properties of the complex in DNA binding were tested via fluorescence spectra,and the results indicated that the complex behaves with fairly strong binding affinities in DNA interaction.In addition,the complex can be predicted to have potential enzymatic activity by molecular docking.Throughin vitrocytotoxicity experiment,it was testified that the complex possessed a varied level of toxicity to inhibit HepG2,SKOV-3,and cell lines.In particular,the complex expressed higher cytotoxicity than cisplatin and carboplatin.Finally,the complex also showed certain antibacterial activity.It is suggested that the silver complex could be further exploited as a potential anticancer pharmaceutical product.

Acknowledgments:This work is financially supported by the Youth Talent Growth Project of the Department of Education of Guizhou Province(Grant No.KY[2020]123)and Guizhou Key Laboratory of Coal Clean Utilization(Qiankehe Platform Talents[2020]2001).

Supporting information is available at http://www.wjhxxb.cn