Integrating single cell sequencing with a spatial quantitative systems pharmacology model spQSP for personalized prediction of triple-negative breast cancer immunotherapy response

A comprehensive review of computational cell cycle models in guiding cancer treatment strategies

This article reviews the current knowledge and recent advancements in computational modeling of the cell cycle. It offers a comparative analysis of various modeling paradigms, highlighting their unique strengths, limitations, and applications. Specifically, the article compares deterministic and stochastic models, single-cell versus population models, and mechanistic versus abstract models. This detailed analysis helps determine the most suitable modeling framework for various research needs. Additionally, the discussion extends to the utilization of these computational models to illuminate cell cycle dynamics, with a particular focus on cell cycle viability, crosstalk with signaling pathways, tumor microenvironment, DNA replication, and repair mechanisms, underscoring their critical roles in tumor progression and the optimization of cancer therapies. By applying these models to crucial aspects of cancer therapy planning for better outcomes, including drug efficacy quantification, drug discovery, drug resistance analysis, and dose optimization, the review highlights the significant potential of computational insights in enhancing the precision and effectiveness of cancer treatments. This emphasis on the intricate relationship between computational modeling and therapeutic strategy development underscores the pivotal role of advanced modeling techniques in navigating the complexities of cell cycle dynamics and their implications for cancer therapy.

Biomarkers and computational models for predicting efficacy to tumor ICI immunotherapy

Numerous studies have shown that immune checkpoint inhibitor (ICI) immunotherapy has great potential as a cancer treatment, leading to significant clinical improvements in numerous cases. However, it benefits a minority of patients, underscoring the importance of discovering reliable biomarkers that can be used to screen for potential beneficiaries and ultimately reduce the risk of overtreatment. Our comprehensive review focuses on the latest advancements in predictive biomarkers for ICI therapy, particularly emphasizing those that enhance the efficacy of programmed cell death protein 1 (PD-1)/programmed cell death-ligand 1 (PD-L1) inhibitors and cytotoxic T-lymphocyte antigen-4 (CTLA-4) inhibitors immunotherapies. We explore biomarkers derived from various sources, including tumor cells, the tumor immune microenvironment (TIME), body fluids, gut microbes, and metabolites. Among them, tumor cells-derived biomarkers include tumor mutational burden (TMB) biomarker, tumor neoantigen burden (TNB) biomarker, microsatellite instability (MSI) biomarker, PD-L1 expression biomarker, mutated gene biomarkers in pathways, and epigenetic biomarkers. TIME-derived biomarkers include immune landscape of TIME biomarkers, inhibitory checkpoints biomarkers, and immune repertoire biomarkers. We also discuss various techniques used to detect and assess these biomarkers, detailing their respective datasets, strengths, weaknesses, and evaluative metrics. Furthermore, we present a comprehensive review of computer models for predicting the response to ICI therapy. The computer models include knowledge-based mechanistic models and data-based machine learning (ML) models. Among the knowledge-based mechanistic models are pharmacokinetic/pharmacodynamic (PK/PD) models, partial differential equation (PDE) models, signal networks-based models, quantitative systems pharmacology (QSP) models, and agent-based models (ABMs). ML models include linear regression models, logistic regression models, support vector machine (SVM)/random forest/extra trees/k-nearest neighbors (KNN) models, artificial neural network (ANN) and deep learning models. Additionally, there are hybrid models of systems biology and ML. We summarized the details of these models, outlining the datasets they utilize, their evaluation methods/metrics, and their respective strengths and limitations. By summarizing the major advances in the research on predictive biomarkers and computer models for the therapeutic effect and clinical utility of tumor ICI, we aim to assist researchers in choosing appropriate biomarkers or computer models for research exploration and help clinicians conduct precision medicine by selecting the best biomarkers.

Identification of indole-based natural compounds as inhibitors of PARP-1 against triple-negative breast cancer: a computational study

