Application of systems biology-based in silico tools to optimize treatment strategy identification in Still's disease

The emerging paradigm in pediatric rheumatology: harnessing the power of artificial intelligence

Artificial intelligence algorithms, with roots extending into the past but experiencing a resurgence and evolution in recent years due to their superiority over traditional methods and contributions to human capabilities, have begun to make their presence felt in the field of pediatric rheumatology. In the ever-evolving realm of pediatric rheumatology, there have been incremental advancements supported by artificial intelligence in understanding and stratifying diseases, developing biomarkers, refining visual analyses, and facilitating individualized treatment approaches. However, like in many other domains, these strides have yet to gain clinical applicability and validation, and ethical issues remain unresolved. Furthermore, mastering different and novel terminologies appears challenging for clinicians. This review aims to provide a comprehensive overview of the current literature, categorizing algorithms and their applications, thus offering a fresh perspective on the nascent relationship between pediatric rheumatology and artificial intelligence, highlighting both its advancements and constraints.

Aflibercept Off-Target Effects in Diabetic Macular Edema: An In Silico Modeling Approach

Intravitreal aflibercept injection (IAI) is a treatment for diabetic macular edema (DME), but its mechanism of action (MoA) has not been completely elucidated. Here, we aimed to explore IAI's MoA and its multi-target nature in DME pathophysiology with an in silico (computer simulation) disease model. We used the Therapeutic Performance Mapping System (Anaxomics Biotech property) to generate mathematical models based on the available scientific knowledge at the time of the study, describing the relationship between the modulation of vascular endothelial growth factor receptors (VEGFRs) by IAI and DME pathophysiological processes. We also undertook an enrichment analysis to explore the processes modulated by IAI, visualized the effectors' predicted protein activity, and specifically evaluated the role of VEGFR1 pathway inhibition on DME treatment. The models simulated the potential pathophysiology of DME and the likely IAI's MoA by inhibiting VEGFR1 and VEGFR2 signaling. The action of IAI through both signaling pathways modulated the identified pathophysiological processes associated with DME, with the strongest effects in angiogenesis, blood-retinal barrier alteration and permeability, and inflammation. VEGFR1 inhibition was essential to modulate inflammatory protein effectors. Given the role of VEGFR1 signaling on the modulation of inflammatory-related pathways, IAI may offer therapeutic advantages for DME through sustained VEGFR1 pathway inhibition.

Expert consensus on the treatment of patients with adult-onset still's disease with the goal of achieving an early and long-term remission

We performed a comprehensive systematic targeted literature review and used the Delphi method to formulate expert consensus statements to guide the treatment of adult-onset Still's disease (AOSD) to achieve an early and long-term remission. Seven candidate statements were generated and reached consensus in the first round of voting by the panel of experts. We postulate: (i) In patients with AOSD with predominant arthritis at onset who achieved no disease control with glucocorticoids (GCs), the use of methotrexate can be considered, whereas the use of cyclosporin A and low-dose GCs should not (Statements 1-3); (ii) In patients with AOSD with poor prognostic factors at diagnosis, an IL-1 inhibitor (IL-1i) in addition to GCs should be taken into consideration as early as possible (Statement 4); (iii) A switch to an IL-6 inhibitor (IL-6i) may be considered in patients with AOSD with prevalent joint involvement, who are unresponsive or intolerant to IL-1i (Statement 5); (iv) Drug tapering or discontinuation may be considered in patients who achieved a sustained clinical and laboratory remission with IL-1i (Statement 6); (v) In patients with AOSD who failed to attain a good clinical response with an IL-1i, switching to an IL-6i may be considered in alternative to a different IL-1i. TNF-inhibitors may be considered as a further choice in patients with a prominent joint involvement (Statement 7). These statements will help clinicians in treatment decision making in patients with AOSD.

A quantitative systems pharmacology model for certolizumab pegol treatment in moderate-to-severe psoriasis

Background

Psoriasis is a chronic immune-mediated inflammatory systemic disease with skin manifestations characterized by erythematous, scaly, itchy and/or painful plaques resulting from hyperproliferation of keratinocytes. Certolizumab pegol [CZP], a PEGylated antigen binding fragment of a humanized monoclonal antibody against TNF-alpha, is approved for the treatment of moderate-to-severe plaque psoriasis. Patients with psoriasis present clinical and molecular variability, affecting response to treatment. Herein, we utilized an in silico approach to model the effects of CZP in a virtual population (vPop) with moderate-to-severe psoriasis. Our proof-of-concept study aims to assess the performance of our model in generating a vPop and defining CZP response variability based on patient profiles.

