含苯并咪唑基配體鈷(Ⅱ)配合物的合成、晶體結(jié)構(gòu)和兒茶酚酶活性

張前　韓燕　焦元紅＊,(湖北理工學(xué)院化學(xué)與化工學(xué)院，黃石435003)(新鄉(xiāng)學(xué)院化學(xué)化工學(xué)院，新鄉(xiāng)453003)

含苯并咪唑基配體鈷(Ⅱ)配合物的合成、晶體結(jié)構(gòu)和兒茶酚酶活性

張前1韓燕2焦元紅＊,1
(1湖北理工學(xué)院化學(xué)與化工學(xué)院，黃石435003)
(2新鄉(xiāng)學(xué)院化學(xué)化工學(xué)院，新鄉(xiāng)453003)

合成和表征了一種鈷(Ⅱ)配合物，[CoLCl]2[CoLCl0.5(H2O)0.5]2ClO4·10H2O(1)(LH=(2-(二((1H-苯并咪唑-2-基)甲基)氨基)乙酸))。以3，5-二叔丁基兒茶酚(3,5-DTBC)為反應(yīng)底物，用紫外光譜測試了1的兒茶酚酶催化活性。研究結(jié)果表明：在配合物的晶體結(jié)構(gòu)中，不對稱單元的2個Co(Ⅱ)都形成變形的三角雙錐構(gòu)型。在pH=5～11范圍內(nèi)，配合物1對3，5-DTBC的氧化顯示了pH值依賴性，它的兒茶酚酶催化活性隨著溫度的升高而升高，并且其催化氧化3，5-DTBC的動力學(xué)符合米氏方程模型。

鈷(Ⅱ)配合物；晶體結(jié)構(gòu)；兒茶酚酶；催化活性

Catecholase(EC 1.10.3.1)is a type-3 copper enzyme usually encountered in plant tissues and in some insects and crustaceans.It catalyzes the oxidation of a broad range of o-diphenols to the corresponding o-quinones coupled with the reduction of oxygen to water[1].Catecholase is of economic importance due to its help to protect damaged plants against both bacterial and fungal disease.In the crystal structure of catecholase isolated from Ipomoea batatas(sweet potato),the active site contains a binuclear copper(Ⅱ)center,and each copper(Ⅱ)ion is coordinated by three imidazole nitrogen atoms from three histidine units[2]. It′s inspired numerous scientists to study the structure-activity relationships of catecholase. Benzimidazole could mimic the histidine residue owing to its same imidazole group,and recently many catecholase model complexes containing benzimidazole and imidazole ligands have been reported based on their active sites[3-8].The synthesis of artificial catecholase could provide a deeper insight into the structure, catalytic activity and mechanism of the natural catecholase,and it helps us to design more effective catecholase models.Herein,we report the synthesis and crystal structure of a cobalt(Ⅱ)complex(1) containing benzimidazole ligand.Furthermore,its catecholase catalytic activity was investigated using 3,5-di-tert-butylcatechol(3,5-DTBC)as the substrate.

1　Experimental

1.1Materials and physical measurements

All reagents were purchased from commercial companies and directly used unless stated otherwise. The melting point was determined with an XT4A micromelting point apparatus and was uncorrected. The pH value was adjusted by PHS-3C digital pH meter.UV-Vis spectra were measured on an Analytik jena Specord 210 spectrophotometer.The IR spectra were carried out on a Perkin-Elmer Spectrum BX FTIR instrument in tablets with potassium bromide. Elemental analyses were examined on a Perkin-Elmer 2400 instrument.Electrospray ionization mass spectra (ESI-MS)were performed on an Applied Biosystems API 2000 LC/MS/MS system.1H NMR spectra were recorded on a Varian Mercury 400 spectrometer at 400 MHz.

1.2Synthesis of the ligand LH

The two-step synthesis of the ligand LH is outlined in Scheme 1.The intermediate bis(benzimidazol-2-ylmethyl)amine(a)were prepared according to the literature method[9].The ligand LH was synthesized using a modified literature procedure[10].

Scheme 1 Synthesis of the ligand LH

