A comprehensive review of integrative pharmacology-based investigation: A paradigm shift in traditional Chinese medicine.

Over the past decade, traditional Chinese medicine (TCM) has widely embraced systems biology and its various data integration approaches to promote its modernization. Thus, integrative pharmacology-based traditional Chinese medicine (TCMIP) was proposed as a paradigm shift in TCM. This review focuses on the presentation of this novel concept and the main research contents, methodologies and applications of TCMIP. First, TCMIP is an interdisciplinary science that can establish qualitative and quantitative pharmacokinetics-pharmacodynamics (PK-PD) correlations through the integration of knowledge from multiple disciplines and techniques and from different PK-PD processes in vivo. Then, the main research contents of TCMIP are introduced as follows: chemical and ADME/PK profiles of TCM formulas; confirming the three forms of active substances and the three action modes; establishing the qualitative PK-PD correlation; and building the quantitative PK-PD correlations, etc. After that, we summarize the existing data resources, computational models and experimental methods of TCMIP and highlight the urgent establishment of mathematical modeling and experimental methods. Finally, we further discuss the applications of TCMIP for the improvement of TCM quality control, clarification of the molecular mechanisms underlying the actions of TCMs and discovery of potential new drugs, especially TCM-related combination drug discovery.

Introduction

Traditional Chinese medicine (TCM) is one of the oldest traditional medical systems and is characterized by personalized, holistic, and multicomponent therapeutic strategies. TCM plays a crucial role in modern medical and health care with its unique theory and philosophy and has attracted worldwide attention in recent years. TCM and natural plants are important sources of modern drugs, and approximately 45% of today's bestselling drugs originated from natural products or their derivatives, including some modern “blockbuster” drugs such as artemisinin, aspirin, digitoxin, and morphine. Youyou Tu was inspired by the detailed methods of artemisia usage described in Ge Hong's ancient book to discover artemisinin (a famous antimalarial drug) and won the 2015 Nobel Prize in Physiology or Medicine. Characterized by its holistic concept, the TCM formula is the main form of its clinical application. Growing clinical evidence has shown that TCM has favorable clinical efficacy based on randomized double-blind experiments and placebo-controlled studies, such as with MaXingShiGan-YinQiaoSan, arsenic, and retinoic acid, and exerts definite curative effects against H1N1 influenza and acute promyelocytic leukemia,.

Over the past decades, TCM has been widely embraced by systems biology and its various data integration approaches to promote the modernization and internationalization of TCM. Recently, some novel concepts and methods have been proposed to systemically identify the active constituents and reveal the molecular basis of TCM in the treatment of various diseases or syndromes, such as network pharmacology of TCM (TCMNP),, system pharmacology of TCM (TCMSP), Fangjiomics, and chinmedomics. However, the following accumulating problems have emerged. (1) The chemical profile of TCM formulas is greatly influenced by the different herb origins, production areas, harvest times, processing methods and preparation technologies, etc. However, a large amount of chemical information has been collected from chemical databases, which lack information on the qualitative identification and quantitative determination of TCM formulas. (2) It has widely been accepted that absorption, distribution, metabolism, excretion and pharmacokinetics (ADME/PK) are critical features for determining which chemical components are likely to be active and what kind of mode of action they may adopt to achieve their therapeutic effects. TCMNP and TCMSP only focus on direct action based on the absorbed prototype constituents interacting with therapeutic-related targets but often ignore regulation by the intestinal flora and constituent–constituent interactions. (3) A variety of ingredients in different amounts play a synergistic role in the treatment of various diseases. However, there is no appropriate method to determine the synergistic effects based on the multidimensional associations of multiple components, targets, pathways and bioactivities at the cellular, tissue and organ levels both in vitro and in vivo. Notably, there is poor consistency active substances between different in vitro and in vivo pharmacological studies, and there is also a lack of a dose–effect relationship between the multiple components of TCM and their bioactivities, which may hinder access to reliable and high-level experimental evidence to verify the pharmacological findings in vivo.

To address these problems, an integrated pharmacology-based research strategy, a paradigm shift in traditional Chinese medicine (TCMIP), was proposed and developed, which may inject new vitality into the modernization and internationalization of TCM,. In this review, we introduce the concept and routine research strategy of TCMIP, as well as summarize its recent developments. We also focus on the main research contents and data resources and methodologies of TCMIP in addition to its applications in the research field of TCM, its major challenges and future directions, which may create a paradigm shift in TCM studies.

The concept of TCMIP

TCMIP, under the guidance of TCM theories, is an interdisciplinary science that comprehensively explores the interactions between the multiple constituents of TCM and the body at multiple levels, such as the molecular, cellular, tissue, organ and animal levels. TCMIP focuses on the construction and evaluation of the multidimensional associations among the chemical and ADME/PK profiles of TCM, the disease–syndrome–formula association network, and pharmacological actions to qualitatively and quantitatively assess the PK–PD correlations in vivo, which can be subsequently applied to identify the bioactive constituents, reveal the pharmacological mechanisms and elucidate the compatibility of TCM. The basic research framework of TCMIP is illustrated in Fig. 1. First, high-throughput and high-sensitivity analytical techniques in combination with chemical databases and virtual prediction methods are employed to fully characterize the chemical and ADME/PK profiles of TCM formulas. According to the ADME/PK properties, the active substances of TCM are divided into three forms: the absorbed prototypes (AP), the absorbed metabolites (AM) and the unabsorbed constituents (UAC). Second, TCMIP utilizes omics-based data and imbalanced biomolecular network analysis to comprehensively reveal the molecular basis associated with the initiation and progression of diseases and TCM syndromes. Then, three active TCM substances are associated with the therapeutic mechanisms through three modes of action, including the direct effects of the AP and AM interacting with the target network (PMTN), the indirect effects of the AP, AM and UAC regulating the gut microbiota (PMCGM), and the auxiliary effects of constituent–constituent interactions based on the actions of the AP and AM with ADME/PK-related enzymes or transporters (CCIs-PMPK). Third, a multidimensional association network is constructed for the qualitative PK–PD correlation to screen out key active constituents (KACs) of TCM, critical molecular targets (CMTs) and critical pharmacological effects (CPEs). Moreover, TCMIP characterizes the quantitative PK–PD correlations of the KACs of TCM and CPEs using in vitro PK–PD dynamic complex models (IV–PK/PD–DCMs) and AI-related algorithms. Finally, knock-in/out of constituents and genes will be utilized to validate the pharmacological actions of KACs and CMTs, as well as the synergetic mechanism of multiple KACs and CMTs.

Figure 1

The basic research framework of TCMIP. Step 1: System characterization of the chemical profile of a TCM using LC–MS in combination with chemical databases. Step 2: Systemic identification of the ADME/PK profile of a TCM and uncovering the three action forms of the active substances of the TCM (AP, AM, UAC). Step 3: Comprehensive investigation of the therapeutic mechanisms based on the three action modes (PMTN, PMCGM, and CCIs-PMPK). Step 4: Establishment of the qualitative and quantitative PK–PD correlation by multidimensional association network and mathematical modeling. Step 5: Verification of the KACs and CMTs by knock-in/out of constituents and genes. AP: absorbed prototypes; AM: absorbed metabolites; UAC: unabsorbed constituents; PMTN: the direct interactions of AP and AM with the therapeutics target network; PMCGM: indirect effects of AP, AM and UAC regulating gut microbiota; CCIs-PMPK: auxiliary effects of constituent–constituent interactions based on AP and AM Action with ADME/PK-related enzymes or transporters; KACs: key active constituents; CMTs: critical molecular targets.

