Biological and Pharmaceutical Bulletin
Online ISSN : 1347-5215
Print ISSN : 0918-6158
ISSN-L : 0918-6158
Regular Articles
Metabolomic Profiling for Identification of Potential Biomarkers in the Protective Effects of Modified Sinisan against Liver Injury in Dimethylnitrosamine Treated Rats
Cheng-Gang LiuXiao-Li WangXiao-Wei DuDe-You JiangNai-Zhi GengShuo-Xin ZhangYuan-Yuan ZhouHai-Xue Kuang
Author information
JOURNAL FREE ACCESS FULL-TEXT HTML

2013 Volume 36 Issue 11 Pages 1700-1707

Details
Abstract

Metabolomics is a new platform based on the comprehensive analysis of low molecular weight metabolites and provides a powerful approach to discover biomarkers in biological systems. Modified Sinisan (MSNS), a traditional Chinese medicine formula, displayed bright prospects in the prevention and therapy of liver injury. However, its molecular mechanism of hepatoprotective effects remains unclear. This paper was designed to explore the effects and potential mechanisms of MSNS against dimethylnitrosamine-induced liver injury. Global metabolic profiling was performed by ultra-performance liquid chromatography/electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC/ESI-Q-TOF-MS) in conjunction with multivariate data analysis and pathway analysis. Eleven serum biomarkers were identified and pathway analysis results showed that phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, tryptophan metabolism, retinol metabolism, tyrosine metabolism were perturbed by liver injury. More importantly, MSNS has showed satisfactory pharmacological effect on liver injury through partially regulating the perturbed pathways, correlates well to the biochemical and histopathological detection results. The present study proved that the robust metabolomics approach is promising for unraveling hepatoprotective effects of MSNS and these findings provide new insights into mechanisms of the liver injury, and its pathophysiologic processes.

Sinisan is a widely used traditional Chinese medicine (TCM) formula and was first described in “Treatise on Cold-induced Febride Diseases (Shanghan Lun),” a treatise on exogenous febrile disease written by the famous Chinese physician Zhang Zhongjing (150 to 219 A.D. in the Chinese Eastern Han Dynasty). It is composed of four herbs including bupleuri radix, paeoniae radix alba, aurantii fructus immaturus and glycyrrhizae radix et rhizoma. Sinisan has been popularly used in folk medicine as a hepatic protectant.1) To enhance therapeutic efficacy of liver injury, Sinisan was modified by adding astragali radix, ostreae concha and xiongdan to achieve the property of soothing liver and invigorating Qi, softening hard lumps and dispeling nodes, clearing heat and detoxicating. Modified Sinisan (MSNS) has been clinically used as a common formula and produced evident effect on liver injury.2) However, its potential remedy mechanisms are unclear.

Metabolomics, a promising ‘omics’ platform, is the comprehensive analysis of low molecular weight metabolites in a biological sample.3) In recent years, metabolomics has shown great potential in many fields such as disease diagnosis,4) cancer,58) natural product discovery,9) toxicology,1012) and traditional medicine.1317) As a significant part of systems biology, metabolomics can be expected to provide new insights into the potential mechanisms of TCM.

In this study, the mechanisms of MSNS on dimethylnitrosamine (DMN)-induced liver injury rats were studied by examining plasma extracts obtained from the rats repeatedly exposed to DMN by the LC/MS-based coupled with pattern recognition approach. At the same time, serum biochemical analysis and histopathologic examinations were performed. The goal of the study is to acquire an overall understanding on the underlying molecular mechanisms of MSNS on liver injury rats. In addition, the identification of potential biomarkers in the protective effects of MSNS against liver injury in serum of dimethylnitrosamine treated rats was investigated by metabolomic approach for the first time.

MATERIALS AND METHODS

Reagents

Acetonitrile, formic acid were purchased from Merck (Darmstadt, Germany); Leucine enkephalin was purchased from Sigma-Aldrich (St. Louis, MO, U.S.A.); Dimethylnitrosamine was purchased by Tianjin Chemistry Reagent Research Institution (Tianjin, China); Deionized water was puried “in-house” using a Milli-Q ultra pure water system (Millipore, Bedford, MA, U.S.A.).

Preparation of MSNS