25 mL of sodium hydroxide aqueous solution (0.02 mol·L-1)was added dropwise to a mixture of chloroacetic acid(0.95 g,10 mmol)and the intermediate a(2.77 g,10 mmol)in ethanol(100 mL).The resulting mixture was stirred rapidly and refluxed for 6 hours,then the solvent was evaporated to give the crude product,which was dissolved in water and filtered.The solution was acidified to pH 3～4 with hydrochloric acid(0.1 mol·L-1)and further filtered to obtain the white product.Yield:2.18 g(65%).1H NMR(DMSO-d6,400 MHz):3.42(s,2H),4.28(s, 4H),5.06(s,2H),7.25～7.61(m,8H),10.78(s,1H). ESI-MS m/z:335.20(M+).Anal.Calcd.for C18H17N5O2(%):C,64.47;H,5.11;N,21.88.Found(%):C, 64.28;H,5.41;N,21.42.

1.3Synthesis of the cobalt(Ⅱ)complex(1)

The ligand LH(0.34 g,1 mmol)was dissolved in methanol(40 mL),and the resulting solution was adjusted pH value to 8.0 using 0.5 mol·L-1NaOH. Then CoCl2·6H2O(0.24 g,1 mmol)and NaClO4(0.12 g,1 mmol)were added to the stirred solution at 333 K for 6 hours.The solution was cooled to room temperature and filtrated to remove the insoluble substance.Red crystals suitable for X-ray diffraction studies were obtained after five days.m.p.214～216℃.UV-Vis spectra(MeOH solution):278 nm(ε= 24 400 L·mol-1·cm-1).IR(KBr,cm-1):3 025(m),2 942(s),1 716(s),1 545(s),1 460(s),1 405(s),1 056 (m),755(m).Anal.Calcd.for C72H86N20O23Co4Cl4(%):C, 43.70;H,4.35;N,14.16.Found(%):C,43.34;H, 4.08;N,14.42.

1.4Synthesis of the other complexes(2～7)

The other preparations of transition metal(Mn(Ⅱ), 2;Fe(Ⅱ),3;Ni(Ⅱ),4;Cu(Ⅱ),5;Zn(Ⅱ),6)complexes containing the ligand LH were similar to the cobalt(Ⅱ)complex(1)using metal chlorides(namely,MnCl2· 4H2O,FeCl3·6H2O,NiCl2·6H2O,CuCl2·2H2O,ZnCl2), but NaClO4was not added.The dicobalt(Ⅱ)complex [Co2(EGTB)(NO3)2(DMF)2](NO3)2·2DMF(7)was prepared according to the literature[11].

[MnLCl](2):Grey crystalline powder.Yield:53%, m.p.282～284℃.UV-Vis spectra(MeOH solution): 282 nm(ε=22 840 L·mol-1·cm-1).Molar conductance, ΛM:(DMF,S·cm2·mol-1)12.5.IR(KBr,cm-1):3 044 (m),2 935(s),1 722(s),1 578(s),1 452(s),1 412(s), 1 066(m).ESI-MS m/z:424.54(M+).Anal.Calcd.for C18H16N5O2MnCl(%):C,50.88;H,3.77;N,16.49. Found(%):C,51.14;H,3.52;N,16.02.

[FeLCl]Cl(3):Orange crystalline powder.Yield: 42%,m.p.＞300℃.UV-Vis spectra(MeOH solution): 276 nm(ε=24 510 L·mol-1·cm-1).Molar conductance, ΛM:(DMF,S·cm2·mol-1)72.8.IR(KBr,cm-1):3 012 (m),2 936(s),1 755(s),1 568(s),1 444(s),1 408(s), 1 036(m).ESI-MS m/z:425.56(M+).Anal.Calcd.for C18H16N5O2FeCl2(%):C,46.85;H,3.47;N,15.18. Found(%):C,47.08;H,3.89;N,14.86.

[NiLCl]·2H2O(4):Pale-green crystalline powder. Yield:62%,m.p.252～253℃.UV-Vis spectra(MeOH solution):272 nm(ε=20 980 L·mol-1·cm-1).Molar conductance,ΛM:(DMF,S·cm2·mol-1)9.6.IR(KBr, cm-1):3 048(m),2 955(s),1 762(s),1 526(s),1 438(s), 1 402(s),1 048(m).ESI-MS m/z:428.51(M-2H2O)+. Anal.Calcd.for C18H20N5O4NiCl(%):C,46.53;H, 4.31;N,15.08.Found(%):C,46.16;H,4.22;N,15.56.

