In Silico Modeling of the Induction of Apoptosis by Cryptotanshinone in Osteosarcoma Cell Lines.

Osteosarcoma (OS) is the most common primary malignant bone tumor of both children and pet canines. Its characteristic genomic instability and complexity coupled with the dearth of knowledge about its etiology has made improvement in the current treatment difficult. We use the existing literature about the biological pathways active in OS and combine it with the current research involving natural compounds to identify new targets and design more effective drug therapies. The key components of these pathways are modeled as a Boolean network with multiple inputs and multiple outputs. The combinatorial circuit is employed to theoretically predict the efficacies of various drugs in combination with Cryptotanshinone. We show that the action of the herbal drug, Cryptotanshinone on OS cell lines induces apoptosis by increasing sensitivity to TNF-related apoptosis-inducing ligand (TRAIL) through its multi-pronged action on STAT3, DRP1 and DR5. The Boolean framework is used to detect additional drug intervention points in the pathway that could amplify the action of Cryptotanshinone.

Comprehensive live-cell imaging analysis of cryptotanshinone and synergistic drug-screening effects in various human and canine cancer cell lines.

BACKGROUND: Several studies have highlighted both the extreme anticancer effects of Cryptotanshinone (CT), a Stat3 crippling component from Salvia miltiorrhiza, as well as other STAT3 inhibitors to fight cancer. METHODS: Data presented in this experiment incorporates 2 years of in vitro studies applying a comprehensive live-cell drug-screening analysis of human and canine cancer cells exposed to CT at 20 µM concentration, as well as to other drug combinations. As previously observed in other studies, dogs are natural cancer models, given to their similarity in cancer genetics, epidemiology and disease progression compared to humans. RESULTS: Results obtained from several types of human and canine cancer cells exposed to CT and varied drug combinations, verified CT efficacy at combating cancer by achieving an extremely high percentage of apoptosis within 24 hours of drug exposure. CONCLUSIONS: CT anticancer efficacy in various human and canine cancer cell lines denotes its ability to interact across different biological processes and cancer regulatory cell networks, driving inhibition of cancer cell survival.

A Gaussian Mixture-Model Exploiting Pathway Knowledge for Dissecting Cancer Heterogeneity.

In this work, we develop a systematic approach for applying pathway knowledge to a multivariate Gaussian mixture model for dissecting a heterogeneous cancer tissue. The downstream transcription factors are selected as observables from available partial pathway knowledge in such a way that the subpopulations produce some differential behavior in response to the drugs selected in the upstream. For each subpopulation, each unique (drug, observable) pair is considered as a unique dimension of a multivariate Gaussian distribution. Expectation-maximization (EM) algorithm with hill-climbing is then used to rank the most probable estimates of the mixture composition based on the log-likelihood value. A major contribution of this work is to examine the efficacy of the EM based approach in estimating the composition of experimental mixture sets from cell-by-cell measurements collected on a dynamic cell imaging platform. Towards this end, we apply the algorithm on hourly data collected for two different mixture compositions of A2058, HCT116, and SW480 cell lines for three scenarios: untreated, Lapatinib-treated, and Temsirolimus-treated. Additionally, we show how this methodology can provide a basis for comparing the killing rate of different drugs for a heterogeneous cancer tissue. This obviously has important implications for designing efficient drugs for treating heterogeneous malignant tumors.

Phenotype Classification Using Moment Features of Single-Cell Data.

Features for standard expression microarray and RNA-Seq classification are expression averages over collections of cells. Single cell provides expression measurements for individual cells in a collection of cells from a particular tissue sample. Hence, it can yield feature vectors consisting of higher order and mixed moments. This article demonstrates the advantage of using these expression moments in cancer-related classification. We use synthetic data generated from 2 real networks, the mammalian cell cycle network and a melanoma-related pathway network, and real single-cell data generated via fluorescent protein reporters from 2 cell lines, HT-29 and HCT-116. The networks consist of hidden binary regulatory networks with Gaussian observations. The steady-state distributions of both the original and mutated networks are found, and data are drawn from these for moment-based classification using the mean, variance, skewness, and mixed moments. For the real data, we only observe 1 gene at a time, so that only the mean, variance, and skewness are considered, the analysis being done for 2 genes, EGFR and ERRB2. For the synthetic data, classification improves as we move from just the mean to mean, variance, and skewness and then to these plus the mixed moments. Comparisons are done with 3, 4, or 5 features, using feature selection. Sample size effects are considered. For the real data, we only consider mean, variance, and skewness, with results improving when the higher order moments are used as features.

