Evidence-Based Complementary and Alternative Medicine

Evidence-Based Complementary and Alternative Medicine / 2013 / Article

Research Article | Open Access

Volume 2013 |Article ID 136732 | https://doi.org/10.1155/2013/136732

Zhao Fang, Meixia Zhang, Zhenghui Yi, Chengping Wen, Min Qian, Tieliu Shi, "Replacements of Rare Herbs and Simplifications of Traditional Chinese Medicine Formulae Based on Attribute Similarities and Pathway Enrichment Analysis", Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 136732, 9 pages, 2013. https://doi.org/10.1155/2013/136732

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

Academic Editor: Shu-Ming Wang
Received14 Sep 2012
Revised24 Dec 2012
Accepted26 Dec 2012
Published21 Feb 2013

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 [810]. 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.


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

10.80.60.40.2
10.750.50.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.


Cold CoolWarmHot

Herb 0.250.250.250.25
Herb 0.8000

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.


Herb Function 1Function 2Function 3Function 4Function 5

Herb Function 2Function 4Function 6Function 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.


Formula Herbs

SJYJJS Mulberry leafChrysanthemumForsythia suspense Reed rhizome
Common hogfennel rootBitter apricot kernelPlatycodon rootLiquorice
SKT Mulberry leaf Fermented soybean Thunberg fritillary bulb Radix adenophorae
White mulberry root-barkCape jasmine fruitBitter apricot kernelLiquorice
SJZBY Mulberry leaf Chrysanthemum Thunberg fritillary bulb Reed rhizome
Common hogfennel root Bitter apricot kernelPlatycodon rootLiquorice

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.

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).


Formula Herbs

DHJSTJJF Root of doubleteeth pubescent angelica MistletoeRadix gentianae macrophyllae
Radix saposhnikoviae Manchurian wildginger Chinese angelica
Szechwan lovage rhizomeChinese herbaceous peonyDrying rehmannia root
Bark of eucommia Radix achyranthis bidentataeGinseng root
TuckahoeCassia barkLiquorice
FFT Radix saposhnikoviae Chinese angelica Radix gentianae macrophyllae
Cassia twig Notopterygium rootBitter apricot kernel
TuckahoeLiquorice
FFDHT Radix saposhnikoviae Chinese angelicaRadix gentianae macrophyllae
TuckahoeRoot of doubleteeth pubescent angelicaLiquorice

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.


NumberPathway name value

1Leishmaniasis
2Pathways in cancer
3Pertussis
4HTLV-I infection
5Chagas disease (American trypanosomiasis)
6T-cell receptor signaling pathway
7Cytokines and inflammatory response
8Measles
9Legionellosis
10Amoebiasis
11Free radical-induced apoptosis
12Cytokine network
13African trypanosomiasis
14IL-5 signaling pathway
15Colorectal cancer
16Influence of Ras and Rho proteins on G1 to S transition
17NF- B activation by nontypeable Haemophilus influenzae
18Rheumatoid arthritis
19Toll-like receptor signaling pathway
20Signal 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.


NumberPathway name value

1Pathways in cancer
2Cytokines and inflammatory response
3Colorectal cancer
4Cytokine network
5Malaria
6Chagas disease (American trypanosomiasis)
7Pancreatic cancer
8Amoebiasis
9Bladder cancer
10Leishmaniasis
11Pertussis
12Tuberculosis
13Legionellosis
14Rheumatoid arthritis
15Small cell lung cancer
16Chronic myeloid leukemia
17Prostate cancer
18HTLV-I infection
19Influenza A
20African trypanosomiasis


NumberPathway name value

1Pathways in cancer
2Colorectal cancer
3Cytokines and inflammatory response
4Cytokine network
5Prostate cancer
6Amoebiasis
7Pancreatic cancer
8Chagas disease (American trypanosomiasis)
9Chronic myeloid leukemia
10Pertussis
11Leishmaniasis
12Apoptotic signaling in response to DNA damage
13Tuberculosis
14Small cell lung cancer
15Influence of Ras and Rho proteins on G1 to S Transition
16p53 signaling pathway
17Toxoplasmosis
18Measles
19Malaria
20HTLV-I infection