Methods

We built a quantitative systems pharmacology (QSP) model of a clinical trial-like vPop with moderate-to-severe psoriasis treated with two dosing schemes of CZP (200 mg and 400 mg, both every two weeks for 16 weeks, starting with a loading dose of CZP 400 mg at weeks 0, 2, and 4). We applied different modelling approaches: (i) an algorithm to generate vPop according to reference population values and comorbidity frequencies in real-world populations; (ii) physiologically based pharmacokinetic (PBPK) models of CZP dosing schemes in each virtual patient; and (iii) systems biology-based models of the mechanism of action (MoA) of the drug.

Results

The combination of our different modelling approaches yielded a vPop distribution and a PBPK model that aligned with existing literature. Our systems biology and QSP models reproduced known biological and clinical activity, presenting outcomes correlating with clinical efficacy measures. We identified distinct clusters of virtual patients based on their psoriasis-related protein predicted activity when treated with CZP, which could help unravel differences in drug efficacy in diverse subpopulations. Moreover, our models revealed clusters of MoA solutions irrespective of the dosing regimen employed.

Conclusion

Our study provided patient specific QSP models that reproduced clinical and molecular efficacy features, supporting the use of computational methods as modelling strategy to explore drug response variability. This might shed light on the differences in drug efficacy in diverse subpopulations, especially useful in complex diseases such as psoriasis, through the generation of mechanistically based hypotheses.

In silico evaluation of the role of lisdexamfetamine on attention-deficit/hyperactivity disorder common psychiatric comorbidities: mechanistic insights on binge eating disorder and depression

Attention-deficit/hyperactivity disorder (ADHD) is a psychiatric condition well recognized in the pediatric population that can persist into adulthood. The vast majority of patients with ADHD present psychiatric comorbidities that have been suggested to share, to some extent, the pathophysiological mechanism of ADHD. Lisdexamfetamine (LDX) is a stimulant prodrug approved for treating ADHD and, in the US, also for binge eating disorder (BED). Herein, we evaluated, through a systems biology-based in silico method, the efficacy of a virtual model of LDX (vLDX) as ADHD treatment to improve five common ADHD psychiatric comorbidities in adults and children, and we explored the molecular mechanisms behind LDX's predicted efficacy. After the molecular characterization of vLDX and the comorbidities (anxiety, BED, bipolar disorder, depression, and tics disorder), we created a protein-protein interaction human network to which we applied artificial neural networks (ANN) algorithms. We also generated virtual populations of adults and children-adolescents totaling 2,600 individuals and obtained the predicted protein activity from Therapeutic Performance Mapping System models. The latter showed that ADHD molecular description shared 53% of its protein effectors with at least one studied psychiatric comorbidity. According to the ANN analysis, proteins targeted by vLDX are predicted to have a high probability of being related to BED and depression. In BED, vLDX was modeled to act upon neurotransmission and neuroplasticity regulators, and, in depression, vLDX regulated the hypothalamic-pituitary-adrenal axis, neuroinflammation, oxidative stress, and glutamatergic excitotoxicity. In conclusion, our modeling results, despite their limitations and although requiring in vitro or in vivo validation, could supplement the design of preclinical and potentially clinical studies that investigate treatment for patients with ADHD with psychiatric comorbidities, especially from a molecular point of view.

The circadian rhythm: an influential soundtrack in the diabetes story

Type 2 Diabetes Mellitus (T2DM) has been the main category of metabolic diseases in recent years due to changes in lifestyle and environmental conditions such as diet and physical activity. On the other hand, the circadian rhythm is one of the most significant biological pathways in humans and other mammals, which is affected by light, sleep, and human activity. However, this cycle is controlled via complicated cellular pathways with feedback loops. It is widely known that changes in the circadian rhythm can alter some metabolic pathways of body cells and could affect the treatment process, particularly for metabolic diseases like T2DM. The aim of this study is to explore the importance of the circadian rhythm in the occurrence of T2DM via reviewing the metabolic pathways involved, their relationship with the circadian rhythm from two perspectives, lifestyle and molecular pathways, and their effect on T2DM pathophysiology. These impacts have been demonstrated in a variety of studies and led to the development of approaches such as time-restricted feeding, chronotherapy (time-specific therapies), and circadian molecule stabilizers.

