In silico Approach for Anti-Thrombosis Drug Discovery: P2Y1R Structure-Based TCMs Screening

Cardiovascular diseases (CVDs), including thrombosis, which is induced by platelet aggregation, are the leading cause of mortality worldwide. The P2Y₁ receptor (P2Y₁R) facilitates platelet aggregation and is thus an important potential anti-thrombotic drug target. The P2Y₁R protein structure contains a binding site for receptor antagonist MRS2500 within its seven-transmembrane bundle, which also provides suitable pockets for numerous other ligands to act as nucleotide antagonists of P2Y₁R. The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) comprises 499 Chinese Pharmacopoeia-registered herbs and the structure information for 29,384 ingredients. In silico docking of these compounds into the P2Y₁R protein structure within the MRS2500 pocket can identify potential antithrombotic drugs from natural medicinal plants. Docking studies were performed and scored to evaluate ligand-binding affinities. In this study, a total of 8987 compounds from Traditional Chinese Medicine (TCM) were filtered by Lipinski's rule of five, and their ideal oral-intake properties were evaluated. Of these, 1656 compounds distributed in 443 herbs docked into the P2Y₁R-MRS2500 structure in 16,317 poses. A total of 38 compounds were ranked with a DockScore above 70, and these may have significant potential for development into anti-thrombosis drugs. These computational results suggested that licorice (Glycyrrhiza uralensis Fisch), cimicifugae (Cimicifuga foetida L.), and ganoderma (Ganoderma lucidum Karst) and their chemical constituents, which have not previously been widely used for anti-thrombosis, may have unexpected effects on platelet aggregation. Moreover, two types of triterpene scaffolds summarized from 10 compounds were distributed in these three herbs and also docked into P2Y₁R. These scaffold structures may be utilized for the development of drugs to inhibit platelet aggregation.

Introduction

Cardiovascular disease (CVD) is the leading cause of mortality worldwide. CVD is multifactorial, and its risk factors include stroke, hypertension, arrhythmias, and thrombosis (Mozaffarian et al., 2016). Platelet aggregation-induced thrombosis obstructs blood circulation, playing a central role in acute, and chronic arterial vascular diseases (Radomski et al., 2005). Antiplatelet drugs decrease thrombus formation, and their estimated market is worth 24 billion United States dollars (USD).

G protein-coupled P2Y receptors belong to the nucleotide receptor G protein-coupled receptor (GPCR) family and have eight mammalian subtypes (P2Y_{1, 2, 4, 6, 11-14}) (Kim et al., 2003). P2Y₁ and P2Y₁₂ belong to the human purinergic GPCRs and can be activated by adenosine 5′-diphosphate (ADP) to induce platelet activation (Gurbel et al., 2015). ADP is the first small-molecular weight platelet agonist, and its receptors, such as P2Y₁ receptors, can couple to activated phospholipase C. The activation of serotonin receptor supplements signaling through the P2Y₁ receptor, demonstrating that it is a specific antagonist able to block ADP-induced platelet aggregation (Jin and Kunapuli, 1998). The human P2Y₁ receptor protein structure and its two ligand-binding sites for the nucleotide-like antagonist MRS2500 and allosteric antagonist 1-(2-(2-tert-butylphenoxy)pyridin-3-yl)-3-(4-(trifluoromethoxy)phenyl)urea (BPTU) were reported in 2015 (Protein Database [PDB] ID: 4XNW, 4XNV) (Zhang et al., 2015). (1′R,2′S,4′S,5′S)-4-(2-Iodo-6-methylaminopurin-9-yl)-1-[(phosphato) methyl]-2(phosphato)bicycle[3.1.0]-hexane (MRS2500) is a reported antagonist candidate that exerts its effect via its unique chemical structure. This compound binds the recombinant human P2Y₁ receptor and inhibits the platelet aggregation caused by ADP with an 50% inhibitory concentration (IC50) value in the nanomolar range. It also effectively reduces arterial thrombosis and prolongs bleeding time and has been evaluated as a prototypical antithrombotic agent both ex vivo and in vivo (Hechler et al., 2006). Unlike P2Y₁₂R, P2Y₁R has a highly conserved in class A GPCR residue P229. The pocket for MRS2500 binding to P2Y₁R mainly defined by residues from the N terminus, ECL2, and its helices structures. In P2Y₁R, the antagonist MRS2500 potentially prevents the movements of these helices and stabilized the receptor in an inactive state by interacting with helices In the P2Y₁R–MRS2500 structure, each terminal oxygen of the two phosphates forms at least one contact with the receptor. The hydrogen bonds from 3′-phosphate with Arg195 and Thr201, meanwhile, it is engaged in two salt-bridge interactions with Lys46 at the N terminus. The 5′-phosphate forms a salt-bridge with Thr205 and makes hydrogen bonds with Asp204 and Arg310. P2Y₁R and P2Y₁₂R structures reveal very different features in binding their nucleotide-like ligands even though recognized by the same endogenous ligand ADP. Most significantly, the binding site of MRS2500 in P2Y₁R locates much closer to the extracellular surface than the other known GPCR structures affiliated small-molecule ligand-binding sites. Due to its more safety advantage over the P2Y₁R inhibitors of reducing bleeding liabilities than P2Y₁₂R, it has been suggested to discovery as a whole new drug targets (Gachet, 2008). Moreover, the P2Y1R is also enrolled in other procedure in human body, such as activation of extracellular signal-regulated kinase in astrocytes and vascular inflammation (Zerr et al., 2011). In recent years, many study were performed on P2Y₁₂R, however, the P2Y₁R protein and its ligand pocket crystal structure was firstly reported in 2015. And P2Y₁R has the specificity bind-model and its diversity of regulation mechanisms compared with P2Y₁₂R. This give us a hit that it is easier for other unknown small-molecule can binding to the pocket and thus be a potential antagonist of P2Y₁R.

