Evidence-Based Complementary and Alternative Medicine
Volume 2018 (2018), Article ID 1854972, 11 pages
https://doi.org/10.1155/2018/1854972

Research Article

A Component Formula of Chinese Medicine for Hypercholesterolemia Based on Virtual Screening and Biology Network

Xiaoqian Huo, Fang Lu, Liansheng Qiao, Gongyu Li, and Yanling Zhang

School of Chinese Pharmacy, Beijing University of Chinese Medicine, Beijing 100102, China

Correspondence should be addressed to Yanling Zhang; zhangyanling@bucm.edu.cn

Received 20 November 2017; Revised 18 May 2018; Accepted 4 June 2018; Published 28 June 2018

Academic Editor: Juntra Karbwang

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

Abstract.

Hypercholesterolemia is a risk factor to atherosclerosis and coronary heart disease II. The abnormal rise of cholesterol in plasma is the main symptom. Cholesterol synthesis pathway is an important pathway of the origin of cholesterol, which is an essential pathway for the therapy of hypercholesterolemia. The 3-hydroxy-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase), squalene synthase (SQS), and sterol regulatory element binding protein-2 (SREBP-2) are closely connected with the synthesis of cholesterol. The inhibition of these targets can reduce the cholesterol in plasma. This study aimed to build a component formula including three Traditional Chinese Medicines (TCM) components with the inhibition activity of these targets by using virtual screening and biological network. Structure-based pharmacophore models of HMG-CoA reductase and SQS and ligand-based pharmacophore model of SREBP-2 were constructed to screen the Traditional Chinese Medicine Database (TCMD). Molecular docking was used for further screening of components of HMG-CoA reductase and SQS. Then, metabolic network was constructed to elucidate the comprehensive interaction of three targets for lipid metabolism. Finally, three potential active compounds were obtained, which are poncimarin, hexahydrocurcumin, and forsythoside C. The source plants of the compounds were also taken into account, which should have known action of lowering hyperlipidemia. The lipid-lowering effect of hexahydrocurcumin was verified by experiment in vitro. The components that originated from TCMs with lipid-lowering efficacy made up a formula with a synergistic effect through the computer aid drug design methods. The research provides a fast and efficient method to build TCM component formula and it may inspire the study of the explanation of TCM formula mechanism.

1. Introduction

Hypercholesterolemia, a kind of hyperlipidemia, is the abnormal rise of serum cholesterol, which is an induced factor to atherosclerosis and coronary heart disease [1, 2]. The cholesterol synthesis pathway generates most of cholesterol in the human body, which provides a significant approach for antihypercholesterolemia. So the targets in the cholesterol synthesis pathway are always hot spots for new lipid-lowering drug development. There are multiple targets in the synthesis pathway. Besides, some targets have been important targets for the drug design because of better pesticide effect. 3-Hydroxyl-3-methylglutaryl-coenzyme A reductase (HMG-CoA reductase) is one of the most commonly used targets in the therapy of the hyperlipidemia, which is the limit enzyme of cholesterol synthesis. As the inhibitor of HMG-CoA reductase, statins provided many benefits for patients including effective reduction of TC, TG, low-density lipoprotein cholesterol (LDL-C) levels, and neovascularization and lowering cerebral injury with curative effect and security [3–6]. Squalene synthase (SQS) was another key enzyme in the cholesterol synthesis, which can catalyze the biosynthesis of the key cholesterol precursor, squalene [7]. The design of SQS inhibitor has never stopped. Some SQS inhibitors have favorable therapeutic effect in reducing cholesterol [8, 9]. Because of the protein-protein interaction between HMG-CoA reductase and SQS, the inhibitor of the SQS can be a complementary drug of the HMG-CoA reductase. The research provided evidence of SQS inhibitors as a monotherapy and adjunctive therapy for the patients who cannot achieve a satisfied cholesterol target goal [10]. The coadministration with statins may decrease statin-induced myotoxicity [11, 12]. New SQS inhibitors have been in continuous research [13]. Sterol regulatory element binding protein-2 (SREBP-2), a member of the SREBP family (SREBPs), is an important cholesterol regulator in cells that can directly activate the expression of more than 30 genes including HMG-CoA reductase, SQS, and other important enzymes and regulatory factors in the lipid metabolism pathway [14–16]. It has been discovered that the inhibitors of SREBP-2 can decrease cholesterol level and have well application prospect in the treatment of hypercholesterolemia. The targets were chosen including two important enzymes in the synthesis pathway and an important regulator. The protein-protein interaction between the targets can get more satisfied antihyperlipidemia curative effect [17, 18]. The Traditional Chinese Medicines (TCM) were considered as a source of inspiration of new drugs [19]. The TCM were a complex system that includes multiple components. Components take synergistic effect by acting on different targets for enhancing curative effect and reducing the side effects [20]. So, microconstituents in TCM can play the same or a better role of the megadose of single component. The formula was a traditional drug combination method that combines Chinese medicines with different effect under the guidance of TCM [21]. Component formula of the Chinese medicine was a combination of the TCM components or extracts. Component formula covers the shortage of single-ingredient drug and provides the benefit such as the clear action targets and mechanism of herbal pharmacology. The component formula of the Chinese medicine is a theoretical innovation in the Chinese materia medica theory.