A Bayesian approach to determine the composition of heterogeneous cancer tissue.

BACKGROUND: Cancer Tissue Heterogeneity is an important consideration in cancer research as it can give insights into the causes and progression of cancer. It is known to play a significant role in cancer cell survival, growth and metastasis. Determining the compositional breakup of a heterogeneous cancer tissue can also help address the therapeutic challenges posed by heterogeneity. This necessitates a low cost, scalable algorithm to address the challenge of accurate estimation of the composition of a heterogeneous cancer tissue. METHODS: In this paper, we propose an algorithm to tackle this problem by utilizing the data of accurate, but high cost, single cell line cell-by-cell observation methods in low cost aggregate observation method for heterogeneous cancer cell mixtures to obtain their composition in a Bayesian framework. RESULTS: The algorithm is analyzed and validated using synthetic data and experimental data. The experimental data is obtained from mixtures of three separate human cancer cell lines, HCT116 (Colorectal carcinoma), A2058 (Melanoma) and SW480 (Colorectal carcinoma). CONCLUSION: The algorithm provides a low cost framework to determine the composition of heterogeneous cancer tissue which is a crucial aspect in cancer research.

Systems biology for organotypic cell cultures.

Translating in vitro biological data into actionable information related to human health holds the potential to improve disease treatment and risk assessment of chemical exposures. While genomics has identified regulatory pathways at the cellular level, translation to the organism level requires a multiscale approach accounting for intra-cellular regulation, inter-cellular interaction, and tissue/organ-level effects. Tissue-level effects can now be probed in vitro thanks to recently developed systems of three-dimensional (3D), multicellular, "organotypic" cell cultures, which mimic functional responses of living tissue. However, there remains a knowledge gap regarding interactions across different biological scales, complicating accurate prediction of health outcomes from molecular/genomic data and tissue responses. Systems biology aims at mathematical modeling of complex, non-linear biological systems. We propose to apply a systems biology approach to achieve a computational representation of tissue-level physiological responses by integrating empirical data derived from organotypic culture systems with computational models of intracellular pathways to better predict human responses. Successful implementation of this integrated approach will provide a powerful tool for faster, more accurate and cost-effective screening of potential toxicants and therapeutics. On September 11, 2015, an interdisciplinary group of scientists, engineers, and clinicians gathered for a workshop in Research Triangle Park, North Carolina, to discuss this ambitious goal. Participants represented laboratory-based and computational modeling approaches to pharmacology and toxicology, as well as the pharmaceutical industry, government, non-profits, and academia. Discussions focused on identifying critical system perturbations to model, the computational tools required, and the experimental approaches best suited to generating key data.

Assessing Combinational Drug Efficacy in Cancer Cells by Using Image-based Dynamic Response Analysis.

The landscape of translational research has been shifting toward drug combination therapies. Pairing of drugs allows for more types of drug interaction with cells. In order to accurately and comprehensively assess combinational drug efficacy, analytical methods capable of recognizing these alternative reactions will be required to prioritize those drug candidates having better chances of delivering appreciable therapeutic benefits. Traditional efficacy measures are primarily based on the "extent" of drug inhibition, which is the percentage of cells being killed after drug exposure. Here, we introduce a second dimension of evaluation criterion, speed of killing, based on a live cell imaging assay. This dynamic response trajectory approach takes advantage of both "extent" and "speed" information and uncovers synergisms that would otherwise be missed, while also generating hypotheses regarding important mechanistic modes of drug action.

Evaluating gene set enrichment analysis via a hybrid data model.