Chinese medicines were purchased from Harbin Shiyitang Drug Store (Harbin, China) for preparation of the MSNS extract. MSNS extract is a mixture of 7 crude ingredients: Bupleurum chinense DC., Citrus aurantium L., Paeonia lactiflora PALL., Processed Glycyrrhiza uralensis FISCH., Astragalus membranaceus (FISCH.) BGE. var. mongholicus (BGE.) HSIAO, Ostrea gigas THUNBERG, xiongdan. All the crude drugs were of high quality and were authenticated by Prof. Xiaowei Du, Department of Pharmacognosy of Heilongjiang University of Chinese Medicine.

MSNS was prepared at the pharmacognosy laboratory according to the method we optimized before. Ostrea gigas THUNBERG was immersed in 1900 mL water and refluxed for 30 min, then the other 6 drugs except xiongdan were added and refluxed for another 60 min. The solution was filtered with eight-layer gauze, and the residue was further extracted one time in the same way. The solutions of twice extractions were mixed and condensed under decompression to roughly 190 mL, and then 2 g xiongdan was added. Finally, the solution was prepared as freeze-dried power. The dried powder was stored at 4°C before use.

Treatment of Rats

In this study, a DMN-induced liver injury model was used because it has been well documented and demonstrated to have good repeatability and reliability.1820) One experiment cycle continued for 4 weeks, and the following procedure was performed immediately prior to the experiments. Male Wistar-derived rats weighing 180–220 g were used. They were provided by The Good Laboratory Practice Centre of Heilongjiang University of Chinese Medicine (Harbin, China), and certificate number was SCXK (Beijing) 2007-0001. The animals had free access to food and water under the standard conditions of humidity (50%±10%), temperature (25±2°C) and 12 h light-dark cycle for 1 week prior to treatment. All protocols were approved by the Animal Use and Care Committee of Heilongjiang University of Chinese Medicine. In this experiment, rats were divided into a treatment group (n=8, treated with 42.45 g crude drug/kg body weight of MSNS each day for 28 d and 10 mg/kg body weight of DMN for 3 consecutive days per week for four weeks); a model group (n=8, treated with oral saline solution each day and 10 mg/kg body weight of DMN for 3 consecutive days per week for four weeks); and a control group (n=8, treated only with oral saline solution each day). After the end of the experimental period, all the animals were sacrificed for tissue sampling and blood collection. The blood samples were allowed to coagulate at 4°C for 30 min. The serum was then separated by centrifugation at 1400×g for 10 min at 4°C. The supernatants were immediately frozen and stored at −80°C until biochemical indicator detection and ultra-performance liquid chromatography/electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC/ESI-Q-TOF-MS) analysis. The livers were removed immediately and used for histopathologic examinations.

Histopathologic Examination and Biochemical Analysis

Necropsies were performed on all the animals at the end of the experiment. The liver specimens were fixed in 10% neutral buffered formalin, dehydrated, and paraffin-embedded. After being sectioned, the thin sections of the tissues were stained with hematoxylin and eosin for histological observation.

Clinical chemistry analysis of serum samples was carried out with an Automatic Analyzer using appropriate kits. The following parameters were tested: alanine aminotransferase (ALT), aspartate aminotransferase (AST), total bilirubin (TBIL), total protein (TP) and albumin (ALB).

Serum Sample Preparation for UPLC/MS Analysis

For UPLC/MS and UPLC/MSn analysis, 500 µL of serum sample was loaded into the centrifuge tube followed by precipitation with 500 µL acetonitrile. Extraction was conducted by vortexing vigorously for 1 min and centrifuged at 16000×g for 10 min, then all the supernatant was transferred to another tube and evaporated at 40°C to dryness. Finally, the residue was reconstituted in 150 µL acetonitrile–water (3 : 7, v/v), then filtered through 0.22 µm filter membrane and clear solution was transferred to autosampler vial.

UPLC/MS Analysis