This research hopes to build a component formula based on the computer aid drug design (CADD). Pharmacophore and molecular docking models were used for the virtual screening of the TCM component with potential inhibition activity of HMG-CoA reductase, SQS, and SREBP-2. A network of these three targets was built to show the synergistic effect between three targets, which provided the mechanism of action for the component formula. Besides, the TCM components were chosen under the consideration of TCM theory.

Pharmacophore and molecular docking models were built for the virtual screening. The components with potential activity of chosen targets were chosen. To expound the synergistic effect between the components, a network based on three targets was built. Pharmacophore and molecular docking were two important methods in the virtual screening. Pharmacophore model, which was an important method of molecular modeling for virtual screening, has been widely used to search the potential active compounds [22, 23]. Ligand-based pharmacophore model is generated by collecting important atoms and atoms groups of active ligands, which shows the necessary features of active ligands. Galahad is one of the most frequently used methods of ligand-based pharmacophore [24]. Receptor-based pharmacophore models are based on the known three-dimensional (3D) structure of targets, which can reflect the spatial relationships between protein structure and ligand in active site. Structure-based pharmacophore (SBP) method is a kind of receptor-based pharmacophore model. In the present study, pharmacophore models of HMG-CoA reductase inhibitors and SQS inhibitors were built by SBP, respectively, while pharmacophore models of SREBP-2 inhibitors were built by Galahad because of lacking protein structures. The best pharmacophore models of each target were selected and used to search the Traditional Chinese Medicine Database (TCMD, vision 2009) [25, 26]. Then, virtual screening results of HMG-CoA reductase and SQS were further screened by molecular docking. Molecular docking is another important method in CADD, which can demonstrate the interaction between ligand and the receptor. The information of the interaction was valuable in the drug design. Recently, biological network provides a new method for complex system research. Metabolic network, an important type of biological network, can reflect the interaction among compounds, enzymes, and regulatory factors involved in the metabolic process for the analysis. The application of metabolic network provides a bridge to study on the mode of action of TCM, which has been widely used for thousands of years by taking advantage of “multiple ingredients and multiple targets” [27–29]. In this paper, one metabolic network model was constructed for analyzing the mechanism of interaction among HMG-CoA reductase, SQS, and SREBP-2 to explain how effectiveness was enhanced by drug combination through the interaction among the targets.

This research hopes to build a component formula, which has obvious advantages such as exact dose and relatively clear mechanisms of action. This work combined the CADD methods and TCM theory for a component formula. The CADD can effectively screen component in a short time and provided the information of targets and interaction of the component for further study. Network elucidates the synergistic effect between the targets, or it may aid the discovery of potential targets and the mechanism analysis. Then, the source plant can diminish the range of suitable components that are more likely to have activity.

In summary, the virtual screening is a powerful method for searching active compounds from TCM and the mechanism of complex effect of multicompounds can be clearly illuminated based on biological network. The component formula with clear material foundation and certified mechanism can be formed through the combination of potential active compounds with different targets’ inhibition activity. The article aims to provide a clue to new component formula design and the study of mechanism of TCM component formula based on virtual screen and biological network.

2. Materials and Methods

2.1. Construction of Pharmacophore

