Evidence-Based Complementary and Alternative Medicine

Volume 2013 (2013), Article ID 136732, 9 pages

http://dx.doi.org/10.1155/2013/136732

Research Article

Replacements of Rare Herbs and Simplifications of Traditional Chinese Medicine Formulae Based on Attribute Similarities and Pathway Enrichment Analysis

Zhao Fang,¹ Meixia Zhang,² Zhenghui Yi,³ Chengping Wen,⁴ Min Qian,¹ and Tieliu Shi¹

¹Center for Bioinformatics and Computational Biology, Shanghai Key Laboratory of Regulatory Biology, Institute of Biomedical Sciences and School of Life Sciences, East China Normal University, Shanghai 200241, China

²Department of Ophthalmology, West China Hospital, Sichuan University, Chengdu, Sichuan 610041, China

³Schizophrenia Program, Shanghai Mental Health Center, School of Medicine, Shanghai Jiao Tong University, Shanghai 200030, China

⁴TCM Clinical Basis Institute, Zhejiang University of Chinese Medicine, Hangzhou, Zhejiang 310053, China

Received 14 September 2012; Revised 24 December 2012; Accepted 26 December 2012

Academic Editor: Shu-Ming Wang

Copyright © 2013 Zhao Fang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A Traditional Chinese Medicine (TCM) formula is a collection of several herbs. TCM formulae have been used to treat various diseases for several thousand years. However, wide usage of TCM formulae has results in rapid decline of some rare herbs. So it is urgent to find common available replacements for those rare herbs with the similar effects. In addition, a formula can be simplified by reducing herbs with unchanged effects. Based on this consideration, we propose a method, called “formula pair,” to replace the rare herbs and simplify TCM formulae. We show its reasonableness from a perspective of pathway enrichment analysis. Both the replacements of rare herbs and simplifications of formulae provide new approaches for a new formula discovery. We demonstrate our approach by replacing a rare herb “Forsythia suspensa” in the formula “the seventh of Sang Ju Yin plus/minus herbs (SSJY)” with a common herb “Thunberg Fritillary Bulb” and simplifying two formulae, “the fifth of Du Huo Ji Sheng Tang plus/minus herbs (FDHJST)” and “Fang Feng Tang” (FFT) to a new formula “Fang Feng Du Huo Tang” (FFDHT).

1. Introduction

Traditional Chinese Medicine is an ancient system used in disease treatments for several thousand years already [1, 2]. Currently, TCM is popular only in Asia, but also in United States, Europe, and other countries as a complementary or alternative medicine [3, 4]. Nearly 100,000 TCM formulae have been recovered [5, 6], each of which normally contains several herbs. Here, “herbs” refer to not only plants, but also animals and minerals with effects of treatments [7].

An herb normally has five attributes: they have nature, taste, channel tropism, functions, and indications [8–10]. Nature refers to a basic characteristic based on the patient’s reaction to an herb, and it has four properties: cold, cool, warm, and hot; taste indicates the actions an herb has on human body and it includes five characters: spicy, sour, sweet, bitter, and salty; channel tropism denotes an herb’s selective therapeutic effects on a certain part of the body, and it includes twelve different organs: heart, liver, spleen, lung, kidney, stomach, large intestine, small intestine, bladder, gallbladder, pericardium, three warmers; functions signify the mechanisms of an herb’s therapeutic effects. For example, a common herb and “Chinese angelica,” has the functions of “harmonizing the blood,” “regulating menses,” and “moistening intestines”; indications describe the diseases or symptoms treated by an herb. Several herbs compose a formula in the order of Master, Adviser, Soldier, and Guild [11], while a later order indicates less importance [12].

A formula usually contains many active compounds. These compounds target many molecules in the cell and work together to increase therapeutic efficacy and reduce adverse effects [6, 11]. Although great efforts have been made, the mechanisms of most formulae are still unknown [6, 13]. Recently, systems biology approaches [14] and TCM informatics [15] have been applied to explore the mechanisms for different formulae, but with a little significant progress. In this study, we explore the pharmacological effects of formulae from a perspective of pathway enrichment analysis.

The availability of some herbs has a rapid decline with their wide usage [16]; for example “Shi-He Ming Yan Wan” [17], a powerful formula, has been used to treat many different diseases. It has been believed that the treatment power of the formula partially comes from a rare herb “rhinoceros horn,” which is widely used in many TCM formulae. Due to poaching and illegal trade, the number of black rhinos has shrunk by 96% since 1970 in Africa. The same situation has happened for “Pholidota,” “Forsythia suspensce,” “Cordyceps sinensis,” and so forth. So to continually use those formulae to treat the diseases, it is urgent to find the common available replacements for those rare herbs with the same effects [18].