[CuLCl](5):Green crystalline powder.Yield:58%, m.p.＞300℃.UV-Vis spectra(MeOH solution):276 nm(ε=24 380 L·mol-1·cm-1).Molar conductance,ΛM: (DMF,S·cm2·mol-1)11.2.IR(KBr,cm-1):3 052(m), 2 938(s),1 762(s),1 555(s),1 428(s),1 412(s),1 026 (m).ESI-MS m/z:433.14(M+).Anal.Calcd.for C18H16N5O2CuCl(%):C,49.88;H,3.67;N,16.17.Found(%): C,50.14;H,4.02;N,16.62.

[ZnLCl](6):Pale-yellow crystalline powder.Yield: 54%,m.p.＞300℃.UV-Vis spectra(MeOH solution): 278 nm(ε=20 850 L·mol-1·cm-1).Molar conductance, ΛM:(DMF,S·cm2·mol-1)11.8.IR(KBr,cm-1):3 048 (m),2 944(s),1 732(s),1 568(s),1 466(s),1 422(s), 1 014(m).ESI-MS m/z:434.55(M+).Anal.Calcd.for C18H16N5O2ZnCl(%):C,49.71;H,3.68;N,16.11.Found (%):C,49.48;H,3.26;N,16.58.

1.5Structure determ ination of com p lex 1

Red crystal of compound 1 having approximate dimensions of 0.30 mm×0.30 mm×0.20 mm was mounted on a glass fiber in a random orientation at 298(2)K.The determination of unit cell and the data collection were performed with Mo Kα radiation(λ= 0.071 073 nm)on a Brucker Smart-2000 CCD diffactometer.A total of 16 976 reflections were collected in the range of 1.63°<θ<25.68°at room temperature. The original data were integrated and corrected for Lorentz polarization effects using SAINT program and were corrected for absorption effects[12-13].Structures were solved by direct methods using SHELXS-97,and all non-hydrogen atoms were directly located from difference Fourier maps and successfully refined with anisotropic displacement parameters[14].The coordinates of hydrogen atoms only participated in the structure factor calculation and not included in the refinement. Hydrogen atoms in the complexes were placed in idealized positions,and those in water molecules were located in the difference Fourier maps.The Cl(2)and O(5)are both half occupancy.Due to the Cl(2)is too close to O(5),it was necessary to apply SIMU restraints to the O(5)to avoid the high ellipsoid in O(5). Because the guest solvent molecules are highly disordered and impossible to refine using conventional discrete-atom models,the QUEEZE subroutine of the PLATON software suite was applied to remove the scattering from the free water molecules[15].A SQUEEZE analysis showed that the void space was occupied by 104 electrons per formula,corresponding to 10 water molecules per formula unit.The number of free water molecules were also determined by element analysis. The final R=0.045 9 and wR=0.150 2(I＞2σ(I)),S= 1.122,(Δ/σ)max=0.001,(Δρ)max=852 e·nm-3and(Δρ)min=-375 e·nm-3.Crystal data and structural refinement were listed in Table 1.The selected bond lengths and bond angles were summarized in Table 2,and the hydrogen bond lengths and bond angles were given in Table 3.

CCDC:950329,1.

Table 1 C rystal data and structure refinements of com p lex 1

Table 2 Selected bond lengths(nm)and bond angles(°)of com p lex 1

Table 3 Hydrogen bond lengths(nm)and bond angles(°)

1.6Catecholase m im ic activity

3,5-di-tert-butylcatechol(3,5-DTBC)was used as the substrate to investigate the catalytic activity of the catecholase model complex[7].Kinetic experiments for the oxidation of 3,5-DTBC were performed spectrophotometrically on a Analytik jena Specord 210 spectrophotometer by the characteristic absorption band of the expected product,3,5-di-tert-butyl-obenzoquinone,at 400 nm.Freshly prepared stock solutions of the complex(0.10 mmol·L-1)and 3,5-DTBC(1.0 mmol·L-1)were used under aerobic conditions at 25℃.Oxidative reactions of 3,5-DTBC (0.010 mmol·L-1～1.0 mmol·L-1)were initiated by the addition of the complex(0.010 mmol·L-1)in methanol-tris-HCl(tris=trihydroxymethylaminomethane,0.01 mol·L-1)buffer solution.The pH values of the reaction system changed in the range of 5～11.The absorbance versus wavelength(wavelength scan)of the solutions was recorded at a regular time interval of 5 min.A kinetic treatment on the basis of the Michaelis-Menten approach was applied and the results were evaluated from Lineweaver-Burk double-reciprocal plots.In each kinetic trial,8 points were recorded to obtain the kinetic data.The 3,5-DTBC solution in the absence ofthe complex was taken as a control.