The pharmacophore models of HMG-CoA reductase and SQS were constructed by “Receptor-Ligand Pharmacophore Generation Protocol” based on the crystal structures of proteins that were obtained from the Protein Data Bank (http://www.rcsb.org/pdb/home/home.do) in the Discovery Studio (DS, vision 4.0) [30]. By entering the “HMG-CoA reductase” and “Squalene synthase” as the search terms in the PDB, the structures with high resolution and known active ligands were chosen to be applied in the research. The information of structures was shown in Table 1.

Table 1: The information of the structures of HMG-CoA reductase and SQS.

PDB ID
Resolution (Å)
Ligand ID
(nM)
Target

3ASX
2.00
D99
20.0
SQS

1HWL
2.10
FBI
0.9
HMG-CoA reductase

The conformations for all compounds were generated by the BEST conformation generation method. Active site was the site of the ligand combined with the structure and the location of the pharmacophore. Active site was searched automatically. All crucial features, including Hydrogen Bond Acceptor (A), Hydrogen Bond Donor (D), and Hydrophobic (H), would be mapped in the active sites. Furthermore, A, D, and H features were clustered automatically and the most crucial features were chosen to build pharmacophore model by employing Cluster Current Features protocol. Besides, Ev (Exclusion volume) spheres were added to the inaccessible areas in active site.

Because of lacking the whole structure of SREBP-2, ligand-based pharmacophore method, Galahad, was applied to construct SREBP-2 inhibitor pharmacophore models [31]. Seven compounds with the known activity of SREBP-2 were chosen from the Binding database (http://www.bindingdb.org) as the training set to construct the pharmacophore models. The structure and the IC₅₀values of the compounds in the training set were shown in Figure 1. The force field was set to Tripos and the charge was set to Gasteiger-Huckel atomic partial charges. The energy was minimized by the Powell method. The optimal set of the conformations was fixed as rigid bodies in Cartesian space. Then twenty models were generated by the genetic algorithm with different features, conformations, and overlay of the molecules [32, 33].

Figure 1: The structures and IC₅₀ values of the compounds in the training set for SREBP-2 pharmacophore modeling.

2.2. Pharmacophore Validation

Test set, including active compounds and inactive compounds, is used to assess the discrimination capacity of pharmacophore models. The numbers of the compounds in the test sets were shown in Table 2. Active compounds were obtained from the Binding database, while the inactive compounds were chosen from the MDDR (MDL Drug Data Report: Version 2007.2).

Table 2: The validation results of each pharmacophore model.

Name
Feature
Specificity
Energy
An
Dn
Ha
Ht
HRA
IEI
CAI

M-HMG-CoA reductase-1
AH₃Ev₁₄
—
—
92
375
52
72
0.56
2.94
1.66

M-HMG-CoA reductase-2
AH₃Ev₁₄
—
—
92
375
60
203
0.25
1.20
0.79

M-SQS-1
A₃DH₂Ev₁₈
—
—
96
329
41
54
0.43
2.60
1.11

M-SQS-2
A₃DH₂Ev₁₈
—
—
96
329
54
140
0.56
0.56
0.74

M-SREBP-1
A₂DH₃
3.7
6.26
20
60
20
40
100
2.00
2

M-SREBP-2
A₂DH₃
3.9
3.79
20
60
20
24
100
3.33
3.33

M-SREBP-3
A₂DH₃
4.6
5.17
20
60
20
22
100
3.64
3.64

aAn: active compounds number. bDn: all compounds in test number. cHa: hit active compounds number. dHt: hit compounds number. eHRA: the hit ratio of active compounds. fIEI: identification of effective index. gCAI: comprehensive appraisal index.

Hit ratio of active compounds (HRA), identification of effective index (IEI), and comprehensive appraisal index (CAI) values were employed to evaluate the SBP models. The high HRA shows a great ability of pharmacophore model to identify active compounds in the test set. The high IEI means an outstanding ability to identify active compounds from inactive compounds. CAI is calculated to evaluate the models comprehensively [34, 35].

The SBP model with highest CAI and IEI would be chosen as the best pharmacophore model to screen the TCMD. The best pharmacophore models built by Galahad could be evaluated by both test set result and the internal evaluation, which was selected by the following three criteria: (1) the model should be in high IEI and CAI values; (2) the model needs to have lower energy; (3) the model needs to have high specificity. The best pharmacophore model would be used to screen the TCMD.

2.3. TCMD Screening

The best pharmacophore models of three targets were utilized as queries to search the TCMD, which include 233,033 compounds from the Traditional Chinese Medicine. The “flexible database search” was carried out in virtual screening process. FIT value of each compound was calculated to assess the mapping quality of the compounds. The higher FIT value suggested a better mapping between the compounds and the pharmacophore model. Then the hit compounds would satisfy Lipinski’s rules. Finally, the compounds with drug-like character and potential inhibited activity were listed for further docking research.

2.4. Molecule Docking

The crystal structure 3ASX and 1HWL were used for molecular docking of HMG-CoA reductase and SQS. The binding site was the same to pharmacophore model. LibDock and CDOCKER were preliminarily selected in this work. To choose a more appropriate algorithm, the initial compounds FBI and D99 were extracted from the binding site and redocked into the structure for calculating the RMSD value. The RMSD was used to evaluate a rationality of the algorithms and docking parameter settings. In general, RMSD of less than 2.00Å indicated that the algorithm could reproduce the observed binding mode for the initial compound [36, 37]. A smaller RMSD value means that the docking method better fits the system. The compounds of each target from the pharmacophore results were docking into the protein structure by the same methods of the initial compounds. Scoring functions were used to evaluate the results. The compounds with higher scores may have better interaction.

2.5. Metabolic Network Construction

The pathways of the three targets were obtained from the Reactome database (http://www.reactome.org) [38]. Then the pathways were imported into the Cytoscape (vision 2.8.3) [39, 40]. The construction of the network was completed in the Cytoscape by merging the pathways of three targets. Then the isolated nodes and self-loop edges would be removed and the repeated sides and points were deleted.

The parameters were calculated to certify the stability and fault tolerance of the network. Network connectivity is the character for evaluating the connected condition of network. The network has high antidestroying ability when the connectivity value reaches 100%. Scale-free property of network is evaluated by degree distribution. Network has a scale-free property when the degree distribution of the network obeyed power-law distribution, which has high fault tolerance. Barabasi-Albert model (BA model) put forward by Barabasi and Albert [41] pointed the standard of the scale-free network.

In generation, b<3 is the criterion to evaluate whether the network is in scale-less property. The networks can be applied intuitively to show the interaction among the three targets and the mechanism of formula synergistic effect.

2.6. The Design of TCM Component Formula

The component was chosen by the following two criteria: (1) component of HMG-CoA reductase and SQS should have a high rank of the sum of the FIT value ranks and docking score ranks of the compounds. The sum can evaluate the activity of the compounds. The compound that has a higher rank may have a higher activity; (2) the source plants of HMG-CoA reductase, SQS, and SREBP-2 should have a known lipid-lowering action or form a formula related to lipid-lowering. The source plants showed the traditional use in the TCM, which provide a clue of the activity of component. Network of the targets was built for the analysis of synergistic effect between the components. The related targets were explored in the network, which inspired other effects or the mechanism of action of the formula.

2.7. The Activity of TCM Compounds

The activity of TCM compounds was verified on the cell. HepG2 cell lines were used for the MTT assay and activity assays. Cells were cultured at 1×104 cells/well in 96-well culture plate. The cells were adhered to and treated with different concentrations of the targeted compounds after 24 h and cultured at 37°C in the atmosphere of 5% CO₂. Then, MTT in 0.5 mg/ml was added to each well after the supernatants were discarded and incubated at 37°C in 5% CO₂ for an additional 4 h. After that, 150 μL of formazan in dimethyl sulfoxide (DMSO) was added after the MTT was removed. And then the plates were shaken for 5 min at low speed. Cell proliferation was evaluated by measuring the absorbance at 570 nm using ELISA Plate Reader. The IC₅₀ values were calculated by SPSS Statistics 17.0.

3. Results and Discussion

3.1. Pharmacophore Generation

Two pharmacophore models of HMG-CoA reductase were built in this study. The models contained hydrophobic features and hydrogen bond acceptor feature. SQS pharmacophore models were composed of the hydrophobic features, hydrogen bond acceptor features, and hydrogen bond donor feature. Ev features were added to each pharmacophore where in the protein structure the compounds cannot approach. The number of the Ev features in HMG-CoA reductase was fourteen and the number in the SQS was eighteen. Twenty models of SREBP-2 were generated [42, 43]. The results of three of them were shown in Table 2. The pharmacophore models of SREBP-2 contained hydrophobic features, hydrogen bond acceptor features, and hydrogen bond donor feature.

3.2. Pharmacophore Validation

The pharmacophore validation results were shown in Table 2. M-HMG-CoA reductase-1 and M-SQS-1 with the highest IEI and CAI values were chosen as the best models to screen the TCMD. M-SREBP-2 has the least energy and high specificity, IEI, and CAI value. So it was chosen to screen the TCMD. The best pharmacophore was shown in Figure 2.

Figure 2: The best pharmacophores of each target.

The M-HMG-CoA reductase-1 included one hydrogen bond acceptor feature, three hydrophobic features, and fourteen exclusion volume spheres. To check the rationality of the HMG-CoA reductase pharmacophore model, the initial compound FBI was mapped with the best pharmacophore model. The mapping result was shown in Figure 3. The initial compound FBI was mapped well with all the features of the pharmacophore model of HMG-CoA reductase.

Figure 3: The mapping results of each initial compound.

Three hydrogen bond acceptor features, one hydrogen bond donor feature, two hydrophobic features, and eighteen exclusion volume spheres composed the M-SQS-1. The initial compound D99 was mapped with the M-SQS-1. The methyl of the D99 was mapped with the hydrophobic feature. The hydroxyls in the aromatic ring were mapped with hydrogen accept features. The mapping result was shown in Figure 3. Three hydrophobic features, two hydrogen bond acceptor features, and one hydrogen bond donor feature formed the M-SREBP-2.

Hexahydrocurcumin was the potential compound with HMG-CoA reductase inhibition activity to build the component formula. Hexahydrocurcumin was mapped well with all the features of the HMG-CoA reductase pharmacophore model. The aromatic ring and the methyl in the aromatic ring were mapped with the hydrophobic features. The hydroxyl was mapped with the hydrogen bond acceptor feature. Forsythoside C was the potential inhibitor of the SQS to form the component formula. The mapping situation between the compound and the pharmacophore was shown in Figure 4. Forsythoside C has a satisfactory mapping situation with the SREBP-2 pharmacophore model.

Figure 4: The mapping results of each TCM compound.

3.3. Virtual Screening of the TCM Compounds

The best models were used to screen the TCMD. The numbers of compounds with potential activity of HMG-CoA reductase and drug-like properties were obtained was 2101. The number of SQS compounds was 1430 and 29 compounds of the SREBP-2 were screened. Then the compounds of HMG-CoA reductase and SQS would be further studied by molecular docking.

3.4. Molecular Docking

The binding site of HMG-CoA reductase was defined with a radius of 8.54 Å, which was the same as the feature generation of pharmacophore. The CDOCKER method got a lower RMSD of 0.37 Å than 2.28 Å of the LibDock. Therefore, the CDOCKER method was chosen for the docking study of HMG-CoA reductase. The interaction of the initial compounds was shown in Figure 5. The initial compounds FBI generated the hydrophobic interactions with the residues CYS561, LEU562, HIS752, LEU853, and LEU857 and formed hydrogen bond interactions with the residues ARG590, ASP690, LYS691, LYS692, SER565, LYS735, ASN686, GLU559, and SER684. The docking of TCM compounds used the same parameter of initial compound.

Figure 5: Interaction between the initial components and the HMG-CoA reductase and SQS.

The binding site of SQS was defined with a radius of 9.32 Å as the feature generation site of pharmacophore model. The RMSD of the CDOCKER was 0.77 Å, which was chosen for the docking study of SQS. The initial compound D99 generated hydrogen bond interactions in SER51, ARG52, SER53, and PHE54. Hydrophobic interactions existed in PHE54, TYR73, VAL179, LEU183, LEU211, PHE288, and PRO292 between D99 and protein. Besides, a Pi-Sulfur bond was generated in CYS289.

Hexahydrocurcumin formed hydrogen bond interactions with SER661, SER684, and GLU559 and formed hydrophobic interactions with LEU853, which has a similar binding mode to the initial compound. It had been reported that the fermented Curcuma longa has regulation of adipogenesis in high-fat diet-induced obese rats [45]. Besides, the extra of Curcuma longa lessens high-fat diet-induced inflammation in subcutaneous adipose tissue with the white pepper [46]. In the molecular docking, Forsythoside C had the hydrophobic bond interactions with PHE288, VAL179, LEU211, LEU183, PHE54, TYR73, and PRO292. LEU183 and PRO292 were also the key amino acids of the initial compound D99. Besides, Forsythoside C and D99 both formed the Pi-Sulfur bond in CYS289, which suggested that the Pi-Sulfur bond may be a necessary interaction between the compounds and SQS. Poncimarin from Aurantii Fructus Immaturus was chosen from the SREBP-2 result. The compound was mapping well with the pharmacophore. Aurantii Fructus Immaturus was a part of the formula with hypercholesterolemia treatment function [47]. The interaction of the TCM compounds was shown in Figure 6.

Figure 6: Interaction between the TCM components and the HMG-CoA reductase and SQS.

3.5. Biological Network Construction

The metabolic network was constructed based on the pathways of HMG-CoA reductase, SQS, and SREBP-2 to elucidate the synergistic action mechanism of compounds. The network was shown in Figure 7. The parameters of the network were calculated to validate the stability and fault tolerance of the network and shown in Table 3.

Table 3: Topological parameters of the network.

Parameters
Network
Parameters
Network

Nodes
705
Edges
659

Connected components
1
Network diameter
26

Network radius
15
Network centralization
0.065

Shortest paths
496320 (100%)
Characteristic path length
10.635

Average number of neighbors
2.357
Network heterogeneity
1.181

Figure 7: Metabolic network of the HMG-CoA reductase, SQS, and SREBP-2.

The connectivity of the network was 100%, which proved the network with a fine connected condition. The network built in this research was in scale-less property because the value Y in the Barabási-Albert model was -1.557.

Through the analysis of parameter, stability and fault tolerance of network can be proven. Then the network was applied intuitively to show the interaction among the three targets and the mechanism of formula synergistic effect.

The network has shown a synergistic effect of the three targets. Cholesterol is generated by 25-step catalytic synthesis as end-product in the pathway. HMG-CoA reductase, an important rate-limiting enzyme in the synthesis process, plays an important role in catalysis to generate a mevalonic acid (MVA). MVA would generate farnesyl pyrophosphate (FPP) through several steps of catalytic synthesis reaction. FPP is the precursor of SQ. The generation of SQ takes place in the catalytic action of SQS. After several steps of enzymatic reaction, SQ eventually generates cholesterol. The inhibition of the enzyme activity can reduce the yield of catalysate. In the pathway, cholesterol levels would significantly reduce when the HMG-CoA reductase and SQS are both inhibited. While the genes expression of HMG-CoA reductase and SQS is regulated by SREBP-2, when SREBP-2 is restrained, HMG-CoA reductase and SQS transcription levels would decline. The reduction of HMG-CoA reductase and SQS transcription would reduce the enzymatic reaction rate. Combining the inhibition effect of three targets, the content of cholesterol would be reduced even more.

Besides, the network has shown that three targets would produce a better effect by affecting other important targets. Figure 7 had shown that the PPARs family was one of the important parts of the network. In this family, PPAR-alpha is one of the most important targets in the treatment of hypercholesterolemia. Its agonist, Fibrates, can reduce the levels of triglycerides by 30-50% [48]. They can also reduce LDL-C potentially. Besides, the levels of HDL-C increased by 5-15% depending on lipid phenotype and baseline concentration [49]. HDL is a significant protective factor in the plasma against coronary heart disease [50]. So the potential action of PPAR-alpha through the network may play an aided role in the overall effect of the component formula for a better effect.

3.6. The Design of TCM Component Formula

The information of the components was shown in Table 4.

Table 4: Information of the potential active TCM components.

Target
Name
Source Plant

SREBP-2
Poncimarin
Aurantii Fructus Immaturus [44]

HMG-CoA reductase
Hexahydrocurcumin
Zingiber officinale, Alpinia officinarum, Curcuma longa

SQS
Forsythoside C
Forsythia suspensa

The three compounds were chosen according to the criteria in Materials and Methods. The source plants of the compounds were commonly used plants in the TCM. The extra of Curcuma longa lessens high-fat diet-induced inflammation in subcutaneous adipose tissue with the white pepper [51]. Aurantii Fructus Immaturus was a part of the formula with hypercholesterolemia treatment function [47].

The combination of the three compounds had synergistic effect to decrease cholesterol levels, which was clearly illuminated by biological network. The mechanism action was separated from network model and shown in Figure 8. Hexahydrocurcumin and Forsythoside C reduce the cholesterol through inhibition of HMG-CoA reductase and SQS. Then poncimarin helps controlling the blood lipid by the PPARs family. Through synergistic effect of three compounds, the inhibition of this pathway can be enhanced by lowering enzyme activity and reducing the content of the enzyme so as to achieve the production of reinforcement in cholesterol synthesis inhibition.

Figure 8: The synergistic effect of TCM component formula.

3.7. The Activity of the TCM Compounds

The test of lipid-lowering effect of hexahydrocurcumin (CAS number: 36062-05-2) with various doses was treated on cancer cell lines and HepG2 cell lines. Cell viability was examined by the MTT assay for 24 h. Then the total triglyceride (TG) was measured. The results indicated that the total triglyceride displayed a dose-dependent inhibition (10, 20, 30, 40, and 50 μmol/L) with 20.11% decrease of TG level in 30 μmol/L (Figure 9). The verification of the TCM compounds on targets level will be further studied.

Figure 9: The TG-decrease effect of hexahydrocurcumin.

4. Conclusion

This research built a TCM component formula formed by the TCM components acting on given targets chosen by virtual screen. A network was built for the interaction analysis. The component formula that contains poncimarin, hexahydrocurcumin, and Forsythoside C can get an enhanced lipid-lowering effect by the combined action of the targets. The component formula can enhance the inhibiting effect of the cholesterol synthesis.

This research hopes to provide a method to build a component formula based on components acting on the known targets. Through selection of targets, component formula can be designed. The components combine the advantages of the virtual screening and the TCM traditional action. Besides, other types of interaction can be shown in the network, which provided more design direction of component formula. Through the analysis of the component formula, the methods may provide a clue to the Traditional Chinese Medicine formula. The traditional formula may generate a complex interaction by the target which the component acted on.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors gratefully acknowledge the support given to this work by the National Natural Science Foundation of China (no. 81173522, no. 81573831, and no. 81430094).

References

J. L. Goldstein, W. R. Hazzard, H. G. Schrott, E. L. Bierman, and A. G. Motulsky, “Hyperlipidemia in coronary heart disease. I. Lipid levels in 500 survivors of myocardial infarction,” The Journal of Clinical Investigation, vol. 52, no. 7, pp. 1533–1543, 1973.

J. L. Goldstein, H. G. Schrott, W. R. Hazzard, E. L. Bierman, and A. G. Motulsky, “Hyperlipidemia in coronary heart disease. II. Genetic analysis of lipid levels in 176 families and delineation of a new inherited disorder, combined hyperlipidemia,” The Journal of Clinical Investigation, vol. 52, no. 7, pp. 1544–1568, 1973.

E. S. Istvan and J. Deisenhofer, “Structural mechanism for statin inhibition of HMG-CoA reductase,” Science, vol. 292, no. 5519, pp. 1160–1164, 2001.

C. P. Cannon, E. Braunwald, C. H. McCabe et al., “Intensive versus moderate lipid lowering with statins after acute coronary syndromes,” The New England Journal of Medicine, vol. 350, no. 15, pp. 1495–1504, 2004.

Y. Kureishi, Z. Luo, I. Shiojima et al., “The HMG-CoA reductase inhibitor simvastatin activates the protein kinase Akt and promotes angiogenesis in normocholesterolemic animals,” Nature Medicine, vol. 6, no. 9, pp. 1004–1010, 2000.

M. Endres, U. Laufs, Z. Huang et al., “Stroke protection by 3-hydroxy-3-methylglutaryl (HMG)-CoA reductase inhibitors mediated by endothelial nitric oxide synthase,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 95, no. 15, pp. 8880–8885, 1998.

J. Pandit, D. E. Danley, G. K. Schulte et al., “Crystal structure of human squalene synthase. A key enzyme in cholesterol biosynthesis,” The Journal of Biological Chemistry, vol. 275, no. 39, pp. 30610–30617, 2000.

M. H. Davidson, “Squalene synthase inhibition: A novel target for the management of dyslipidemia,” Current Atherosclerosis Reports, vol. 9, no. 1, pp. 78–80, 2007.

T. Nishimoto, Y. Amano, R. Tozawa et al., “Lipid-lowering properties of TAK-475, a squalene synthase inhibitor, in vivo and in vitro,” British Journal of Pharmacology, vol. 139, no. 5, pp. 911–918, 2003.

H. Hiyoshi, M. Yanagimachi, M. Ito et al., “Effect of ER-27856, a novel squalene synthase inhibitor, on plasma cholesterol in rhesus monkeys: Comparison with 3-hydroxy-3-methylglutaryl-CoA reductase inhibitors,” Journal of Lipid Research, vol. 41, no. 7, pp. 1136–1144, 2000.

O. P. Flint, B. A. Masters, R. E. Gregg, and S. K. Durham, “Inhibition of cholesterol synthesis by squalene synthase inhibitors does not induce myotoxicity in vitro,” Toxicology and Applied Pharmacology, vol. 145, no. 1, pp. 91–98, 1997.

T. Nishimoto, E. Ishikawa, H. Anayama et al., “Protective effects of a squalene synthase inhibitor, lapaquistat acetate (TAK-475), on statin-induced myotoxicity in guinea pigs,” Toxicology and Applied Pharmacology, vol. 223, no. 1, pp. 39–45, 2007.

J. R. Burnett, “Drug evaluation: TAK-475 - An oral inhibitor of squalene synthase for hyperlipidemia,” Current Opinion in Infectious Diseases, vol. 7, no. 9, pp. 850–856, 2006.

J. D. Horton, J. L. Goldstein, and M. S. Brown, “SREBPs: activators of the complete program of cholesterol and fatty acid synthesis in the liver,” The Journal of Clinical Investigation, vol. 109, no. 9, pp. 1125–1131, 2002.

M. S. Brown and J. L. Goldstein, “The SREBP pathway: regulation of cholesterol metabolism by proteolysis of a membrane-bound transcription factor,” Cell, vol. 89, no. 3, pp. 331–340, 1997.

W. Kong, J. Wei, P. Abidi et al., “Berberine is a novel cholesterol-lowering drug working through a unique mechanism distinct from statins,” Nature Medicine, vol. 10, no. 12, pp. 1344–1351, 2004.

S. M. Grundy, G. L. Vega, Z. Yuan, W. P. Battisti, W. E. Brady, and J. Palmisano, “Effectiveness and tolerability of simvastatin plus fenofibrate for combined hyperlipidemia (the SAFARI trial),” American Journal of Cardiology, vol. 95, no. 4, pp. 462–468, 2005.

E. A. Stein, M. H. Davidson, C. A. Dujovne et al., “Efficacy and tolerability of low-dose simvastatin and niacin, alone and in combination, in patients with combined hyperlipidemia: A prospective trial,” Journal of Cardiovascular Pharmacology and Therapeutics, vol. 1, no. 2, pp. 107–116, 1996.

R. Verpoorte, D. Crommelin, M. Danhof et al., “Commentary: ‘A systems view on the future of medicine: inspiration from Chinese medicine?’,” Journal of Ethnopharmacology, vol. 121, no. 3, pp. 479–481, 2009.

X. Liu, Q. Wang, G. Song, G. Zhang, Z. Ye, and E. M. Williamson, “The classification and application of toxic chinese materia medica,” Phytotherapy Research, vol. 28, no. 3, pp. 334–347, 2014.

D. L. Klayman, “Qinghaosu (Artemisinin): An antimalarial drug from China,” Science, vol. 228, no. 4703, pp. 1049–1055, 1985.

M. Kołaczkowski, M. Nowak, M. Pawłowski, and A. J. Bojarski, “Receptor-based pharmacophores for serotonin 5-HT7R antagonists - Implications to selectivity,” Journal of Medicinal Chemistry, vol. 49, no. 23, pp. 6732–6741, 2006.

K.-H. Kim, N. D. Kim, and B.-L. Seong, “Pharmacophore-based virtual screening: A review of recent applications,” Expert Opinion on Drug Discovery, vol. 5, no. 3, pp. 205–222, 2010.

N. J. Richmond, C. A. Abrams, P. R. N. Wolohan, E. Abrahamian, P. Willett, and R. D. Clark, “GALAHAD: 1. Pharmacophore identification by hypermolecular alignment of ligands in 3D,” Journal of Computer-Aided Molecular Design, vol. 20, no. 9, pp. 567–587, 2006.

C. Y.-C. Chen, “Computational screening and design of traditional Chinese medicine (TCM) to block phosphodiesterase-5,” Journal of Molecular Graphics and Modelling, vol. 28, no. 3, pp. 261–269, 2009.

C.-Y. Chen, Y.-H. Chang, D.-T. Bau et al., “Discovery of potent inhibitors for phosphodiesterase 5 by virtual screening and pharmacophore analysis,” Acta Pharmacologica Sinica, vol. 30, no. 8, pp. 1186–1194, 2009.

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

T. Nakayama, S. Suzuki, H. Kudo, S. Sassa, M. Nomura, and S. Sakamoto, “Effects of three Chinese herbal medicines on plasma and liver lipids in mice fed a high-fat diet,” Journal of Ethnopharmacology, vol. 109, no. 2, pp. 236–240, 2007.

G. B. Zhang, Q. Y. Li, Q. L. Chen, and S. B. Su, “Network pharmacology: a new approach for chinese herbal medicine research,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 621423, 9 pages, 2013.

C. Barillari, G. Marcou, and D. Rognan, “Hot-spots-guided receptor-based pharmacophores (HS-pharm): A knowledge-based approach to identify ligand-anchoring atoms in protein cavities and prioritize structure-based pharmacophores,” Journal of Chemical Information and Modeling, vol. 48, no. 7, pp. 1396–1410, 2008.

S.-Y. Yang, “Pharmacophore modeling and applications in drug discovery: challenges and recent advances,” Drug Discovery Therapy, vol. 15, no. 11-12, pp. 444–450, 2010.

H. S. Kim, M. Ohno, B. Xu et al., “2-Substitution of adenine nucleotide analogues containing a bicyclo[3.1.0]hexane ring system locked in a northern conformation: enhanced potency as P2Y₁ receptor antagonists,” Journal of Medicinal Chemistry, vol. 46, no. 23, pp. 4974–4987, 2003.

G. Mustata, A. V. Follis, D. I. Hammoudeh et al., “Discovery of novel myc-max heterodimer disruptors with a three-dimensional pharmacophore model,” Journal of Medicinal Chemistry, vol. 52, no. 5, pp. 1247–1250, 2009.

O. F. Guner, Pharmacophore Perception, Development, and Use in Drug Design, International University Line, La Jolla, CA, USA, 2000.

Z. Yang, Y. Zhang, X. Wang, and Y. Qiao, “Pharmacophore model generation of P2Y 12 inhibitor,” in Proceedings of the 2012 International Conference on Biomedical Engineering and Biotechnology, iCBEB 2012, pp. 396–399, chn, May 2012.

Z. Zhou, A. K. Felts, R. A. Friesner, and R. M. Levy, “Comparative performance of several flexible docking programs and scoring functions: Enrichment studies for a diverse set of pharmaceutically relevant targets,” Journal of Chemical Information and Modeling, vol. 47, no. 4, pp. 1599–1608, 2007.

B. D. Bursulaya, M. Totrov, R. Abagyan, and C. L. Brooks III, “Comparative study of several algorithms for flexible ligand docking,” Journal of Computer-Aided Molecular Design, vol. 17, no. 11, pp. 755–763, 2003.

D. Croft, G. O'Kelly, G. Wu et al., “Reactome: a database of reactions, pathways and biological processes,” Nucleic Acids Research, vol. 39, no. 1, pp. D691–D697, 2011.

P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software Environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.

M. E. Smoot, K. Ono, J. Ruscheinski, P. L. Wang, and T. Ideker, “Cytoscape 2.8: new features for data integration and network visualization,” Bioinformatics, vol. 27, no. 3, pp. 431-432, 2011.

A. Barabási and Z. N. Oltvai, “Network biology: understanding the cell's functional organization,” Nature Reviews Genetics, vol. 5, no. 2, pp. 101–113, 2004.

J. Caballero, “3D-QSAR (CoMFA and CoMSIA) and pharmacophore (GALAHAD) studies on the differential inhibition of aldose reductase by flavonoid compounds,” Journal of Molecular Graphics and Modelling, vol. 29, no. 3, pp. 363–371, 2010.

G. Wolber, T. Seidel, F. Bendix, and T. Langer, “Molecule-pharmacophore superpositioning and pattern matching in computational drug design,” Drug Discovery Therapy, vol. 13, no. 1-2, pp. 23–29, 2008.

A. Guiotto, P. Rodighiero, and U. Quintily, “Poncimarin, a new coumarin from Poncirus trifoliata L.,” Zeitschrift fur Naturforschung Section C Biosciences, vol. 30, no. 3, p. 420, 1975.

J. H. Kim, O.-K. Kim, H.-G. Yoon et al., “Anti-obesity effect of extract from fermented Curcuma longa L. through regulation of adipogenesis and lipolysis pathway in high-fat diet-induced obese rats,” Food and Nutrition Research, vol. 60, 2016.

A. M. Neyrinck, M. Alligier, P. B. Memvanga et al., “Curcuma longa extract associated with white pepper lessens high fat diet-induced inflammation in subcutaneous adipose tissue,” PLoS ONE, vol. 8, no. 11, Article ID e81252, 2013.

X. Song, J. Wang, P. Wang, N. Tian, M. Yang, and L. Kong, “¹H NMR-based metabolomics approach to evaluate the effect of Xue-Fu-Zhu-Yu decoction on hyperlipidemia rats induced by high-fat diet,” Journal of Pharmaceutical and Biomedical Analysis, vol. 78-79, pp. 202–210, 2013.

B. Staels, J. Dallongeville, J. Auwerx, K. Schoonjans, E. Leitersdorf, and J. Fruchart, “Mechanism of action of fibrates on lipid and lipoprotein metabolism,” Circulation, vol. 98, no. 19, pp. 2088–2093, 1998.

W. M. Freitas, L. A. Quaglia, S. N. Santos et al., “Low HDL cholesterol but not high LDL cholesterol is independently associated with subclinical coronary atherosclerosis in healthy octogenarians,” Aging Clinical and Experimental Research, vol. 27, no. 1, pp. 61–67, 2015.

T. Gordon, W. P. Castelli, M. C. Hjortland, W. B. Kannel, and T. R. Dawber, “High density lipoprotein as a protective factor against coronary heart disease. The Framingham study,” American Journal of Medicine, vol. 62, no. 5, pp. 707–714, 1977.

M. Rana, S. S. Reddy, P. Maurya et al., “Turmerone enriched standardized Curcuma longa extract alleviates LPS induced inflammation and cytokine production by regulating TLR4-IRAK1-ROS-MAPK-NFκB axis,” Journal of Functional Foods, vol. 16, pp. 152–163, 2015.