Modeling drug response using network-based personalized treatment prediction (NetPTP) with applications to inflammatory bowel disease

For many prevalent complex diseases, treatment regimens are frequently ineffective. For example, despite multiple available immunomodulators and immunosuppressants, inflammatory bowel disease (IBD) remains difficult to treat. Heterogeneity in the disease across patients makes it challenging to select the optimal treatment regimens, and some patients do not respond to any of the existing treatment choices. Drug repurposing strategies for IBD have had limited clinical success and have not typically offered individualized patient-level treatment recommendations. In this work, we present NetPTP, a Network-based Personalized Treatment Prediction framework which models measured drug effects from gene expression data and applies them to patient samples to generate personalized ranked treatment lists. To accomplish this, we combine publicly available network, drug target, and drug effect data to generate treatment rankings using patient data. These ranked lists can then be used to prioritize existing treatments and discover new therapies for individual patients. We demonstrate how NetPTP captures and models drug effects, and we apply our framework to individual IBD samples to provide novel insights into IBD treatment.

For many prevalent complex diseases, treatment regimens are frequently ineffective. For example, despite multiple available immunomodulators and immunosuppressants, inflammatory bowel disease (IBD) remains difficult to treat. Heterogeneity in the disease across patients makes it challenging to select the optimal treatment regimens, and some patients do not respond to any of the existing treatment choices. Drug repurposing strategies for IBD have had limited clinical success and have not typically offered individualized patient-level treatment recommendations. In this work, we present NetPTP, a Network-based Personalized Treatment Prediction framework which models measured drug effects from gene expression data and applies them to patient samples to generate personalized ranked treatment lists. To accomplish this, we combine publicly available network, drug target, and drug effect data to generate treatment rankings using patient data. These ranked lists can then be used to prioritize existing treatments and discover new therapies for individual patients. We demonstrate how NetPTP captures and models drug effects, and we apply our framework to individual IBD samples to provide novel insights into IBD treatment.

Author summary
Offering personalized treatment results is an important tenant of precision medicine, particularly in complex diseases which have high variability in disease manifestation and treatment response. We have developed a novel framework, NetPTP (Network-based Personalized Treatment Prediction), for making personalized drug ranking lists for patient samples. Our method uses networks to model drug effects from gene expression data and applies these captured effects to individual samples to produce tailored drug treatment rankings. We applied NetPTP to inflammatory bowel disease, yielding insights into the treatment of this particular disease. Our method is modular and generalizable, and thus

Introduction
Drug development is an expensive and lengthy endeavor, on average costing approximately a billion dollars to successfully bring a drug to market [1]. As such, drug repurposing, also known as drug repositioning, has become an important avenue for discovering existing treatments for new indications, saving time and money in the quest for new therapies. With increasing data available on drugs and diseases, computational approaches for drug repositioning have shown great potential by integrating multiple sources of information to discover novel matchings of drugs and diseases. Using transcriptomic data, multiple existing computational approaches for drug repurposing are based on constructing representations of diseases and drugs and assessing their similarity. For example, Li and Greene et al used differentially expressed genes to construct and compare disease and drug signatures and van Noort et al applied a similar approach using 500 probe sets in colorectal cancer [2,3]. However, by representing the disease as an aggregate, these methods can be limited in their ability to capture patient and disease heterogeneity. Furthermore, by treating each gene or probe set individually, these methods frequently fail to capture different combinations of perturbations that cause similar disease phenotypes, which contributes to disease heterogeneity. For complex, heterogeneous diseases, there are frequently multiple avenues of treatment targeting different aspects of the disease, and many patients do not respond to the same set of therapies. Such diseases could benefit from a generative method that produces more personalized therapeutic strategies that target an individual's disease state.
One such condition is inflammatory bowel disease (IBD), which consists of two main subtypes, ulcerative colitis (UC) and Crohn's disease (CD). Both are chronic inflammatory conditions of the gastrointestinal system which together affect over 1.5 million people in the United States [4]. As a heterogeneous disease, different IBD patients frequently respond to different treatment drugs that target specific pathways unique to the disease pathogenesis seen in that particular patient. As such, there currently exist multiple different treatments for IBD which have different mechanisms of action, such as sulfasalazine, infliximab, azathioprine, and steroids [5]. However, it is frequently unclear which patients would derive the most benefit from each of these classes of drugs. Furthermore, many patients do not respond or develop nonresponse to these therapies, resulting in escalation of their treatment regimens or surgery.
There exist a few previous computational repurposing methods that have been applied to IBD. For example, Dudley et al compared drugged gene expression signatures from the Connectivity Map (CMap) to IBD gene expression data identified topiramate as a potential therapeutic candidate [6]. Another approach overlapped IBD genes implicated in genome wide association studies with known drug targets for IBD [7]. More recently, newer approaches have incorporated gene interactions by examining sets of genes in the same pathway. For example, Grenier et al employed a pathway-based approach using genetic loci from IBD gene wide association studies and pathway set enrichment analysis to identify new candidate drugs [8]. While these methods have yielded some new potential therapies, there is still a great need for identifying responders and for additional therapeutic strategies for nonresponders.
We present Network-based Personalized Treatment Prediction (NetPTP), a novel systems pharmacological approach for modeling drug effects, which incorporates the interactions of genes and proteins with drug targets in order repurpose and prioritize drugs in individual patients. Using publicly available human and mouse gene expression data, we show how our approach can be used to identify drugs based on each patient's disease profile. We further aggregate these individual results to the disease level to gain new insights into the treatment of CD and UC.