Traditional Chinese Medicines (TCMs) is a relevant source of biological compounds with innovative mechanisms of action. Although TCMs concepts diverge from Western medicine, concepts, some correlations exist between the two systems. For example, Huoxuehuayu promotes blood circulation and alleviates blood stasis. Several compounds, including ligustrazine, saffionin A, resveratrol, propylgallate, notoginseng triterpenes, crudione, lerulic acid, salvianolic acid A and B, tanshinone IIA, and safflower flavin, which are extracted from the TCM herbs: Salvia miltiorrhiza Bge., Rheum palmatum L., and Ligusticum chuanxiong hort., can inhibit platelet aggregation induced by arachidonic acid (AA), ADP, platelet-activating factor (PAF) and collagen or thrombin both in vivo and in vitro (Chen, 1981). Moreover, because the mechanism of thrombosis involves multiple pathways and requires a multi-target approach for treatment, various small molecule compounds that are not limited to the TCMs activity for the treatment of platelet aggregation can be used to target different proteins, tissues and organs. Therefore, other TCM plants and their active ingredients may act on the P2Y₁ receptor protein to inhibit platelet aggregation via an innovative approach.

Computational approaches have been widely applied to molecular biology and medicinal research. Structure-based methods, such as ligand-protein docking, are efficient and reliable tools for novel drug discovery and design. In silico docking can elucidate the interactions and binding mechanisms between the protein target and its suitable ligands. The Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) comprises all 499 Chinese herbs registered in the Chinese Pharmacopoeia and 29,384 ingredients (Ru et al., 2014). The TCMSP provides a total of twelve absorption, distribution, metabolism, and excretion (ADME)-related properties, including oral bioavailability and half-life.

The objective of this study was to use an in silico docking approach to investigate the potential of compounds in the TCMSP and their distributed herbs as P2Y₁R antagonists and potential new antiplatelet drugs. In this study, we collected all the information of herbs and their contained ingredients from the TCMSP. Then performing the docking procedure of these total 29,384 compounds into the pocket structure of receptor P2Y₁R. Analyzed results of these high efficient docked compounds, summarized there general characters such as scaffold and there distributed herbs will provide a new direction in antiplatelet drugs innovation.

Materials and Methods

Software