Herb switches and simplification in a formula have been tested. For example, “Liu-wei-di-huang” (LWDH) [12], a famous formula with six herbs, was simplified from “Ba-wei-di-huang” (BWDH) by discarding two herbs based on herbs’ attributes. In our study, we adopted the similar concept and simplified TCM formula with a “formula pair” method based on large-scale computation. We then showed the reasonableness of our approach with a perspective of pathway enrichment analysis, same as to the replacements of rare herbs.

2. Materials and Methods

2.1. Data Sources

We collected 4,343 formulae and 6,171 herbs from SIRC/TCM database (http://www.tcm120.com/1w2k/q_pres.asp). If several formulae have a same name but different herbs, we took the classical one (e.g., formula A) as a leader and named others as “the th of formula A plus/minus herbs” (). Disease ontology identities (DOIDs) of diseases treated by formulae were obtained from human disease ontology.

2.2. Score System of Combinational Degree (CD) for Formula Pairs

Combinational degree refers to the similarity between two formulae. The CD was calculated to evaluate combinational degree between two formulae based on the herbs they share and the weights of shared herbs. Assuming there are two formulae and . and are the number of herbs in and , respectively, is the number of shared herbs, herb indicates one of the shared herbs, and and signify the sequential numbers of herb ’s positions in and . The weights of herb in and , named and , are calculated as

The CD of and is calculated as

An example is shown in Table 1.

Table 1: Weight of each herb in and .

Formula
herb 1
herb 2
herb 3
herb 4
herb 5
herb 6
herb 7

1
0.8
0.6
0.4
0.2

1

0.75

0.5
0.25

    includes 5 herbs: herb 1, herb 2, herb 3, herb 4, and herb 5, while includes 4 herbs: herb 2, herb 4, herb 6, and herb 7. Herb 2 and herb 4 are assumed as shared herbs of and . Weight of each herb is shown in this table. So CD of and is calculated as = = 1.475. Then was divided by 7.5 to 0.197 for normalization.

We calculated all formula pairs of 4,343 formulae. The generated values for those CDs are from 0 to 7.5; we then normalized them to for convenient purpose. Combinational degree was used as an evaluation index of the similarity between two formulae. For a pair of formulae with a CD value of 0, we considered these two formulae without combination, which also means that they have no shared herbs. indicates the two formulae as low combination, while as middle combination, and as high combination based on the number of shared herbs in two formulae and the position of each herb in the formulae.

2.3. Disease Similarity Calculation for Each Formula Pair

DOSim [19] is an R package used to calculate the similarity between diseases using DOID based on disease ontology. For 4,343 formulae-treated diseases, only l,769 formulae-treated diseases can be mapped to DOIDs. Assuming that and represent the number of diseases treated by formulae and , respectively, DOSim will generate similarity values for this formula pair. Each similarity value can be denoted as . The similarity of treated diseases between and , named as Sim, can be calculated as follows:

For a pair of formulae with a , we consider that those diseases the two formulae treat have no similarity, and a indicates those diseases with significant similarity for a formula pair, a with high similarity.

2.4. Score System of Attribute Similarities (AS) for Herb Pairs

To evaluate whether the attributes of two herbs are similar or not, a sore system was constructed for the purpose. For the five attributes in a formula, we define each one with a weight of 1. Based on a previous study, the detailed algorithms are as follows.

(1) The four natures of an herb are represented as , , , and , respectively. Different herbs have different natures. If an herb has only one nature, we reward this nature a value of 1 and other natures as 0. If the nature of an herb is “ping,” which means neutral, a value of 0.25 will be assigned to each of the four natures. A slight nature, for example, “slight cold,” will get a value of 0.8; for a severe nature, such as “severe cold,” we assign it a value of 1.2. For two herbs and , the similarity of their natures, named , is calculated as

An example is shown in Table 2.

Table 2: Values of natures for herb and herb .

Cold
Cool
Warm
Hot

Herb
0.25
0.25
0.25
0.25

Herb
0.8
0
0
0

  herb has a nature of “ping” and herb has a nature of “slight cold.” Values of natures are showed in this table. So is calculated as = .

(2) The five tastes are represented as , respectively. According to the similar principal, if an herb has several tastes, we reward each of those tastes a value of 1 and 0 to other unidentified tastes. For a slight taste, such as “slight sour,” a value of 0.8 will be assigned. If the taste of an herb is described as “tasteless,” the value of “sweet” will be added 0.5. If the taste of an herb is depicted as “pucker,” the value of “sour” will get 0.5 bonus. For the two herbs and , the similarity of tastes, named , is calculated as

(3) The twelve channel tropisms are represented as , respectively. According to the similar principal, if an herb has several channel tropisms, we reward each of those channel tropisms a value of 1 and 0 to other unidentified channel tropisms. For the two herbs and , their similarity of tastes, named , is calculated as

(4) The numbers of functions for herbs and are named as and . is the number of shared functions between two herbs. The similarity of the functions between two herbs, named , is calculated as

An example is shown in Table 3.

Table 3: Functions and of herb and herb .

Herb
Function 1
Function 2
Function 3
Function 4
Function 5

Herb
Function 2
Function 4
Function 6
Function 7

  herb has 5 functions and herb has 4 functions. Function 2 and Function 4 are shared functions. So is calculated as = .

(5) The numbers of indications for herbs and are named as and . is the number of shared indications between two herbs. The similarity of the indications between two herbs, named , is calculated as

So, attribute similarity (AS) of herb and herb is calculated as

If herbs and have an AS value greater than 3, they are considered to have high attributes’ similarity.

3. Results

It has been a long history to use herbs to treat diseases. Herbs usually are any part of plants or certain animals with their medicinal effects. Since rare herbs decline rapidly with their wide usage, they face great threat of extinction. To protect those invaluable plants or animals, it is urgent to find the replacement for the formulae with rare herbs. Meanwhile, formulae are not changeless and can be simplified for the cost efficiency and availability. However, the prerequisites for replacements of rare herbs and simplifications of existed formulae are that the new formulae should not change the medical effects compared to original ones. According to this concept, we designed related strategy and used it to computationally detect the possibility of the replacement for all of the formulae we collected.

3.1. Replacements of Rare Herbs

Combinational degree refers to the similarity of two formulae; the smaller the degree, the less similarity between two formulae. Accordingly, the formula pairs with were used for herb replacement. With this criterion, two formulae in the pair share few herbs, thus, avoiding rare herbs shared by the two formulae, since if a rare herb is the shared herb for a formula pair, the herb used to replace the rare herb may be the rare herb itself according to our method. Assuming herb is a rare herb in formulae and herb is a common one in , to replace herb with herb , we built a strict model as follows:

(1),(2), (3).

The replacement of herb with should meet all of the three conditions above. Condition (1) guarantees that herbs and contribute highly similar effects. Condition (3) assures that herbs and have similar weights. Under such circumstance, when herb replaces herb in , the order of Master, Adviser, Soldier, and Guild in can be kept, and the impact of the replacement will be minimized and can be ignored.

Forsythia suspensa is a rare herb protected by law due to its severe decline. According to our model, we replaced it in the formula “the seventh of Sang Ju Yin plus/minus herbs” (SSJY) with a common herb “Thunberg Fritillary Bulb,” which is a component (herb) in another formula “Sang Ke Tang” (SKT). Both SSJY and SKT can treat chronic bronchitis. Herbs in SSJY and SKT are shown in Table 4.

Table 4: Herbs in SSJY, SKT, and SJZBY.

Formula

Herbs

SJYJJS
Mulberry leaf
Chrysanthemum
Forsythia suspense
Reed rhizome

Common hogfennel root
Bitter apricot kernel
Platycodon root
Liquorice

SKT
Mulberry leaf
Fermented soybean
Thunberg fritillary bulb
Radix adenophorae

White mulberry root-bark
Cape jasmine fruit
Bitter apricot kernel
Liquorice

SJZBY
Mulberry leaf
Chrysanthemum
Thunberg fritillary bulb
Reed rhizome

Common hogfennel root
Bitter apricot kernel
Platycodon root
Liquorice

SSJY and SKT, meet the requirements in the replacement model based on the following results:

(1); (2); (3)in SSJY, ; in SKT, ; /; (4).

So, we replaced “Forsythia suspensa” in SSJY with “Thunberg Fritillary Bulb” in SKT. The replacement resulted in a new formula and we named it as “Sang Ju Zhe Bei Yin” (SJZBY, Table 4).

3.2. Formulae Simplifications

To simplify a formula, we also built a model. For a formula pair with and , both of them can treat the same or highly similar diseases, we believe that the high similarity of the two formulae is caused by the shared herbs and the rest herbs with . Under such assumption, two formulae can be simplified to a new formula which only consists of shared herbs and herbs with high attributes’ similarity.

Assuming formulae and are a formula pair with and , has herbs, and has herbs; the number of shared herbs is . For other herbs in and herbs in , we calculate AS for each herb pair, and totally there are pairs. As a result, we then find out those pairs of herbs with their . If there are () herbs in those formula pairs with and there are herbs in and herbs in , the new formula can be generated as follows:

(1)if , new formula will consist of shared herbs and herbs according to the principle of “least herbs.” This formula is named as “new ”; (2)if , new formula will consist of shared herbs and herbs. This formula is named as “new ”; (3)if , assuming is the total weights of herbs in and is the total weights of herbs in ,

(1) if , new formula will be a “new ” according to “higher weights, more importance” concept, the “new ” formula will consist of shared herbs plus herbs;(2) if , new formula will be a “new ”; (3) if , both “new ” and “new ” will be new formulae after discarding those unshared herbs with their . The workflow is shown in Figure 1.

Figure 1: Workflow of simplification.

For example, both formulae “the fifth of Du Huo Ji Sheng Tang plus/minus herbs” (FDHJST) and “Fang Feng Tang” (FFT) can treat rheumatoid arthritis (RA). FDHJST includes 15 herbs (Table 5) and FFT has 8 herbs (Table 5). They have 5 shared herbs: “Radix saposhnikoviae,” “Chinese angelica,” “Radix gentianae macrophyllae,” “Tuckahoe,” “Liquorice”. We calculated AS scores of herb-pairs for the rest unshared herbs and found that was 3.44 and greater than 3. In FDHJST, ; in FFT, . Since the weight of “root of doubleteeth pubescent angelica” in is greater than the weight (0.375) of notopterygium root in , the new formula, named as “Fang Feng Du Huo Tang (FFDHT)”, will contain 5 shared herbs and “Root of doubleteeth pubescent angelica” (Table 5).

Table 5: Herbs in FDHJST, FFT, and FFDHT.

Formula
Herbs

DHJSTJJF
Root of doubleteeth pubescent angelica
Mistletoe
Radix gentianae macrophyllae

Radix saposhnikoviae
Manchurian wildginger
Chinese angelica

Szechwan lovage rhizome
Chinese herbaceous peony
Drying rehmannia root

Bark of eucommia
Radix achyranthis bidentatae
Ginseng root

Tuckahoe
Cassia bark
Liquorice

FFT
Radix saposhnikoviae
Chinese angelica
Radix gentianae macrophyllae

Cassia twig
Notopterygium root
Bitter apricot kernel

Tuckahoe
Liquorice

FFDHT
Radix saposhnikoviae
Chinese angelica
Radix gentianae macrophyllae

Tuckahoe
Root of doubleteeth pubescent angelica
Liquorice

3.3. Pathway Enrichment Analysis for SSJY, SKT, and SJZBY

To further verify that our formula replacement is relevant from biomedical view, we carried out pathway enrichment analysis with those target proteins for each formula. Protein targets of herbs in each formula were obtained from TCMID [10]. For SSJY and SKT, they have 21 shared compounds targeting 78 proteins. We used ClueGO [20] to conduct pathway enrichment analysis for those targets with 0.01 as the threshold of value. The results showed that those targets were enriched in 39 pathways with values less than 0.01. We ranked those pathway according to the value of each enriched pathways in an ascending order and selected top 20 pathways (Table 6) for further mechanism analyses.

Table 6: Top 20 pathways enriched by shared targets of SSJY and SKT with values less than 0.01.

Number
Pathway name
value

1
Leishmaniasis

2
Pathways in cancer

3
Pertussis

4
HTLV-I infection

5
Chagas disease (American trypanosomiasis)

6
T-cell receptor signaling pathway

7
Cytokines and inflammatory response

8
Measles

9
Legionellosis

10
Amoebiasis

11
Free radical-induced apoptosis

12
Cytokine network

13
African trypanosomiasis

14
IL-5 signaling pathway

15
Colorectal cancer

16
Influence of Ras and Rho proteins on G1 to S transition

17
NF-B activation by nontypeable Haemophilus influenzae

18
Rheumatoid arthritis

19
Toll-like receptor signaling pathway

20
Signal transduction through IL1R

Among the top 20 pathways enriched by those 78 shared targets, we found that the pathway of cytokines and inflammatory response ranked 7th and the pathway of free radical-induced apoptosis ranked 11th. Both the two pathways were closely related to chronic bronchitis.

The results showed that there were six targets enriched in the pathway of cytokines and inflammatory response. They are granulocyte-macrophage colony-stimulating factor (GM-CSF), tumor necrosis factor (TNF), interleukin-2 (IL-2), interleukin-4 (IL-4), interleukin-6 (IL-6), and interleukin-10 (IL-10).

Inflammation has been proved to be a central factor to the development and progression of chronic bronchitis [21]. GM-CSF is capable of generating both granulocyte and macrophage colonies from precursor cells, and it has important functions in host responses to external stimuli and in inflammatory conditions [21]. Increased levels of GM-CSF have been found in the epithelium and bronchoalveolar lavage fluid from patients with chronic bronchitis [22]. For other targets, TNF shows statistically significant evidence of association with the susceptibility of chronic obstructive pulmonary disease such as chronic bronchitis [23]. IL-4, IL-6, and IL-10 are cytokines with anti-inflammatory effect [24]. The ingredients of both SSJY and SKT have significant impact on those targets, resulting in effective treatment to chronic bronchitis. Therefore, this pharmacological action of SSJY and SKT could be one of the main mechanisms for chronic bronchitis treatment.

Another pathway closely connected with chronic bronchitis is free radical-induced apoptosis. It has been reported that apoptosis of structural cells in the lung may contribute to the pathogenesis of chronic bronchitis [25, 26]. There is an increase in endothelial and epithelial apoptosis in the lungs of patients with chronic bronchitis which cannot be counterbalanced by an increase in proliferation, resulting in a destruction of lung tissue. Therefore, apoptosis has been proposed to be a therapeutic target [26]. Since targets of SSJY and SKT enrich in this pathway, the two formulae may inhibit the apoptosis in endothelial and epithelial cells by a multitarget effect. This pharmacological action of SSJY and SKT may be another main mechanism for chronic bronchitis treatment. Since those shared targets of SSJY and SKT are enriched in the two pathways, this could be the reason why both SSJY and SKT can treat chronic bronchitis.

The new formula, SJZBY, also includes the same 78 targeted proteins. Therefore, it is reasonable to say that SJZBY should have the similar effect on the treatment of chronic bronchitis. Pathway enrichment analysis for shared targets of those formulae shows the reasonableness of this replacement.

3.4. Pathway Enrichment Analysis for FDHJST, FFT, and FFDHT

Pathway enrichment analysis was also applied to explore the potential mechanism for formula simplification. We collected the potential targets for formulae—FDHJST, FFT, and FFDHTl; they are 182, 133, and 95 proteins, respectively. The results show that targets of FDHJST, FFT, and FFDHT are enriched in 73, 64, and 53 pathways with values less than 0.01, respectively. Same as the above mentioned method, we ranked those pathways according to the value of each enriched pathway in an ascending order and selected top 20 pathways (Tables 7, 8, and 9) for further mechanism analyses.

Table 7: Top 20 pathways enriched by targets of FDHJST with values less than 0.01.

Number
Pathway name
value

1
Pathways in cancer

2
Cytokines and inflammatory response

3
Colorectal cancer

4
Cytokine network

5
Malaria

6
Chagas disease (American trypanosomiasis)

7
Pancreatic cancer

8
Amoebiasis

9
Bladder cancer

10
Leishmaniasis

11
Pertussis

12
Tuberculosis

13
Legionellosis

14
Rheumatoid arthritis

15
Small cell lung cancer

16
Chronic myeloid leukemia

17
Prostate cancer

18
HTLV-I infection

19
Influenza A

20
African trypanosomiasis

Table 8: Top 20 pathways enriched by targets of FFT with values less than 0.01.

Number
Pathway name
value

1
Pathways in cancer

2
Colorectal cancer

3
Cytokines and inflammatory response

4
Cytokine network

5
Prostate cancer

6
Amoebiasis

7
Pancreatic cancer

8
Chagas disease (American trypanosomiasis)

9
Chronic myeloid leukemia

10
Pertussis

11
Leishmaniasis

12
Apoptotic signaling in response to DNA damage

13
Tuberculosis

14
Small cell lung cancer

15
Influence of Ras and Rho proteins on G1 to S Transition

16
p53 signaling pathway

17
Toxoplasmosis

18
Measles

19
Malaria

20
HTLV-I infection

Table 9: Pathways enriched by targets of FFDHT with values less than 0.01.

Number
Pathway name
value

1
Pathways in cancer

2
Colorectal cancer

3
Prostate cancer

4
p53 signaling pathway

5
Pertussis

6
Small cell lung cancer

7
Influence of Ras and Rho proteins on G1 to S transition

8
Endometrial cancer

9
Pancreatic cancer

10
Amyotrophic lateral sclerosis (ALS)

11
Cytokines and inflammatory response

12
Bladder cancer

13
Tuberculosis

14
Amoebiasis

15
Apoptosis

16
HTLV-I infection

17
Cytokine network

18
RB tumor suppressor/checkpoint signaling in response to DNA damage

19
Apoptotic signaling in response to DNA damage

20
Chronic myeloid leukemia

In those top 20 pathways enriched by targets of FDHJST, we found that the pathway of cytokines and inflammatory response were closely connected with RA. The result showed that the value of this pathway ranked second in the top 20 pathways. Further research showed that there were thirteen targets of FDHJST enriched in this pathway. They are GM-CSF, interferon beta (IFNB), interferon-gamma (IFNG), transforming growth factor beta-1 (TGFB1), TNF, interleukin-1 alpha (IL-1A), IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IL-13.

A previous report has confirmed that antagonism of GM-CSF represents a novel therapeutic approach for a variety of autoimmune-mediated inflammatory diseases, including RA [27]. IFNG can stimulate the production of chemokines and is a powerful activator of mononuclear phagocytes; IFNG has also been tried in immune-mediated diseases such as RA [28]. IFNB has shown antirheumatic potential [28]. TNF is proven to be expressed at high levels in rheumatoid joint tissue, where they contribute significantly to inflammation and articular destruction. TNF is the first cytokine to be fully validated as a therapeutic target for RA [29]. TGFB1 is highly expressed in joints in RA and is considered to be a regulator of anti-inflammation in RA [30]. Enhanced expression of TGFB1 protein has been detected in RA synovia and it may be related to the active pathological changes in RA synovia including synoviocyte hyperplasia, inflammatory cell infiltration, sublining angiogenesis, and granuloma formation [31].

For those interleukins, IL-6 and IL-8 can be found in RA pathway in KEGG pathway annotation [32]. IL-10, also known as human cytokine synthesis inhibitory factor, is an anti-inflammatory cytokine and has been proposed to treat RA in clinical practice because of its capacity to inhibit cellular immunity and deactivate macrophages [33]. In summary, considering the proteins of compounds from FDHJST target, the pathway enrichment analysis demonstrates the potential mechanisms of the formula on RA treatment.

Among those top 20 enriched pathways by targets of FFT, the pathway of cytokines and inflammatory response was closely connected with RA. The result showed that this pathway ranked third among the top 20 pathways according the value. There were ten targets of FFT enriched in this pathway. They are GM-CSF, IFNB, TNF, IL-2, IL-4, IL-5, IL-6, IL-8, IL-10, IL-13, all of which are the same targets of FDHJST. This result indicates that the potential mechanism of FFT to treat RA is highly similar to that of FDHJST.

After simplification, the resulting new formula, FFDHT, was also enriched in the pathway of cytokines and inflammatory response with the value ranked 11th among the top 20 pathways. Seven targets were enriched in this pathway. They are GM-CSF, TNF, IL-2, IL-4, IL-6, IL-8, IL-10, all of which are also the targets of FDHJST or FFT. Accordingly, it is reasonable to state that FFDHT inherits the pharmacological effects of original ones and has a highly potential effect the treatment of RA, which shows the reasonableness of this simplification.

4. Discussion

Many herbs used in Traditional Chinese Medicine are endangered, such as tiger bone used to treat rheumatism. Its widely usage results in the rapid decline of tigers with the poaching and illegal trade, which push tigers to extinction [17]. Another example is Cordyceps sinensis; the huge commercial demand of Cordyceps sinensis for its powerful tonic function has led to the excessive harvest and its dramatic decline [34]. Therefore, those invaluable herbs are hard to be obtained currently, and the practices in TCM face a great challenge to use those related formulae to treat the diseases. Finding the common replacements for those invaluable rare herbs is in high demand. In this study, we proposed a method to replace rare herbs with common available ones and showed its reasonableness from a perspective of pathway enrichment. The case study suggests the applicability for the replacement of rare herbs, which opens the gate for wide implementation in the field and could have the profound impact on this field.

Moreover, we also proposed a method to simplify formulae based on the similar rationale. A new formula can be formulated with “less herbs but same effect” concept to the original one. Pathway enrichment analysis also shows the reasonableness of the simplification. Our approaches provide an alternative way to reformulate those traditional prescriptions.

Although herbs have been widely used for thousands of years, most of their targets are still unclear and the mechanisms underling their effects remain unknown. And that has strongly prevented the modernization of traditional Chinese Medicine. For example, in the method of score system of attributes’ similarities for herb pairs, we found that “rhinoceros horn” and “Buffalo Horn” have high attributes’ similarity (). Replacement of “rhinoceros horn” with “Buffalo Horn” had been used in clinic, we would like to explore the rationale from a perspective of pathway enrichment. However, no targets of the two herbs have been inferred. And similar situation between “walnut kernel” and “Cordyceps sinensis,” we also found that “walnut kernel” has a potential to replace “Cordyceps sinensis” since the AS between them is 3.2. However, no targets have been identified yet. With the progress of pharmacological research on herbs, it can be anticipated that more and more targets will be identified; we believe that more replacements can be carried out based on our approaches and the potential mechanisms behind those replacements can be explored from a perspective of pathway enrichment.

In this work, both replacement of rare herbs and simplification of formulae were computationally tested; our approaches provide an alternative way for new TCM formulation and mechanism inference. To fully verify our method and test the effects of those new formulae, more preclinical experiments need to be conducted. By the combination of in silico and web lab approaches, we expect that the modernization of Traditional Chinese Medicine will be speeded up; thus, people will benefit from this progress.

Acknowledgment

The authors are grateful to Ms. Bingxing Lu for her help in paper preparation. This work was supported by the National 973 Key Basic Research Program (Grant no. 2010CB945401 and 2012CB910400), the National Natural Science Foundation of China (Grant no. 31171264, 31071162, 31000590, and 81171272), and the Science and Technology Commission of Shanghai Municipality (11DZ2260300).

References

J. Zhao, P. Jiang, and W. D. Zhang, “Molecular networks for the study of TCM pharmacology,” Briefings in Bioinformatics, vol. 11, no. 4, pp. 417–430, 2009.

S. Li, Z. Q. Zhang, L. J. Wu, X. G. Zhang, Y. D. Li, and Y. Y. Wang, “Understanding ZHENG in traditional Chinese medicine in the context of neuro-endocrine-immune network,” IET Systems Biology, vol. 1, no. 1, pp. 51–60, 2007.

E. Chan, M. Tan, J. N Xin, S. Sudarsanam, and D. E. Johnson, “Interactions between traditional Chinese medicines and Western therapeutics,” Current Opinion in Drug Discovery and Development, vol. 13, no. 1, pp. 50–65, 2010.

F. Cheung, “TCM: made in China,” Nature, vol. 480, no. 7378, pp. S82–S83, 2011.

L. Wang, G. B. Zhou, P. Liu et al., “Dissection of mechanisms of Chinese medicinal formula Realgar-Indigo naturalis as an effective treatment for promyelocytic leukemia,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 12, pp. 4826–4831, 2008.

Z. N. Wen, Z. J. Wang, S. Wang et al., “Discovery of molecular mechanisms of traditional Chinese medicinal formula Si-Wu-Tang using gene expression microarray and connectivity map,” PLoS ONE, vol. 6, no. 3, Article ID e18278, 2011.

J. Shi, Y.-L. Liu, Y.-X. Fang, G.-Z. Xu, H.-F. Zhai, and L. Lu, “Traditional Chinese medicine in treatment of opiate addiction,” Acta Pharmacologica Sinica, vol. 27, no. 10, pp. 1303–1308, 2006.

B. Xiao, O. Tao, J. Luo, and Y. Wang, “Relationships between properties and functional targets of Chinese herbs,” Journal of Chinese Integrative Medicine, vol. 9, no. 7, pp. 789–793, 2011.

S. Lukman, Y. He, and S. C. Hui, “Computational methods for traditional chinese medicine: a survey,” Computer Methods and Programs in Biomedicine, vol. 88, no. 3, pp. 283–294, 2007.

R. C. Xue, Z. Fang, M. X. Zhang, et al., “TCMID: traditional Chinese medicine integrative database for herb molecular mechanism analysis,” Nucleic Acids Research, vol. 41, no. 1, pp. D1089–D1095, 2012.

S. K. Yan, J. Zhao, S. S. Dou, et al., “Methodology of modernization research in traditional Chinese medicine based on systems biology and network biology,” Chinese Journal of Natural Medicines, vol. 7, no. 4, pp. 249–259, 2009.

S. Li, B. Zhang, D. Jiang, Y. Wei, and N. Zhang, “Herb network construction and co-module analysis for uncovering the combination rule of traditional Chinese herbal formulae,” BMC Bioinformatics, vol. 11, no. 11, article S6, 2010.

P. Tian, “Convergence: where west meets east,” Nature, vol. 480, no. 7378, pp. S84–S86, 2011.

H. Kitano, “Systems biology: a brief overview,” Science, vol. 295, no. 5560, pp. 1662–1664, 2002.

F. Rui, “Progress in TCM Informatics,” Chinese Journal of Information on TCM, vol. 16, no. 1, pp. 2–7, 2009.

J. D. Li, “System study of endangerment medicinal animals in China,” journal of Chinese Materia Medica, vol. 26, no. 11, pp. 728–730, 2001.

D. Graham-Rowe, “Biodiversity: endangered and in demand,” Nature, vol. 480, no. 7378, pp. S101–S103, 2011.

M. K. Ang-Lee, J. Moss, and C.-S. Yuan, “Herbal medicines and perioperative care,” Journal of the American Medical Association, vol. 286, no. 2, pp. 208–216, 2001.

J. Li, B. S. Gong, X. Chen et al., “DOSim: an R package for similarity between diseases based on disease ontology,” BMC Bioinformatics, vol. 12, no. 1, article 266, 2011.

G. Bindea, B. Mlecnik, H. Hackl et al., “ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks,” Bioinformatics, vol. 25, no. 8, pp. 1091–1093, 2009.

M. Shen, R. Vermeulen, R. S. Chapman et al., “A report of cytokine polymorphisms and COPD risk in Xuan Wei, China,” International Journal of Hygiene and Environmental Health, vol. 211, no. 3-4, pp. 352–356, 2008.

B. Balbi, C. Bason, E. Balleari et al., “Increased bronchoalveolar granulocytes and granulocyte/macrophage colony-stimulating factor during exacerbations of chronic bronchitis,” European Respiratory Journal, vol. 10, no. 4, pp. 846–850, 1997.

P. J. Castaldi, M. H. Cho, M. Cohn et al., “The COPD genetic association compendium: a comprehensive online database of COPD genetic associations,” Human Molecular Genetics, vol. 19, no. 3, pp. 526–534, 2010.

S. M. Opal and V. A. DePalo, “Anti-inflammatory cytokines,” Chest, vol. 117, no. 4, pp. 1162–1172, 2000.

L. Segura-Valdez, A. Pardo, M. Gaxiola, B. D. Uhal, C. Becerril, and M. Selman, “Upregulation of gelatinases A and B, collagenases 1 and 2, and increased parenchymal cell death in COPD,” Chest, vol. 117, no. 3, pp. 684–694, 2000.

I. K. Demedts, T. Demoor, K. R. Bracke, G. F. Joos, and G. G. Brusselle, “Role of apoptosis in the pathogenesis of COPD and pulmonary emphysema,” Respiratory Research, vol. 7, no. 1, article 53, 2006.

A. L. Cornish, I. K. Campbell, B. S. McKenzie, S. Chatfield, and I. P. Wicks, “G-CSF and GM-CSF as therapeutic targets in rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 5, no. 10, pp. 554–559, 2009.

J. Van Holten, C. Plater-Zyberk, and P. P. Tak, “Interferon-β for treatment of rheumatoid arthritis?” Arthritis Research, vol. 4, no. 6, pp. 346–352, 2002.

P. C. Taylor and M. Feldmann, “Anti-TNF biologic agents: still the therapy of choice for rheumatoid arthritis,” Nature Reviews Rheumatology, vol. 5, no. 10, pp. 578–582, 2009.

Y. Sugiura, T. Niimi, S. Sato et al., “Transforming growth factor beta-1 gene polymorphism in rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 61, no. 9, pp. 826–828, 2002.

L. Dai, Z. Q. Ye, M. A. Tang, and P. D. Yin, “Expression of TGFB1 in rheumatoid synovia and its relationship with synovial pathological changes,” Chinese Journal of Rheumatology, vol. 4, no. 6, pp. 357–360, 2000.

H. Ogata, S. Goto, K. Sato, W. Fujibuchi, H. Bono, and M. Kanehisa, “KEGG: Kyoto encyclopedia of genes and genomes,” Nucleic Acids Research, vol. 27, no. 1, pp. 29–34, 1999.

E. W. St Clair, “Interleukin 10 treatment for rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 58, no. 1, pp. I99–I102, 1999.

D. R. Yang, “Regional and ecological distribution of Cordyceps sinensis,” Journal of Anhui Agricultral University, vol. 26, pp. 143–148, 1999.