Gene set enrichment analysis (GSA) methods have been widely adopted by biological labs to analyze data and generate hypotheses for validation. Most of the existing comparison studies focus on whether the existing GSA methods can produce accurate P-values; however, practitioners are often more concerned with the correct gene-set ranking generated by the methods. The ranking performance is closely related to two critical goals associated with GSA methods: the ability to reveal biological themes and ensuring reproducibility, especially for small-sample studies. We have conducted a comprehensive simulation study focusing on the ranking performance of seven representative GSA methods. We overcome the limitation on the availability of real data sets by creating hybrid data models from existing large data sets. To build the data model, we pick a master gene from the data set to form the ground truth and artificially generate the phenotype labels. Multiple hybrid data models can be constructed from one data set and multiple data sets of smaller sizes can be generated by resampling the original data set. This approach enables us to generate a large batch of data sets to check the ranking performance of GSA methods. Our simulation study reveals that for the proposed data model, the Q2 type GSA methods have in general better performance than other GSA methods and the global test has the most robust results. The properties of a data set play a critical role in the performance. For the data sets with highly connected genes, all GSA methods suffer significantly in performance.

Newton, laplace, and the epistemology of systems biology.

For science, theoretical or applied, to significantly advance, researchers must use the most appropriate mathematical methods. A century and a half elapsed between Newton's development of the calculus and Laplace's development of celestial mechanics. One cannot imagine the latter without the former. Today, more than three-quarters of a century has elapsed since the birth of stochastic systems theory. This article provides a perspective on the utilization of systems theory as the proper vehicle for the development of systems biology and its application to complex regulatory diseases such as cancer.

Assessing the efficacy of molecularly targeted agents on cell line-based platforms by using system identification.

BACKGROUND: Molecularly targeted agents (MTAs) are increasingly used for cancer treatment, the goal being to improve the efficacy and selectivity of cancer treatment by developing agents that block the growth of cancer cells by interfering with specific targeted molecules needed for carcinogenesis and tumor growth. This approach differs from traditional cytotoxic anticancer drugs. The lack of specificity of cytotoxic drugs allows a relatively straightforward approach in preclinical and clinical studies, where the optimal dose has usually been defined as the "maximum tolerated dose" (MTD). This toxicity-based dosing approach is founded on the assumption that the therapeutic anticancer effect and toxic effects of the drug increase in parallel as the dose is escalated. On the contrary, most MTAs are expected to be more selective and less toxic than cytotoxic drugs. Consequently, the maximum therapeutic effect may be achieved at a "biologically effective dose" (BED) well below the MTD. Hence, dosing study for MTAs should be different from cytotoxic drugs. Enhanced efforts to molecularly characterize the drug efficacy for MTAs in preclinical models will be valuable for successfully designing dosing regimens for clinical trials. RESULTS: A novel preclinical model combining experimental methods and theoretical analysis is proposed to investigate the mechanism of action and identify pharmacodynamic characteristics of the drug. Instead of fixed time point analysis of the drug exposure to drug effect, the time course of drug effect for different doses is quantitatively studied on cell line-based platforms using system identification, where tumor cells' responses to drugs through the use of fluorescent reporters are sampled over a time course. Results show that drug effect is time-varying and higher dosages induce faster and stronger responses as expected. However, the drug efficacy change along different dosages is not linear; on the contrary, there exist certain thresholds. This kind of preclinical study can provide valuable suggestions about dosing regimens for the in vivo experimental stage to increase productivity.

Tracking transcriptional activities with high-content epifluorescent imaging.

High-content cell imaging based on fluorescent protein reporters has recently been used to track the transcriptional activities of multiple genes under different external stimuli for extended periods. This technology enhances our ability to discover treatment-induced regulatory mechanisms, temporally order their onsets and recognize their relationships. To fully realize these possibilities and explore their potential in biological and pharmaceutical applications, we introduce a new data processing procedure to extract information about the dynamics of cell processes based on this technology. The proposed procedure contains two parts: (1) image processing, where the fluorescent images are processed to identify individual cells and allow their transcriptional activity levels to be quantified; and (2) data representation, where the extracted time course data are summarized and represented in a way that facilitates efficient evaluation. Experiments show that the proposed procedure achieves fast and robust image segmentation with sufficient accuracy. The extracted cellular dynamics are highly reproducible and sensitive enough to detect subtle activity differences and identify mechanisms responding to selected perturbations. This method should be able to help biologists identify the alterations of cellular mechanisms that allow drug candidates to change cell behavior and thereby improve the efficiency of drug discovery and treatment design.