Triple-negative breast cancer (TNBC) is the most aggressive kind of breast cancer known to mankind. It is a heterogeneous disease that is formed due to the missing estrogen, progesterone and human epidermal growth factor 2 receptors. Poly(ADP-ribose) polymerase-1 (PARP-1) protein helps in the development of TNBC by repairing the cancer cells, which proliferate and spread metastatically. To determine the potential PARP-1 inhibitors (PARPi), 0.2 million natural products from Universal Natural Product Database were screened using molecular docking and six hit compounds were selected based on their binding affinity towards PARP-1. The bio-availability and drug-like properties of these natural products were evaluated using ADMET analysis. Molecular dynamics simulations were conducted for these complexes for 200 ns to examine their structural stability and dynamic behaviour and further compared with the complex of talazoparib (TALA), an FDA-approved PARPi. Using MM/PBSA calculations, we conclude that the complexes HIT-3 and HIT-5 (-25.64 and -23.14 kcal/mol, respectively) show stronger binding energies with PARP-1 than TALA with PARP-1 (-10.74 kcal/mol). Strong interactions were observed between the compounds and hotspot residues, Asp770, Ala880, Tyr889, Tyr896, Ala898, Asp899 and Tyr907, of PARP-1 due to the existence of various types of non-covalent interactions between the compounds and PARP-1. This research offers critical information about PARPi, which could potentially be incorporated into the treatment of TNBC. Moreover, these findings were validated by comparing them with an FDA-approved PARPi.

Informing virtual clinical trials of hepatocellular carcinoma with spatial multi-omics analysis of a human neoadjuvant immunotherapy clinical trial

Human clinical trials are important tools to advance novel systemic therapies improve treatment outcomes for cancer patients. The few durable treatment options have led to a critical need to advance new therapeutics in hepatocellular carcinoma (HCC). Recent human clinical trials have shown that new combination immunotherapeutic regimens provide unprecedented clinical response in a subset of patients. Computational methods that can simulate tumors from mathematical equations describing cellular and molecular interactions are emerging as promising tools to simulate the impact of therapy entirely in silico. To facilitate designing dosing regimen and identifying potential biomarkers, we developed a new computational model to track tumor progression at organ scale while reflecting the spatial heterogeneity in the tumor at tissue scale in HCC. This computational model is called a spatial quantitative systems pharmacology (spQSP) platform and it is also designed to simulate the effects of combination immunotherapy. We then validate the results from the spQSP system by leveraging real-world spatial multi-omics data from a neoadjuvant HCC clinical trial combining anti-PD-1 immunotherapy and a multitargeted tyrosine kinase inhibitor (TKI) cabozantinib. The model output is compared with spatial data from Imaging Mass Cytometry (IMC). Both IMC data and simulation results suggest closer proximity between CD8 T cell and macrophages among non-responders while the reverse trend was observed for responders. The analyses also imply wider dispersion of immune cells and less scattered cancer cells in responders' samples. We also compared the model output with Visium spatial transcriptomics analyses of samples from post-treatment tumor resections in the original clinical trial. Both spatial transcriptomic data and simulation results identify the role of spatial patterns of tumor vasculature and TGFβ in tumor and immune cell interactions. To our knowledge, this is the first spatial tumor model for virtual clinical trials at a molecular scale that is grounded in high-throughput spatial multi-omics data from a human clinical trial.

A transcriptome-informed QSP model of metastatic triple-negative breast cancer identifies predictive biomarkers for PD-1 inhibition

Triple-negative breast cancer (TNBC), a highly metastatic breast cancer subtype, has limited treatment options. While a small number of patients attain clinical benefit with single-agent checkpoint inhibitors, identifying these patients before the therapy remains challenging. Here, we developed a transcriptome-informed quantitative systems pharmacology model of metastatic TNBC by integrating heterogenous metastatic tumors. In silico clinical trial with an anti-PD-1 drug, pembrolizumab, predicted that several features, such as the density of antigen-presenting cells, the fraction of cytotoxic T cells in lymph nodes, and the richness of cancer clones in tumors, could serve individually as biomarkers but had a higher predictive power as combinations of two biomarkers. We showed that PD-1 inhibition neither consistently enhanced all antitumorigenic factors nor suppressed all protumorigenic factors but ultimately reduced the tumor carrying capacity. Collectively, our predictions suggest several candidate biomarkers that might effectively predict the response to pembrolizumab monotherapy and potential therapeutic targets to develop treatment strategies for metastatic TNBC.

Biomarkers in Breast Cancer: An Old Story with a New End