Application of systems biology to identify pharmacological mechanisms of thrombotic microangiopathy evoked by combined activated prothrombin complex concentrate and emicizumab

Emicizumab is a bispecific monoclonal antibody that substitutes for the function of missing or deficient factor VIII (FVIII) in people with hemophilia A (PwHA). Long-term safety and efficacy of emicizumab have been demonstrated in several clinical trials. Nevertheless, in the first of these, three cases of thrombotic microangiopathy (TMA) occurred in PwHA treated with emicizumab receiving high doses of activated prothrombin complex concentrate (aPCC), a bypassing agent used for treating breakthrough bleeds when FVIII neutralizing antibodies (inhibitors) make FVIII replacement ineffective. The aim of the present work is to offer a method to elucidate the pathophysiological and pharmacological mechanisms involved in this treatment-induced TMA. Systems biology and machine learning-based Therapeutic Performance Mapping System is a validated in silico technology that allowed us to construct models of potential mechanisms behind induced TMA. Two drug combinations were modeled and assessed: emicizumab plus aPCC and emicizumab plus recombinant activated factor VII (another bypassing agent). Our models showed that both combinations were related to activation of the coagulation cascade. However, mechanisms involved mainly in platelet activation and possibly in complement activation were detected only for emicizumab plus aPCC, potentially explaining the occurrence of TMA only in this combination.

In silico clinical trial evaluating lisdexamfetamine's and methylphenidate's mechanism of action computational models in an attention-deficit/hyperactivity disorder virtual patients' population

Introduction

Attention-deficit/hyperactivity disorder (ADHD) is an impairing psychiatric condition with the stimulants, lisdexamfetamine (LDX), and methylphenidate (MPH), as the first lines pharmacological treatment.

Methods

Herein, we applied a novel in silico method to evaluate virtual LDX (vLDX) and vMPH as treatments for ADHD applying quantitative systems pharmacology (QSP) models. The objectives were to evaluate the model's output, considering the model characteristics and the information used to build them, to compare both virtual drugs' efficacy mechanisms, and to assess how demographic (age, body mass index, and sex) and clinical characteristics may affect vLDX's and vMPH's relative efficacies.

Results and Discussion

We molecularly characterized the drugs and pathologies based on a bibliographic search, and generated virtual populations of adults and children-adolescents totaling 2,600 individuals. For each virtual patient and virtual drug, we created physiologically based pharmacokinetic and QSP models applying the systems biology-based Therapeutic Performance Mapping System technology. The resulting models' predicted protein activity indicated that both virtual drugs modulated ADHD through similar mechanisms, albeit with some differences. vMPH induced several general synaptic, neurotransmitter, and nerve impulse-related processes, whereas vLDX seemed to modulate neural processes more specific to ADHD, such as GABAergic inhibitory synapses and regulation of the reward system. While both drugs' models were linked to an effect over neuroinflammation and altered neural viability, vLDX had a significant impact on neurotransmitter imbalance and vMPH on circadian system deregulation. Among demographic characteristics, age and body mass index affected the efficacy of both virtual treatments, although the effect was more marked for vLDX. Regarding comorbidities, only depression negatively impacted both virtual drugs' efficacy mechanisms and, while that of vLDX were more affected by the co-treatment of tic disorders, the efficacy mechanisms of vMPH were disturbed by wide-spectrum psychiatric drugs. Our in silico results suggested that both drugs could have similar efficacy mechanisms as ADHD treatment in adult and pediatric populations and allowed raising hypotheses for their differential impact in specific patient groups, although these results require prospective validation for clinical translatability.

Artificial Intelligence and laboratory data in rheumatic diseases

Artificial intelligence (AI)-based medical technologies are rapidly evolving into actionable solutions for clinical practice. Machine learning (ML) algorithms can process increasing amounts of laboratory data such as gene expression immunophenotyping data and biomarkers. In recent years, the analysis of ML has become particularly useful for the study of complex chronic diseases, such as rheumatic diseases, heterogenous conditions with multiple triggers. Numerous studies have used ML to classify patients and improve diagnosis, to stratify the risk and determine disease subtypes, as well as to discover biomarkers and gene signatures. This review aims to provide examples of ML models for specific rheumatic diseases using laboratory data and some insights into relevant strengths and limitations. A better understanding and future application of these analytical strategies could facilitate the development of precision medicine for rheumatic patients.

