Abstract

With the rapid development of cheminformatics, computational biology, and systems biology, great progress has been made recently in the computational research of Chinese herbal medicine with in-depth understanding towards pharmacognosy. This paper summarized these studies in the aspects of computational methods, traditional Chinese medicine (TCM) compound databases, and TCM network pharmacology. Furthermore, we chose arachidonic acid metabolic network as a case study to demonstrate the regulatory function of herbal medicine in the treatment of inflammation at network level. Finally, a computational workflow for the network-based TCM study, derived from our previous successful applications, was proposed.

1. Introduction

The past century has witnessed tremendous achievements in orthodox medicine and pharmaceutical science as thousands of drugs were created to save human lives throughout the world. Meanwhile, traditional Chinese medicine (TCM), one of the major alternative medicines, actively embraces science and technology, resulting in fruitful clinical outcomes [1]. As an essential part of TCM, Chinese herbal medicine plays an important role in treating patients under the framework of TCM theory, which is characterized by being holistic, systematic, and individualistic, while it is nowadays echoed by the topics of personalized and precise medicine [2–4].

Chinese herbal medicine, originated from the first literature Shennong’s Materia Medica (~220 CE), has accumulated an abundance of healing knowledge in both theory and practice through thousands of years. To be specific, medical materials are classified according to four natures (hot, warm, cool, and cold), five flavours (acrid, sweet, bitter, sour, and salty), and different meridians, which require careful combination in a formula for each individual patient [5]. Here also the issue of medical materials’ compatibility is raised, which further determines the four major components in a formula: principal, associate, assistant, and coordinator [6]. Enriched by tremendous practical data, formulas are developing and evolving in the history with better clinical performance.

To better utilize the hidden wisdom in Chinese herbal medicine, one needs to interpret its property and function using analytical tools well developed in modern science [7]. Such methods include but are not limited to analytical chemistry, molecular pharmacology, animal model, and computational method [8–10]. In this paper, we will discuss from the computational perspective the topic of Chinese herbal medicine addressing the biological networks of complex diseases and try to propose a computational workflow derived from our previous research [11–13].

2. Computational Method

Generally speaking, the computational methods applied in herbal medicine are largely associated with computational chemistry and computational biology because these two disciplines provide researchers with intermediate tools to process chemicals from herbs and targets in pathways. For example, given a small molecule discovered from a herb of interests, one can use chemical similarity searching to find compounds with similar structure and presumably similar bioactivity [14]. Furthermore, docking is a powerful method to quantify the binding energy of a ligand to its target [15]. Based on this method, scientists have developed a web server to identify potential drug targets from multiple therapeutic areas for natural compounds in herbal medicine, which is also known as reverse docking [16–18]. Molecular dynamics is another widely used method to simulate the interaction of a small molecule and a macromolecule, shedding light on the binding kinetics of a biological system at atomic level [19, 20]. For example, Chang et al. integrated this method with docking in the screening of a TCM database against hemagglutinin, a target in influenza virus [21]. On the other hand, we can apply systems biology to study the behaviour of a biological network with and without the treatment of a certain herb in order to simulate the efficacy and side effects of the herb at network level [11, 22]. For example, given the topology and its related parameters, a set of ordinary differential equations (ODEs) can faithfully describe the dynamic property of the network with the perturbation of herbal medicine [11, 23]. However, since ODEs are highly dependent on detailed biological pathways, lack of either topology or parameter information would result in failure of this method. In this case, Boolean network modelling could be a useful alternative [24]. In conclusion, all of these methods are developed to study a biological system (e.g., herb and cell) at certain levels while the combination of them can cover a wide range of scales on both time and space (Figure 1).





3. TCM Compound Database

The major goal of TCM research is to understand the mechanisms of herb-human interplay, which in one modern way that can be translated into the interaction of chemical components and macromolecular targets at molecular level. With the help of cheminformatics and data science, several TCM compound databases of high qualities were curated: Traditional Chinese Medicine Database [25], Chinese Traditional Medicinal Herbs Database [26], TCM Database@Taiwan [27], Traditional Chinese Medicine Systems Pharmacology Database [28], Traditional Chinese Medicines Integrated Database [29], and Herbal Ingredients’ Targets Database [30].