Therefore, TCMIP embodies the unity of system theory and reductionism in vivo and in vitro. TCMIP is also a database- and algorithm-dependent research strategy. Big data and artificial intelligence technologies will drive the rapid development of TCMIP. Moreover, there are several differences between TCMIP and similar concepts (TCMNP, TCMSP, etc.), as shown in Fig. 2. First, TCMIP focuses more on qualitative and quantitative analysis to characterize the reliable and precise chemical profiles of TCM, while TCMNP and TCMSP often obtain TCM chemical information from chemical databases. Second, TCMIP comprehensively elucidates the mechanism of TCM, including the three forms of active substances (AP, AM and UAC) and three modes of action of TCM (PMTN, PMCGM and CCIs-PMPK). In particular, gut microbiota and ADME/PK-related enzymes/transporters, in addition to the PMTN of AP and AM, are also important action modes of TCM based on PMCGM and CCIs-PMPK,. However, TCMNP and TCMSP mainly focus on active substances (AP) in their single form and single action modes (PMTN). In addition, the multidimensional association network in TCMIP has more complex interaction types, such as constituent–target interactions, constituent–gut microbiota interactions, constituent–constituent interactions, gut microbiota–target interactions, and target–target interactions. Additionally, in TCMIP, IV–PK/PD–DCMs and the knock-in/out of constituents are characteristically used as key techniques to establish the quantitative PK–PD correlation and verify the synergistic actions of multiple KACs of TCM.

Figure 2

Comparison of TCMIP and other pharmacologies from the following aspects: chemical material basis, the forms of active substances, the action modes of TCM, multidimensional integration for PK–PD correlations, and experimental validation.

Main research contents of TCMIP

Systemic characterization of the chemical and ADME/PK profiles of TCM formulas and identification of the three forms of active substances

The chemical profile is the pharmacodynamic material basis of TCM and is a very complex chemical system including hundreds or thousands of known and unknown components, ranging from trace to very large amounts of ingredients and from small molecules to large molecules. After oral administration, the chemical constituents contained in TCM undergo a series of ADME/PK processes in vivo, which is critical to understand the pharmacodynamic material basis and the therapeutic mechanisms of TCM formulas. Recently, with the rapid development of chromatographic techniques and mass spectrometry, ultra-performance liquid chromatography coupled with high-resolution mass spectrometry (UPLC–HRMS) has been presented as a high-throughput and high-sensitivity platform to systemically characterize the chemical and ADME/PK profiles of TCMs. After oral administration, the chemical constituents contained in plasma, urine, feces and target organs are qualitatively identified and quantitatively determined at different times to obtain AP that are absorbed into the body, AM that are metabolized by the gut microbiota and cytochrome P450 (CYP450) enzymes, and UAC through excretion of feces. In particular, tissue and cellular pharmacokinetics are the key points that focus on ADME combinational models in vitro to prepare samples for ADME/PK profiling in accordance with the varieties and concentrations of AP and AM in vivo to reveal the molecular mechanism related to the ADME/PK process and constituent–constituent interactions of TCM.

However, current ADME/PK studies have paid insufficient attention to UAC and drug concentrations in target tissues and cells, including the varieties and concentrations of AP and AM. The complexity and microconstituents of TCMs present a major technological and economic challenge in ADME/PK studies due to the lack of chemical standards, inadaptability and time consumption of analytical technologies. Fortunately, in silico PK prediction provides a quick, inexpensive, and high-throughput technique that could be used to obtain the in vivo ADME/PK profiles of TCMs more conveniently. However, the existing in silico models have been established based on training datasets to predict the influence of selected physicochemical properties relevant to PK on the basis of the properties of a large number of individual components determined using high-throughput experimental assays. The current in silico PK models are unsuitable for TCM-related research because of the multiple complex ingredient–ingredient interactions. Although we developed a TCM-ADMEpred method as a novel strategy for poly-pharmacokinetic prediction of TCM according to the theory that chemical components with similar structures often show similar pharmacokinetic properties, many open questions and unexplored limitations still need to be considered in future studies.

Comprehensive investigation of the three action modes of TCM in TCMIP

TCMIP utilizes omics-based data and imbalanced biomolecular network analysis to comprehensively reveal the molecular basis associated with the initiation and progression of diseases and TCM syndromes. In addition, TCMIP determines the efficacy of TCM formulas and the corresponding pharmacological mechanisms via multilevel network analysis and experimental validation. TCMIP will establish the in vivo relationship between the three forms of active substances contained in TCM (AP, AM, UAC) and the therapeutic molecular basis of TCM through the corresponding action modes of PMTN, PMCGM, and CCIs-PMPK. Various TCM formulas have different chemical bases, active substances, action modes and therapeutic molecular networks, which are characterized by the synergistic effects of TCM based on multiple components, modes and targets.

TCMIP reveals the mechanisms underlying the direct effects of PMTN in target tissues and cells

PMTN is the main mode of action of TCM, of which AP and AM directly interact with the corresponding therapeutic targets, including receptors, enzymes, ion channels, transporters, nucleic acids, etc., which further form a constituent–target regulatory network, as shown in Fig. 3. The main modes of constituent–target interactions are as follows22, 23, 24, 25, 26, 27: AP or AM combine with the binding sites of transporters, AP or AM bind with receptors, AP or AM are converted into active ingredients to exert their effects by enzymatic catalysis, AP or AM control the “gating” of ion channels, etc. In order to systemically verify the interaction between the constituents of TCMs (AP, AM) and the molecular targets of a disease/syndrome, the molecular characteristics of TCM therapeutics are revealed by omics, big data and artificial intelligence technologies. Moreover, using virtual prediction and experimental verification, constituent–target interactions can be obtained, and the composition, target, pathway and pharmacological action of the multidimensional complex association network can be further constructed. On this basis, TCMIP emphasizes the construction and analysis of the multidimensional complex association network of PMTN to identify candidate bioactive compounds and elucidate the mechanisms of action of TCM formulas.

Figure 3

TCMIP uncovers the action modes of PMTN between the constituents of TCM (AP, AM) and the molecular targets of disease. ①: AP or AM as transport substrates combine with the binding sites of transporters to cause conformational changes and complete the process of absorption or expulsion. ②: AP or AM as ligands bind with receptors based on the following three steps: primary recognition of the receptor, orientation and change in the structural conformation, and physical binding. ③: AP or AM are converted into active ingredients to exert their effects by the substrate binding to the enzyme to convert the substrate into a product for release. ④: AP or AM control the “gating” of ion channels to open or close and influence the effects of inorganic ions. ⑤: AP or AM interfere with or block the synthesis of nucleic acids by bacteria, viruses and tumor cells to destroy their proliferation. PMTN: the direct interactions of AP and AM with the therapeutics target network; AP: absorbed prototypes; AM: absorbed metabolites.

However, there are still very few high-quality and high-level PMTN studies on TCM because of the following points. First, there are no in-depth studies on the ADME/PK of TCMs. Although TCMs have certain pharmacological effects in vivo after oral administration, the varieties and drug concentrations of the AP and AM in target tissues or cells are still unclear. Therefore, the existing PMTN studies are inconsistent with the active ingredients in vivo. In particular, in traditional TCM pharmacological research, crude herbal extracts are usually directly added to cell or organ culture systems for in vitro pharmacological evaluation, while the ADME/PK process is notably ignored, resulting in the active ingredients in the crude drugs not being the true effective components in vivo. In addition, constituent–target interactions are mainly based on virtual prediction without experimental verification. Furthermore, there is no dose–effect relationship between multiple components and the bioactivity of TCMs. Thus, we need to investigate the ADME/PK of TCMs in depth. Cellular pharmacokinetics has been proposed as a novel concept to reveal the intracellular localization and dynamic process of AP and AM to judge the possibility of the constituent–target interactions of TCMs in the future. Additionally, the constituent–target interactions and regulatory mechanisms should be validated by high-throughput target fishing techniques or knock-in/out constituents or genes, etc. Moreover, in order to study the constituent–target interactions and dose–effect relationships between multiple constituents of TCMs and their bioactivities in depth, it is urgent to establish IV–PK/PD–DCMs in tissues and cells.

TCMIP elucidates the mechanisms underlying the indirect effects of PMCGM in the gut microbiota