Breast cancer is the second most frequent cancer in the world. It is a heterogeneous disease and the leading cause of cancer mortality in women. Advances in molecular technologies allowed for the identification of new and more specifics biomarkers for breast cancer diagnosis, prognosis, and risk prediction, enabling personalized treatments, improving therapy, and preventing overtreatment, undertreatment, and incorrect treatment. Several breast cancer biomarkers have been identified and, along with traditional biomarkers, they can assist physicians throughout treatment plan and increase therapy success. Despite the need of more data to improve specificity and determine the real clinical utility of some biomarkers, others are already established and can be used as a guide to make treatment decisions. In this review, we summarize the available traditional, novel, and potential biomarkers while also including gene expression profiles, breast cancer single-cell and polyploid giant cancer cells. We hope to help physicians understand tumor specific characteristics and support decision-making in patient-personalized clinical management, consequently improving treatment outcome.

Quantifying Intratumoral Heterogeneity and Immunoarchitecture Generated In-Silico by a Spatial Quantitative Systems Pharmacology Model

Spatial heterogeneity is a hallmark of cancer. Tumor heterogeneity can vary with time and location. The tumor microenvironment (TME) encompasses various cell types and their interactions that impart response to therapies. Therefore, a quantitative evaluation of tumor heterogeneity is crucial for the development of effective treatments. Different approaches, such as multiregional sequencing, spatial transcriptomics, analysis of autopsy samples, and longitudinal analysis of biopsy samples, can be used to analyze the intratumoral heterogeneity (ITH) and temporal evolution and to reveal the mechanisms of therapeutic response. However, because of the limitations of these data and the uncertainty associated with the time points of sample collection, having a complete understanding of intratumoral heterogeneity role is challenging. Here, we used a hybrid model that integrates a whole-patient compartmental quantitative-systems-pharmacology (QSP) model with a spatial agent-based model (ABM) describing the TME; we applied four spatial metrics to quantify model-simulated intratumoral heterogeneity and classified the TME immunoarchitecture for representative cases of effective and ineffective anti-PD-1 therapy. The four metrics, adopted from computational digital pathology, included mixing score, average neighbor frequency, Shannon's entropy and area under the curve (AUC) of the G-cross function. A fifth non-spatial metric was used to supplement the analysis, which was the ratio of the number of cancer cells to immune cells. These metrics were utilized to classify the TME as "cold", "compartmentalized" and "mixed", which were related to treatment efficacy. The trends in these metrics for effective and ineffective treatments are in qualitative agreement with the clinical literature, indicating that compartmentalized immunoarchitecture is likely to result in more efficacious treatment outcomes.

Quantitative Spatial Profiling of Immune Populations in Pancreatic Ductal Adenocarcinoma Reveals Tumor Microenvironment Heterogeneity and Prognostic Biomarkers

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive disease with poor 5-year survival rates, necessitating identification of novel therapeutic targets. Elucidating the biology of the tumor immune microenvironment (TiME) can provide vital insights into mechanisms of tumor progression. In this study, we developed a quantitative image processing platform to analyze sequential multiplexed IHC data from archival PDAC tissue resection specimens. A 27-plex marker panel was employed to simultaneously phenotype cell populations and their functional states, followed by a computational workflow to interrogate the immune contextures of the TiME in search of potential biomarkers. The PDAC TiME reflected a low-immunogenic ecosystem with both high intratumoral and intertumoral heterogeneity. Spatial analysis revealed that the relative distance between IL10+ myelomonocytes, PD-1+ CD4+ T cells, and granzyme B+ CD8+ T cells correlated significantly with survival, from which a spatial proximity signature termed imRS was derived that correlated with PDAC patient survival. Furthermore, spatial enrichment of CD8+ T cells in lymphoid aggregates was also linked to improved survival. Altogether, these findings indicate that the PDAC TiME, generally considered immuno-dormant or immunosuppressive, is a spatially nuanced ecosystem orchestrated by ordered immune hierarchies. This new understanding of spatial complexity may guide novel treatment strategies for PDAC.

SIGNIFICANCE

Quantitative image analysis of PDAC specimens reveals intertumoral and intratumoral heterogeneity of immune populations and identifies spatial immune architectures that are significantly associated with disease prognosis.

The tumor microenvironment and triple-negative breast cancer aggressiveness: shedding light on mechanisms and targeting

INTRODUCTION