Identifying mechanistic similarities in drug responses.

MOTIVATION: In early drug development, it would be beneficial to be able to identify those dynamic patterns of gene response that indicate that drugs targeting a particular gene will be likely or not to elicit the desired response. One approach would be to quantitate the degree of similarity between the responses that cells show when exposed to drugs, so that consistencies in the regulation of cellular response processes that produce success or failure can be more readily identified. RESULTS: We track drug response using fluorescent proteins as transcription activity reporters. Our basic assumption is that drugs inducing very similar alteration in transcriptional regulation will produce similar temporal trajectories on many of the reporter proteins and hence be identified as having similarities in their mechanisms of action (MOA). The main body of this work is devoted to characterizing similarity in temporal trajectories/signals. To do so, we must first identify the key points that determine mechanistic similarity between two drug responses. Directly comparing points on the two signals is unrealistic, as it cannot handle delays and speed variations on the time axis. Hence, to capture the similarities between reporter responses, we develop an alignment algorithm that is robust to noise, time delays and is able to find all the contiguous parts of signals centered about a core alignment (reflecting a core mechanism in drug response). Applying the proposed algorithm to a range of real drug experiments shows that the result agrees well with the prior drug MOA knowledge. AVAILABILITY: The R code for the RLCSS algorithm is available at http://gsp.tamu.edu/Publications/supplementary/zhao12a.

Characterization of drug efficacy regions based on dosage and frequency schedules.

This paper proposes a framework to study the drug effect at the molecular level in order to address the following question of current interest in the drug community: Given a fixed total delivered drug, which is better, frequent small or infrequent large drug dosages? A hybrid system model is proposed to link the drug's pharmacokinetic and pharmacodynamic information, and allows the drug effects for different dosages and treatment schedules to be compared. A hybrid model facilitates the modeling of continuous quantitative changes that leads to discrete transitions. An optimal dosage-frequency regimen and the necessary and sufficient conditions for the drug to be effective are obtained analytically when the drug is designed to control a target gene. Then, we extend the analysis to the case where the target gene is part of a genetic regulatory network. A crucial observation is that there exists a "sweet spot," defined as the "drug efficacy region (DER)" in this paper, for certain dosage and frequency arrangements given the total delivered drug. This paper quantifies the therapeutic benefits of dosage regimen lying within the DER. Simulations are performed using MATLAB/SIMULINK to validate the analytical results.

Causality, randomness, intelligibility, and the epistemology of the cell.

Because the basic unit of biology is the cell, biological knowledge is rooted in the epistemology of the cell, and because life is the salient characteristic of the cell, its epistemology must be centered on its livingness, not its constituent components. The organization and regulation of these components in the pursuit of life constitute the fundamental nature of the cell. Thus, regulation sits at the heart of biological knowledge of the cell and the extraordinary complexity of this regulation conditions the kind of knowledge that can be obtained, in particular, the representation and intelligibility of that knowledge. This paper is essentially split into two parts. The first part discusses the inadequacy of everyday intelligibility and intuition in science and the consequent need for scientific theories to be expressed mathematically without appeal to commonsense categories of understanding, such as causality. Having set the backdrop, the second part addresses biological knowledge. It briefly reviews modern scientific epistemology from a general perspective and then turns to the epistemology of the cell. In analogy with a multi-faceted factory, the cell utilizes a highly parallel distributed control system to maintain its organization and regulate its dynamical operation in the face of both internal and external changes. Hence, scientific knowledge is constituted by the mathematics of stochastic dynamical systems, which model the overall relational structure of the cell and how these structures evolve over time, stochasticity being a consequence of the need to ignore a large number of factors while modeling relatively few in an extremely complex environment.