Here are some highlights of several essential TCM databases. Traditional Chinese Medicine Database is one of the pioneers in this field. It offers detailed information for more than 20000 TCM natural compounds, especially for the annotation of related experiments and references [25]. TCM Database@Taiwan claims to be the largest TCM database containing more than 30000 compounds from 352 TCM ingredients. These ingredients were organized in different categories according to the TCM theory, while the molecules under each ingredient were curated in both cdx (2D) and mol2 (3D) formats [27]. Traditional Chinese Medicine Systems Pharmacology Database aims to be a full-stack assistant in computational TCM research, offering not only 3D molecular structures but also physicochemical and pharmacokinetic properties of 12144 compounds. In addition, it provides 3311 potential drug targets and user-tailored networks to show the interaction of drug-target-disease [28].

While most of the researches focus on virtual screening of the databases by various computational tools like chemical similarity search, docking, and molecular dynamics [31–34], drug-likeness analysis was also applied to the TCM databases in order to investigate the pharmaceutical potential of natural compounds. For example, Shen et al. published a series of papers on the topic of drug-likeness analysis of traditional Chinese medicines by comparing the molecular properties and scaffold architectures for drug-like compounds, non-drug-like compounds, and natural compounds from traditional Chinese medicines. Subsequently, they applied machine learning approaches on a TCM database, identifying almost 60% of the molecules as drug-like [35–37]. Furthermore, Xue et al. developed a method with the combination of network topology analysis and cheminformatics measurements to predict the safety for natural compounds of an in-house TCM database. With that, they were able to discover that a promising lead compound Silibinin was surprisingly very similar to a withdrawn drug called Plicamycin [38].

4. TCM Network Pharmacology

The concept of network pharmacology was proposed by Hopkins in 2008 [39], although the idea of multitarget drugs in a network approach has been prevailing since an earlier date [40]. Because the previous paradigm (one gene, one drug, one disease) in drug discovery no longer satisfied the treatment of complex diseases like AIDS, cancer, and neurological disorders, scientists realized the necessity to shift this paradigm into a network-targeted combination therapy with robust efficacy and low toxicity [39, 41].

In the field of TCM network pharmacology, people have begun to address essential questions by utilizing various computational methods. For example, Li et al. proposed a TCM research framework as “phenotype network, biological network, herb network”, with which they developed several useful tools to elucidate the network pharmacology alongside the systematic interpretation of TCM theory and practice [42–45]. On the other hand, using network analysis and molecular docking, Gu et al. introduced multiple computational approaches in the research of herbal pharmacology and efficacy evaluation [46, 47]. Other novel studies in both methods and applications have been accumulating recently with rapid advances in the understanding of herbal medicine [48–54].

Among the studies mentioned above, Zhao and Li developed a systematic approach focusing on three parts of TCM network pharmacology. The first one was drug target identification by relating pharmacological and genomic spaces. The computational framework drugCIPHER was used to infer drug-target interactions in a genome-wide scale [43]. The second one was building on the research of drug-gene-disease relationship. They applied a comodule approach to elucidate drug-disease associations, which was automated by the software called comCIPHER [45]. The third part included several methods that contribute to the analysis of TCM combination in the biological network scenario: for example, the Distance-Based Mutual Information Model (DMIM) to uncover the combination rule of TCM formulas and Network Target Based Identification of Multicomponent Synergy (NIMS) to screen synergistic TCM drug combination [55, 56].

5. TCM Interacting with Biological Networks of Complex Diseases

For the last decade, we have been exploring network-based drug design to address complex diseases like inflammation, HIV, and cancer [57–59]. In order to shift a biological network from a disease state to the normal state, we developed a Multitarget Optimum Intervention (MTOI) method to identify key targets in a network for intervention [60, 61]. Here, the arachidonic acid (AA) metabolic network is taken to illustrate our strategy.