2　Results and discussion

2.1Crystal structure of complex 1

In the crystal structure of complex 1(Fig.1),the asymmetric unit has two cobalt(Ⅱ)ions.One cobalt(Ⅱ)ion is five-coordinated by one chloride ion Cl(1),two benzimidazole nitrogen atoms(N(2),N(4)),one carboxy oxygen atom O(1)and one amine nitrogen atom(N(1)) of the ligand LH,forming a distorted trigonal biyramid coordination geometry(angular structural parameter(τ) equals to 0.96)[16].Cl(1)and N(1)occupy the axial position,and N(2),N(4)and O(1)located in the equatorial plane.The Cl(1)-Co(1)-N(1)bond angle is 177.35(6)°,which is slightly deviated from ideal 180° (Table 1).However，the other cobalt(Ⅱ)ion is partly coordinated by one water oxygen atom(O(5))(or one chloride ion Cl(2)),which is unusual for such complex. Similar example has also been reported for[Cu(4-atrz)4(Cl)0.5(H2O)0.5](ClO4)1.5with equal Cl/O ratio in the same position,and this may offer an opportunity to fine control over the Cl/O composition and properties for the complex[17].Since the Cl(2)is half occupancy,the ClO4-ion is located to balance the charge.N(6)and O(5)(or Cl(2))reside in the axial position,and N(7), N(9)and O(3)make up the equatorial plane.The Co-N bond distances range from 2.035(2)to 2.320(2)nm, and the distance of Co-N(amino)is obviously longer than that of Co-N(benzimidazole),thus forming a distorted trigonal biyramid coordination geometry(τ= 0.80).It′s also found that the carboxyl group is deprotonized,which leads to a single negative charge ligand L-.The dihedral angles between two benzimidazole rings of the deprotonized ligand L-in the unit cell are 56.68°and 56.01°,respectively.As shown in Table 2 and Fig.2,the uncoordinated anion ClO4-and water molecules participate in the formation of hydrogen bonds.The molecules are stabilized by intermolecular N-H…O,C-H…O and C-H…Cl hydrogen bonds,leading to the formation of a three dimension network.

Fig.1 An ORTEP view of Structure of comp lex 1

Fig.2 Packing of complex 1 in unit cell

2.2Studies of catecholase activity

3,5-DTBC has been widely chosen as a substrate in catecholase model complex studies.It′s readily oxidized to the corresponding 3,5-di-tert-butyl-oquinone(3,5-DTBQ)owing to its low redox potential, and its bulky substituents prevent further reactions such as ring opening[18].3,5-DTBQ is considerably stable and has a strong absorption at λmax=400 nm(ε= 1 900 L·mol-1·cm-1),so that the activity and reaction rate could be determined using UV spectroscopy[19]. Fig.3 showed the variation of the spectral behavior of complex 1(0.01 mmol·L-1)in the presence of 3,5-DTBC(100 equiv.)under aerobic conditions.The band that corresponded to 3,5-DTBQ was observed at 400 nm after addition of the substrate(3,5-DTBC)to the solution of complex 1.

Fig.3 UV-Vis spectra of the methanol solution of complex 1(0.1 mmol·L-1)after addition of 3,5-DTBC(10 mmol·L-1)for 4 h(a),complex 1(b)and 3,5-DTBC(c)