In contrast to other breast cancer subtypes, there are currently limited options of targeted therapies for triple-negative breast cancer (TNBC). Immense research has demonstrated that not only cancer cells but also stromal cells and immune cells in the tumor microenvironment (TME) play significant roles in the progression of TNBC. It is thus critical to understand the components of the TME of TNBC and the interactions between the various cell populations.

AREAS COVERED

The components of the TME of TNBC identified by single-cell technologies are reviewed. Furthermore, the molecular interactions between the cells and the potential therapeutic targets contributing to the progression of TNBC are discussed.

EXPERT OPINION

Single-cell omics studies have contributed to the classification of cells in the TME and the identification of important cell types involved in the progression and the treatment of the tumor. The interactions between cancer cells and stromal cells/immune cells in the TME have led to the discovery of potential therapeutic targets. Experimental data with spatial and temporal resolution will further boost the understanding of the TME of TNBC.

Virtual clinical trials of anti-PD-1 and anti-CTLA-4 immunotherapy in advanced hepatocellular carcinoma using a quantitative systems pharmacology model

BACKGROUND

Hepatocellular carcinoma (HCC) is the most common form of primary liver cancer and is the third-leading cause of cancer-related death worldwide. Most patients with HCC are diagnosed at an advanced stage, and the median survival for patients with advanced HCC treated with modern systemic therapy is less than 2 years. This leaves the advanced stage patients with limited treatment options. Immune checkpoint inhibitors (ICIs) targeting programmed cell death protein 1 (PD-1) or its ligand, are widely used in the treatment of HCC and are associated with durable responses in a subset of patients. ICIs targeting cytotoxic T-lymphocyte-associated protein 4 (CTLA-4) also have clinical activity in HCC. Combination therapy of nivolumab (anti-PD-1) and ipilimumab (anti-CTLA-4) is the first treatment option for HCC to be approved by Food and Drug Administration that targets more than one immune checkpoints.

METHODS

In this study, we used the framework of quantitative systems pharmacology (QSP) to perform a virtual clinical trial for nivolumab and ipilimumab in HCC patients. Our model incorporates detailed biological mechanisms of interactions of immune cells and cancer cells leading to antitumor response. To conduct virtual clinical trial, we generate virtual patient from a cohort of 5,000 proposed patients by extending recent algorithms from literature. The model was calibrated using the data of the clinical trial CheckMate 040 (ClinicalTrials.gov number, NCT01658878).

RESULTS

Retrospective analyses were performed for different immune checkpoint therapies as performed in CheckMate 040. Using machine learning approach, we predict the importance of potential biomarkers for immune blockade therapies.

CONCLUSIONS

This is the first QSP model for HCC with ICIs and the predictions are consistent with clinically observed outcomes. This study demonstrates that using a mechanistic understanding of the underlying pathophysiology, QSP models can facilitate patient selection and design clinical trials with improved success.

Dynamics of tumor-associated macrophages in a quantitative systems pharmacology model of immunotherapy in triple-negative breast cancer

Quantitative systems pharmacology (QSP) modeling is an emerging mechanistic computational approach that couples drug pharmacokinetics/pharmacodynamics and the course of disease progression. It has begun to play important roles in drug development for complex diseases such as cancer, including triple-negative breast cancer (TNBC). The combination of the anti-PD-L1 antibody atezolizumab and nab-paclitaxel has shown clinical activity in advanced TNBC with PD-L1-positive tumor-infiltrating immune cells. As tumor-associated macrophages (TAMs) serve as major contributors to the immuno-suppressive tumor microenvironment, we incorporated the dynamics of TAMs into our previously published QSP model to investigate their impact on cancer treatment. We show that through proper calibration, the model captures the macrophage heterogeneity in the tumor microenvironment while maintaining its predictive power of the trial results at the population level. Despite its high mechanistic complexity, the modularized QSP platform can be readily reproduced, expanded for new species of interest, and applied in clinical trial simulation.

•

A mechanistic model of quantitative systems pharmacology in immuno-oncology

•

Dynamics of tumor-associated macrophages are integrated into our previous work

•

Conducting in silico clinical trials to predict clinical response to cancer therapy

Simulations of tumor growth and response to immunotherapy by coupling a spatial agent-based model with a whole-patient quantitative systems pharmacology model