Chromatographic separations were performed on an ACQUITY UPLC BEH C18 column (1.7 µm, 2.1 mm×100 mm, Waters Corp., Milford, U.S.A.), using a ACQUITY™ UPLC system. Column was maintained at 40°C and eluted with a linear gradient of 2–100%B, where A=water with 0.1% formic acid and B=acetonitrile with 0.1% formic acid. The duration of the gradient program was 0–1.5 min, 2%B; 1.5–3 min, 2–30%B; 3–22 min, 30–80%B; 22–24 min, 80–100%B; 24–30 min wash with 100%B; 30–32, 100–2%B and a 5 min recycle time, at a flow rate of 400 µL/min. Five microliter aliquots of sample were loaded. The column eluent was analyzed by Q-TOF mass spectrometry, using positive ion electrospray ionization (ESI). The column eluent was directed to the mass spectrometer without split.

All MS and MSn mass spectra were acquired on a Micromass Q-TOF Micro™ Mass Spectrometer operating in the positive ion mode. The desolvation gas was set to 750 L/h at a temperature of 350°C. The cone gas was set to 50 L/h and the source temperature was set to 120°C. The capillary and cone voltages were set to 2700 V and 25 V, respectively. The Q-TOF Micro was operated with an alternating collision energy of 6 eV. The data acquisition rate was set to 0.2 s, with a 0.02 s interscan delay. All analyses were acquired using the lockspray to ensure accuracy and reproducibility. Leucine enkephalin was used as the lockmass (m/z 556.2771) at a concentration of 50 fmol/L and a flow rate of 5 µL/min. The mass spectrometric full-scan data were acquired from m/z 80 to m/z 1500 Da with a lockspray frequency of 10 s, and data averaging over 10 scans.

Statistical Analysis

The statistical analyses on UPLC/MS data were performed using Markerlynx Applications Manager Version 4.1 (Waters Corp., Manchester, U.K.). Pareto scaling was applied to the MS integration data prior to principal component analysis (PCA) and orthogonal projection to latent structures discriminate analysis (OPLS-DA). Markers of liver injury were assessed using Student’s t-test.

Clinical chemistry analysis results were expressed as the means±standard deviation. All data were analyzed using the Statistical Package for the Social Sciences for Windows, version 11.0 (SPSS Inc., Chicago, IL, U.S.A.). The groups were compared using ANOVA and multivariate statistical analysis, and statements of statistical significance were based on p<0.05(*).

RESULTS AND DISCUSSION

Serum Biochemical Measurement

Serum AST, ALT activities, ALB, TBIL and TP levels were all changed significantly from 4 weeks after the treatment of DMN (Table 1). MSNS supplementation attenuated the elevation of AST and ALT activities in the rats treated with DMN. A statistically significant restoration in ALB, TBIL and TP levels was observed in rats treated with MSNS compared with rats in model group. These results demonstrated that MSNS had a good performance in recovering the clinical biochemistry in the liver injury model.

Table 1. Effects of Modified Sinisan Supplementation on Selected Parameters of Liver Injury in Rats Induced by Dimethylnitrosamine
ControlDMNMSNS+DMN
ALT (units/mL serum)34.60±5.46102.18±15.28*56.00±14.77#
AST (units/mL serum)155.50±29.65277.82±47.81*196.36±37.60#
Total protein (mg/mL)57.27±3.8437.09±4.85*42.31±5.11#
Total bilirubin (mg/dL serum)1.48±0.3721.89±5.48*17.26±2.93#
Albumin (mg/mL)37.49±1.6022.14±3.80*29.97±4.35#

Each value represents the mean±S.D. for eight rats. *Significant difference compared with normal control group (p<0.05); #Significant difference compared with DMN treatment group (p<0.05).

Histopathology Examinations

Histopathological changes of rat livers were shown in Fig. 1, the liver tissue from DMN-treated rats (Fig. 1B) showed centrilobular necrosis and had more fibrosis, necrotic hepatocytes and degenerative hepatocytes than that of normal control rats (Fig. 1A). DMN-injured rat livers showed histological pattern of fibrotic septa and incorporation of degenerated hepatocytes into pseudolobules. Central veins were surrounded by collagen fibers. MSNS supplementation markedly alleviated the degree of liver fibrosis and significantly reduced the collagen deposition (Fig. 1C).

Fig. 1. Micrographs of Rat Liver Specimens for (A) Normal, (B) Model, and (C) MSNS Treatment Groups

Establishment of Metabolic Fingerprints