Artificial intelligence assessment of the potential of tocilizumab along with corticosteroids therapy for the management of COVID-19 evoked acute respiratory distress syndrome

Acute respiratory distress syndrome (ARDS), associated with high mortality rate, affects up to 67% of hospitalized COVID-19 patients. Early evidence indicated that the pathogenesis of COVID-19 evoked ARDS is, at least partially, mediated by hyperinflammatory cytokine storm in which interleukin 6 (IL-6) plays an essential role. The corticosteroid dexamethasone is an effective treatment for severe COVID-19 related ARDS. However, trials of other immunomodulatory therapies, including anti-IL6 agents such as tocilizumab and sarilumab, have shown limited evidence of benefit as monotherapy. But recently published large trials have reported added benefit of tocilizumab in combination with dexamethasone in severe COVID-19 related ARDS. In silico tools can be useful to shed light on the mechanisms evoked by SARS-CoV-2 infection and of the potential therapeutic approaches. Therapeutic performance mapping system (TPMS), based on systems biology and artificial intelligence, integrate available biological, pharmacological and medical knowledge to create mathematical models of the disease. This technology was used to identify the pharmacological mechanism of dexamethasone, with or without tocilizumab, in the management of COVID-19 evoked ARDS. The results showed that while dexamethasone would be addressing a wider range of pathological processes with low intensity, tocilizumab might provide a more direct and intense effect upon the cytokine storm. Based on this in silico study, we conclude that the use of tocilizumab alongside dexamethasone is predicted to induce a synergistic effect in dampening inflammation and subsequent pathological processes, supporting the beneficial effect of the combined therapy in critically ill patients. Future research will allow identifying the ideal subpopulation of patients that would benefit better from this combined treatment.

Unveiling chronic spontaneous urticaria pathophysiology through systems biology

BACKGROUND

Chronic spontaneous urticaria (CSU) is a rare, heterogeneous, severely debilitating, and often poorly controlled skin disease resulting in an itchy eruption that can be persistent. Antihistamines and omalizumab, an anti-IgE mAb, are the only licensed therapies. Although CSU pathogenesis is not yet fully understood, mast cell activation through the IgE:high-affinity IgE receptor (FcεRI) axis appears central to the disease process.

OBJECTIVE

We sought to model CSU pathophysiology and identify in silico the mechanism of action of different CSU therapeutic strategies currently in use or under development.

METHODS

Therapeutic performance mapping system technology, based on systems biology and machine learning, was used to create a CSU interactome validated with gene expression data from patients with CSU and a CSU model that was used to evaluate CSU pathophysiology and the mechanism of action of different therapeutic strategies.

RESULTS

Our models reflect the known role of mast cell activation as a central process of CSU pathophysiology, as well as recognized roles for different therapeutic strategies in this and other innate and adaptive immune processes. They also allow determining similarities and differences between them; anti-IgE and Bruton tyrosine kinase inhibitors play a more direct role in mast cell biology through abrogation of FcεRI signaling activity, whereas anti-interleukins and anti-Siglec-8 have a role in adaptive immunity modulation.

CONCLUSION

In silico CSU models reproduced known CSU and therapeutic strategies features. Our results could help advance understanding of therapeutic mechanisms of action and further advance treatment research by patient profile.

Systems biology and artificial intelligence analysis highlights the pleiotropic effect of IVIg therapy in autoimmune diseases with a predominant role on B cells and complement system

Intravenous immunoglobulin (IVIg) is used as treatment for several autoimmune and inflammatory conditions, but its specific mechanisms are not fully understood. Herein, we aimed to evaluate, using systems biology and artificial intelligence techniques, the differences in the pathophysiological pathways of autoimmune and inflammatory conditions that show diverse responses to IVIg treatment. We also intended to determine the targets of IVIg involved in the best treatment response of the evaluated diseases. Our selection and classification of diseases was based on a previously published systematic review, and we performed the disease characterization through manual curation of the literature. Furthermore, we undertook the mechanistic evaluation with artificial neural networks and pathway enrichment analyses. A set of 26 diseases was selected, classified, and compared. Our results indicated that diseases clearly benefiting from IVIg treatment were mainly characterized by deregulated processes in B cells and the complement system. Indeed, our results show that proteins related to B-cell and complement system pathways, which are targeted by IVIg, are involved in the clinical response. In addition, targets related to other immune processes may also play an important role in the IVIg response, supporting its wide range of actions through several mechanisms. Although B-cell responses and complement system have a key role in diseases benefiting from IVIg, protein targets involved in such processes are not necessarily the same in those diseases. Therefore, IVIg appeared to have a pleiotropic effect that may involve the collaborative participation of several proteins. This broad spectrum of targets and 'non-specificity' of IVIg could be key to its efficacy in very different diseases.