Ethics statement
The experimental portion of this study (protocol ID 32963) was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) at Stanford University.

Data preprocessing
The connectivity map. We downloaded all instances in CMap as raw CEL files. We first preprocessed the instances using robust multi-array average [9]. Probes were then mapped to genes using the corresponding platform annotation files. For genes corresponding to multiple probes, the expression values were averaged across all probes corresponding to the same gene. We then corrected for batch effects using ComBat [10] with the sva package [11], correcting for batch number as given by CMap and for cell line. Only instances which mapped to a drug in DrugBank and had at least one gene target measured were included in downstream analysis.
Human IBD disease datasets. We processed four publicly available IBD colonic sample datasets from GEO: GSE16879 [12], GSE9686 [13], GSE10616 [14], and GSE36807 [15] (S1 Table). Each study contains baseline colonic sample prior to treatment of CD and UC patients, and also contains healthy control colonic samples. Of note, GSE16879 consists of patients refractory to corticosteroids and/or immunosuppression, and also contains samples before and after infliximab treatment. We likewise preprocessed these datasets using the same framework as described for the CMap data.
Methotrexate response data in rheumatoid arthritis. Aside from IBD, we processed a rheumatoid arthritis dataset, GSE45867 [18], which contains human synovial biopsy samples before and after methotrexate treatment (N = 8) and before and after tocilizumab treatment (N = 12). As tocilizumab is not a drug present CMap, we analyzed only the methotrexate samples, comparing our simulated drugged samples generated from the untreated samples to the measured treated samples.

Network construction
We next downloaded all drug data available online from DrugBank [19]. We extracted drug targets and converted all gene names to the Entrez gene identifier. Only drugs that were included in CMap and have a gene target that was measured were included in downstream analysis.
We extracted all human pathways from Reactome [20] and connected all pathways into one large network. Only genes that were measured in CMap were included. Undirected edges from Reactome are represented as two directed edges in opposite directions.
As shown in Fig 1, step 1, each circular node represents a gene, and each directed arrow represents an interaction from Reactome. We next add drug nodes, depicted as triangles, to the overall network or genes (Fig 1, step 1). A drug may have multiple targets, which would be represented as having multiple outgoing edges. Similarly, a gene may be the target of multiple drugs, depicted by multiple incoming edges.

Drug effect modeling
The network is then parameterized by drugged gene expression data, where each circular gene node is a linear regression of the incoming nodes, and triangular drug nodes are binary PLOS COMPUTATIONAL BIOLOGY variables (Fig 1, step 2). Regularized regression was not used as the median number of input genes for each node, 13, was small compared to the 3,400 training instances, lowering the concern for overfitting. The network of regressions can be intuitively thought of as diffusion of the drug effect through the gene network. Thus, we have now created a model where drugs can be turned "on" or "off", where the network is capturing the changes induced by these effects.
With the parameterized network modeling drug effects, we then apply this network to IBD disease data in order to discover which drugs may be most effective for each individual patient. For each IBD sample, we take the gene expression values and overlay them onto the nodes of the network. We then turn drugs on or off (Fig 1, step 3), and propagate these effects through the network (Fig 1 step 4). Thus, for each IBD sample, we create a "drugged IBD sample" gene expression sample. Using the healthy data, we average all healthy samples, to create an "average healthy" gene expression sample (Fig 1, step 5). We then compare the average healthy sample to each drugged IBD sample, using Euclidean distance, to create a personalized, ranked drug list (Fig 1, step 6). In essence, we are comparing the effects of different drugs to see which set can bring the original IBD sample closest to the average healthy sample. This way, we discover drugs that not only treat the disease symptoms, but may have fewer side effects than drugs that may have a beneficial effect on the disease but do not result in an expression profile closer to the average healthy sample.