Samples of serum were analyzed by LC/MS to determine what changes in metabolome are associated with a deleterious effect of DMN on liver. To optimize the experimental conditions, a preinvestigation was conducted before the full study. Fingerprints of a small batch of serum samples were acquired in positive and negative mode, respectively. Relatively, adequate information of metabolites could be detected in the ESI positive mode. Therefore, full-scan detection was eventually set at ESI positive mode. We investigated the response of analytes after the addition of different concentrations (0.01, 0.1 and 1%) of formic acid into the mobile phase and found that better separation was obtained when 0.1% formic acid was used. Considering resolution and sensitivity meanwhile, the flow rate was optimized at 0.4 mL/min. The flow and temperature of desolvation gas were set at 750 L/h and 350°C, respectively to remove redundant solvent resulting from the high flow rate for mass spectrometer. Ten pooled QC samples were run to “condition” and provide assurance that the system was suitable for metabonomic analysis. During the analytical run, one of these pooled QC samples was interspersed between every five biological samples. After analyzing by UPLC/MS, more than 6000 ions could be detected in serum.

OPLS-DA and PCA Processing of UPLC/MS Data

In our study, OPLS-DA was used to obtain better discrimination between the control and model group and to identify potential biomarkers. The OPLS-DA score plots (Fig. 2A) separated serum samples into two blocks, the control block and the model block, indicating that the liver injury model was successfully reproduced. The corresponding loading plot (Fig. 2B) indicated that differentiating metabolites were attributable to the clustering observed in the score plot. Then S-plot (Fig. 2C) was used to exhibiting the responsibility of each ion for these variations more intuitively. As we can observe in the S-plot, most of the ions were clustered around the origin point, only a few of them scattered in the margin region. Interestingly, just these few ions contributed to the clustering observed in the score plot.

Fig. 2. OPLS Score Plot (A), Loading Plot (B) and S-Plot (C) of the Control and Model Groups

PCA score plot (Fig. 3) separated serum samples into different blocks. Samples subjected to the same treatment were located on the same trajectory. With regard to information analyst of PCA, the control and liver injury groups were significantly divided into two classes. MSNS group developed the distant tendency to the diverse extent of the liver injury group, indicating that MSNS could have some therapeutic effect on the live injury, which was consistent with blood biochemistry examinations. Moreover, some of the samples in the control and MSNS groups clustered, suggesting that the MSNS and the control groups could have similar metabolic profiles.

Fig. 3. Trajectory Analysis of PCA Score Plots for the Liver Injury Rats after Regiments MSNS Treatment

Identification of Liver Injury-Related Metabolites

Metabolites that significantly contributed to the clustering were identified according to a threshold of variable importance in the projection (VIP) values, which could be generated after OPLS-DA processing. Ions with VIP values exceeding two were selected. Then Student’s t-test was used to select potential biomarkers worthy of preferential study in the next step. The critical p-value was set to 0.05 for significantly differential variables. Following the criterion above, 11 significantly differential endogenous metabolites were selected for further study. The possible elemental compositions of the selected metabolites were generated by using the MarkerLynx. The calculated mass, mass deviation (ppm and mDa), speculated formula, double bond equivalent and i-fit value (the isotopic pattern) were calculated according to the selected m/z ions. The smaller mass deviation and lower i-fit value indicate a more accurate elemental composition. The structure information was predicted by comparing the accurate MS and MS/MS fragments with the metabolites searching in databases (http://www.genome.jp, http://www.hmdb.ca, http://www.scripps.edu). The results of identification were shown in Table 2.

Table 2. Identification of Biomarkers in Serum of Liver Injury Rats Detected by UPLC/ESI-Q-TOF-MS
No.VIPObserved [M+H]+RtActual-MProposed compoundProposed structureTrend in model groupa)Trend in MSNS groupb)
19.30188.06863.20187.0633Indoleacrylic acid↓***↑***
27.11166.08612.89165.0790L-Phenylalanine↑***↓***
34.53160.07282.77159.0684Indoleacetaldehyde↓***
44.44303.228815.96302.2246Retinyl ester↓***↑*
54.42165.05260.86164.0473Enol-phenylpyruvate↑**
63.57205.09363.20204.0899L-Tryptophan↓***↑**
73.55182.07960.86181.0739L-Tyrosine↑***↓*
83.35162.10890.69162.1130Carnitine↑***↓***
93.28165.05151.55164.0473Phenylpyruvic acid↑*↓*
102.44289.215711.58288.2089Testosterone↑***↓***
112.17400.338016.06399.3348L-Palmitoylcarnitine↑*