Identification of the Potential Molecular Mechanisms Linking RUNX1 Activity with Nonalcoholic Fatty Liver Disease, by Means of Systems Biology

Nonalcoholic fatty liver disease (NAFLD) is the most prevalent chronic hepatic disease; nevertheless, no definitive diagnostic method exists yet, apart from invasive liver biopsy, and nor is there a specific approved treatment. Runt-related transcription factor 1 (RUNX1) plays a major role in angiogenesis and inflammation; however, its link with NAFLD is unclear as controversial results have been reported. Thus, the objective of this work was to determine the proteins involved in the molecular mechanisms between RUNX1 and NAFLD, by means of systems biology. First, a mathematical model that simulates NAFLD pathophysiology was generated by analyzing Anaxomics databases and reviewing available scientific literature. Artificial neural networks established NAFLD pathophysiological processes functionally related to RUNX1: hepatic insulin resistance, lipotoxicity, and hepatic injury-liver fibrosis. Our study indicated that RUNX1 might have a high relationship with hepatic injury-liver fibrosis, and a medium relationship with lipotoxicity and insulin resistance motives. Additionally, we found five RUNX1-regulated proteins with a direct involvement in NAFLD motives, which were NFκB1, NFκB2, TNF, ADIPOQ, and IL-6. In conclusion, we suggested a relationship between RUNX1 and NAFLD since RUNX1 seems to regulate NAFLD molecular pathways, posing it as a potential therapeutic target of NAFLD, although more studies in this field are needed.

Identification of Potential Muscle Biomarkers in McArdle Disease: Insights from Muscle Proteome Analysis

Glycogen storage disease type V (GSDV, McArdle disease) is a rare genetic myopathy caused by deficiency of the muscle isoform of glycogen phosphorylase (PYGM). This results in a block in the use of muscle glycogen as an energetic substrate, with subsequent exercise intolerance. The pathobiology of GSDV is still not fully understood, especially with regard to some features such as persistent muscle damage (i.e., even without prior exercise). We aimed at identifying potential muscle protein biomarkers of GSDV by analyzing the muscle proteome and the molecular networks associated with muscle dysfunction in these patients. Muscle biopsies from eight patients and eight healthy controls showing none of the features of McArdle disease, such as frequent contractures and persistent muscle damage, were studied by quantitative protein expression using isobaric tags for relative and absolute quantitation (iTRAQ) followed by artificial neuronal networks (ANNs) and topology analysis. Protein candidate validation was performed by Western blot. Several proteins predominantly involved in the process of muscle contraction and/or calcium homeostasis, such as myosin, sarcoplasmic/endoplasmic reticulum calcium ATPase 1, tropomyosin alpha-1 chain, troponin isoforms, and alpha-actinin-3, showed significantly lower expression levels in the muscle of GSDV patients. These proteins could be potential biomarkers of the persistent muscle damage in the absence of prior exertion reported in GSDV patients. Further studies are needed to elucidate the molecular mechanisms by which PYGM controls the expression of these proteins.

Strategy and Challenges of Paraclinical Examinations in Adult-Onset Still’s Disease

Adult-onset Still’s disease is a complex autoinflammatory disease with a multifactorial etiology. Its presentation is less stereotypical than that of a monogenic autoinflammatory disease and is actually relatively common with few specific signs. To avoid under- or over-prescription of complementary examinations, it is useful to advance in a structured manner, taking into consideration the actual added value of each supplemental examination. In this review, we detail the different complementary tests used in adult Still’s disease. We consider them from three different angles: positive diagnostic approach, the differential diagnosis, and the screening for complications of the disease. After discussing the various tests at our disposal, we look at the classical diagnostic strategy in order to propose a structured algorithm that can be used in clinical practice. We conclude with the prospects of new complementary examinations, which could in the future modify the management of patients.

Ceruletide and Alpha-1 Antitrypsin as a Novel Combination Therapy for Ischemic Stroke