The kinetic studies on the oxidation of 3,5-DTBC were carried out by the method of initial rates.It′s found that the reaction between complex 1 and substrate(3,5-DTBC)was a first order reaction,so their relationship could be described using a formula: ln[(Af-Ai)/(Af-At)]=kt[20].Here,k is the rate constant, Af,Aiand Atwere the final,initial and the time t absorbance of the reaction solution,respectively.The maximum absorbance of the growth of the quinone (3,5-DTBQ)at 400 nm was measured as a function of time,and Fig.4 showed a linear relationship in the ln[(Af-Ai)/(Af-At)]versus reaction time(t)plot in the presence of complex 1.The effect of the complexes (1～7)concentration on the oxidation of 3,5-DTBQ was also observed at pH 9.0 and 25℃(Fig.5).The reaction rate was initially increased linearly with increasing concentration of these complexes,indicating a firstorder dependence on each complex concentration.As shown in Fig.6,the pH dependence of the catecholase activity was promoted by the complexes(1～7)in the range of 5.0～11.0.The plots of v0(initial rate)vs pH value were observed to exhibit sigmoid-shape profiles, and the catalytic activities of all complexes were nearly increased in the pH value range from 7.0 to 9.0.

Fig.4 Plot of ln[(Af-Ai)/(Af-At)]versus the reaction time at pH 9.0 and 25℃

Fig.5 Complexes(1～7)dependence of oxidation of 3,5-DTBC(1 mmol·L-1)at pH 9.0 and 25℃

Fig.6 pH value dependence for the oxidation of 3,5-DTBC catalyzed by the complexes(1～7)

Fig.7 Plot of the reaction rates vs concentration of 3,5-DTBC for complex 1

The data depicted in Fig.7 revealed saturation kinetics of complex 1 with Michaelis-Menten-likebehavior.So an analysis based on the Michaelis-Menten model,originally developed for enzyme kinetics,was applied.The Michaelis-Menten equation is usually expressed as following:

Here,vmaxis the maximum initial rate,Kmis Michaelis constant,kcatis the turnover number,namely,the amounts(mol)of substrate converted to product per minute per mole of mimetic enzyme.cSand cMare the concentrations of substrate and model complex, respectively.The Michaelis-Menten equation can be turned into Lineweaver-Burk equation by taking its double reciprocal:

The inset in Fig.7 showed 1/v0was linear with 1/c3,5-DTBCfor complex 1,implying a first order reaction dependence.The slope of the straight line was Km/vmax, and the intercept on the vertical axis was 1/vmax. Several kinetic parameters vmax,Kmand kcatwere calculated from the Lineweaver-Burk plot model (Table 4).It is found that kcatof complex 1 toward 3,5-DTBC increased from 4.35 min-1(20℃)to 13.55 min-1(45℃).So the catalytic activity of complex 1 increased with increasing of temperatures.The kinetic data of the complexes(1～7)were also measured at pH 9.0 and 25℃.As shown in Table 5,the catalytic activity of these complexes were easily arranged in descending order:3＞7＞5＞1＞4＞2＞6.The catecholase model complexes have been reported to relate to several factors such as metal-metal distance,lability of exogenous ligands,coordination configuration and electrochemical properties[21].Of all the complex,3 showed the best catalytic activity,its difference with the other complexes(1,2,4,5,6)could be attributed to the redox potential effect of various metal ions on thecatecholasemodelcomplexes,and the results were in accordance with the literatures reported previously[11,22-23]. The catalytic activity of complex 7 was better than complex 1,because the former was a dicobalt(Ⅱ)complex,which may be dependent upon the interaction between each active center and substrate due to its long dinuclear distance[11].A probable explanation for the catecholase catalytic reactivity of complex 1 was mainly related to its coordination environment.The five-coordinated cobalt(Ⅱ)ion has an unsaturated vacancy,and it could combine with the substrates[7,11].

Table 4 K inetic data of com p lex 1 at pH 9.0

Tab le 5 K inetic data of the com p lexes 1～7 at pH 9.0 and 25℃

3　Conclusions

In summary,we have prepared a cobalt(Ⅱ)complex,which was characterized by single crystal X-ray diffraction.The complex showed pH value dependence for the oxidation of 3,5-DTBC in the range of 5～11,and the kinetics of 3,5-DTBC catalyzed by the complex obeyed the Michaelis-Mentent equation,its catecholase catalytic activity increased with the rise of the temperature.

[1]Koval I A,Gamez P,Belle C,et al.Chem.Soc.Rev.,2006, 35:814-840

[2]Klabunde T,Eicken C,Sacchettini J C,et al.Nat.Struct. Biol.,1998,5:1084-1090

[3]Patra A,Giri G C,Sen T K,et al.Polyhedron,2014,67:495-504

[4]Caglar S,Aydemir I E,Adiguzel E,et al.Inorg.Chim.Acta, 2013,408:131-138