Curation of IBD drugs in literature
We curated a list of drugs that have been previously studied in the context IBD based on work by Percha, et al. [21,22] Briefly, entities and dependency paths are extracted from over 16 million MEDLINE abstracts. These paths are then clustered using the Ensemble Biclustering for Classification method to produce clusters of entity relationships. Specifically, we used the chemical-disease relationship output, which consists of 6 themes: 1) (T) Treatment/therapy, 2) (C) Inhibits cell growth, 3) (Sa) side effect/adverse event, 4) (Pr) prevents, suppresses, 5) (Pa) alleviates, reduces, 6) (J) role in pathogenesis. For our known drugs, we extracted any drug that had a T, Pr, and/or Pa relationship with CD, UC, and/or IBD.

Calculating alternative drug-disease score rankings
We compared NetPTP to the drug rankings produced by the method used in Dudley et al [6], as our work uses similar data sources and is also applied to the same disease, IBD. In their method, a drug-disease score (DDS) is derived for each drug based on how anti-correlated the drug signature is as compared to the disease signature. We calculated the DDS as described in Sirota and Dudley et al [23] for all drugs and CD, and all drugs and UC using the samples in GEO9686, GEO36807, and GEO10616. GEO16879 was not included as these patients were refractory to the mainstay treatments for IBD. Following their method, we first used significance analysis of microarrays (SAM) [24] to derive lists of up-regulated and down-regulated genes, comparing diseased samples to the healthy samples in each study. We then calculated the up-regulated and down-regulated enrichment score, and subsequently the drug-disease score (DDS) for each drug-CD and drug-UC pair in each of the three studies.
We then ranked all drugs by their DDS, from most anti-correlated to least anti-correlated. The original method calls for DDS to be set to 0 if the up-regulated enrichment score and down-regulated enrichment score are the same direction, otherwise DDS is set to the difference between the two. This results in a large number of drugs all having a score of 0, and we differentiated drugs with a score of 0 by ranking them by the difference between the up-regulated and down-regulated enrichment scores.
In order to compare the rankings from the two methods, we curated the FDA approved treatments for IBD that are in CMap. These were budesonide, prednisone, prednisolone, methylprednisolone, azathioprine, mercaptopurine, sulfasalazine, mesalazine, and methotrexate. We then compared the rankings of these nine known treatments based on the DDS to our NetPTP rankings across all studies.

Experimental protocol for drug evaluation in TNBS mice
We conducted a pilot experimental study to assess the top drug prediction in TNBS mice, amrinone. The study was reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) at Stanford University. We purchased twenty 6-week-old C57BL/6 mice from The Jackson Laboratory (Bar Harbor, ME). This species was chosen to match the species used from the source data used for drug prediction from Dohi et al [17], though this species is known to be more resistant to the development of TNBS colitis as compared to others [25,26]. The mice were divided evenly into drug treatment and control groups. TNBS 5% w/v in methanol was purchased from Fisher Scientific Company, LLC (Hampton, NH). TNBS at a dosage of 100mg/ kg in 50% ethanol in a volume of 150μL was administered intrarectally to each mouse on day 0. For control mice, 150μL of 50% ethanol solution was used. Amrinone was purchased from Sigma-Aldrich, Inc (St. Louis, MO) and was administered at 10mg/kg in 1% v/v DMSO with saline at a dilution of 1mg/mL via intraperitoneal injection. Control mice received a 10mg/kg injection of saline. The treatment injections were administered daily from day 1 to day 4.
Mice were sacrificed on day 4, which was chosen to further match the experiments from the source data in Dohi et al [17], where mice were sacrificed on day 4. The colon was then dissected and harvested from each mouse. The tissue was fixed in formalin and paraffin-embedded. A longitudinal cross section slice of the colon was mounted on a slide and stained with hematoxylin and eosin.
Each slide was assessed in a blinded and randomized fashion by an independent veterinary pathologist. Each slide was graded based on ulceration, inflammation, edema, and fibroplasia. Ulceration, edema, and fibroplasia were considered as present or absent. Inflammation was graded as follows: 0 = no inflammation; 1 = small, focal areas limited to the lamina propria; 2 = multifocal or coalescing areas extending into the submucosa; 3 = transmural.
Given the heterogeneous and patchy nature of TNBS colitis, we assessed for the presence of fibroplasia, or wound healing, while taking into account the degree of induced colitis. Using R 3.6.1, we constructed a logistic regression model to assess the presence of fibroplasia. We used the degree of inflammation, the presence of ulceration, the presence of edema, and the treatment group as the covariates in our model.

