pH-Controlled Drug Release of Nitrogen Doped Carbon Dots Delivering Doxorubicine for Synergetic Photo-Thermal Therapy and Chemotherapy

Zheng YANG An-Jian XIE Yu-Hua SHEN

Citation:  YANG Zheng, XIE An-Jian, SHEN Yu-Hua. pH-Controlled Drug Release of Nitrogen Doped Carbon Dots Delivering Doxorubicine for Synergetic Photo-Thermal Therapy and Chemotherapy[J]. Chinese Journal of Inorganic Chemistry, 2018, 34(10): 1775-1782. doi: 10.11862/CJIC.2018.236 shu

pH药物控释的氮掺杂碳点载运阿霉素用于光热与化疗的协同治疗

    通讯作者: 沈玉华, s_yuhua@163.com
  • 基金项目:

    国家自然科学基金(No.21571002和No.21671001)和安徽省环境友好高分子材料重点实验室项目

    国家自然科学基金 21671001

    国家自然科学基金 21571002

摘要: 提出一种利用氮掺杂碳点(N-CDs)的光热性能与化疗药物阿霉素(Dox)相结合的联合治疗肿瘤模式。实验结果标明,所制备的N-CDs的分散液经光照后升温可达10℃,是性能优异的光热剂。同时,载运Dox后的N-CDs-Dox纳米复合物还具有pH触发的药物控释作用。因此,这种多功能的N-CDs-Dox纳米复合物能够实现光热和化疗的协同作用,有效杀伤肿瘤细胞。