The AA metabolic network contains two major pathways in inflammatory processes with five key enzymes: cyclooxygenase-1/2 (COX1/2), 5-lipoxygenase (5LOX), microsomal prostaglandin E synthase-1 (PGES), and leukotriene A4 hydrolase (LTA4H). Overproduction of two metabolites prostaglandin E2 (PGE2) and leukotriene B4 (LTB4) leads to inflammatory syndromes; for example, PGE2 is highly associated with arthritis while LTB4 is associated with coughs and asthma. Furthermore, the concentration ratio of prostacyclin (PGI2) and thromboxane A2 (TXA2) is linked to drugs’ side effects (Figure 2) [60, 61]. With a dynamic model of 26 ODEs, we faithfully described its metabolic process and successfully screened several multitarget inhibitors [62–64].





Using a combined approach of molecular docking and network simulation, we were able to reveal the regulatory roles of herbal medicine in AA metabolic network [11]. Firstly, we collected 28 anti-inflammatory herbs and selected all the available molecules from Traditional Chinese Medicine Database under some criteria. The structures of key enzymes in AA network were either obtained from Protein Data Bank (PDB) or modelled. Then, we used AutoDock to investigate the binding energy in each complex among 5 proteins and 237 ligands. The docking results were further imported in the network simulation described by 26 ODEs. Finally, the reduction of inflammatory mediators (PGE2 and LTB4) and the ratio of PGI2 / TXA2 as the efficacy and side effects, respectively, were calculated for each herb and formula in the treatment of inflammation.

The network output provided us with the molecular interpretation of the regulatory function of each herb or formula. Concretely, we found that herbs decreased LTB4 production more than PGE2 production, while most nonsteroidal anti-inflammatory drugs (NSAIDs) mainly reduce the production of PGE2. In addition, [PGI2]/[TXA2] ratio of each herb indicated mild side effects. From the simulation of herbal formula, we discovered an interesting phenomenon that herbs in a formula tended to inhibit different enzymes in this network, thus covering most targets to outperform a single herb. Moreover, our method helped to understand the dosage effect of a formula, which could reach the same therapeutic effect at a low dose rather than an individual herb at a high dose. Finally, the simulation results enabled us to design several new herbal combinations for better PGE2 reduction, which was not achieved in the previous collection of herbs.

6. A Computational Workflow for Network-Based TCM Study

On the basis of our previous applications, here we would like to propose a computational workflow for herbal medicine study (Figure 3). In this workflow, we start a project from extracting medical and biological data computationally or manually. These data are usually from various resources like books, papers, and clinical reports. As one of the building blocks in computational herbal medicine, construction of biological networks associated with diseases is vital and inevitable. Fortunately, Kyoto Encyclopedia of Genes and Genomes (KEGG) provides us with a fundamental framework consisting of many essential biological networks [65]. Nevertheless, literature survey is still necessary before one constructs a detailed network with sufficient kinetic parameters. With enough information obtained, we can move on to collect the compounds from TCM databases and the protein structures from PDB. For those proteins without experimental data available, methods like homology modelling can be used to build the 3D structure. In order to quantify the interactions between compounds and proteins, we may apply docking or other methods like molecular dynamics. The predicted dissociation constant of each protein-ligand complex will be imported into the ODE-based simulation to model the dynamics of the biological network. With the outputted data, the efficacy and toxicity of a herb or a formula could be interpreted at a system level. Last but not least, based on the simulation and further validations from wet experiments, one can design new TCM formulas for better performance.





7. Conclusions

Computational studies on Chinese herbal medicine in the scenario of biological networks associated with diseases are a new yet promising direction, because they enable researchers to investigate the pharmacological mechanism of TCM at the molecular and network levels, which in return illuminate scientists to reuse the ancient medical knowledge in a systematic and predictable manner. Although not mature in its initial stage, this field has already gained fruitful outcomes in theory, shedding light on the rigorous applications of herbal medicine and a new paradigm of drug discovery. As a foundation of computational herbal medicine, the interplay of TCM database and biological network generated numerous insightful results during the past decade, providing us with a modern understanding of this complex system. Furthermore, our proposed workflow, including many cutting-edge computational techniques, aims to offer a novel approach to elucidate herbal medicine in silico, providing guides for the in vitro and in vivo experiments.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work has been supported by the Ministry of Science and Technology of China (Grant no. 2016YFA05023032) and the National Natural Science Foundation of China (Grants nos. 81273436 and 21673010).