The 2-dimensional (2D) and 3D structures of TCM compounds were drawn using ChemBioOffice 2010 (PerkinElmer Inc., Cambridge, MA, USA) and MarvinSketch V15.8 (ChemAxon Ltd., Budapest, Hungary). The P2Y₁ protein and MRS2500 structure (PDB ID: 4XNW) were downloaded from the PDB database (http://www.rcsb.org/pdb) and optimized with Molecular Operating Environment (MOE) 2014.09 (Chemical Computing Group Inc., Montreal, QC, Canada). Candidate compounds were downloaded from the TCMSP database (http://lsp.nwsuaf.edu.cn/tcmsp.php). The chemical structure format was transformed by OpenBabelGUI (O'Boyle et al., 2011) version 2.3.2 (OpenBableGUI 2006 by Chris Morley). Virtual screening was performed using Discovery Studio Client V4.5 (DS 4.5; Accelrys Inc., San Diego, CA, USA). The network and functional analyses were generated using QIAGEN's Ingenuity Pathway Analysis (IPA^®, QIAGEN Redwood City, CA, USA; http://www.qiagen.com/ingenuity).

Lipinski's Rule of Five and ADME/Toxicity (T) Predictions

Those compounds collected in the TCMSP database were evaluated for Lipinski's rule of five (RO5; Lipinski et al., 2012) before screening. Poor pharmacokinetics and toxicity are important causes of costly failures during the downstream process of drug development. This evaluation procedure was performed in DS 4.5, the “filter by Lipinski and Veber Rules” which located in “Small Molecules” module. All setting as default values which individual property thresholds according to the original publications. Their ADME/T properties (van de Waterbeemd and Gifford, 2003) must be predicted and are now computable with in silico approaches, accelerating the drug discovery process. Therefore, we can also filter and design drug-like molecules. The DS 4.5 ADME/T Description module was used to predict the pharmacokinetic properties of all candidate compounds. This program offers reliable predictions for pharmaceutically relevant properties, including aqueous solubility, blood-brain barrier penetration (BBP), cytochrome P450 2D6 inhibition, hepatotoxicity, human intestinal absorption (HIA), and plasma protein binding. The ADME/T Description predicts the physically significant descriptors of pharmaceutically relevant properties of organic molecules. The BBP-95&99 and Absorption-95&99 confidence intervals are demonstrated as four elliptical rings in the ADME/T_AlogP98 & ADME/T_PSA_2D (Y–X axis) coordinate system.

Docking Analysis

The virtual screening docking process was performed with the LigandFit module. LigandFit gives the best poses at the binding site using a stochastic conformational search and the energy of the ligand-protein complex (Sato et al., 2006). Compounds from the TCMSP database were docked to the P2Y₁ protein active site fit for MRS2500 reported by Zhang (Zhang et al., 2015). All procedures were completed using Chemistry at HARvard Molecular Mechanics (CHARMm; Brooks et al., 2009). The 3D multi-conformational compound molecular structures were generated using Monte Carlo algorithm-based conformational analysis and rigid body minimization completed using the Smart Minimizer module. The DockScore was adopted as the LigandFit scoring to rank the TCMSP compounds that docked into the P2Y₁R ligand-binding pocket for MRS2500. The LigandFit docking procedure consists of the following two steps: (1) identify and select the region of the protein as the active site for docking by cavity detection and (2) dock the candidate ligands to the selected site. The docking cavity was defined using the DS site search module. For all potential drugs, the docking site was derived from the position of the MRS2500 and P2Y₁ co-crystallized construction (PDB ID: 4XNW), and the grid resolution was set to 0.5 Å (default). The ligand-accessible grid was defined as the minimum distance between a grid point and the protein and was 2.0 Å for hydrogen and 2.5 Å for heavy atoms. This confirmed binding site was used to calculate the non-bonded interactions between all potential compounds and the P2Y₁ receptor protein residues. The docking procedure was initiated with the generation of random ligand conformations. The following procedures were performed after a new conformation was generated: The shapes of the potential ligands were compared with the active site, and if the result was acceptable, the dock energy (DockScore) was computed between the protein and ligand trial conformation. Variable numbers of Monte Carlo steps were used to generate different ligand conformations. Scoring was performed with six scoring functions: DockScore, LigScore 1 and 2, Piecewise Linear Potential (PLP) 1 and 2, and Potential of Mean Force (PMF). We assumed that the bioactive orientations ranked by DockScore and other scoring functions were used to retain each molecule and compute the enrichment factors.

Results and Discussion

Screening and Structure Analysis

Figure 1 summarizes the technology roadmap and the screening results after each step. In 1997, Lipinski led to the well-known “Rule of Five” (RO5) for selecting drug-like molecules, which demonstrated orally administered drugs should have good oral bioavailability in order to be effective. According to the “Rule of Five,” a drug-like molecule should have no more than one of the following criterias: No more than five hydrogen bond donors, no more than 10 hydrogen bond acceptors, a molecular mass <500 daltons, and an octanol-water partition coefficient log P not >5. To ensure this study with a reliable and available results. All compounds need to filtered by the RO5. From the 29,384 native TCM compounds collected from the TCMSP database, 8987 compounds were filtered by Lipinski's rule of five and prepared for further study. Successful docking to P2Y₁R was observed for 1656 compounds (16,317 poses), and their ADME/T descriptors were calculated (Figure 2, Supplementary Datasheet 1). According to the DS 4.5 suggested theory, the Human Intestinal Absorption (HIA) model can embodied in ADMET plot the validate the predictions. The ADMET—Human Intestinal Absorption model predicts human intestinal absorption (HIA) after oral administration. Intestinal absorption is defined as a percentage absorbed rather than as a ratio of concentrations. A well-absorbed compound is one that is absorbed at least 90% into the bloodstream in humans. These levels are defined by the 95% (red line) and 99% (green line) confidence ellipsoids in Figure 2. For the ADME/T_PSA_2D, ADME/T_AlogP98 plane, two confidence ellipsoids of 95 and 99% for oral absorption surrounded most of these natural compounds. Two ellipsoids of 95 and 99% for BBP with high penetrations further refined the range of compounds suitable for clinical development as oral-intake drugs. Total view of these 1656 compounds, majority number within an acceptable range which indicated that most candidates compounds have an ideal human intestinal absorption which can further refined the range of compounds suitable for clinical development as oral-intake drugs. Candidate ligand poses can be evaluated and prioritized according to the DockScore function. Each docked compound has an optimal pose indicated by its highest DockScore. A total of 38 compounds had a DockScore above 70.0, and these presented high potential for further studies (Table 1).

FIGURE 1

Figure 1. Methodology roadmap and major results for P2Y₁R-related anti-thrombosis drug design from the TCMSP database.

FIGURE 2

Figure 2. ADME/T prediction plots for compounds that successfully docked in P2Y₁R (1656) from among the compounds in the TCMSP database.

TABLE 1

Table 1. Docking results of the top 38 (DockScore > 70) and their plant sources.

Traditional Chinese Herb Analysis

Although many TCMs are considered useful for treating CVDs, each herb contains numerous compounds with different P2Y₁R DockScores. A total DockScore was obtained for each herb by summing the DockScores of its compounds as recorded in the TCMSP (Figure S1). Figure 3 shows the 10 herbs that possessed the highest total DockScores. The IPA indicated that 1890 compounds present in these 10 herbs may be considered drug candidates for hematological diseases and CVDs (Figure S2).

FIGURE 3

Figure 3. Top 10 herbs selected by estimating the number of hit compounds (orange area between bars), max DockScore (pink line) and total DockScore (blue bars).

Additional GPCRs were mapped by the software. Hydroxycarboxylic acid receptor 2 (HCAR2), a common target for the widely prescribed dyslipidemia drug niacin, is associated with platelet aggregation inhibition (Pike, 2005). In addition, data indicate that the proto-oncogene tyrosine-protein kinase (SRC), phosphoinositide 3-kinase (PI3K), mitogen-activated protein kinase (MAPK), N-methyl-D-aspartate receptor (NMDAR), receptor tyrosine-protein kinase erbB-2 (ERBB2), nuclear factor kappa-light-chain-enhancer of activated B cells (NF-kB), and protein kinase a (Pka) pathways (Figure 4) are activated by the top 10 herbs.

FIGURE 4

Figure 4. IPA-mapped CVD pathway for the 1890 compounds distributed in the top 10 herbs.

Candidate Compound Analysis

Sixteen compounds were associated with platelet aggregation, and all of the relevant herbs were validated in the literature, as shown in Table 2. Among them, ascorbic acid, also known as vitamin C, is a well-known dietary antioxidant that inactivates oxygen free radicals and exerts a variety of protective effects on blood vessels and platelets (Yang et al., 2015). Vitexin and caffeic acid showed significant antiplatelet aggregation activity and were relatively selective inhibitors of platelet aggregation induced by AA and collagen (Afifi and Abu-Dahab, 2012; Lee and Bae, 2015). Rutin, naringin and quercetin, three flavonoid glycosides, were metabolized to phenolic acids via aglycones by human intestinal microflora and have been considered natural prodrugs with anti-platelet activity since the 1990s (Ndhlala et al., 2015). Esculetin is a derivative of coumarin and exists as glycosides and caffeic acid conjugates in Cimicifugae Rhizoma (Cimicifuga foetida L.) and many other medicinal plants. The systemic usage of esculetin-containing preparations has an anticoagulant effect and may interact with anticoagulant drugs, such as warfarin (Karnewar et al., 2016). Curcumin, a major component of turmeric, is also found in Cimicifuga foetida L. and inhibits arachidonate-, adrenaline-, and collagen-induced platelet aggregation. Curcumin inhibited exogenous arachidonate-produced thromboxane B₂ (TXB₂) in washed platelets and the incorporation of [¹⁴C]AA into platelet phospholipids. Moreover, it also inhibited the deacylation of AA-labeled phospholipids (Du et al., 2016). Citric acid and L-malic acid, two main organic acids in Fructus Choerospondiatis [Choerospondias axillaris (Roxb.) Burtt et Hill] and Lycii Fructus (Lycium chinense Mill.), possess antioxidant, anti-inflammatory, and antiplatelet aggregation activities (Tang et al., 2013).

TABLE 2

Table 2. Compounds related to platelet aggregation extracted from the top 10 herbs.

Cimicifugae, Ganoderma, and Licorice Analysis

According to the dual analyses of results related to both these compounds and plants, we focused on three particular TCM herbs, cimicifugae (Cimicifuga foetida L.), ganoderma (Ganoderma lucidum Karst), and licorice (Glycyrrhiza uralensis Fisch).

Cimicifugae (Cimicifuga foetida L.) (also known as black cohosh, black snakeroot, and rattlesnake root) has a long tradition of use as a beneficial herbal remedy for joints, muscles, nerve aches, and gynecological disorders in both China and North America. Cimicifugae is used as an anti-inflammatory, analgesic and antipyretic remedy in TCM. However, there are few reports of its clinical use in anti-thrombosis. Sodium ferulate (SF) is an active principle compound from cimicifugae and has been used in TCM and approved by the State Drug Administration of China for the treatment of cardiovascular and cerebrovascular diseases. SF has antithrombotic, platelet aggregation inhibitory, and antioxidant activities in both animals and humans.

Ganoderma (Ganoderma lucidum Karst) is a popular medicinal mushroom and has been used in TCM in Asia for the past two millennia to treat numerous diseases. Its regular consumption is believed to preserve human vitality and promote longevity. Ganoderma has been used worldwide, and its constituents have been extensively identified; however, the biological potency of this fungus has not yet been sufficiently investigated. According to clinical observations and pharmacologic experiments, ganoderma can improve microcirculation in organs and increase the amount of blood. Moreover, in recent years, a series studies focused on ganoderma and anti-thrombosis (Kawagishi et al., 1993). The administration of this herb in vitro and in vivo was found to significantly inhibit ADP-induced platelet aggregation. These inhibitory effects may involve the release of intrinsic ADP from platelets or the metabolites of AA and PAF. Other research showed that purified ganoderma fibrinolytic protease protected mice against thrombotic death or paralysis induced by collagen and epinephrine and also activated the partial thromboplastin time (APTT) and thrombin time (TT) in rat platelets (Kumaran et al., 2011). The polysaccharide of ganoderma was also found to prolong the blood clotting and bleeding time in rats and decrease the thrombosis and fibrinogen level in blood plasma (Klupp et al., 2016). Evidence indicates that ganoderma has potential as an antithrombotic herb because of its antiplatelet activity.

Licorice (Glycyrrhiza uralensis Fisch; Chinese name Gan-cao) is one of the most indispensable crude drugs in TCM preparations (constituting ~60%), and its wide clinical use widely began in ancient Egyptian, Greek and China in the Second century B.C. Licorice comprises the roots of three Glycyrrhiza species (Glycyrrhiza uralensis Fisch, Glycyrrhiza glabra L., and Glycyrrhiza inflata Batalin). In traditional clinical treatment, this botanical drug is typically used for allergic-inflammatory disease, gastrointestinal problems, cancer, CVD, and bladder and kidney ailments. Approximately 300 diverse compounds have been isolated from licorice, including triterpene saponins and flavonoids, which are responsible for the antiviral, anti-inflammatory, antitumor, antitussive, anti-oxidant, antispasmodic, and metabolic syndrome preventive activities. However, evidence suggests that licorice has anti-thrombosis activity and can be used to treat CVDs.

These three herbs contained 14 compounds with a DockScore above 70 (Table 3). By calculating the binding energies between P2Y₁R and these 14 compounds, the binding efficiency and loss of conformational entropy and energy of a bound ligand can be obtained. These compounds have different ranges of ligand energies, Van der Waals energies and bond energies compared with the ligand MRS2500. However, according to CHARMm command-based calculations, all 14 have similar rotational and translational entropy values, indicating that these compounds all have the same effective binding efficiency as the original ligand MRS2500 (Figure 5).

TABLE 3

Table 3. Chemical information of candidate compounds from Cimicifuga foetida L., Ganoderma lucidum Karst, and Glycyrrhiza uralensis Fisch.

FIGURE 5

Figure 5. Comparison of 14 compounds in the binding model structure with MRS2500.

A drug-likeness (DL) analysis was performed to determine whether these compounds have potential uses as drugs. The DL index values of these compounds were calculated using the Tanimoto coefficient and are presented in the TCMSP database. A molecule with DL > 0.18 was considered a drug-like compound and was selected for the following processes. Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of a wide range of substrates in the liver, and its inhibition is related to most drug-drug interactions. The CYP2D6 model can predict candidate compounds for CYP2D6 enzyme inhibition; therefore, investigating CYP2D6 inhibition is a required part of the regulatory procedures involved in the drug discovery and development process. Discovery Studio ADME/T Descriptors contain a computational model for compounds inhibiting the CYP2D6 enzyme, which was developed from a training set of 151 structurally diverse compounds with known CYP2D6 inhibition constants. This experiment was performed with modified Bayesian learning, and all 14 compounds scored below the cutoff Bayesian score of 0.161 and were not classified as CYP2D6 inhibitors. Therefore, these compounds are unlikely to have drug-drug interactions and can improve the efficiency and safety of drugs, which has significant clinical importance. The hepatotoxicity model can predict potential organ toxicity and was developed from the available literature data for 436 compounds known to exhibit liver toxicity or trigger dose-related elevated aminotransferase levels in more than 10% of the human population. For these 14 compounds, eight had no hepatotoxic activity, including four compounds from licorice (Table 4).

TABLE 4

Table 4. P2Y₁R binding energies for compounds PD-1~PD-14.

The scaffold and pharmacophore are two basic elements of drug molecules. The pharmacophore consists of discrete atoms, groups, and fragments responsible for the physic-chemical features, and their arrangement controls the specific activities. However, the pharmacophore must incorporate an integrated scaffold to form an actual molecule, and thus, implanting the same pharmacophore on various scaffolds results in structurally diverse compounds that act on the same biological target. The promiscuity of receptors indicates that scaffolds can be rationally transformed based on the flexibility and plasticity of the receptors, thereby imbuing receptor binding sites with diversity and variability. Compared to traditional high-throughput screening methods, the scaffold-based approach can facilitate higher initial synthetic efforts guided by the co-crystal structure.

Among these 14 compounds, two triterpene scaffolds are found in 10 compounds. Scaffold 1 was found in compounds PD-4, 5, 6, 7, 8, and 9, whereas scaffold 2 was found in compounds PD-11, 12, 13, and 14 (Figure 6). These two scaffolds could be further developed into more potent and selective drug-like ligands by an iterative process including co-crystallography of the complex constructed by candidate compounds with the P2Y₁R protein. These candidate compounds could subsequently be used for chemical syntheses followed by structure-guided computational design. An ideal DL scaffold should not only form key interactions with the protein target but also possess a conserved binding mode that can tolerate small substitutions.

FIGURE 6

Figure 6. Two types of scaffolds in Cimicifuga foetida L., Ganoderma lucidum Karst, and Glycyrrhiza uralensis Fisch.

Conclusion

Platelet aggregation occurs when platelets adhere to each other at sites of vascular injury. This process has long been recognized as critical for hemostatic plug formation and thrombosis, which is the final product of the aggregated platelets that formed a platelet plug. Vitamin K antagonists and warfarin are often taken orally to reduce thromboembolic occurrence; however, both of these drugs have an unwanted side effect of an increased risk of bleeding.

Uncontrolled platelet aggregation is critical in thrombosis, which may cause heart attacks, unstable angina and other systematic disorders. The inhibition of platelet aggregation is fundamentally important in the treatment and prevention of CVDs. Platelet aggregation can regulated by G protein-coupled P2Y receptor families, particularly P2Y₁₂ and P2Y₁ (Aslam et al., 2013). P2Y₁ is required for platelet shape change in response to ADP and is also a principal receptor in mediating both physiological and pathological ADP-induced platelet aggregation. At higher ADP concentrations, P2Y₁-deficient platelets become partially aggregated. In vivo, the bleeding time was increased by the lack of P2Y₁ expression, which could protect against collagen- and ADP-induced thrombosis. It is sufficient to block ADP-induced platelet aggregation by utilize specific antagonists inhibit the signaling through P2Y1 receptor (Jin and Kunapuli, 1998).

MRS2500 is a new antiplatelet drug that has strong antithrombotic activity in systemic thromboembolism induced by adrenaline infusion. In P2Y₁R, the antagonist MRS2500 may prevent the movement of P2Y₁R helix structures and stabilize the receptor in an inactive state. The nucleotide-like antagonist MRS2500 of P2Y₁R has numerous anchor points on different receptor domains in hydrophilic and charged sites. The interactions of both P2Y₁ and P2Y₁₂R with their nucleotide ligands reveal disparate binding modes between these two P2YR subfamilies, broadening our understanding of the diversity of signal-recognition mechanisms in GPCRs and providing additional potential ligands and diversity binding models to regulate the proteins of these subfamilies. MRS2500 has is an ideal antagonist of P2Y₁R that can prevent platelet aggregation caused by thrombosis.

The pharmacological actions of TCMs activate blood and resolve stasis to treat thrombosis. Many Traditional Chinese herbs and Chinese patent medicines are widely applied in clinical use to treat thrombosis and have good prospects in Chinese Medicine. These drugs are used for promoting blood circulation for removing blood stasis (PBCRBS) and are involved in the clinical treatment of cardiovascular-related diseases, particularly as anti-thrombosis agents. Compared to Western medicines, TCM herbs have advantages of mild action, multiple pathways, multiple targets, and fewer adverse reactions. According to the network pharmacology analysis of PBCRBS herbs, the most frequently used TCMs contain Danshen (Salvia miltiorrhiza Bge.), Chuanxiong (Ligusticum chuanxiong hort), Honghua (Carthamus tinctorius L.), Sanqi (Panax notoginseng Burk), and Danggui [Angelica sinensis (Oliv.) Diels]. However, analyzing these herbs and their compounds revealed that they mainly exerted anti-thrombosis activity by affecting the nitric oxide synthase 2 (NOS2), prostaglandin-endoperoxide synthase 2 (PTGS2), tumor necrosis factor (TNF), and interleukin 1 beta (IL-1β) pathways. Currently, no research into the application of these herbs/compounds to P2Y₁ receptors is available. Because the construction of MRS2500 and the combination model between MRS2500 and P2Y₁R has been established, numerous TCMs can be utilized to treat thrombosis in a novel and mild manner. A combination of antiplatelet dugs is more efficient than a single drug at inhibiting the multiple pathways of platelet activation. Furthermore, natural compounds extracted from herbs have fewer side effects than synthetic compounds. We investigated the herbs and their compounds with anti-platelet aggregation activities similar to those of MRS2500.

Overall, in this study, a total of 8987 compounds from TCMs filtered by Lipinski's rule of five were predicted have ideal oral-intake properties. Among them, 1656 compounds distributed in 443 herbs were able to dock into the P2Y₁R-MRS2500 structure in 16,317 poses. A total of 38 compounds had a DockScore above 70.0 and may have significant potential for development into anti-thrombosis drugs. We also investigated Cimicifuga foetida L., Ganoderma lucidum Karst and Glycyrrhiza uralensis Fisch. These herbs and their chemical constituents have not been widely used in anti-thrombosis treatment and may have potential as new drugs for the inhibition of platelet aggregation.

Author Contributions

HL and PX designed the experiment. FY performed the experimental computing, data analysis, figure drawing, and paper writing. LS, YP, CH, and LX collected the compound information.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This study was supported by the CAMS Innovation Fund for Medical Sciences (CIFMS) ID: 2016-I2M-3-015 and 2016-I2M-1-012. We also very grateful to Dr. Jun Xu (Research Center for Drug Discovery, School of Pharmaceutical Sciences, Sun Yat-sen University, Guangzhou) for his valuable suggestions. We would like to thank Mr. Ke-hao Wu (College of Life Science, Hunan Normal University, Changsha) for his assistance with data processing and Miss Xi-chun Wei (C+D Design Research Center, DC ALLIANCE Company, Shanghai) for her assistance with figure drawing.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphar.2016.00531/full#supplementary-material

Supplementary Data Sheet 1. Total original DockScores of compounds collected from TCMSP.

Figure S1. Total plants DockScore results analysis. (Blue columns represent total DockScore and Orange polyline represents the tendency of the highest DockScore of each plant).

Figure S2. Canonical pathway analysis calculated in Ingenuity Pathway Analysis (IPA).

References

Afifi, F. U., and Abu-Dahab, R. (2012). Phytochemical screening and biological activities of Eminium spiculatum (Blume) Kuntze (family Araceae). Nat. Prod. Res. 26, 878–882. doi: 10.1080/14786419.2011.565558

PubMed Abstract | CrossRef Full Text | Google Scholar

Aslam, M., Sedding, D., Koshty, A., Santoso, S., Schulz, R., Hamm, C., et al. (2013). Nucleoside triphosphates inhibit ADP, collagen, and epinephrine-induced platelet aggregation: role of P2Y1 and P2Y12 receptors. Phytochemistry 32, 239–241. doi: 10.1016/j.thromres.2013.08.021

PubMed Abstract | CrossRef Full Text

Brooks, B. R., Brooks, C. L. III, Mackerell, A. D. Jr., Nilsson, L., Petrella, R. J., Roux, B., et al. (2009). CHARMM: the biomolecular simulation program. J. Comput. Chem. 30, 1545–1614. doi: 10.1002/jcc.21287

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, K. J. (1981). Certain progress in the treatment of coronary heart disease with traditional medicinal plants in China. Am. J. Chin. Med. 9, 193–196. doi: 10.1142/S0192415X81000251

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, L., Feng, X., Xiang, X., and Jin, Y. (2016). Wound healing effect of an in situ forming hydrogel loading curcumin-phospholipid complex. Curr. Drug Deliv. 13, 76–82. doi: 10.2174/1567201813666151202195437

PubMed Abstract | CrossRef Full Text | Google Scholar

Gachet, C. (2008). P2 receptors, platelet function and pharmacological implications. Thromb. Haemost. 99, 466–472. doi: 10.1160/th07-11-0673

PubMed Abstract | CrossRef Full Text | Google Scholar

Gurbel, P. A., Kuliopulos, A., and Tantry, U. S. (2015). G-protein–coupled receptors signaling pathways in new antiplatelet drug development. Arterioscler. Thromb. Vasc. Biol. 35, 500–512. doi: 10.1161/ATVBAHA.114.303412

PubMed Abstract | CrossRef Full Text | Google Scholar

Hechler, B., Nonne, C., Roh, E. J., Cattaneo, M., Cazenave, J. P., Lanza, F., et al. (2006). MRS2500 [2-iodo-N6-methyl-(N)-methanocarba-2′-deoxyadenosine-3′,5′-bisphosphate], a potent, selective, and stable antagonist of the platelet P2Y1 receptor with strong antithrombotic activity in mice. J. Pharmacol. Exp. Ther. 316, 556–563. doi: 10.1124/jpet.105.094037

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, J., and Kunapuli, S. P. (1998). Coactivation of two different G protein-coupled receptors is essential for ADP-induced platelet aggregation. Proc. Natl. Acad. Sci. U.S.A. 95, 8070–8074. doi: 10.1073/pnas.95.14.8070

PubMed Abstract | CrossRef Full Text | Google Scholar

Karnewar, S., Vasamsetti, S. B., Gopoju, R., Kanugula, A. K., Ganji, S. K., Prabhakar, S., et al. (2016). Mitochondria-targeted esculetin alleviates mitochondrial dysfunction by AMPK-mediated nitric oxide and SIRT3 regulation in endothelial cells: potential implications in atherosclerosis. Sci. Rep. 6:24108. doi: 10.1038/srep24108

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawagishi, H., Fukuhara, F., Sazuka, M., Kawashima, A., Mitsubori, T., and Tomita, T. (1993). 5′-deoxy-5′-methylsulphinyladenosine, a platelet-aggregation inhibitor from ganoderma-lucidum. Phytochemistry 32, 239–241. doi: 10.1016/S0031-9422(00)94974-4

CrossRef Full Text | Google Scholar

Kim, H. S., Ohno, M., Xu, B., Kim, H. O., Choi, Y., Ji, X. D., et al. (2003). 2-Substitution of adenine nucleotide analogues containing a bicyclo [3.1.0] hexane ring system locked in a northern conformation: enhanced potency as P2Y1 receptor antagonists. J. Med. Chem. 46, 4974–4987. doi: 10.1021/jm030127+

PubMed Abstract | CrossRef Full Text | Google Scholar

Klupp, N. L., Kiat, H., Bensoussan, A., Steiner, G. Z., and Chang, D. H. (2016). A double-blind, randomised, placebo-controlled trial of Ganoderma lucidum for the treatment of cardiovascular risk factors of metabolic syndrome. Sci. Rep. 6:29540. doi: 10.1038/srep29540

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumaran, S., Palani, P., Nishanthi, R., and Kaviyarasan, V. (2011). Studies on screening, isolation and purification of a fibrinolytic protease from an isolate (VK12) of Ganoderma lucidum and evaluation of its antithrombotic activity. Med. Mycol. J. 52, 153–162. doi: 10.3314/jjmm.52.153

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, W., and Bae, J. S. (2015). Antithrombotic and antiplatelet activities of orientin in vitro and in vivo. J. Funct. Food 17, 388–398. doi: 10.1016/j.jff.2015.05.037

CrossRef Full Text | Google Scholar

Lipinski, C. A., Lombardo, F., Dominy, B. W., and Feeney, P. J. (2012). Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 64, 4–17. doi: 10.1016/j.addr.2012.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Mozaffarian, D., Benjamin, E. J., Go, A. S., Arnett, D. K., Blaha, M. J., Cushman, M., et al. (2016). Heart disease and stroke statistics–2016 update a report from the american heart association. Circulation 133, e599. doi: 10.1161/CIR.0000000000000366

PubMed Abstract | CrossRef Full Text

Ndhlala, A. R., Ghebrehiwot, H. M., Ncube, B., Aremu, A. O., Gruz, J., Šubrtová, M., et al. (2015). Antimicrobial, anthelmintic activities and characterisation of functional phenolic acids of Achyranthes aspera linn.: a medicinal plant used for the treatment of wounds and ringworm in east africa. Front. Pharmacol. 6:274. doi: 10.3389/fphar.2015.00274

PubMed Abstract | CrossRef Full Text | Google Scholar

O'Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., and Hutchison, G. R. (2011). Open Babel: an open chemical toolbox. J. Cheminform. 3, 1–14. doi: 10.1186/1758-2946-3-33

PubMed Abstract | CrossRef Full Text | Google Scholar

Pike, N. B. (2005). Flushing out the role of GPR109A (HM74A) in the clinical efficacy of nicotinic acid. J. Clin. Invest. 115, 3400–3403. doi: 10.1172/JCI27160

PubMed Abstract | CrossRef Full Text | Google Scholar

Radomski, A., Jurasz, P., Alonso-Escolano, D., Drews, M., Morandi, M., Malinski, T., et al. (2005). Nanoparticle-induced platelet aggregation and vascular thrombosis. Brit. J. Pharmacol. 146, 882–893. doi: 10.1038/sj.bjp.0706386

PubMed Abstract | CrossRef Full Text | Google Scholar

Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 1–6. doi: 10.1186/1758-2946-6-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Sato, H., Shewchuk, L. M., and Tang, J. (2006). Prediction of multiple binding modes of the CDK2 inhibitors, anilinopyrazoles, using the automated docking programs GOLD, FlexX, and LigandFit: an evaluation of performance. J. Chem. Inf. Model. 46, 2552–2562. doi: 10.1021/ci600186b

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, X., Liu, J., Dong, W., Li, P., Li, L., Lin, C., et al. (2013). The cardioprotective effects of citric acid and L-malic acid on myocardial ischemia/reperfusion injury. Evid. Based Complement Alternat. Med. 2013, 1–11. doi: 10.1155/2013/820695

PubMed Abstract | CrossRef Full Text | Google Scholar

van de Waterbeemd, H., and Gifford, E. (2003). ADMET in silico modelling: towards prediction paradise? Nat. Rev. Drug Discov. 2, 192–204. doi: 10.1038/nrd1032

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, M., Bhopale, V. M., and Thom, S. R. (2015). Ascorbic acid abrogates microparticle generation and vascular injuries associated with high-pressure exposure. J. Appl. Physiol. 119, 77–82. doi: 10.1152/japplphysiol.00183.2015

PubMed Abstract | CrossRef Full Text | Google Scholar

Zerr, M., Hechler, B., Freund, M., Magnenat, S., Lanois, I., Cazenave, J. P., et al. (2011). Major contribution of the P2Y1 receptor in purinergic regulation of TNF α-induced vascular inflammation. Circulation 123, 2404–2413. doi: 10.1161/CIRCULATIONAHA.110.002139

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, D., Gao, Z. G., Zhang, K., Kiselev, E., Crane, S., Wang, J., et al. (2015). Two disparate ligand-binding sites in the human P2Y1 receptor. Nature 520, 317–321. doi: 10.1038/nature14287

PubMed Abstract | CrossRef Full Text | Google Scholar