Characterization of the effectiveness of reporting lists of small feature sets relative to the accuracy of the prior biological knowledge.

When confronted with a small sample, feature-selection algorithms often fail to find good feature sets, a problem exacerbated for high-dimensional data and large feature sets. The problem is compounded by the fact that, if one obtains a feature set with a low error estimate, the estimate is unreliable because training-data-based error estimators typically perform poorly on small samples, exhibiting optimistic bias or high variance. One way around the problem is limit the number of features being considered, restrict features sets to sizes such that all feature sets can be examined by exhaustive search, and report a list of the best performing feature sets. If the list is short, then it greatly restricts the possible feature sets to be considered as candidates; however, one can expect the lowest error estimates obtained to be optimistically biased so that there may not be a close-to-optimal feature set on the list. This paper provides a power analysis of this methodology; in particular, it examines the kind of results one should expect to obtain relative to the length of the list and the number of discriminating features among those considered. Two measures are employed. The first is the probability that there is at least one feature set on the list whose true classification error is within some given tolerance of the best feature set and the second is the expected number of feature sets on the list whose true errors are within the given tolerance of the best feature set. These values are plotted as functions of the list length to generate power curves. The results show that, if the number of discriminating features is not too small-that is, the prior biological knowledge is not too poor-then one should expect, with high probability, to find good feature sets.

Conditioning-based modeling of contextual genomic regulation.

A more complete understanding of the alterations in cellular regulatory and control mechanisms that occur in the various forms of cancer has been one of the central targets of the genomic and proteomic methods that allow surveys of the abundance and/or state of cellular macromolecules. This preference is driven both by the intractability of cancer to generic therapies, assumed to be due to the highly varied molecular etiologies observed in cancer, and by the opportunity to discern and dissect the regulatory and control interactions presented by the highly diverse assortment of perturbations of regulation and control that arise in cancer. Exploiting the opportunities for inference on the regulatory and control connections offered by these revealing system perturbations is fraught with the practical problems that arise from the way biological systems operate. Two classes of regulatory action in biological systems are particularly inimical to inference, convergent regulation, where a variety of regulatory actions result in a common set of control responses (crosstalk), and divergent regulation, where a single regulatory action produces entirely different sets of control responses, depending on cellular context (conditioning). We have constructed a coarse mathematical model of the propagation of regulatory influence in such distributed, context-sensitive regulatory networks that allows a quantitative estimation of the amount of crosstalk and conditioning associated with a candidate regulatory gene taken from a set of genes that have been profiled over a series of samples where the candidate's activity varies.

Integrating global gene expression and radiation survival parameters across the 60 cell lines of the National Cancer Institute Anticancer Drug Screen.

The 60 cell lines of the National Cancer Institute Anticancer Drug Screen (NCI-60) constitute the most extensively characterized in vitro cancer cell model. They have been tested for sensitivity to more than 100,000 potential chemotherapy agents and have been profiled extensively at the DNA, RNA, protein, functional, and pharmacologic levels. We have used the NCI-60 cell lines and three additional lines to develop a database of responses of cancer cells to ionizing radiation. We compared clonogenic survival, apoptosis, and gene expression response by microarray. Although several studies have profiled relative basal gene expression in the NCI-60, this is the first comparison of large-scale gene expression changes in response to genotoxic stress. Twenty-two genes were differentially regulated in cells with low survival after 2-Gy gamma-rays; 14 genes identified lines more sensitive to 8 Gy. Unlike reported basal gene expression patterns, changes in expression in response to radiation showed little tissue-of-origin effect, except for differentiating the lymphoblastoid cell lines from other cell types. Basal expression patterns, however, discriminated well between radiosensitive and more resistant lines, possibly being more informative than radiation response signatures. The most striking patterns in the radiation data were a set of genes up-regulated preferentially in the p53 wild-type lines and a set of cell cycle regulatory genes down-regulated across the entire NCI-60 panel. The response of those genes to gamma-rays seems to be unaffected by the myriad of genetic differences across this diverse cell set; it represents the most penetrant gene expression response to ionizing radiation yet observed.