Integrating the connectivity map, drugbank, and reactome
Our network consists of 6,982 genes, that corresponds to the intersection of genes measured in CMap and genes present in Reactome. From DrugBank and CMap, we curated 453 drugs, which corresponded to 3,400 instances from CMap, including controls. These 453 drugs targeted 496 genes, with a median of 4 gene targets per drug.

Similar drug mechanisms cluster together
We first assessed the effects of the drugs on individual patients, to examine if similar drugs result in similar effects. As an example, using a drugged IBD sample derived from the first CD sample (GSM244753) in GSE9686, we clustered the drugged profiles (S1 Fig). The dendrogram reveals that drugs with similar mechanisms of action cluster together. Starting from the right, we see that the steroids cluster together, which includes topical and systemic steroids (Fig 2A). The neighborhood of drugs clustered near the steroids include additional anti-inflammatory drugs and immunosuppressants. For example, the cluster right below the steroids includes leflunomide, chloroquine, sirolimus, and tacrolimus (Fig 2A). Amrinone is present in this cluster, with chloroquine as its nearest neighbor. Though the steroids belong to four different categories of ATC codes due to being used in different disease processes, they all have a similar mechanism of action.
As we move counterclockwise, we come across a group of antibiotics mixed with chemotherapeutic agents (Fig 2B). These drugs block various forms of topoisomerase, with the antibiotics blocking bacterial topoisomerase and the chemotherapeutic agents blocking human topoisomerase.
Continuing along, the next large cluster along the top contains drugs that act on various receptors within the body, such as beta-adrenergic and dopamine receptors (Fig 2C). On the left, we see another section of anti-inflammatory drugs, including cyclooxygenase inhibitors and immunomodulatory drugs such as azathioprine (Fig 2D). The bottom portion of the circle contains drugs that are used for neuropsychiatric diseases, liver, and kidney issues, such as antipsychotics, diabetes and cholesterol medications, and diuretics (S1 Fig).

Assessing drug response in GSE16879
Using the samples from GSE16879, we visualized treatment responders, treatment nonresponders, and controls for CD ( Fig 3A) and UC patients (Fig 3B). Before treatment,

PLOS COMPUTATIONAL BIOLOGY
responders and nonresponders tend to be located in close proximity in both CD and UC. For treatment nonresponders, there was minimal movement of the samples after treatment. For treatment responders, there is a shift of the after-response samples toward the cluster of healthy control samples. Overall, treatment responders after treatment were located significantly closer to healthy controls than treatment nonresponders for CD (p < 0.001) and for UC (p < 0.001). This indicates that the treatment appears to have an effect on the tissue that brings the sample closer to the healthy tissue, in part reversing some of the effects caused by the disease. The treated samples ultimately appear to exhibit an expression profile partway between the original disease state and healthy state. Our approach leverages this when creating our ranked drug list, comparing simulated drugged samples to healthy control samples.

Comparing measured versus predicted methotrexate response
We visualized our simulated methotrexate treatment samples versus untreated and actual methotrexate treatment samples using principal components analysis (Fig 4). Notably, untreated samples and methotrexate treated samples did not cluster together but were scattered across both principal component 1 (PC1) and principal component 2 (PC2) (Fig 4A). PC1 and PC2 accounted for 29.1% and 13.2% of the variance respectively. The simulated methotrexate samples were generally located on the trajectory between the untreated and treated samples, with some deviance towards the center of the plot. When analyzing PC1 and PC2 separately (Fig 4B and 4C), the simulated methotrexate samples always fell in between the untreated and actual treated samples for PC2 (Fig 4C).
GSE16879 consists of patients that are refractory to corticosteroids and/or immunosuppression, with 7 UC patients and 6 CD patients on corticosteroids at baseline, per Table 1 of Arijs et al [12]. We evaluated the rank of prednisone in GSE16879 versus the other studies. For patients in GSE16879, the rank of prednisone was significantly higher than the patients from the other 3 studies (p = 0.028), indicating prednisone was predicted to be less effective in patients in GSE16879. The median in GSE16879 was 261 with an interquartile range (IQR) of 154.5-338.5 versus a median of 178 and an IQR of 68. 25-297.25 in the other patients.
We next assessed the rankings of drugs previously associated with IBD in the literature for CD (Fig 5A) and UC samples (Fig 5B). In CD patients (Fig 5A), three different subgroups of patients appear, and samples do not cluster by study. Patients on the left are predicted to have good response to anti-inflammatory medication and steroids, but not to immunomodulators such as azathioprine. Those in the middle show the opposite pattern, with good predicted response to immunomodulators but not anti-inflammatory medications, and those on the right tend to have high ranks for all classes.
In UC patients (Fig 5B), we again see that samples do not cluster by study, but instead in two main groups. Patients on the left exhibit a high ranking for steroids, anti-inflammatory medications, and immunomodulators. In contrast, patients on the right are predicted to have poor efficacy to immunomodulators but maintain favorable rankings for steroids and antiinflammatory medications.