Quantitative systems pharmacology (QSP) models and spatial agent-based models (ABM) are powerful and efficient approaches for the analysis of biological systems and for clinical applications. Although QSP models are becoming essential in discovering predictive biomarkers and developing combination therapies through in silico virtual trials, they are inadequate to capture the spatial heterogeneity and randomness that characterize complex biological systems, and specifically the tumor microenvironment. Here, we extend our recently developed spatial QSP (spQSP) model to analyze tumor growth dynamics and its response to immunotherapy at different spatio-temporal scales. In the model, the tumor spatial dynamics is governed by the ABM, coupled to the QSP model, which includes the following compartments: central (blood system), tumor, tumor-draining lymph node, and peripheral (the rest of the organs and tissues). A dynamic recruitment of T cells and myeloid-derived suppressor cells (MDSC) from the QSP central compartment has been implemented as a function of the spatial distribution of cancer cells. The proposed QSP-ABM coupling methodology enables the spQSP model to perform as a coarse-grained model at the whole-tumor scale and as an agent-based model at the regions of interest (ROIs) scale. Thus, we exploit the spQSP model potential to characterize tumor growth, identify T cell hotspots, and perform qualitative and quantitative descriptions of cell density profiles at the invasive front of the tumor. Additionally, we analyze the effects of immunotherapy at both whole-tumor and ROI scales under different tumor growth and immune response conditions. A digital pathology computational analysis of triple-negative breast cancer specimens is used as a guide for modeling the immuno-architecture of the invasive front.

Prognostic/predictive markers in systemic therapy resistance and metastasis in breast cancer

Breast cancer is a highly heterogeneous group of diseases posing a significant challenge in biomarker-driven research and the development of effective targeted therapies. Especially the treatment of metastatic breast cancer poses even more challenges, as we still lose more than 42,000 women and men each year in the United States alone. New biological insight helps to improve breast cancer treatment through early detection, adaptation to chemotherapy resistance, and tailoring to find the right size of care. This review focuses on existing and new areas of predictive biomarkers under development to tailor the management of breast cancer and the application of integrative approaches that have resulted in the promising candidate biomarker discovery. Furthermore, we review new methods to detect metastatic progression using imaging, and blood-based assays. We hope to increase the attention and awareness of a new generation of therapeutic development strategies in metastatic breast cancer.

From Fitting the Average to Fitting the Individual: A Cautionary Tale for Mathematical Modelers

An outstanding challenge in the clinical care of cancer is moving from a one-size-fits-all approach that relies on population-level statistics towards personalized therapeutic design. Mathematical modeling is a powerful tool in treatment personalization, as it allows for the incorporation of patient-specific data so that treatment can be tailor-designed to the individual. Herein, we work with a mathematical model of murine cancer immunotherapy that has been previously-validated against the average of an experimental dataset. We ask the question: what happens if we try to use this same model to perform personalized fits, and therefore make individualized treatment recommendations? Typically, this would be done by choosing a single fitting methodology, and a single cost function, identifying the individualized best-fit parameters, and extrapolating from there to make personalized treatment recommendations. Our analyses show the potentially problematic nature of this approach, as predicted personalized treatment response proved to be sensitive to the fitting methodology utilized. We also demonstrate how a small amount of the right additional experimental measurements could go a long way to improve consistency in personalized fits. Finally, we show how quantifying the robustness of the average response could also help improve confidence in personalized treatment recommendations.

Multi-Scale Spatial Analysis of the Tumor Microenvironment Reveals Features of Cabozantinib and Nivolumab Efficacy in Hepatocellular Carcinoma

Background

Concomitant inhibition of vascular endothelial growth factor (VEGF) and programmed cell death protein 1 (PD-1) or its ligand PD-L1 is a standard of care for patients with advanced hepatocellular carcinoma (HCC), but only a minority of patients respond, and responses are usually transient. Understanding the effects of therapies on the tumor microenvironment (TME) can provide insights into mechanisms of therapeutic resistance.

Methods

14 patients with HCC were treated with the combination of cabozantinib and nivolumab through the Johns Hopkins Sidney Kimmel Comprehensive Cancer Center. Among them, 12 patients (5 responders + 7 non-responders) underwent successful margin negative resection and are subjects to tissue microarray (TMA) construction containing 37 representative tumor region cores. Using the TMAs, we performed imaging mass cytometry (IMC) with a panel of 27-cell lineage and functional markers. All multiplexed images were then segmented to generate a single-cell dataset that enables (1) tumor-immune compartment analysis and (2) cell community analysis based on graph-embedding methodology. Results from these hierarchies are merged into response-associated biological process patterns.