The gut microbiota has gained increasing attention because it not only affects gastrointestinal physiology but also performs some basic functions in the immunological, metabolic, structural and neurological landscapes of the human body. PMCGM promotes the paradigm shift from the direct effects of PMTN to the structural and functional correlations between different tissues and between tissues and their microenvironment, which is of great help to the pathophysiology of complex diseases and the treatment of TCM. Recently, PMCGM has been considered an important mode of TCM treatment through the regulation of the gut microbiota. As shown in Fig. 4, TCMIP elucidates the action approaches for PMCGM of the constituents (AP, AM and UAC) of the TCM–gut microbiota–host axis based on the following aspects30, 31, 32, 33: the constituents of TCM directly or indirectly regulate the composition of gut microbiota, directly or indirectly modulate the metabolism of the gut microbiota by altering the composition of microbiota, the metabolites of the gut microbiota enter into body and regulate host function, as well as regulation of the host intestinal barrier, etc.

Figure 4

TCMIP elucidates the action approaches for PMCGM of the constituents (AP, AM and UAC) of TCM–gut microbiota–host. ①: The constituents of TCM directly regulate the composition of the gut microbiota by promoting the growth of beneficial microbiota or selectively inhibiting the growth of harmful microbiota. ②: The constituents of TCM indirectly regulate the composition of the gut microbiota by altering intermediate factors, such as the pH of the gastrointestinal tract. ③: The constituents of TCM modulate the metabolism of gut microbiota indirectly by altering the composition of the microbiota. ④: The constituents of TCM modulate the metabolism of the gut microbiota indirectly by means of increasing or reducing the activities of enzymes related to the gut microbiota. ⑤: The metabolites of the gut microbiota enter the body and regulate host function. ⑥: The TCM or gut microbiota regulate the host intestinal barrier to prevent gut microbiota or harmful substances from entering the body. PMCGM: indirect effects of AP, AM and UAC regulating gut microbiota; AP: absorbed prototypes; AM: absorbed metabolites; UAC: unabsorbed constituents.

TCMIP focuses on the investigation of the gut microbiota to understand the pharmacological effects and molecular basis of TCM therapeutics, deciphers the mystery of the oral bioavailability of TCMs, and explains the compatibility theory of TCM, etc.. Growing evidence shows that TCMs may interact with various gut microbiota to exhibit their therapeutic effects on a variety of diseases, including diabetes, obesity, and cancer,,. TCMIP is expected to become a new way to study the human gut microbiota by employing sequencing technology and metagenomic and bioinformatics methods. Then, gut microbiota depletion and fecal transplantation strategies are utilized to confirm the relationship between the gut microbiota and the TCM,. However, it is very difficult to clarify the causal relationship between the constituents of the TCM and gut microbiota, and new virtual prediction and high-throughput experimental methods need to be developed. The interaction among TCMs, the gut microbiota and the host is a complex and dynamic process that provides a huge challenge for the establishment of simulation models in vitro. Thus, establishing PK–PD dynamic complex models of PMCGM in vitro may be an urgent task to evaluate the therapeutic effects of TCMs during the regulation of gut microbiota.

TCMIP clarifies the mechanisms underlying the auxiliary effects of CCIs-PMPK

Although growing evidence shows that herb–herb interactions may be widespread in TCM formulas, CCIs-PMPK have received very little attention compared with PMTN and PMCGM,. At the micro level, CCIs-PMPK are the basis of herb–herb interactions and is able to alter the ADME properties of active or toxic constituents, including their oral bioavailability, protein binding capacity, blood–brain barrier permeability and half-life, etc., which may be associated with the risks and benefits of TCM formulas. To the best of our knowledge, there are two types of transporter- and enzyme-based CCIs to indicate the existence of CCIs-PMPK, as shown in Fig. 5. Transporter-based CCIs-PMPK are divided into activation or inhibition of activity or expression, blocking substrate binding site(s), inhibiting ATP hydrolysis, and altering the integrity of the cell membrane, etc. Enzyme-based CCIs-PMPK include irreversible inhibition, competitive inhibition, noncompetitive inhibition, uncompetitive inhibition to form a substrate–enzyme–inhibitor complex, etc.

Figure 5

TCMIP clarifies the mechanisms of CCIs-PMPK based on the following aspects. Transporter-based CCIs-PMPK: ①: Activation or inhibition of the activity or expression of transport proteins; ②: Competitively, noncompetitively or allosterically blocking the substrate binding site(s); ③: Inhibiting ATP hydrolysis; and ④: Altering the integrity of the cell membrane. Enzyme-based CCIs-PMPK: ⑤: Irreversible inhibition by irreversible competition with the substrate for the same binding site(s); ⑥: Competitive inhibition through reversible binding of the inhibitor to the enzyme; ⑦: Noncompetitive inhibition by binding enzymes and substrates with two sequences to form substrate–enzyme–inhibitor complexes; ⑧: Uncompetitive inhibition by the enzyme–substrate complex binding to form a substrate–enzyme–inhibitor complex; and ⑨: Inhibiting the enzymes by uncompetitive inhibition. CCIs-PMPK: auxiliary effects of constituent–constituent interactions based on AP and AM action with ADME/PK-related enzymes or transporters.

In order to carry out CCI-PMPK studies, PK and PD experiments of different constituent–constituent and herb–herb combinations are first investigated to determine the ADME/PK profiles of KACs of the TCM and PD profiles related to the TCM therapeutic, confirm the existence of obvious constituent–constituent and herb–herb interactions, and establish the PK–PD correlation of the CCIs-PMPK. After that, efflux transporter and CYP450 enzymes are further employed to clarify the molecular mechanism of the constituent–constituent and herb–herb interactions. In particular, physiology-based pharmacokinetic (PBPK) modeling is an appropriate method for the high-throughput evaluation of membrane transporters and CYP450-mediated CCIs. Transcriptome analysis has been applied to better understand CCIs and address the difficulties caused by the complex components of TCM formulas. However, CCIs-PMPK in TCM are more complex than the study of the phenotypes and mechanisms of the CCIs, and the relationship between them can explain their efficacy or toxicity. Similar to PMCGM, virtual prediction methods and IV–PK/PD–DCMs of CCIs-PMPK will be developed for the high-throughput identification of CCIs and in-depth study of the relationship between CCIs and the benefits or risks of TCMs.

TCMIP constructs the qualitative association of PK–PD and identifies KACs, CMTs and CPEs by establishing a multidimensional network

The three types of active substances and the three action modes of TCM offer a good starting point for the establishment of the multidimensional association among the chemical and ADME/PK profiling of the TCM formulas, their regulatory targets and molecular profiling, pharmacological actions and clinical phenotypes, TCM theories, etc. There are many types of networks involved in the complex interactions between TCM formulas and the body, including constituent–constituent interaction networks, constituent–target interaction networks, constituent–gut microbiota interaction networks, protein–protein interaction networks, and gut microbiota–metabolite interaction networks, etc. On this basis, a multidimensional association complex network can be constructed from different standpoints, ranging from molecule–molecule interactions to cell–cell interactions, from unipartite networks to bipartite networks (constituent–target interactions) or even tripartite networks (constituent–protein–phenotype interactions). Thanks to the rapid advances of network science and high-throughput technologies, bioinformatics algorithms and computational biology approaches have driven the development of a complex multidimensional association network of TCMIP. Network-based computational algorithms, such as network centrality, network controllability and network propagation, could be used to characterize the topological structure and heterogeneity of the network to better understand and simulate the dynamics and functions of a multidimensional association network. In particular, computational network-based approaches, such as deep-learning-based methods, show a strong capacity to extract hidden information and provide new insights. A multidimensional association network can be used to identify KACs of TCMs, mine CMTS-related diseases and reveal CPEs of TCM therapy.

However, the complexities of TCMs pose a great challenge to TCMIP research on multidimensional associations. The existing methods of multidimensional association networks are mainly suitable for PMTN patterns and rarely involve PMCGM and CCIs-PMPK. Therefore, new network-based computational algorithms and approaches have been developed for PMCGM and CCIs-PMPK. In addition, the multidimensional association network involves multiple links of the PK–PD correlation without effective correlations and cross-validation computational and experimental methods. Thus, some cross-validation methods, such as IV–PK/PD–DCMs, will be developed into functional universal methods for the analysis and validation of multidimensional association networks, systemically revealing the role of TCMs ranging from the macro to micro level.