Comparison to drug-disease score rankings
We compared the ranking of nine known IBD treatments generated by NetPTP against the rankings generated by the DDS from the method in Sirota and Dudley et al [23]. Rankings were generated for each drug-disease pair for 9 drugs, 2 diseases (CD, UC), and in 3 studies, resulting in N = 54. Using the Wilcoxon rank-sum paired test, our NetPTP rankings of the known IBD treatments were significantly more favorable than the DDS rankings (p = 0.0011). Overall, where a lower ranking is more favorable, our rankings of known IBD treatments had a median ranking of 102.5 (interquartile range 46.5-205.2) and the DDS rankings had a median of 243.5 (interquartile range 128.2-331.8). Fig 6 depicts the rankings of previous IBD drugs in literature for the DSS day 2, day 4, day 6, and TNBS samples from GSE22307 [16] and GSE53835 [17]. As we progress from day 2 to day 6 of DSS administration and the inflammation is being induced, the drug ranking pattern changes from immunomodulators being highly ranked to ant-inflammatory drugs and steroids being highly ranked. For example, sulfasalazine significantly improves in rank as the inflammation progresses (Fig 7, p = 0.01). By the time we reach day 6, the pattern generally resembles the TNBS samples, which represent acute inflammation of the colon. For day 6 DSS mice and TNBS mice, the top 10 drugs are shown in Table 4.

Experimental evaluation of top ranked DRUG in TNBS mice
We constructed a logistic regression model to assess the presence of wound healing in TNBS mice given our top prediction, amrinone, versus saline controls (Table 5). Drug status was treated as 1 for amrinone, and 0 for saline. A higher level of inflammation or the presence of ulceration led to a greater likelihood of the presence of recovery, whereas edema had the opposite effect. Those mice that were in the amrinone treatment group had a coefficient of 2.06, corresponding to an odds ratio of 7.75.