A cell proliferation and chromosomal instability signature in anaplastic thyroid carcinoma.

Here, we show that the anaplastic thyroid carcinoma (ATC) features the up-regulation of a set of genes involved in the control of cell cycle progression and chromosome segregation. This phenotype differentiates ATC from normal tissue and from well-differentiated papillary thyroid carcinoma. Transcriptional promoters of the ATC up-regulated genes are characterized by a modular organization featuring binding sites for E2F and NF-Y transcription factors and cell cycle-dependent element (CDE)/cell cycle gene homology region (CHR) cis-regulatory elements. Two protein kinases involved in cell cycle regulation, namely, Polo-like kinase 1 (PLK1) and T cell tyrosine kinase (TTK), are part of the gene set that is up-regulated in ATC. Adoptive overexpression of p53, p21 (CIP1/WAF1), and E2F4 down-regulated transcription from the PLK1 and TTK promoters in ATC cells, suggesting that these genes might be under the negative control of tumor suppressors of the p53 and pRB families. ATC, but not normal thyroid, cells depended on PLK1 for survival. RNAi-mediated PLK1 knockdown caused cell cycle arrest associated with 4N DNA content and massive mitotic cell death. Thus, thyroid cell anaplastic transformation is accompanied by the overexpression of a cell proliferation/genetic instability-related gene cluster that includes PLK1 kinase, which is a potential molecular target for ATC treatment.

Optimization of oligonucleotide microarray fabricated by spotting 65-mer.

DNA microarrays currently provide measurements of sufficiently high quality to allow a wide variety of sound inferences about gene regulation and the coordination of cellular processes to be drawn. Nonetheless, a desire for greater precision in the measurements continues to drive the microarray research community to seek higher measurement quality through improvements in array fabrication and sample labeling and hybridization. We prepared oligonucleotide microarrays by printing 65-mer on aldehyde functional group-derivatized slides as described in a previous study. We could improve the reliability of data by removing enzymatic bias during probe labeling and hybridizing under a more stringent condition. This optimized method was used to profile gene expression patterns for nine different mouse tissues and organs, and multidimensional scaling (MDS) analysis of data showed both strong similarity between like samples and a clear, highly reproducible separation between different tissue samples. Three other microarrays were fabricated on commercial substrates and hybridized following the manufacturer's instructions. The data were then compared with in-house microarray data and reverse transcription-polymerase chain reaction (RT-PCR) data. The microarray printed on the custom aldehyde slide was superior to microarrays printed on commercially available substrate slides in terms of signal intensities, background, and hybridization characteristics. The data from the custom substrate microarray generally showed good agreement in quantitative changes up to 100-fold changes of transcript abundance with RT-PCR data. However, more accurate comparisons will be made as more genomic sequence information is gathered in the public data domain.

Performance characteristics of 65-mer oligonucleotide microarrays.

Microarray fabrication using presynthesized long oligonucleotide is becoming increasingly important, but a study of large-scale array productions has not yet been published. We addressed the issue of fabricating oligonucleotide microarrays by spotting commercial presynthesized 65-mers with 5' amines representing 7500 murine genes. Amine-modified oligonucleotides were immobilized on glass slides having aldehyde groups via transient Schiff base formation followed by reduction to produce a covalent conjugate. When RNA derived from the same source was used for Cy3 and Cy5 labeling and hybridized to the same array, signal intensities spanning three orders of magnitude were observed and the coefficient of variance (CV) between the two channels for all spots was 8 to 10%. To ascertain the reproducibility of ratio determination of these arrays, two triplicate hybridizations (with fluorochrome reversal) comparing RNAs from a fibroblast (NIH3T3) and a breast cancer (JC) cell line were carried out. The 95% confidence interval for all spots in the six hybridizations was 0.60 to 1.66. This level of reproducibility allows use of the full range of pattern finding and discriminant analysis typically applied to complementary DNA (cDNA) microarrays. Further comparative testing was carried out with oligonucleotide microarrays, cDNA microarrays, and reverse transcription (RT)-PCR assays to examine the comparability of results across these different methodologies.