English

  • Cancer, especially the malignant tumor, which is one of the severe threats to anthropic health, is an arduous challenge till to now[1]. Heretofore, massive methods have been developed to cure cancer, such as radiation therapy, surgical therapy, and chemotherapy[2]. Although accompanying with the successful cancer ablation, the insurmountable defects of these treatments such as low selectivity, huge invasion, and serious side effects, abate heavy pains of patients suffering from cancer. To overcome the fatal disease, photo-thermal therapy (PTT) basing on the principle of heating target cells has become a feasible alternative for the cancer treatment recently[3]. Near infrared (NIR) light-induced PTT which could cure tumor in deeper tissues is an emerging therapy in anticancer domain[4]. Two ranges of wavelengths are called bio-windows benefitting for PTT, which are 650~950 nm and 1 000 ~1 350 nm respectively[5]. Comparing to traditional cancer treatments, PTT destroys tumor selectively by extracorporeal laser irradiation, demonstrating that PTT is a noninvasive therapy with high selectivity. In past decades, massive amounts of optical absorbing materials for PTT have been developed, such as gold nanomaterials[7], carbon nanoparticles[8], and other metal-recombined nanocomposites[9]. However, the reasonable combination of PTT and other therapies to cure cancer more effectively is still main barrier in synergistic cancer treatment.

    For solving the above-mentioned challenges, the rapid development of nanotechnology becomes a key avenue. Lately, fluorescent N-CDs possessing high water-solubility, facile surface modification, limited toxicity, excellent biocompatibility and extraordinary photo-stability, have been a kind of quite momentous nanomaterials in biomedical field[10]. On account of aforementioned features, N-CDs could convert the absorbed NIR light into heat, proving to be an ideal PTT agent. However, augmenting photo-thermal effect and drug delivery capacity of N-CDs for enhanced synergetic effect between PTT and chemotherapy is still an insurmountable problem.

    In this paper, we synthesized the N-CDs by a facile one-step hydrothermal process. Even more important, the synthesized N-CDs show a broad and forceful absorption band in the range of visible to NIR light (400~800 nm). Utilizing this property of N-CDs, we successfully prepared a simple and low-cost photo-thermal agent with high-performance for PTT in cancer therapy. Using the N-CDs as anti-cancer drug carrier, the Dox-loaded N-CDs could achieve the combined PTT and chemotherapy synergistically.

    Critic acid, urea, and dimethyl sulfoxide (DMSO) were purchased from Aladdin Industrial Inc. 3-(4, 5-Dimethylthiazol-2-yl)-2, 5-diphenyltetrazolium bromide (MTT), high glucose dulbecco′s modified eagle′s medium (DMEM), phosphate buffer solution (PBS), and fetal bovine serum (FBS), were obtained from Sinopharm Chemical Reagent Co., Ltd. Penicillin, streptomycin, pancreatin solution (25% (w/w)), hoechst 33342, and propidium iodide (PI), were purchased from Sangon Biotech Inc. All the agents were analytical pure and used as received without further purification. The experimental water were double distilled water (18.0 MΩ·cm).

    Morphological details of the obtained samples were recorded by Transmission electron microscopy (TEM) at an accelerating voltage of 100 kV (JEM-2100, Japan). Raman spectrum was studied using a Via-Reflex laser confocal Raman spectroscopy under a 785 nm laser excitation with the power of 340 mW. For investigating the crystallinity of samples, a DX-2700 X-ray diffractometer (Cu source, λ=0.154 056 nm, 40 V, 100 mA) was used to obtain the power X-ray diffraction (XRD) patterns over the 2θ range of 10° to 70°. The fourier transform infrared (FT-IR) spectra was observed by a NEXUS-870 FT-IR spectrometer using the KBr pellet technique. An ESCALAB250 spectrometer with Al source (=1 486.6 eV, 150 W, 1×10-6 Pa) was used to analyze the element distri-bution of samples by X-ray photoelectron spectros-copy. The UV-Vis spectra were obtained using a Shimadzu UV-3600. The photothermal effects of N-CDs dispersions were studied using a Fluke Ti32 infrared thermography camera to obtain the thermographies.

    Using a reported facile one-step hydrothermal treatment procedure, we synthesized the N-CDs[11]. Briefly, 1.0 g of citric acid and 0.7 g of urea were added into 30 mL of DI water. Then the mixed solution was sonicated to be transparent solution. After that, the as-prepared solution was poured to a 50.0 mL Teflon-lined stainless autoclave with sealing tightly and heated to 180 ℃ with maintaining for 5 h by hydro-thermal reaction. After naturally cooling down to room temperature, the hyaline solution with brown and yellow color was dialyzed in a dialysis bag (3.5 kDa) over night. Finally, the dialyzed suspension was frozen in ultra-low temperature freezer (-86 ℃) for 10 min and then was subjected to freeze-drying to obtain the black-green solid powders, i.e., N-CDs.

    For measuring the PTT effect of the N-CDs, 2 mL of the various samples dispersions with different concentrations (0.5, 1.0, 2.0, 5.0 and 10.0 mg·mL-1) were added in a small glass bottle and then irradiated with an 650 nm NIR laser at regular time intervals, respectively. The temperatures of dispersions were recorded by the Fluke Ti32 thermal infrared camera immediately after irradiation every 1 min. The equal volume of water was used as controls. All data were acquired from three independent experiments to ensure the reliability of the data.

    For loading N-CDs with Dox, N-CDs (20 mg) were added to Dox solution (150 μg·mL-1, 20 mL) under sonicating. Stirring by a mechanical agitator for 24 h, then the hybrid nanocomposites were dialyzed in a dialysis bag (3.5 kDa) over night. The final products were frozen in ultra-low temperature freezer (-86 ℃) for 10 min and then were subjected to freeze-drying. To calculate the efficacy of Dox loading, 1 mg of N-CDs-Dox powders were dispersed in 5 mL deionized water. The absorption peak at 510 nm was used for validating the presence of Dox and estimating the Dox loading efficacy by comparison with the standard curve.

    In order to investigate the drug release in different environments, 2 mg of N-CDs-Dox powders were added to 1 mL of PBS solution at pH values of 5.0, 6.5 and 7.4 in dialysis bags respectively. Then the dialysis bags were put into 9 mL of PBS solution with the same pH value respectively. At different time points during 0 to 72 h, 1 mL of dialyzed PBS solution was collected and analyzed by UV-Vis spectrophoto-meter at 510 nm to calculate the drug releasing efficiency according to the standard curve of Dox solution. After determination, the collected solution was put back to original solution for continuous dialysis.

    The as-prepared N-CDs and N-CDs-Dox nanocomposites were dispersed in DMEM medium to form mixed dispersions with different concentrations of 1.0, 2.0, 5.0, 7.5 and 10.0 mg·mL-1. The HeLa cells suspended in DMEM medium was seeded into 96-well plates at a density of 5×103 each well (200 μL) and incubated for 24 h. After washing each well by PBS for 3 times, the cells were incubated with only 100 μL of DMEM medium (control), and with 100 μL of DMEM medium containing the samples of N-CDs and N-CDs-Dox nanocomposites (0.5, 1.0, 2.0, 5.0, 7.5 and 10.0 mg·mL-1, respectively) for 24 h, respectively. The DMEM without samples was used as control group at the same condition. After washing each well by PBS for 3 times, 20 μL of MTT and 100 μL of DMEM medium were added into each well with sequential incubation for 4 h. As contrast, the irradiated groups were taken out with the irradiation of 650 nm laser (1 W·cm-2) for 10 min and then incubated for 24 h. 100 μL of DMSO solution was then added for dissolving the insoluble formazan crystals. To investigate the selectivity of the prepared N-CDs normal cells, 293T cells were used as control incubating at the same conditions, following with the equal post-treatments. The optical density of samples was measured using an Elisa reader at 490 nm to calculate the cell viability by the formula:

    ${\rm{Cell\;\;viability = }}\frac{{{\rm{O}}{{\rm{D}}_{{\rm{exp}}}} - {\rm{O}}{{\rm{D}}_{{\rm{bla}}}}}}{{{\rm{O}}{{\rm{D}}_{{\rm{con}}}} - {\rm{O}}{{\rm{D}}_{{\rm{bla}}}}}} \times 100 $

    Where ODbla, ODexp, and ODcon are the optical densities of blank, experimental and control groups, respectively.

    To further understand the therapeutic abilities of as-prepared N-CDs and N-CDs-Dox nanocomposites against cancer cells, the Hela cells were seeded into 6-well plates at a density of 5×103, incubating at 37 ℃ under 5% (V/V) CO2 atmosphere (95 % (V/V) air) in the CO2 constant temperature incubator for 24 h. Then 2 mL of medium with the sample concentration of 10.0 mg·mL-1 was added into part of wells, using pure medium as control. When the incubation time lasting for 4 h, the irradiated groups were taken out and irradiated using 650 nm laser (1 W·cm-2) for 10 min. After further 24 h incubation, PBS buffer (pH=7.4) was used to replace the medium, then 0.5 mL of Hoechst 33342 (10 μg·mL-1) and 0.5 mL of PI (10 μg·mL-1) were added into each well for 15 min. The final cells were washed twice with PBS and observed by an inverted fluorescence microscopy to obtain fluorescent images.

    According to Fig. 1a, TEM image of N-CDs investigates that the as-prepared N-CDs possess uniform size with the diameter about 7 nm[12]. Further evidences coming from Fig. 1b and 1c, it is seen that the synthesized N-CDs are dispersed uniformly on the substrate and possess regular shapes with a size of about 7.0 nm and a height of ca. 2.2 nm, which is consistent with that of TEM image. As shown in Fig. 1d, Raman spectroscopic analysis shows that the N-CDs exhibit two bands around 1 352 and 1 597 cm-1 relating to the D band of disordered carbon structure and the G band of graphitic layer of the N-CDs, respectively[12]. Fig. 1e is the XRD pattern of N-CDs showing that a broad diffraction peak emerges at 2θ ranging from 10° to 40°, which illustrates the N-CDs are comprised by disorder carbon structure mainly. The result coming from XRD pattern is corresponding to that of Raman spectrum of N-CDs. FT-IR spectrum (Fig. 1f) reveals that the broad peak at 3 430 cm-1 is related to the stretching vibration of C-OH and N-H, while the peak at 1 128 cm-1 is corresponded to the asymmetric stretching vibration of C-NH-C[12]. Mean-while, the peak at 1 572 cm-1 is attributed to the flexure vibration of N-H, along with the peaks at 1 635 and 1 082 cm-1 are belonged to the specific vibration of C=O[12]. Based on upwards results, the surface of N-CDs contain large amounts of hydroxyl groups, carbonyl groups, amino groups, and epoxy groups which are benefitting for solubility of N-CDs in water.

    Figure 1

    Figure 1.  (a) TEM image, (b) AFM image, (c) height profile, (d) Raman spectrum, (e) XRD pattern, (f) FT-IR spectrum, (g) XPS survey spectrum, and (h) high-resolution spectrum of C1s of N-CDs

    According to the XPS survey spectrum of N-CDs in Fig. 1g, the as-synthesized N-CDs are combined by the elements of carbon, oxygen, and nitrogen primarily. As shown in Fig. 1h, the characteristic peaks at 284.75, 286.56 and 287.1 eV in the C1s spectrum, indicate the presence of varietal types of carbon bonds corresponding to sp2 C=C, and C=O as well as C-N, respectively[12]. The mentioned results demonstrated the obtained N-CDs are doped by nitrogen succe-ssfully, within amounts of hydroxyl and carbonyl groups on the surface which are consistent with the FT-IR spectrum completely.

    Further evidences coming from Fig. 2a, the drug loading content achieves about 1 μg·mg-1, illustrating that the N-CDs are promising drug carriers. The loaded drug may be combined with the N-CDs by main electrostatic attraction, also π-π stacking, hydrophobic and van der Waals interactions[13]. To investigate the drug loading mechanism in depth, zeta potentials of N-CDs dispersion, Dox solution, and N-CDs-Dox dispersion were obtained at pH=7.4 respectively (Fig. 2b). ζ potential of N-CDs dispersion is about -24.10 mV, showing the N-CDs are negatively charged materials. The Dox solution possesses the positively charged surface with the ζ potential of ca. 9.54 mV. After drug loading process, the prepared N-CDs-Dox dispersion shows the ζ potential was about -10.25 mV, which was enlarged comparing to that of N-CDs dispersion. The result confirms that Dox-loading of N-CDs is due to the charge interaction between Dox and N-CDs[14]. The release behaviors (Fig. 2c) of Dox from the N-CDs-Dox nanocomposites are studied in PBS buffer solution at different pH values of 5.5, 6.5 and 7.4 to simulate the lysosome, tumor site, and physiological conditions respectively. It is seen that when the pH value is 7.4, only 10% Dox is released from the nanocomposites after 72 h. While the pH value decreases to 6.5, the cumulative amount of released Dox increases to about 30% after the same time. On the contrast, while the pH value lowers to 5.5, the amount of the released Dox from the nanocomposites reached as high as ca. 98% within the same time. The kinetics of Dox release at low pH values may be induced by the increase of H+ in PBS buffer solution, which could weaken and destroy the charge interaction of Dox with N-CDs in the obtained nanocomposites[14]. The above results indicate the fact that the release behavior of Dox from the N-CDs-Dox nanocomposites was pH dependent and the decrease of pH value could enhance the release rate of the Dox from the nanocomposites significantly. All the results mentioned above demonstrate that N-CDs are potential drug carriers which could achieve pH-controlled drug release at tumor site locally.

    Figure 2

    Figure 2.  (a) UV-Vis spectra of Dox before and after loading in the N-CDs; (b) ζ potentials of N-CDs dispersion, Dox solution, and N-CDs-Dox dispersion at pH=7.4, respectively; (c) Release curve of Dox at different pH values from N-CDs-Dox nanocomposites

    The photo-thermal effects of all samples were shown in Fig. 3. Under continuous irradiation of 650 nm laser for 10 min, the temperature elevation of pure water increased about 3 ℃, indicating the used laser could not induce obvious temperature elevation in pure water. However, the irradiated N-CDs dispersions (0.5, 1.0, 2.0, 5.0 and 10.0 mg·mL-1, respectively) emerge remarkable temperature changes depending on concentrations of samples. The highest temperature enhancement of N-CDs dispersion (10.0 mg·mL-1) is about 10 ℃, demonstrating that the prepared N-CDs possess excellent photo-thermal ability. The above evidences indicate the synthesized N-CDs are suitable photo-thermal agents for PTT.

    Figure 3

    Figure 3.  Photo-thermal effects of pure water, and N-CDs dispersions (0.5, 1.0, 2.0, 5.0 and 10.0 mg·mL-1, respectively) with a 650 nm laser irradiation for 10 min

    Fig. 4a shows the results of MTT assay using different concentrations of samples to evaluate the cytotoxicity to HeLa cells. The group which is incubated with N-CDs in dark shows high cell viabil-ities (all above 90%), illustrating the obtained N-CDs possess good biocompatibility. Meanwhile, N-CDs co-incubated group shows the obvious apoptosis of HeLa cells after 10 min laser irradiating, with the cell mortalities about 27%, 40% and 57% at the concen-trations of 5.0, 7.5 and 10.0 mg·mL-1 respectively. The obtained results demonstrate that the prepared N-CDs have excellent property in PTT. Comparing to N-CDs-Dox co-incubated groups in dark, the irradiated N-CDs-Dox group exhibits serious photo-induced cell apoptosis with the cell viabilities of about 63%, 32% and 18% at the concentrations of 5.0, 7.5 and 10.0 mg·mL-1 respectively, which are higher than those of irradiated N-CDs groups. The calculated experimental data performed that the drug loaded N-CDs own outstanding synergistic PTT and chemotherapy to HeLa cells, which are potential multifunctional drug carriers for cancer treatment. For comparison, Fig. 4b illustrates the cell viabilities of 293T cells incubating at the same conditions to HeLa cells. The results of all groups are similar to those of HeLa cells with the uniform condition incubations, demonstrating that the obtained N-CDs do not possess prominent selectivity to normal cells and cancer cells. However, the outstanding photothermal ability of N-CDs is benefitting for synergetic therapy with chemotherapy to superficial tumors by the injection method, such as melanoma etc., which could decrease the damage to normal cells as well as improve the antitumor effect.

    Figure 4

    Figure 4.  Viabilities of (a) HeLa and (b) 293T cells incubated with different concentrations of N-CDs and N-CDs-Dox nanocomposites without or with irradiation (650 nm, 1 W·cm-2, 10 min)

    Base on the evidences shown in Fig. 5, the non-irradiated groups of HeLa cells incubate with control, and N-CDs which are stained by Hoechst 33342 and PI, emitting blue fluorescence (a1, b1) and almost invisible red fluorescence (a2, b2) respectively, as well as blue fluorescence (a3, b3) in the merged images, indicating that all the cells grow well without obvious cell apoptosis. The group in dark incubates with N-CDs-Dox nanocomposites shows partial red fluore-scence (c2), demonstrating that N-CDs-Dox nanocom-posites could achieve chemotherapy due to the loaded drug without irradiation. The irradiated group of control shows blue fluorescence (d1) and negligible red fluorescence (d2), demonstrating that the mere laser irradiation could only induce negligible apoptosis of HeLa cells. Meanwhile, the group of N-CDs exhibits local red fluorescence (e2), illustrating the photo-thermal effect of N-CDs could kill cancer cells effectively. Moreover, when N-CDs-Dox nanocom-posites co-incubated cells are irradiated by laser for 10 min, an intense red fluorescence (f3) could be seen, which shows that N-CDs-Dox nanocomposites caused a large amount of cell apoptosis with laser irradiation. The reason is that the PTT ability of N-CDs and chemotherapy property of Dox in N-CDs-Dox nanocomposites realized excellent synergetic efficacy in cancer treatment. Therefore, the N-CDs are novel PTT agents which could be used as drug carrier to achieve combined PTT and chemotherapy.

    Figure 5

    Figure 5.  Fluorescence microscopy images of HeLa cells incubated with (a, d) control (without samples), (b, e) N-CDs, and(c, f) N-CDs-Dox nanocomposites, (a~c) without laser irradiation and (d~f) with laser irradiation (650 nm, 1 W·cm-2)

    HeLa cells were dyed in blue by Hoechst 33342 (1), red by PI (2) and the merged images of both above (3), respectively

    In summary, the water-soluble N-CDs are synthesized and applied for in vitro photo-thermal ablation of cancer cells successfully. According to the experimental results, the as-prepared N-CDs possess good biocompatibility, and efficacious photo-thermal ability as well as favorable drug delivery property. Meanwhile, the N-CDs-Dox nanocomposites exhibit synergetic efficacy of cancer therapy owing to the chemotherapy induced by the loaded Dox and the photo-thermal effect attributing to N-CDs. This study paves a new way for the design and synthesis of other multi-functional materials.

    1. [1]

      王沛, 刘源岗, 王士斌.科学通报, 2017, 12:1233-1240 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=KXTB201712006&dbname=CJFD&dbcode=CJFQWANG Pei, LIU Yuan-Gang, WANG Shi-Bin. Chin. Sci. Bull., 2017, 12:1233-1240 http://kns.cnki.net/KCMS/detail/detail.aspx?filename=KXTB201712006&dbname=CJFD&dbcode=CJFQ

    2. [2]

      Trotti A. Int. J. Radiat Oncol. Biol. Phys., 2000, 47:1-12 doi: 10.1016/S0360-3016(99)00558-1

    3. [3]

      秦志国, 刘东, 杨芳, 等.药学进展, 2017, 41(11):812-823QIN Zhi-Guo, LIU Dong, YANG Fang, et al. Progress in Pharmaceutical Sciences, 2017, 41(11):812-823

    4. [4]

      Weissleder R. Nat. Biotechnol., 2001, 19:316-317 doi: 10.1038/86684

    5. [5]

      Yu Z F, Shi J P, Li J L, et al. J. Mater. Chem. B, 2018, 6:1238-1243 doi: 10.1039/C7TB03052E

    6. [6]

      Sheng D L, Liu T Z, Deng L M, et al. Biomaterials, 2018, 165:1-13 doi: 10.1016/j.biomaterials.2018.02.041

    7. [7]

      Zheng T T, Li G G, Zhou F, et al. Adv. Mater., 2016, 28:8218-8226 doi: 10.1002/adma.v28.37

    8. [8]

      Ge J C, Lan M H, Zhou B J, et al. Nat. Commun., 2014, 5:17054-17060 https://www.ncbi.nlm.nih.gov/pubmed/25105845?dopt=Abstract

    9. [9]

      Chou S S, Kaehr B, Kim J, et al. Angew. Chem. Int. Ed., 2013, 52:4160-4164 doi: 10.1002/anie.201209229

    10. [10]

      Ge J C, Jia Q Y, Liu W M, et al. Adv. Mater., 2015, 27:4169-4177 doi: 10.1002/adma.v27.28

    11. [11]

      Zhang Y, He Y H, Cui P P, et al. RSC Adv., 2015, 5:40393-40401 doi: 10.1039/C5RA04653J

    12. [12]

      Ahmed M, Byrne J A, McLaughlin J A, et al. e-J. Surf. Sci.Nanotechnol., 2009, 7:217-224 doi: 10.1380/ejssnt.2009.217

    13. [13]

      Wang Q L, Huang X X, Long Y J, et al. Carbon, 2013, 59:192-199 doi: 10.1016/j.carbon.2013.03.009

    14. [14]

      Pandey S, Thakur M, Mewada A, et al. J. Mater. Chem. B, 2013, 1:4972-4982 doi: 10.1039/c3tb20761g

  • Figure 1  (a) TEM image, (b) AFM image, (c) height profile, (d) Raman spectrum, (e) XRD pattern, (f) FT-IR spectrum, (g) XPS survey spectrum, and (h) high-resolution spectrum of C1s of N-CDs

    Figure 2  (a) UV-Vis spectra of Dox before and after loading in the N-CDs; (b) ζ potentials of N-CDs dispersion, Dox solution, and N-CDs-Dox dispersion at pH=7.4, respectively; (c) Release curve of Dox at different pH values from N-CDs-Dox nanocomposites

    Figure 3  Photo-thermal effects of pure water, and N-CDs dispersions (0.5, 1.0, 2.0, 5.0 and 10.0 mg·mL-1, respectively) with a 650 nm laser irradiation for 10 min

    Figure 4  Viabilities of (a) HeLa and (b) 293T cells incubated with different concentrations of N-CDs and N-CDs-Dox nanocomposites without or with irradiation (650 nm, 1 W·cm-2, 10 min)

    Figure 5  Fluorescence microscopy images of HeLa cells incubated with (a, d) control (without samples), (b, e) N-CDs, and(c, f) N-CDs-Dox nanocomposites, (a~c) without laser irradiation and (d~f) with laser irradiation (650 nm, 1 W·cm-2)

    HeLa cells were dyed in blue by Hoechst 33342 (1), red by PI (2) and the merged images of both above (3), respectively

  • 加载中
计量
  • PDF下载量:  8
  • 文章访问数:  273
  • HTML全文浏览量:  65
文章相关
  • 发布日期:  2018-10-10
  • 收稿日期:  2018-06-21
  • 修回日期:  2018-08-22
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

/

返回文章