Results

Image processing on 37 multiplexed-images discriminated 59,453 cells and was then clustered into 17 cell types. Compartment analysis showed that at immune-tumor boundaries from NR, PD-L1 level on tumor cells is significantly higher than remote regions; however, Granzyme B expression shows the opposite pattern. We also identify that the close proximity of CD8⁺ T cells to arginase 1^hi (Arg1^hi) macrophages, rather than CD4⁺ T cells, is a salient feature of the TME in non-responders. Furthermore, cell community analysis extracted 8 types of cell-cell interaction networks termed cellular communities (CCs). We observed that in non-responders, macrophage-enriched CC (MCC) and lymphocyte-enriched CC (LCC) strongly communicate with tumor CC, whereas in responders, such communications were undermined by the engagement between MCC and LCC.

Conclusion

These results demonstrate the feasibility of a novel application of multiplexed image analysis that is broadly applicable to quantitative analysis of pathology specimens in immuno-oncology and provides further evidence that CD163^-Arg1^hi macrophages may be a therapeutic target in HCC. The results also provide critical information for the development of mechanistic quantitative systems pharmacology models aimed at predicting outcomes of clinical trials.

Systems biology of angiogenesis signaling: Computational models and omics

Angiogenesis is a highly regulated multiscale process that involves a plethora of cells, their cellular signal transduction, activation, proliferation, differentiation, as well as their intercellular communication. The coordinated execution and integration of such complex signaling programs is critical for physiological angiogenesis to take place in normal growth, development, exercise, and wound healing, while its dysregulation is critically linked to many major human diseases such as cancer, cardiovascular diseases, and ocular disorders; it is also crucial in regenerative medicine. Although huge efforts have been devoted to drug development for these diseases by investigation of angiogenesis‐targeted therapies, only a few therapeutics and targets have proved effective in humans due to the innate multiscale complexity and nonlinearity in the process of angiogenic signaling. As a promising approach that can help better address this challenge, systems biology modeling allows the integration of knowledge across studies and scales and provides a powerful means to mechanistically elucidate and connect the individual molecular and cellular signaling components that function in concert to regulate angiogenesis. In this review, we summarize and discuss how systems biology modeling studies, at the pathway‐, cell‐, tissue‐, and whole body‐levels, have advanced our understanding of signaling in angiogenesis and thereby delivered new translational insights for human diseases.

This article is categorized under:

Cardiovascular Diseases > Computational Models

Cancer > Computational Models

Plant Derived Bioactive Compounds, Their Anti-Cancer Effects and In Silico Approaches as an Alternative Target Treatment Strategy for Breast Cancer: An Updated Overview

Cancer is one of the most common malignant diseases that occur worldwide, among which breast cancer is the second leading cause of death in women. The subtypes are associated with differences in the outcome and were selected for treatments according to the estrogen receptor, progesterone receptor, and human epidermal growth factor receptor. Triple-negative breast cancer, one of the subtypes of breast cancer, is difficult to treat and can even lead to death. If breast cancer is not treated during the initial stages, it may spread to nearby organs, a process called metastasis, through the blood or lymph system. For in vitro studies, MCF-7, MDA-MB-231, MDA-MB-468, and T47B are the most commonly used breast cancer cell lines. Clinically, chemotherapy and radiotherapy are usually expensive and can also cause side effects. To overcome these issues, medicinal plants could be the best alternative for chemotherapeutic drugs with fewer side effects and cost-effectiveness. Furthermore, the genes involved in breast cancer can be regulated and synergized with signaling molecules to suppress the proliferation of breast cancer cells. In addition, nanoparticles encapsulating (nano-encapsulation) medicinal plant extracts showed a significant reduction in the apoptotic and cytotoxic activities of breast cancer cells. This present review mainly speculates an overview of the native medicinal plant derived anti-cancerous compounds with its efficiency, types and pathways involved in breast cancer along with its genes, the mechanism of breast cancer brain metastasis, chemoresistivity and its mechanism, bioinformatics approaches which could be an effective alternative for drug discovery.