NumberPathway name value

1Pathways in cancer
2Colorectal cancer
3Prostate cancer
4p53 signaling pathway
5Pertussis
6Small cell lung cancer
7Influence of Ras and Rho proteins on G1 to S transition
8Endometrial cancer
9Pancreatic cancer
10Amyotrophic lateral sclerosis (ALS)
11Cytokines and inflammatory response
12Bladder cancer
13Tuberculosis
14Amoebiasis
15Apoptosis
16HTLV-I infection
17Cytokine network
18RB tumor suppressor/checkpoint signaling in response to DNA damage
19Apoptotic signaling in response to DNA damage
20Chronic 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

  1. 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. View at: Publisher Site | Google Scholar
  2. 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. View at: Publisher Site | Google Scholar
  3. 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. View at: Google Scholar
  4. F. Cheung, “TCM: made in China,” Nature, vol. 480, no. 7378, pp. S82–S83, 2011. View at: Publisher Site | Google Scholar
  5. 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. View at: Publisher Site | Google Scholar
  6. 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. View at: Publisher Site | Google Scholar
  7. 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. View at: Publisher Site | Google Scholar
  8. 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. View at: Publisher Site | Google Scholar
  9. 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. View at: Publisher Site | Google Scholar
  10. 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. View at: Publisher Site | Google Scholar
  11. 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. View at: Google Scholar
  12. 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. View at: Publisher Site | Google Scholar
  13. P. Tian, “Convergence: where west meets east,” Nature, vol. 480, no. 7378, pp. S84–S86, 2011. View at: Publisher Site | Google Scholar
  14. H. Kitano, “Systems biology: a brief overview,” Science, vol. 295, no. 5560, pp. 1662–1664, 2002. View at: Publisher Site | Google Scholar
  15. F. Rui, “Progress in TCM Informatics,” Chinese Journal of Information on TCM, vol. 16, no. 1, pp. 2–7, 2009. View at: Google Scholar
  16. J. D. Li, “System study of endangerment medicinal animals in China,” journal of Chinese Materia Medica, vol. 26, no. 11, pp. 728–730, 2001. View at: Google Scholar
  17. D. Graham-Rowe, “Biodiversity: endangered and in demand,” Nature, vol. 480, no. 7378, pp. S101–S103, 2011. View at: Google Scholar
  18. 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. View at: Google Scholar
  19. 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. View at: Google Scholar
  20. 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. View at: Publisher Site | Google Scholar
  21. 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. View at: Publisher Site | Google Scholar
  22. 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. View at: Google Scholar
  23. 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. View at: Google Scholar
  24. S. M. Opal and V. A. DePalo, “Anti-inflammatory cytokines,” Chest, vol. 117, no. 4, pp. 1162–1172, 2000. View at: Google Scholar
  25. 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. View at: Google Scholar
  26. 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. View at: Google Scholar
  27. 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. View at: Publisher Site | Google Scholar
  28. 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. View at: Publisher Site | Google Scholar
  29. 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. View at: Publisher Site | Google Scholar
  30. 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. View at: Publisher Site | Google Scholar
  31. 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. View at: Google Scholar
  32. 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. View at: Publisher Site | Google Scholar
  33. E. W. St Clair, “Interleukin 10 treatment for rheumatoid arthritis,” Annals of the Rheumatic Diseases, vol. 58, no. 1, pp. I99–I102, 1999. View at: Google Scholar
  34. D. R. Yang, “Regional and ecological distribution of Cordyceps sinensis,” Journal of Anhui Agricultral University, vol. 26, pp. 143–148, 1999. View at: Google Scholar

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.


More related articles

2041 Views | 824 Downloads | 7 Citations
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.