Ischemic stroke is a primary cause of morbidity and mortality worldwide. Beyond the approved thrombolytic therapies, there is no effective treatment to mitigate its progression. Drug repositioning combinational therapies are becoming promising approaches to identify new uses of existing drugs to synergically target multiple disease-response mechanisms underlying complex pathologies. Here, we used a systems biology-based approach based on artificial intelligence and pattern recognition tools to generate in silico mathematical models mimicking the ischemic stroke pathology. Combinational treatments were acquired by screening these models with more than 5 million two-by-two combinations of drugs. A drug combination (CA) formed by ceruletide and alpha-1 antitrypsin showing a predicted value of neuroprotection of 92% was evaluated for their synergic neuroprotective effects in a mouse pre-clinical stroke model. The administration of both drugs in combination was safe and effective in reducing by 39.42% the infarct volume 24 h after cerebral ischemia. This neuroprotection was not observed when drugs were given individually. Importantly, potential incompatibilities of the drug combination with tPA thrombolysis were discarded in vitro and in vivo by using a mouse thromboembolic stroke model with t-PA-induced reperfusion, revealing an improvement in the forepaw strength 72 h after stroke in CA-treated mice. Finally, we identified the predicted mechanisms of action of ceruletide and alpha-1 antitrypsin and we demonstrated that CA modulates EGFR and ANGPT-1 levels in circulation within the acute phase after stroke. In conclusion, we have identified a promising combinational treatment with neuroprotective effects for the treatment of ischemic stroke.

Exploring the synergistic effects of cabozantinib and a programmed cell death protein 1 inhibitor in metastatic renal cell carcinoma with machine learning

Clinical evidence supports the combination of cabozantinib with an immune checkpoint inhibitor for the treatment of metastatic clear cell renal cell carcinoma (mccRCC) and suggests a synergistic antitumour activity of this combination. Nevertheless, the biological basis of this synergy is not fully characterized. We studied the mechanisms underpinning the potential synergism of cabozantinib combined with a PD1 inhibitor in mccRCC and delved into cabozantinib monotherapy properties supporting its role to partner these combinations. To model physiological drug action, we used a machine learning-based technology known as Therapeutic Performance Mapping Systems, applying two approaches: Artificial Neural Networks and Sampling Methods. We found that the combined therapy was predicted to exert a wide therapeutic action in the tumour and the microenvironment. Cabozantinib may enhance the effects of PD1 inhibitors on immunosurveillance by modulating the innate and adaptive immune system, through the inhibition of VEGF-VEGFR and Gas6-AXL/TYRO3/MER (TAM) axes, while the PD1 inhibitors may boost the antiangiogenic and pro–apoptotic effects of cabozantinib by modulating angiogenesis and T-cell cytotoxicity. Cabozantinib alone was predicted to restore cellular adhesion and hamper tumour proliferation and invasion. These data provide a biological rationale and further support for cabozantinib plus PD1 inhibitor combination and may guide future nonclinical and clinical research.

Elucidating the Mechanism of Action of the Attributed Immunomodulatory Role of Eltrombopag in Primary Immune Thrombocytopenia: An In Silico Approach

Eltrombopag is a thrombopoietin receptor (MPL) agonist approved for the treatment of primary immune thrombocytopenia (ITP). Recent evidence shows that some patients may sustain platelet counts following eltrombopag discontinuation. The systemic immunomodulatory response that resolves ITP in some patients could result from an increase in platelet mass, caused either by the direct action of eltrombopag on megakaryocytes through MPL stimulation, or potential MPL-independent actions on other cell types. To uncover the possible mechanisms of action of eltrombopag, in silico analyses were performed, including a systems biology-based approach, a therapeutic performance mapping system, and structural analyses. Through manual curation of the available bibliography, 56 key proteins were identified and integrated into the ITP interactome analysis. Mathematical models (94.92% mean accuracy) were obtained to elucidate potential MPL-dependent pathways in non-megakaryocytic cell subtypes. In addition to the effects on megakaryocytes and platelet numbers, the results were consistent with MPL-mediated effects on other cells, which could involve interferon-gamma, transforming growth factor-beta, peroxisome proliferator-activated receptor-gamma, and forkhead box protein P3 pathways. Structural analyses indicated that effects on three apoptosis-related proteins (BCL2L1, BCL2, BAX) from the Bcl-2 family may be off-target effects of eltrombopag. In conclusion, this study proposes new hypotheses regarding the immunomodulatory functions of eltrombopag in patients with ITP.