TCMIP establishes quantitative PK–PD correlations between the KACs of TCM formulas and the CPEs of the corresponding therapeutics

On the basis of a multidimensional association network, quantitative PK–PD correlations are an effective method to characterize the correlation, dynamic and quantitative description between the KACs of the TCM formulas and the CPEs of the therapeutics. Recently, quality markers (Q-markers), pharmacokinetic markers (PK-markers) and biomarkers have been deemed the KACs of TCM formulas and the CPEs of therapeutics and have become the current frontier topics of quality control, therapeutic evaluation, mechanistic exploration and synergistic effects of TCM formulas45, 46, 47. Thus, it is of great significance to establish quantitative PK–PD correlations between KACs and CPEs, which can improve the quality control level of TCMs, improve the curative effects, control the risks, and promote personalized medication, etc.

In the present review, we propose the following workflow to establish a quantitative PK–PD correlation between KACs and CPEs, as shown in Fig. 6. First, IV–PK/PD–DCMs will be developed with good stability, high sensitivity, high intelligence and high consistency to the human body, which reflects the mechanism of internal and external correlation of system. A number of TCM products of different quality can be used to obtain different ADME/PK profiles and different molecular and pharmacological profiles using IV–PK/PD–DCMs. After the standardized pretreatment of these PK and PD data, the multi-index weight of each index contribution is calculated by using all standardized data, and the quantitative PK–PD correlation between the KACs and CPEs is preliminarily constructed based on mathematical modeling methods, such as gray relational analysis, artificial neural networks, support vector machines, and genetic algorithms. In particular, the concept of optimal design should be introduced into quantitative PK–PD correlation models to achieve effective sampling to minimize the parameter estimation error. Finally, samples with knockout/in KACs were obtained to further verify the biological activity and the contribution that can be attributed to the whole activity. Through PK–PD modeling, the overall evaluation of a TCM can be performed, and a number of quantitative indicators closely related to efficacy can be determined to clarify the material basis of the biological effects from an overall perspective of TCM.

Figure 6

Establishment of the quantitative PK–PD correlation from the following aspects. First, a large number of TCM samples with different qualities were designed to obtain PK or PD values in vivo and in vitro with significant differences. Second, all these data were preprocessed (data standardization) in order to eliminate the incompatibility of different indexes due to different dimensions, and then all of these standardized data were used to construct the PK–PD correlation model by using machine learning, such as ANNs, SVMs, genetic algorithms or a combination method. Third, through PK–PD modeling, the overall evaluation of a TCM can be performed, and a number of quantitative and qualitative indicators closely related to the efficacy can be determined to clarify the material basis of the biological effects for the whole drug. Fourth, knockout/in of constituents will be carried out to verify the biological activity and confirm the contribution of the KACs to the overall activity. ANNs: artificial neural networks; SVMs: support vector machines; KACs: key active constituents.

Facing the complexity of the interactions between the TCM formulas and the human body, TCMIP integrates the existing data resources, computational models and experimental methods. In this section, the data resources, computational algorithms and tools as well as experimental methods of TCMIP are summarized based on the following five aspects: (1) chemical profiling of TCMs; (2) ADME/PK profiling of TCMs; (3) the three action modes of TCMs; (4) multidimensional association networks of qualitative PK–PD correlation; and (5) quantitative PK–PD correlations (Table 1).

Table 1

Current data resources, computational models and experimental methods for TCMIP.