a) Change trend compared with control group. b) Change trend compared with model group. The levels of potential biomarkers were labeled with (↓) down-regulated and (↑) up-regulated (* p<0.05; ** p<0.01; *** p<0.001).

Biological Pathway Analysis of Identified Biomarkers

Datailed analysis of pathways and network influenced by liver injury were performed by network analysis. MetPA (Metabolomics Pathway Analysis) is a user friendly, web-based tool that dedicates to the analysis and visualization of metabolomic data within the biological context of metabolic pathways. It combines powerful pathway enrichment analysis with the analysis of pathway topological characteristics to help researchers identify the most relevant metabolic pathways involved in the conditions under study. Biological pathway analysis revealed that metabolites that were identified together are important for the host response to liver injury and are responsible for phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, tryptophan metabolism, retinol metabolism, tyrosine metabolism, steroid hormone biosynthesis, ubiquinone and other terpenoid-quinone biosynthesis, aminoacyl-tRNA biosynthesis, fatty acid metabolism and lysine degradation (Fig. 4 and Table 3). The detailed construction of the phenylalanine, tyrosine and tryptophan biosynthesis pathway (Fig. 5) with the highest score in rats was generated using the reference map by searching KEGG. The impact (pathway impact value calculated from pathway topology analysis) threshold was set to 0.10, above this threshold, were filtered out as potential pathways. We identified 5 unique pathways (listed in Table 3) to the liver injury animals. Results suggested that these pathways showed marked perturbations over the time–course of liver injury and could contribute to the development of liver injury.

Fig. 4. Summary of Pathway Analysis with MetPA

(1) Phenylalanine, tyrosine and tryptophan biosynthesis; (2) Phenylalanine metabolism; (3) Tryptophan metabolism; (4) Retinol metabolism; (5) Tyrosine metabolism; (6) Steroid hormone biosynthesis; (7) Ubiquinone and other terpenoid-quinone biosynthesis; (8) Aminoacyl-tRNA biosynthesis; (9) Fatty acid metabolism; (10) Lysine degradation.

Fig. 5. Construction of the Phenylalanine, Tyrosine and Tryptophan Biosynthesis Pathways in Rats

The map was generated using the reference map by KEGG (http://www.genome.jp/kegg/). The gray boxes: enzymatic activities with putative cases of analogy in rats.

Table 3. Results from Pathway Analysis with MetPA
No.Pathway nameTotalExpectedHitsRaw p−log(p)Impact
1Phenylalanine, tyrosine and tryptophan biosynthesis40.0933.49E-051.03E+011.00
2Phenylalanine metabolism90.1942.00E-051.08E+010.65
3Tryptophan metabolism410.8835.45E-022.91E+000.17
4Retinol metabolism170.3613.09E-011.17E+000.15
5Tyrosine metabolism420.9016.02E-015.07E-010.14
6Steroid hormone biosynthesis701.5045.85E-022.84E+000.06
7Ubiquinone and other terpenoid-quinone biosynthesis30.0616.29E-022.77E+000.00
8Aminoacyl-tRNA biosynthesis671.4331.70E-011.77E+000.00
9Fatty acid metabolism390.8315.75E-015.54E-010.00
10Lysine degradation200.1611.47E-011.92E+000.00

Note: Total is the total number of compounds in the pathway; the Hits is the actually matched number from the user uploaded data; the Raw p is the original p value calculated from the enrichment analysis; the Impact is the pathway impact value calculated from pathway topology analysis.

Dynamic Changes of Biomarkers

Of the 11 metabolites identified, L-phenylalanine, enol-phenylpyruvate, L-tyrosine, L-palmitoylcarnitine, testosterone, phenylpyruvic acid and carnitine were up-regulated; L-tryptophan, retinyl ester, indoleacetaldehyde and indoleacrylic acid were down-regulated after liver injury. Phenylalanine is an essential amino acid and the precursor for the amino acid tyrosine. Like tyrosine, phenylalanine is the precursor of catecholamines in the body (tyramine, dopamine, epinephrine and norepinephrine). Phenylalanine is highly concentrated in the human brain and plasma. Tyrosine is the precursor for hormones, thyroid, catecholestrogens and the major human pigment, melanin. Tyrosine is an important amino acid in many proteins, peptides and even enkephalins. It has been demonstrated by other investigators that large quantities of phenylalanine and tyrosine were detected in patients with liver injury and these amino acids were specifically released in the patients of liver injury.21,22) Significant correlations were observed between phenylalanine and tyrosine concentrations and total bilirubin levels and albumin levels. Serum phenylalanine and tyrosine concentrations correlated well with imaging and laparotomy findings of liver injury.23) It is also well-established that L-phenylalanine is catabolized in the liver and its metabolism depends on the integrity of liver cell function, liver cell injury could potentially cause an elevated L-phenylalanine concentration in the serum.2428) It has been demonstrated by other investigators that central plasma clearance of amino acids accurately predicts hepatocyte function in patients with liver disease and correlates with clinical outcome.29,30) In the present study, we observed that the serum levels of phenylalanine and tyrosine decreased after DMN administration, indicating that phenylalanine, tyrosine and tryptophan biosynthesis pathway was disturbed.