[5]Caglar S,Adiguzel E,Caglar B,et al.Inorg.Chim.Acta,2013, 397:101-109

[6]Safaei E,Sheykhi H,Wojtczak A,et al.Polyhedron,2011, 30:1219-1224

[7]ZHANG Yong(張勇),ZHOU Xia(周霞),HU Jia-Wei(胡家尾),et al.Chinese J.Struct.Chem.(結(jié)構(gòu)化學(xué)),2013,32(9): 1291-1296

[8]Apurba B,Lakshmi K D,Michael G B D,et al.Inorg.Chem., 2012,51:7993-8001

[9]CHEN Yan-Guo(陳彥國),LIAO Zhan-Ru(廖展如),LI Wu-Ke(李武客),et al.Acta Chim.Sinica(化學(xué)學(xué)報),2000,58 (10):1191-1195

[10]Ogawa K,Nakata K,Ichikawa K.Chem.Lett.,1998,27:797-798

[11]Zhang Y,Meng X G,Liao Z R,et al.J.Coord.Chem., 2009,62:876-885

[12]SAINT,Program for Area Detector Absorption Correction, Siemens Analytical X-Ray Instruments Inc.,Madison,W I 53719,USA,1994-1996.

[13]Sheldrick G M.SADABS,Program for Siemens Area Detector Absorption Correction,University of G?ttingen,Germany, 1996.

[14]Sheldrick G M.SHELXS-97,Program for the Solution of Crystal Structures,University of G?ttingen,Germany,1997.

[15]Spek A L.Acta Crystallogr.,2009,D 65:148-155

[16]Addison A W,Rao T N,Reedijk J,et al.J.Chem.Soc. Dalton Trans.,1984:1349-1356

[17]Yi L,Du J Y,Liu S,et al.J.Chem.Res.,2004,1:29-31

[18]Rompel A,Fischer H,Meiwes D,et al.FEBS Lett.,1999, 445:103-110

[19]Seneque O,Campion M,Douziech B,et al.Eur.J.Inorg. Chem.,2002,8:2007-2014

[20]Chen Z F,Liao Z R,Li D F,et al.J.Inorg.Biochem.,2004, 98:1315-1318

[21]Neves A,Rossi L M,Bortoluzzi A J,et al.J.Braz.Chem. Soc.,2001,12:747-754

[22]Qiu J H,Liao Z R,Meng X G.Polyhedron,2005,24:1617-1623

[23]Lu L J,Song Y Y,Liu H,et al.J.Coord.Chem.,2012,65: 1278-1288

Synthesis,Crystal Structure and Catecholase Activity of the Cobalt(Ⅱ)Com plex Containing Benzim idazole Ligand

ZHANG Qian1HAN Yan2JIAO Yuan-Hong＊,1
(1College of Chemistry and Chemical Engineering,Hubei Polytechnic University,Huangshi,Hubei 435003,China)
(2College of Chem istry and Chemical Engineering,Xinxiang University,Xinxiang,Henan 453003,China)

A cobalt(Ⅱ)complex,[CoLCl]2[CoLCl0.5(H2O)0.5]2ClO4·10H2O(1)(LH=(2-(bis((1H-benzo[d]-imidazol-2-yl) methyl)amino)acetic acid))has been synthesized and characterized.Its catecholase catalytic activity has been exploited by UV-Vis spectrophotometric study by using 3,5-di-tert-butylcatechol(3,5-DTBC)as the substrate.The results showed that two cobalt(Ⅱ)ions of the asymmetric unit were both coordinated by the ligands in a distorted trigonal-bipyramidal geometry in the crystal structure.Complex 1 exhibited pH value dependence for the oxidation of 3,5-DTBC in the range of pH 5 to 11,and its catecholase catalytic activity increased with increasing temperature.The kinetics of oxidation of 3,5-DTBC catalyzed by complex 1 accorded with the Michaelis-Mentent equation.CCDC:950329,1.

cobalt(Ⅱ)complex;crystal structure;catecholase;catalytic activity

O614.81+2

A

1001-4861(2016)01-0131-08

10.11862/CJIC.2016.031

2015-07-21。收修改稿日期：2015-11-16。

湖北省教育廳B類項目(No.B20114404，B2013067)資助。

＊通信聯(lián)系人。

E-mail：yhjiao2014@163.com