Research content	Main source and database	Main algorithms and computational software	Main experimental method	Limitation
System characterization of the chemical profile of TCM	• Chemical database 1. ChemSpider (v2020.0.18.0) (http://www.chemspider.com/)2. PubChem (https://pubchem.ncbi.nlm.nih.gov/)• TCM comprehensive database 1. ETCM (http://www.tcmip.cn/ETCM/index.php/Home/Index/)2. ChEMBL (v27) (https://www.ebi.ac.uk/chembl/)3. TCMSP (v2.3) (http://tcmspw.com/tcmsp.php)4. BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/index.php)5. TCMID (v2.03) (http://www.megabionet.org/tcmid/)6. TCM-Mesh (http://mesh.tcm.microbioinformatics.org/)7. CHEM-TCM (http://www.chemtcm.com/)• Mass spectrometry database 1. MassBank (v2.1.4) (http://www.MassBank.jp/)2. Wiley MSforID3. NIST (v2018) (http://chemdata.nist.gov/)	/	HPLC, 2D-LC, GC, MS, NMR, LC–MS, GC–MS, MS–NMR, UV, IR	• Incomplete information on TCMs in the databases; • Lack of cross-references with classical databases; • Lack of integration with experimental data; • Structural confirmation of unknown compounds.
System characterization of the ADME/PK profile of TCM	• TCM comprehensive database 1. ETCM2. TCMSP	• ADME prediction software 1. GastroPlus 2. QikProp3. VolSurf+ 4. ADMET Predictor5. SwissADME 6. DataWarrior7. ACD/Percepta, 8. MetaDrug.• ADME fitting software 1. NONMEM2. Phoenix WinNonlin3. ADAPT 4. Kinetica5. MONOLIX	• ADME model: Caco-2 cell model, primary hepatocytes, drug–plasma protein binding and BBB permeability models, MDCK-MDR1 cell model • Microfluidic chip technology • Organ-on-chip system • HPLC, GC, LC–MS, GC–MS • Isotope tracer technology	• Lack of PK model for the multiconstituent interactions of TCM; • Lack of information on the metabolism of TCMs; • Difficult to detect and identify unknown metabolites.
Analyzing and validating the three action modes of TCM	• Drug target databases 1. DrugBank (v5.1.5) (http://redpoll.pharmacy.ualberta.ca/drugbank/)2. STITCH (v5.0) (http://stitch.embl.de/)3. TTD (v2020) (http://db.idrblab.net/ttd/)4. BindingDB (v2020) (http://www.bindingdb.org/bind/index.jsp)• TCM comprehensive database 1. ETCM 2. TCMIP 3. TCMSP 4. BATMAN-TCM 5. CHEM-TCM6. TCM-Mesh• Omics-related database 1. MetaboAnalyst (v4.0) (https://www.metaboanalyst.ca/)2. MaxQuant (v1.6.3.0) (https://www.maxquant.org/)3. UniProt (https://sparql.uniprot.org/)4. TCGA (https://cancergenome.nih.gov/abouttcga/overview)5. ICGC (v2020) (https://icgc.org/)6. KEGG (https://www.genome.jp/kegg/)7. Reactome (https://reactome.org)• Intestinal flora database 1. VMH (https://www.vmh.life/)2. Gmrepo (https://gmrepo.humangut.info)3. GgutMEGA (http://gutmega.omicsbio.info)4. BIO-ML• Small molecule chemistry database 1. PubChem2. Zinc (http://zinc.docking.org/)• Protein database PDB (http://www.rcsb.org/)	• Molecular docking related software 1. AutoDock 2. eHiTS3. FlexX 4. Glide5. GOLD 6. LigandFit7. Surflex8. Similarity ensemble approach9. SuperPred 10. SwissTargetPrediction11. ChemMapper 12. INVDOCK13. TarfisDock 14. DDI-CPI15. IDTarget• Chemical software 1. ChemDraw2. ChemBio3D• Drug–drug interaction software 1. Micromedex® Drug-Reax2. Micromedex® Healthcare Series3. Drug Interactions Facts®4. Lexi-Interact®5. Pharmavista®6. EpocratesRx®7. MediQ®8. Drug interaction checker®• KEGG pathway analysis DAVID • GO analysis	• Omics technologies • Small-molecule probe technology • Biochip technology • Microfluidic technology • Intestinal flora research methods: 16S rDNA sequencing, metagenomic sequencing, metabolomics• CCI research methods: extracorporeal liver system, drug transport cell model	• Target prediction software has limited accuracy; • The update speed of the database is slow; • Complexity and variety of docking software; • Lack of large-scale, high-throughput screening and analysis combinatorial chips; • Lack of a TCM database of intestinal bacteria metabolism; • Lack of database and simulation software for TCM constituent–constituent interactions.
Construction and analysis of a multidimensional association network	• Disease databases 1. HPO (https://hpo.jax.org/app/)2. OMIM (https://omim.org/)3. MalaCards (https://www.malacards.org/)4. DisGeNET (https://www.disgenet.org/)• Protein interaction databases 1. String (v11.0) (https://string-db.org/)2. HAPPI (http://bio.informatics.iupui.edu/HAPPI/)3. HINT (v2020) (http://hint.yulab.org)4. OPHID (v2.9) (http://ophid.utoronto.ca)5. MINT (v2012) (https://mint.bio.uniroma2.it/)6. HPRD (http://hprd.org/index_html)7. DIP (http://dip.doe-mbi.ucla.edu)8. PDZBase (http://icb.med.cornell.edu/services/pdz/start)9. IntAct (https://www.ebi.ac.uk/intact/)• TCM comprehensive database 1. ETCM 2. TCMSP 3. TCMID4. BATMAN-TCM 5. TCM-Mesh• Pathway analysis 1. KEGG 2. Reactome	• Network visualization software 1. Cytoscape2. NAVIGaTOR3. Matascape4. Pajek• Network-related computing methods 1. Restricted Neighborhood Search Clustering2. Molecular Complex Detection3. Random walks with restart4. Hyperlink-Induced Topic Search5. Ingenuity Pathway Analysis6. ScaffoldGraph7. Network propagation8. Shortest path9. Dijkstra algorithm• Pathway analysis DAVID • GO analysis • Network analysis methods 1. Node connectivity analysis2. Node compactness analysis3. Node mediation analysis	/	• Lack of information about interaction type and directionality in the network; • The algorithm of the network is biased; • The biological network of TCM is intricate and difficult to analyze.
Establishment and validation of a quantitative PK–PD correlation	/	• PK–PD correlation analysis 1. Artificial neural networks2. Least angle regression3. Gray relational analysis4. Canonical correlation analysis5. Least squares support vector machines6. Hierarchical clustering analysis	• Microfluidic technology to build PK–PD model • Biochip technology to build PK–PD model • HPLC, GC, MS, LC–MS, GC–MS, NMR • Gene knock-out technology • Constituent knockout technology	• Difficulties in determining PK marker constituents and PD effect indicators; • It is difficult to truly simulate the process of drugs in human body; • Incomplete constituent knockout technology of TCM formulas.

BATMAN-TCM, bioinformatics analysis tool for molecular mechanism of traditional Chinese medicine; BBB, blood–brain barrier; BIO-ML, Broad Institute-Open Biome Microbiome Library; CHEM-TCM, chemical database-traditional Chinese medicine; DAVID, Database For Annotation, Visualization And Integrated Discovery; DIP, Database of Interacting Proteins; ETCM, Encyclopedia of Traditional Chinese Medicine; GC, gas chromatography; Gmrepo, Data Repository for Gut Microbiota; GO, Gene Ontology; gutMEGA, Gut MEtaGenome Atlas; HAPPI, human annotated protein–protein interaction; HINT, high-quality INTeractomes; HPLC, high-performance liquid chromatography; HPO, Human Phenotype Ontology; HPRD, Human Protein Reference Database; ICGC, International Cancer Genome Consortium; IR, infrared spectroscopy; KEGG, Kyoto Encyclopedia of Genes and Genomes; MINT, Molecular Interaction Database; MS, mass spectrometry; NIST, National Institute of Standards and Technology; NMR, nuclear magnetic resonance; OMIM, Online Mendelian Inheritance in Man; OPHID, Online Predicted Human Interaction Database; PDB, Protein Data Bank; TCGA, The Cancer Genome Atlas; TCMID, Traditional Chinese Medicine Integrated Database; TCMSP, Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform; TD-LC, two-dimensional LC; TTD, Therapeutic Target Database; UV, ultraviolet spectroscopy; VMH, Virtual Metabolic Human.

Chemical profiling of TCMs

Data resources

ChemSpider and PubChem are commonly used chemical component databases that contain the molecular structures, physicochemical properties and spectral information of chemical components. However, the chemical components of TCMs collected from foreign databases are limited. Currently, there are some comprehensive TCM databases that can be used to obtain chemical information on TCM components, such as encyclopedia of traditional Chinese medicine (ETCM), ChEMBL, TCMSP, bioinformatics analysis tool for molecular mechanism of traditional Chinese medicine (BATMAN-TCM), traditional Chinese medicine integrated database (TCMID), TCM-Mesh and chemical database-traditional Chinese medicine (CHEM-TCM). ETCM is a comprehensive database established by our team that contains 402 herbal medicines, 3959 formulas and 7284 ingredients contained in the Pharmacopoeia of the People's Republic of China (2015 version) as well as information from the Fourth National Survey on Chinese Materia Medica Resources. The chemical information could cross search with that of CHEMBL and PubChem. In particular, the ETCM database evaluates the druggability of each constituent by calculating the pharmacokinetic parameters according to the models in the Pipeline Pilot absorption, distribution, metabolism, excretion and toxicity (ADMET) collection. However, these databases also have some limitations. In these databases, some chemical components do not have 3D structures, which makes further experiments, such as molecular docking, inconvenient. Moreover, these databases rarely provide information on the metabolites of the prototype components after metabolism in vivo. In addition, we need to enrich the database of TCM ingredients and establish good cross-references among the databases.

Analytical methods

In recent years, a variety of technologies, such as liquid chromatography-mass spectrometry (LC–MS), high-performance liquid chromatography (HPLC), gas chromatography (GC), mass spectrometry (MS), nuclear magnetic resonance (NMR), ultraviolet spectroscopy (UV) and infrared spectroscopy (IR), have been developed and applied to describe the chemical profiles of TCMs. In particular, LC–MS seems to be the most promising technology. LC–MS is suitable for the identification of alkaloids, flavonoids, terpenes, fatty acids, organic acids, polysaccharides and other TCM constituents54, 55, 56, 57. The combination of ultrahigh-performance liquid chromatography (UHPLC) and high-resolution mass spectrometry (HRMS), such as UHPLC coupled with quadrupole time-of-flight tandem mass spectrometry (UHPLC–Q-TOF-MS/MS) and UHPLC–Q Exactive hybrid quadrupole-orbitrap high-resolution accurate mass spectrometry (UHPLC–Q-Orbitrap HRMS), facilitates the high-throughput detection of multiple components in TCMs. In addition, the combined strategy of NMR and LC–MS is a powerful tool for structural identification and promotes the discovery of new active components in TCMs57, 58, 59. Volatile components, such as the volatiles of olive oil and the volatile oil of turmeric, are more suitable for detection by GC–MS,. The structures of complex polysaccharides can be analyzed and characterized by chemical and instrumental methods, such as Fourier transform infrared spectroscopy (FT-IR), HPLC, gel permeation chromatography, monosaccharide composition, methylation analysis, and NMR62, 63, 64.

Nevertheless, the analysis of the chemical profiles of TCMs are not an easy task, especially in the study of TCM prescriptions, and the diversity of ingredients and mutual influence between multiple components lead to difficulty in component analysis. In particular, the identification of unknown compounds has baffled researchers. Therefore, the penetration of multiple disciplines and technologies is required to ensure the establishment of the chemical profile of TCMs.

ADME/PK profiling of TCMs

ETCM is a comprehensive TCM database that contains information on ADMET BB, ADMET solubility, ADMET hepatotoxicity probability, and ADMET CYP2D6 probability of its ingredients. The TCMSP database can also be used to query oral bioavailability (OB), half-life (HL), Caco-2, blood–brain barrier (BBB), etc. These databases make it easier for researchers to understand the AMDE properties of TCMs. Moreover, we should be aware that these databases contain a limited number of components and ADME parameters. Moreover, the ADME parameters provided by these databases are calculated based on certain models that need to be verified by in vivo experiments. Therefore, the experimental information needs to be integrated into the database to ensure that the information provided by the database is consistent with the process of a drug in vivo.

Computational algorithms and tools

In the study of PK, ADME parameters can be simulated in silico according to drug likeness. The commonly used prediction software packages include GastroPlus™,, QikProp, VolSurf+, etc. GastroPlus™ is pharmacokinetic and pharmacodynamic (PK/PD) simulation software based on physiological models that can simulate intravenous, oral, eye, nasal and pulmonary administration routes and is known as the gold standard among similar software programs. For constituents with clear structures and efficacies, HPLC and LC–MS are usually used to detect the drug concentration in the blood at different time points in vivo, and then the ADME parameters are fitted according to the PK fitting software NONMEM, Phoenix Winnonlin, adapt, etc. In addition to calculating the pharmacokinetic parameters, NONMEM can also estimate errors, especially to separate individual errors from total errors. However, what draws more of our attention is that TCMs have the characteristics of complex multiple components and targets, which is quite different from the traditional PK model of a single component, and most of the PK curves of TCM components do not conform to the conventional compartment model. Therefore, it is necessary to establish a model that is more suitable for the multicomponent interaction of TCMs.

Experimental methods

Numerous tissue- and cell-based in vitro models of ADME already exist, including the Caco-2 cell model, primary hepatocyte model, drug–plasma protein binding model, BBB osmotic model,, and MDCK-MDR1 cell model,. Additionally, the organ-on-a-chip system is based on microfluidic technology and aims to simulate the main physiological characteristics or functions of human organs and faithfully reflect human body functions in vitro. To date, organ chips such as intestinal, liver, kidney, heart, and lung chips have been successfully constructed. The drug concentration in vivo is determined by the content changes of the active ingredients in the TCM in vivo, and analytical techniques usually include HPLC, GC, LC–MS, and GC–MS. However, the prescription is the main clinical application form of TCMs, and biological effect methods are often used to study prescriptions with unclear ingredients. Moreover, isotope tracer technology is often used to study tissue distribution and PK, which is characterized by high sensitivity, location observation of drug distribution and the determination of drug target organs.

Although there are many experimental models and analytical techniques, ADME research in TCM is not an easy task. First, it is necessary to establish a model that fits the human body system to study the multiple constituent ADME. Second, studies on PK and PD should be combined and PK markers that can represent the overall efficacy of TCM should be chosen.

Analysis and validation of the three action modes of TCM

PMTN

Target databases include DrugBank, STITCH, Therapeutic Target Database (TTD), and BindingDB, and the comprehensive TCM databases include ETCM, CHEM-TCM, BARMAN-TCM and TCM-Mesh. In the ETCM database, target information was collected from the predictions of MedChem Studio (version 3.0), and these targets were obtained through DrugBank. At the biological level, omics techniques (genomics, transcriptomics, proteomics, metabolomics, etc.) are often used to research drug targets. MetaboAnalyst, MaxQuant, UniProt The Cancer Genome Atlas (TCGA), and the International Cancer Genome Consortium (ICGC) are websites and database resources related to omics research. Additionally, Databases For Annotation, Visualization And Integrated Discovery (DAVID), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome are often used for pathway and enrichment analyses.

PMCGM

Virtual Metabolic Human (VMH) is a metabolic database of human and intestinal flora that can be combined with human and intestinal flora to reconstruct a metabolic model for large-scale, multi-omics research on diet, metabolism and disease. The data repository for gut microbiota (GMrepo) database is a database of human intestinal metagenomics that achieves data sorting by manually managing the corresponding human host macro data (age, gender, country, etc.). This database covers 58,903 human gut runs/samples from 253 projects, which are related to 92 human phenotypes, and users can easily access the collected data. The gut MEGA database collects published human quantitative gut microbial metagenomic data. In addition, the Broad Institute-OpenBiome Microbiome Library (BIO-ML) is a comprehensive database consisting of 7758 isolates of intestinal bacteria, 3632 genome sequences and longitudinal multiomics data.

Currently, commonly used molecular docking software includes AutoDock, eHiTS, FlexX, Glide, GOLD, LigandFit, and Surflex, among which AutoDock is the most cited free software. The structures of most small molecules can be obtained from the PubChem and Zinc databases. Receptor protein structures can be obtained through the Protein Data Bank (PDB) database. In addition, reverse virtual screening based on reverse molecular docking can also be used to predict drug targets. Importantly, we should note that computer simulation docking data cannot completely replace experimental data. A high docking score and high affinity do not necessarily mean that a ligand has a better effect, as this needs to be verified by biological experiments.

CCIs-PMPK

The current drug interaction software packages used at home and abroad include Micromedex® Drug-Reax, Micromedex® Healthcare Series, Drug Interactions Facts®, Lexi-Interact®, Pharmavista®, EpocratesRx®, MediQ®, and Drug interaction checker®. However, inconsistent standards exist between these databases, and the analysis results obtained by using different database software programs are different. Therefore, it is essential to be cautious when using CCI results as clinical decisions. In addition, it is necessary to establish a multiconstituent CCI model for the complex system of TCMs.

Target fishing has the characteristics of high throughput, high speed and high specificity, and it has a wide range of target screening for TCMs. Small molecules are immobilized on a chip, interact with the protein lysate, and the eluent is detected by MS, thereby fishing out the targets of small molecules. However, some proteins with weak binding forces are unavoidably and often missed in vitro, resulting in false negatives. Moreover, in vitro research is different from the real environment of an organism, which may lead to error in reporting. Chemical proteomics uses synthetic chemical methods to generate probes, which can be used to analyze the protein targets of small molecules in cells. Importantly, it can also be administered to cells in situ, which has significant biological advantages. However, this technology is not perfect, and nonspecific binding may occur. How to deal with promiscuity and background binders is a problem that needs to be solved. In the experiment, multiple methods are needed to achieve complementary effects. In addition, the results to be verified in organisms in order to determine the targets of the small molecules.

The most commonly used methods for detecting intestinal flora include 16S rDNA sequencing, metagenomic blood sequencing and metabolomics. 16S rDNA sequencing can accurately quantify all bacterial species in the intestinal microbes, metagenomic sequencing can discover the enrichment of important coding genes or pathways, and metabolomics is related to the phenotype, which can directly reflect the changes in intestinal microbes. Therefore, the combined use of 16S rDNA sequencing, metagenomic sequencing, and metabonomics can overcome the limitations of single omics research to a certain extent. In animal intestinal flora research experiments, the consumption of intestinal flora can be realized by using spectrum antibiotics, and a “sterile” animal model can be established. Additionally, the most commonly used methods of administration, such as gavage and enema, can be used for bacterial transplantation. Moreover, molecular probe technology and gene chip technology can also be used to carry out research on intestinal flora.

Some extracorporeal liver systems can be used to evaluate the interaction between constituents, such as reconstituted microsome systems, recombinant human CYP enzymes, and hepatocytes prepared from human liver tissues. In addition, the in vitro system of membrane vesicles, bidirectional transport assays with cell-based systems, and uptake assays with cell-based systems can be used to study transporter-mediated CCI. Commonly used cell models for studying drug transporters (P-gp, BCRP, OAT1/3, OCT2, etc.) are Caco-2 cells, Madin–Darby canine kidney cells, Lilly Laboratory cancer porcine kidney 1 cells, Chinese hamster ovary cells, human embryonic kidney 293 cells, etc..

Multidimensional association network of qualitative PK–PD correlations

Disease databases, including Human Phenotype Ontology (HPO), Online Mendelian Inheritance In Man (OMIM), MalaCards and DisGeNET, protein interaction databases, including String, Human Annotated Protein–Protein Interaction (HAPPI), High-Quality INTeractomes (HINT), Online Predicted Human Interaction Database (OPHID), Molecular Interaction Database (MINT), Human Protein Reference Database (HPRD), Database Of Interacting Proteins (DIP), PDZBase and IntAct, and other databases, including GO, DAVID, KEGG, and Reactome, are commonly used databases for the construction of biological networks. Some TCM comprehensive resource databases, such as ETCM, TCMSP, and TCMID, can also be used to construct TCM networks, especially the ETCM database. The protein–protein interaction data in ETCM are collected from five molecular interaction databases, including Reactome, HPRD, MINT, IntAct and DIP. The ETCM database provides a platform for establishing a biological network of formulas, herbs, ingredients, targets, diseases, and pathways, and users can easily build biological networks according to their own needs.

Network construction

In the construction of the network, the involved algorithms include the network module search algorithm and the node sorting algorithm in the biological network. Among them, the restricted neighborhood search clustering (RNSC) algorithm is a clustering algorithm based on the cost function. Random walks with restart (RWR) is an algorithm for the random selection of nodes. In addition, there are other algorithms, such as Ingenuity Pathway Analysis (IPA) and Network Propagation (NP). Cytoscape, Matascape, NAVIGaTOR, and Pajek are often used for the visual analysis of complex networks. Users can integrate biological networks with various molecular data, such as gene expression in a visual environment, and link these networks to functional annotation databases.

Network analysis

Common network analysis technologies include network node centrality analysis, module analysis, global topological attribute analysis, comparison and similarity analysis, dynamic analysis, etc.,. Node connectivity, node compactness, and node mediation are all topological characteristic indicators used to evaluate the centrality of the node network. The nodes with three characteristic values that are all greater than the median of the corresponding characteristic values of all nodes in the network are selected as the core nodes in the network. These indicators can be obtained through TCMIP v2.0 (http://www.tcmip.cn/TCMIP/index.php/Home/Login/login.html), NAViGaTOR, TCMSP and TCMIP. Additionally, we can build a custom multidimensional network based on the following principles in TCMIP V2.0: retain the TCM components with good drug properties; retain the TCM components with a large number of corresponding core network targets; retain the core network targets to be verified and the path nodes that are significantly involved; and retain the core targets closely related to the key pathogenic links of the target disease. All of the above screening is based on the topological importance of the core network targets. Through these methods, the complex network is simplified, and the readability of the network is improved.

However, we cannot ignore some existing problems. The accuracy and comprehensiveness of most biological networks and network analysis algorithms have some limitations. In addition, most existing biological networks are undirected networks, from which it is impossible to know the mode of action between neighboring molecules. Information on the type and direction of the interaction needs to be further improved.

Quantitative PK–PD correlations

IV–PK/PD–DCMs in tissues and cells

There have been in vitro PK–PD models: a dialysis/diffusion closed model to study the antibacterial activity and pharmacodynamics of the cerebrospinal fluid drug concentration and that in the central nervous system for neurological drug evaluation and the intestinal absorption–vascular activity combination model established by our group to study the vasodilator activity. Nevertheless, conventional models are limited due to their lack of physiological relevance, and in vitro models based on microfluidics have become a potential solution. For example, the toxicity of anticancer drugs can be tested by the connection between liver, tumor and bone marrow cell lines using microfluidic technology. Based on the concept of TCMIP, we propose an in vitro PK–PD study with microfluidics as the key technology (Fig. 7). In the entire microfluidic device, the liquid flow can be controlled by a circulating pump, which is set between the chips to control the liquid flow and direction between the chips. With this device, a system of “drugs–gut microbiota metabolism–intestinal absorption–liver metabolism–drug transport–target tissue” is established to simulate the process of the drug in the human body. In the whole system, the key part is that the type and concentration of drug components obtained by the microfluidic device is consistent with the corresponding animal samples. Moreover, the efficacy index and effect intensity of the microfluidic samples are also consistent with those of the in vivo samples. Finally, based on an in vitro microfluidic system, the effects of the drug on the target are studied. However, we should know that although we already have such an idea and the corresponding technology, there is still a long way to go to truly simulate the drug process in the human body, which requires joint efforts of multiple disciplines and continuous technological innovation.

Figure 7

Construction of the IV–PK/PD–DCM via microfluidic-based chip technology in combination with the intestinal flora, intestinal cells, liver drug enzyme system and tissue distribution that simulate the drug ADME/PK and pharmacodynamic process. IV–PK/PD–DCM is utilized to obtain TCM administration samples in vitro with a consistent composition and concentration of target tissues or cells in vivo and to evaluate the bioactivities of a TCM in vitro more scientifically. IV–PK/PD–DCM: in vitro PK–PD dynamic complex models.

Common PK–PD models include the linear model, log-linear model, maximum effect model (Emax), sigmoid Emax model, and β-function model. Moreover, the connection mode of PK–PD can be divided into four attributes. NONMEM is a powerful tool for nonlinear mixed effects modeling and pharmacokinetic and pharmacodynamic data simulation and has been used in the PK–PD study of intranasal and intravenous dexmedetomidine. In addition, the ANN modeling method can be used to establish a PK–PD prediction model,. These models are more suitable for the study of a single active constituent. In the TCM system, there is a complex relationship between multiple constituents and multiple effects, where the constituents also interact with each other. Therefore, it is necessary to develop models in accordance with the characteristics of TCM.

Validation of the key active constituents

Gene knockout is a mature technology that can be used to study key targets. Under the influence of gene knockout, constituent knockout technology has been gradually introduced into the study of the effective material basis of TCM. Its research methods include chromatography knockout and antibody knockout,. The chromatographic knockout method mainly relies on the collection and separation technology of the preparative liquid chromatography system. It has the advantages of simplicity, rapidity, and a wide application range, but it has higher requirements for the separation effects of preparation instruments. The principle of the antibody knockout method is similar to that of gene knockout, which can specifically knock out an active ingredient in the extract. This method is suitable for polar and macromolecular compounds and has high selectivity, but its operation is complicated, and its scope of application is narrow. In the future, more experiments are needed to improve constituent knockout technology.

Applications of TCMIP

Quality marker identification for improving TCM quality control

TCM quality control guarantees safety and efficacy for clinical use. However, there are many challenges in quality control, including the use of the same marker ingredients to evaluate the quality of different herbs, the lack of an effective quality traceability system from medicinal material to product, and the lack of correlation between multiple marker ingredients and the safety and efficacy of a TCM. In recent years, we proposed the Q-marker as a novel concept for quality evaluation and standard elaboration of TCM,. A constituent can be defined as a Q-marker when it possesses specificity, abundance, inherent chemical compounds in the products, appropriate PK properties, and ameliorative effects. Notably, Q-markers may contribute greatly to the therapeutic effects or safety of TCM prescriptions; they are often closely related to the pharmacological mechanisms used to establish the qualitative and quantitative correlation between the quality standard of a TCM and its clinical efficacy. Thus, Q-markers may become a new direction in the quality development of TCMs and a breakthrough in the modernization and internationalization of TCMs.

TCMIP is a systematic approach to screen out Q-markers of TCM based on the following steps. (1) Comprehensive comparison of the chemical profiles of different herbs and TCM formulas to identify the characteristic components. (2) Establishment of the quality traceability system to screen out the stabilizing and inherent ingredients from crude herbal drugs or their products. (3) Characterization of the ADME/PK profiles of herbal and TCM formulas to discover the potential KACs in vivo. (4) Establishment and validation of the qualitative and quantitative PK-PD correlations between KACs and the corresponding bioactivities. (5) Optimization of the minimum combination of the KACs of TCM and establishment of the quality control methods and quality standards of TCMs. According to the above research strategy, several TCMIP-based investigations have been carried out to identify the potential Q-markers of various TCM formulas, such as the Yuanhu Zhitong prescription (YZP) and Xin-Su-Ning capsules.

YZP is a classic Chinese patent medicine (CPM) that contains Angelicae Dahuricae Radix and Corydalis Rhizome for the treatment of gastralgia, costalgia, headache, and dysmenorrhea. The quality of YZP produced by different pharmaceutical companies is quite different due to the different origins, harvest times, pretreatments, manufacturing processes and dosage forms of the medicinal materials used. Based on our research over the last several years, we utilized the TCMIP approach to systematically identify Q-markers and establish a quality standard for YZP. First, the chemical constituents of YZP were characterized by chemical fingerprinting and multicomponent quantitative determination,. Second, the intestinal absorption of the YZP constituents was analyzed in vitro, as well as serum pharmacochemistry and pharmacokinetic evaluation, which revealed the ADME/PK profile and identified potential KACs in YZP in vivo147, 148, 149. Third, a computational constituent–target network was constructed, and its molecular basis mainly involved opioid receptors, dopamine receptors, cation channels, GABA-A receptors, benzodiazepine receptors, etc., which indicated that YZP had analgesic, antianxiety, antidepression, vasodilatation and other activities. Based on this, seven constituents were qualitatively identified as Q-markers of YZT based on their high abundance, specific presence in the individual herbs and final product, appropriate drug-like properties, and critical role in the bioactivity of the mixture of YZT constituents. Finally, data mining methods, such as grey relational analysis (GRA) and least squares support vector machine (LS-SVM), were utilized to precisely characterize the quantitative correlation between the identified Q-markers of the YZTs and their efficacy. Three Q-markers were chosen as a minimum combination to control the quality of the YZTs, indicating that the determined Q-markers were suitable for quality control use of this TCM formulation and may represent a widely applicable, low-cost, rapid, simple TCM quality control method.

Buchang Naoxintong capsule (BNC) is a well-known TCM prescription for the treatment of cardiovascular and cerebrovascular diseases. We first developed qualitative and quantitative analytical methods for the rapid high-throughput screening of the preliminary chemical profile of BNC. Then, based on text mining and predictions of the intestinal microbial metabolism and oral bioavailability of these components, the interaction network between the BNC components and their therapeutic molecular targets was constructed. According to their high abundance in the products, appropriate drug properties and importance to the biological activity, amygdalin and paeoniflorin were identified as KACs of BNC. Finally, based on the middle cerebral artery occlusion (MCAO) model, the neuroprotective effects of amygdalin and paeoniflorin were evaluated. The results showed that amygdalin and paeoniflorin could significantly reduce the cerebral infarction volume and improve the neurological function score. Interestingly, we also demonstrated that amygdalin exerted a partial anticerebral ischemia effect by interacting with the glucocorticoid receptor NR3C1 and serpin family C member 1 (SERPINC1). Therefore, KACs should be an important basis of Q-markers, which can be qualitatively and quantitatively related to the bioactivities of TCMs.

Clarification of the molecular mechanisms of TCM therapeutics and the promotion of accurate clinical applications

For a long time, TCM theories have provided guidance for the clinical applications of TCM based on syndrome differentiation and treatment. However, with the development of molecular biology and systems biology, under the promotion of molecular bases, accurate clinical application will become the future development trend of TCM. CPMs, as drugs approved by the National Medical Products Administration (NMPA), have a large market share and high clinical value. Increasing evidence has shown that some CPMs have been investigated in randomized double-blind experiments and placebo-controlled studies, which shows that many CPMs have clear clinical efficacy and advantages for the treatment of complex diseases,. In previous studies, we employed a TCMIP approach to reveal the molecular mechanism of CPM therapeutics and to screen out KACs and CMTs, such as Yuanhu Zhitong tablets (YZTs), Guanxinjing capsule, and Quanduzhong capsules (QDZJNs). For example, 18 AP and 13 AM of YZTs in vivo were virtually calculated to potentially bind to the opioid receptor, dopamine receptor, cation channel, GABA-A receptor, and benzodiazepine receptor, which suggested that YZTs might have analgesic and antidepressive action. Then, mouse hot plate and writhing tests proved that YZTs exhibited a significant analgesic effect with a dose–effect relationship. Additionally, the forced swimming test and the tail suspension test demonstrated that YZTs had a significant decrease in diving compared to controls and possessed obvious antidepressive activity. Furthermore, the total alkaloids of YZTs were confirmed as the main active constituents of the antidepressive effect in a chronic unpredictable mild stress (CUMS) rat model using 1H NMR-based metabonomics. This evidence indicated that the optimal indications of YZTs would be chronic pain with depression. As an NMPA-approved drug, QDZJN is widely used to control blood pressure and protect renal function. In our study, the potential targets of QDZJN were specifically expressed in the kidneys and were involved in the inflammatory response, blood pressure regulation, LPS response and hypoxia. Based on network robustness assessment, it was suggested that QDZJN may have a greater impact on the glomerular network after clearing its potential targets, which provided more evidence for the precise clinical applications of QDZJN against hypertensive nephropathy.

Limitations and future prospects

TCMIP has emerged as a novel approach to elucidate qualitative and quantitative PK–PD correlations by integrating the knowledge and key techniques of multiple disciplines. TCMIP puts forward the three forms of active substances and the corresponding action modes of TCMs, the dynamic development of PK and PD profiles of TCMs with various treatment time points and action spaces, as well as the relationship between macroscopic and microscopic, in vivo and in vitro, PK and PD, etc. Although great progress has been made in its methodology and applications, there are several problems and limitations of current TCMIP research. First, the development of TCMIP needs high-quality data and efficient algorithms. Although there have been several TCM-related databases and research platforms, comprehensive and structured TCM-related data ranging from ancient TCM literature to current clinical medical records, from clinical trials to basic research data involving drugs, and from TCM syndromes and symptoms to complex molecular bases are urgently needed. Moreover, efficient algorithms that are in accordance with the complex characteristics of TCM are lacking. Second, there is a lack of mature in vitro and in vivo pharmacological technologies to verify the qualitative and quantitative PK–PD correlations between multiple KACs of TCMs, and CMTs and CPEs of TCM therapeutics. Third, it is difficult to carry out research on PMCGM and CCIs-PMPK due to the lack of proper calculation and experimental techniques.

To overcome the above bottlenecks in the TCMIP research field, it is necessary to establish a high quality, comprehensive, structured and ecological big data TCM platform by systematically sorting ancient books and documents of TCM theories, dynamically collecting current clinical medical records of famous TCM doctors, and arranging overall information on the basic characteristics and related molecular properties for TCM syndromes, TCM formulas, Chinese herbal medicines, etc. In addition, we should develop a series of innovative algorithms and experimental methods to comprehensively investigate the three action modes of TCMs (PMTN, PMCGM, CCIs-PMPK) in TCMIP, such as TCM-related ADME/PK virtual prediction models, high-throughput prediction and determination of the interactions between the multiple constituents of TCM and gut microbiota, and high-throughput prediction and determination of ADME/PK enzyme or transporter CCIs. In particular, IV–PK/PD–DCMs will be developed as an intelligent platform for in vitro pharmacokinetic and pharmacodynamic evaluations and molecular mechanism research. Moreover, systemically uncovering the active constituents and therapeutic basis of TCM formulas at the molecular level should be carried out based on the guidance of TCM theories in order to enhance the level of quality control, personalized treatment and precise clinical applications of TCM. Notably, TCMIP must strengthen the correlation and transformation with clinical trials of TCMs. Taken together, TCMIP is still in its infancy, and advancements in this field will undoubtedly elicit a conceptual change in combination drug discovery and make an important contribution to the modernization and globalization of TCM.

Acknowledgments

This work was supported by grants from the (Grant Nos. 81830111 and 81774201), (2017YFC1702104 and 2017YFC1702303), the (No. 2016SF-378, China), the (ZXKT17058, China), the National Science and Technology Major Project of China (2019ZX09201005-001-003). The funding agencies had no role in the study design, the collection, analysis, or interpretation of data, the writing of the report, or the decision to submit the article for publication.

Author contributions

Haiyu Xu, Yanqiong Zhang, Ping Wang, Junhong Zhang and Hong Chen wrote the article. Yanqiong Zhang and Junhong Zhang are responsible for sorting out the table. Hong Chen, Luoqi Zhang, Xia Du, Chunhui Zhao and Dan Wu participated in the drawing of the illustrations. Changxiao Liu, Haiyu Xu, Feng Liu and Hongjun Yang conceived the article.

Conflicts of interest

The authors declare no conflicts of interest.