Discussion
In this work, we present NetPTP, a novel systems pharmacological approach for drug repurposing and modeling drug effects. In our approach, we model drug effects and apply these effects to diseased samples to assess which samples exhibit the closest return to a healthy profile, thus leading to individualized predictions of drug efficacy. NetPTP was motivated in part by previous studies [6] and our analysis of the samples in GSE16879, which showed that infliximab responder samples were located significantly closer to healthy controls than nonresponder samples. By assessing for a shift toward healthy controls, we believe we may not only predict the efficacy of the drug in treating the disease but may also rank drugs higher that have fewer adverse effects.
We assessed NetPTP using methotrexate treated samples from GSE45867, which studied another autoimmune disease, rheumatoid arthritis. Our simulated drugged samples were generally located on the trajectory between the untreated and the treated samples. NetPTP appears to capture the direction of the effect of the drug and then subsequently apply the learned effect to new data samples, producing simulated drugged samples which represent a conservative step from the untreated sample towards the treated sample. In particular, the model's prediction fell between the untreated and treated sample for all eight samples along principal component 2. Thus, this particular principal component may reflect more of the biological changes induced with methotrexate treatment. Delving deeper, the gene with the largest absolute weight in PC2 is PMAIP1, also known as NOXA. NOXA plays a role in the apoptosis pathway as part of the BCL-2 family and is regulated by p53 [27]. Methotrexate has been shown to mediate apoptosis via upregulation of p53 and its downstream targets, including NOXA [28]. NOXA has been shown to induce apoptosis of fibroblast-like synoviocytes [29] and bone osteoclasts [30], both of which are thought to have a role in the pathogenesis of rheumatoid arthritis [31].
We applied NetPTP to multiple publicly available CD and UC datasets with human colonic samples. Clustering the simulated drugged profiles for one patient revealed that drugs with similar mechanisms cluster together (Fig 2). When clustering our patients and their rankings with IBD drugs in literature, we found that patients grouped into multiple different treatment response profiles (Fig 5). This suggests that it may be possible to guide treatment decision making, particularly when deciding between multiple first line treatment options. In CD samples, we see that of the ten top ranked drugs, mesalazine, a known IBD drug, is included ( Table 3). The others are diuretics, antimicrobials, and two GABA antagonists, including thiocolchicoside which has anti-inflammatory effects [32]. This list suggests that these CD samples tend to have a fluid overloaded state, likely due to edema, with some dysregulation of the enteric nervous system and inflammation. In UC, though we have some of the same diuretic medications, we also see steroids more often, for example as the top ranked drug in GSE36807 ( Table 2, Table 3). The top overall ranked drug, diazoxide, is a potassium channel activator, and has been shown to heal acute gastric ulcers in rats [33]. Carbonic anhydrase inhibitors are also present on both top ranked lists, and previous studies have shown targeting carbonic anhydrase I and IV have ameliorated IBD in mouse models [34,35].
We compared the rankings generated by NetPTP to the drug-disease score (DDS) generated by Sirota and Dudley et al [23]. Overall, our method produced significantly better rankings for the nine known IBD drugs in the Connectivity Map. However, as IBD is a heterogenous disease, some of the study patients may have not been responsive to each and every one of these therapies. We compared our results at the study level, as calculating the DDS uses SAM, which aggregates multiple patients to derive up-regulated and down-regulated genes. NetPTP aims to address some of these limitations by providing predictions at the patient level while incorporating the connections between genes to capture more of the underlying biology driving drug response.
In mouse samples, we analyzed the drug rankings over the course of acute inflammation development in DSS mice and in TNBS mice. We found that as inflammation progressed from day 2 to day 6, known treatments such as sulfasalazine became significantly more highly ranked (Fig 7), and the day 6 DSS drug rankings became more similar to the drug rankings seen with the TNBS mice (Fig 6). The top drug ranked in the TNBS mice was amrinone, a phosphodiesterase (PDE) type 3 inhibitor with some effect on PDE type 4, which was clustered with other immunosuppressants (Fig 2A). Phosphodiesterases control the concentration of  We assessed this top drug, amrinone, in a preliminary study using a TNBS mouse model, using the same mouse strain as used in GSE53835. We assessed the presence of fibroplasia in twenty mice, taking into account the degree of inflammation and the presence of ulceration, as this particular strain is moderately resistant to TNBS and to account for the heterogenous nature of the TNBS model [25,26]. We evaluated the drug effect on day 4 in keeping with the source data; however, these experiments may benefit from a longer time course to better characterize the drug effect. Though the drug status did not reach statistical significance, it appears to be one of the more important covariates for predicting the presence of fibroplasia in these mice and would merit further investigation in a larger study.
Our approach involves curating and combining multiple publicly available resources, including Reactome [20], the Connectivity Map [39], DrugBank [19], and the Gene Expression Omnibus. Though CMap includes hundreds of drugs, some IBD therapies, such as infliximab, are not represented. Furthermore, transcriptomics data are shifting from gene expression to RNA-seq. For IBD, currently most published studies for public use are gene expression data. However, our approach can easily translate to RNA-seq data and be expanded to include more drugs, such as data from the Library of Integrated Network-Based Cellular Signatures (LINCS) project (http://www.lincsproject.org/).
In addition to different sources of drugged data, the other aspects of NetPTP are also modular and can be easily adjusted to accommodate different network architectures, network data from different species, additional drugs and drug targets, and additional diseases of interest. As NetPTP can use different sources of healthy control data, it could also be used to rank drugs based on a patient's own healthy samples before he or she developed a disease. In addition to being able to incorporate different sources of transcriptomic data, NetPTP works with any pre-defined set of edges that has been converted to a directed network, such as Reactome, , or networks derived from experimental data. With time, we hope the increasing amount of publicly available data available and a wider selection of drugged data including newer therapies such as monoclonal antibodies will expand the applicability and utility of our method, offering personalized treatment regimens as well as identifying novel treatment avenues for IBD.
Supporting information S1 Fig. Dendrogram of drugged profiles of a CD patient from GSE9686. Dendrogram branches are colored by the first level of the anatomic therapeutic chemical classification system.
(TIF) S1 Table. Dataset summary for publicly available mouse and human data. (DOCX)