Tryptophan is an essential amino acid which is the precursor of serotonin. Serotonin is a brain neurotransmitter, platelet clotting factor and neurohormone found in organs throughout the body. Metabolism of tryptophan to serotonin requires nutrients such as vitamin B6, niacin and glutathione. Niacin is an important metabolite of tryptophan. Some disorders of excess tryptophan in the blood may contribute to mental retardation. Assessment of tryptophan deficiency is done through studying excretion of tryptophan metabolites in the urine or blood. Blood may be the most sensitive test because the amino acid tryptophan is transported in a unique way. Increased urination of tryptophan fragments correlates with increased tryptophan degradation, which occurs with oral contraception, depression, mental retardation, hypertension and anxiety states. The amount of tryptophan in the serum was increased significantly after administration of DMN, which might be induced by disorders of tryptophan metabolism.

Retinyl ester is a substrate for lecithin retinol acyltransferase and retinal pigment epithelium-specific 65 kDa protein. Retinyl ester is one of the major retinoid species that are present in the body.31) Retinoids (vitamin A and its metabolites) are required for life and act essentially in maintaining normal cell proliferation and differentiation, vision system, normal male and female reproduction, and a healthy immune.32) Approximately 80% of the body vitamin A is stored in hepatic stellate cells (HSC) with in the lipid droplets as retinyl esters. After liver injury caused by exposure to hepatotoxins, the HSC becomes “activated.” A characteristic feature of HSC activation is the loss of the lipid droplets and the retinoid stores they contain. Leo and Lieber found that there is approximately 5-fold decrease in total hepatic retinol levels with the development of alcoholic hepatitis and nearly 4-fold decrease with the development of cirrhosis.33) The serum contents of retinyl ester decreased obviously, which could potentially explain why retinol metabolism was disrupted. However, other biomarkers were warranted the further mechanistic exploration.

Taken together, our work highlights the importance of metabolomics to elucidate metabolic characters of liver injury and therapeutic effects of MSNS at the global metabolomic levels. We identified 11 potential biomarkers and predicted the metabolites network of liver injury using the validated pattern recognition methods and computational systems analysis. These biomarkers were related to the acutely perturbed pathways of phenylalanine, tyrosine and tryptophan biosynthesis, phenylalanine metabolism, tryptophan metabolism, retinol metabolism and tyrosine metabolism. Of note, our findings also show that MSNS exhibited preventive efficacy against liver injury by adjusting these multiple metabolic pathways towards their normal state. We have first constructed the metabolomic feature profiling and metabolite interaction network of MSNS against liver injury using pattern recognition methods combined with pathways analysis. Taken together, our findings demonstrated the importance of metabolomics as a platform to assess the drug efficacy is highlighted once more.

Acknowledgements

This work was supported by the Grants from Natural Science Foundation of Heilongjiang Province (D201037), Heilongjiang Provincial Administration of Chinese Medicine Project (ZHY10-Z47), Heilongjiang Postdoctoral Fund (LBH-Z09010), Heilongjiang University of Chinese Medicine Research Fund (200909).

REFERENCES
 
© 2013 The Pharmaceutical Society of Japan
feedback
Top