Modeling Therapeutic Antibody-Small Molecule Drug-Drug Interactions Using a Three-Dimensional Perfusable Human Liver Coculture Platform.

Traditional in vitro human liver cell culture models lose key hepatic functions such as metabolic activity during short-term culture. Advanced three-dimensional (3D) liver coculture platforms offer the potential for extended hepatocyte functionality and allow for the study of more complex biologic interactions, which can improve and refine human drug safety evaluations. Here, we use a perfusion flow 3D microreactor platform for the coculture of cryopreserved primary human hepatocytes and Kupffer cells to study the regulation of cytochrome P450 3A4 isoform (CYP3A4) activity by chronic interleukin 6 (IL-6)-mediated inflammation over 2 weeks. Hepatocyte cultures remained stable over 2 weeks, with consistent albumin production and basal IL-6 levels. Direct IL-6 stimulation that mimics an inflammatory state induced a dose-dependent suppression of CYP3A4 activity, an increase in C-reactive protein (CRP) secretion, and a decrease in shed soluble interleukin-6 receptor (IL-6R) levels, indicating expected hepatic IL-6 bioactivity. Tocilizumab, an anti-IL-6R monoclonal antibody used to treat rheumatoid arthritis, has been demonstrated clinically to impact small molecule drug pharmacokinetics by modulating cytochrome P450 enzyme activities, an effect not observed in traditional hepatic cultures. We have now recapitulated the clinical observation in a 3D bioreactor system. Tocilizumab was shown to desuppress CYP3A4 activity while reducing the CRP concentration after 72 hours in the continued presence of IL-6. This change in CYP3A4 activity decreased the half-life and area under the curve up to the last measurable concentration (AUClast) of the small molecule CYP3A4 substrate simvastatin hydroxy acid, measured before and after tocilizumab treatment. We conclude that next-generation in vitro liver culture platforms are well suited for these types of long-term treatment studies and show promise for improved drug safety assessment.

Time-dependent cellular and transcriptional changes in the osteoblast lineage associated with sclerostin antibody treatment in ovariectomized rats.

Inhibition of sclerostin with sclerostin antibody (Scl-Ab) has been shown to stimulate bone formation, decrease bone resorption, and increase bone mass in both animals and humans. To obtain insight into the temporal cellular and transcriptional changes in the osteoblast (OB) lineage associated with long-term Scl-Ab treatment, stereological and transcriptional analyses of the OB lineage were performed on lumbar vertebrae from aged ovariectomized rats. Animals were administered Scl-Ab 3 or 50mg/kg/wk or vehicle (VEH) for up to 26weeks (d183), followed by a treatment-free period (TFP). At 50mg/kg/wk, bone volume (BV/total volume [TV]) increased through d183 and declined during the TFP. Bone formation rate (BFR/bone surface [BS]) and total OB number increased through d29, then progressively declined, coincident with a decrease in total osteoprogenitor (OP) numbers from d29 through d183. Analysis of differentially expressed genes (DEGs) from microarray analysis of mRNA isolated from laser capture microdissection samples enriched for OB, lining cells, and osteocytes (OCy) revealed modules of genes that correlated with BFR/BS, BV/TV, and osteoblastic surface (Ob.S)/BS. Expression change of canonical Wnt target genes was similar in all three cell types at d8, including upregulation of Twist1 and Wisp1. At d29, the pattern of Wnt target gene expression changed in the OCy, with Twist1 returning to VEH level, sustained upregulation of Wisp1, and upregulation of several other Wnt targets that continued into the TFP. Predicted activation of pathways recognized to integrate with and regulate canonical Wnt signaling were also activated at d29 in the OCy. The most significantly affected pathways represented transcription factor signaling known to inhibit cell cycle progression (notably p53) and mitogenesis (notably c-Myc). These changes occurred at the time of peak BFR/BS and continued as BFR/BS declined during treatment, then trended toward VEH level in the TFP. Concurrent with this transcriptional switch was a reduction in OP numbers, an effect that would ultimately limit bone formation. This study confirms that the initial transcriptional response in response to Scl-Ab is activation of canonical Wnt signaling and the data demonstrate that there is induction of additional regulatory pathways in OCy with long-term treatment. The interactions between Wnt and p53/c-Myc signaling may be key in limiting OP populations, thus contributing to self-regulation of bone formation with continued Scl-Ab administration.

Utilization of human nuclear receptors as an early counter screen for off-target activity: a case study with a compendium of 615 known drugs.

Off-target effects of drugs on nuclear hormone receptors (NHRs) may result in adverse effects in multiple organs/physiological processes. Reliable assessments of the NHR activities for drug candidates are therefore crucial for drug development. However, the highly permissive structures of NHRs for vastly different ligands make it challenging to predict interactions by examining the chemical structures of the ligands. Here, we report a detailed investigation on the agonistic and antagonistic activities of 615 known drugs or drug candidates against a panel of 6 NHRs: androgen, progesterone, estrogen α/ß, and thyroid hormone α/ß receptors. Our study revealed that 4.7 and 12.4% compounds have agonistic and antagonistic activities, respectively, against this panel of NHRs. Nonetheless, potent, unintended NHR hits are relatively rare among the known drugs, indicating that such interactions are perhaps not tolerated during drug development. However, we uncovered examples of compounds that unintentionally agonize or antagonize NHRs. In addition, a number of compounds showed multi-NHR activities, suggesting that the cross-talk between multiple NHRs co-operate to elicit in vivo effects. These data highlight the merits of counter screening drug candidate against NHRs during drug discovery/development.

Retinal Toxicity Induced by a Novel ß-secretase Inhibitor in the Sprague-Dawley Rat.

ß-Secretase 1 (BACE1) represents an attractive target for the treatment of Alzheimer's disease. In the course of development of a novel small molecule BACE1 inhibitor (AMG-8718), retinal thinning was observed in a 1-month toxicity study in the rat. To further understand the lesion, an investigational study was conducted whereby rats were treated daily with AMG-8718 for 1 month followed by a 2-month treatment-free phase. The earliest detectable change in the retina was an increase in autofluorescent granules in the retinal pigment epithelium (RPE) on day 5; however, there were no treatment-related light microscopic changes observed in the neuroretina and no changes observed by fundus autofluorescence or routine ophthalmoscopic examination after 28 days of dosing. Following 2 months of recovery, there was significant retinal thinning attributed to loss of photoreceptor nuclei from the outer nuclear layer. Electroretinographic changes were observed as early as day 14, before any microscopic evidence of photoreceptor loss. BACE1 knockout rats were generated and found to have normal retinal morphology indicating that the retinal toxicity induced by AMG-8718 was likely off-target. These results suggest that AMG-8718 impairs phagolysosomal function in the rat RPE, which leads to photoreceptor dysfunction and ultimately loss of photoreceptors.

A systematic assessment of mitochondrial function identified novel signatures for drug-induced mitochondrial disruption in cells.

Mitochondrial perturbation has been recognized as a contributing factor to various drug-induced organ toxicities. To address this issue, we developed a high-throughput flow cytometry-based mitochondrial signaling assay to systematically investigate mitochondrial/cellular parameters known to be directly impacted by mitochondrial dysfunction: mitochondrial membrane potential (MMP), mitochondrial reactive oxygen species (ROS), intracellular reduced glutathione (GSH) level, and cell viability. Modulation of these parameters by a training set of compounds, comprised of established mitochondrial poisons and 60 marketed drugs (30 nM to 1mM), was tested in HL-60 cells (a human pro-myelocytic leukemia cell line) cultured in either glucose-supplemented (GSM) or glucose-free (containing galactose/glutamine; GFM) RPMI-1640 media. Post-hoc bio-informatic analyses of IC50 or EC50 values for all parameters tested revealed that MMP depolarization in HL-60 cells cultured in GSM was the most reliable parameter for determining mitochondrial dysfunction in these cells. Disruptors of mitochondrial function depolarized MMP at concentrations lower than those that caused loss of cell viability, especially in cells cultured in GSM; cellular GSH levels correlated more closely to loss of viability in vitro. Some mitochondrial respiratory chain inhibitors increased mitochondrial ROS generation; however, measuring an increase in ROS alone was not sufficient to identify mitochondrial disruptors. Furthermore, hierarchical cluster analysis of all measured parameters provided confirmation that MMP depletion, without loss of cell viability, was the key signature for identifying mitochondrial disruptors. Subsequent classification of compounds based on ratios of IC50s of cell viability:MMP determined that this parameter is the most critical indicator of mitochondrial health in cells and provides a powerful tool to predict whether novel small molecule entities possess this liability.

Common handling procedures conducted in preclinical safety studies result in minimal hepatic gene expression changes in Sprague-Dawley rats.

Gene expression profiling is a tool to gain mechanistic understanding of adverse effects in response to compound exposure. However, little is known about how the common handling procedures of experimental animals during a preclinical study alter baseline gene expression. We report gene expression changes in the livers of female Sprague-Dawley rats following common handling procedures. Baseline gene expression changes identified in this study provide insight on how these changes may affect interpretation of gene expression profiles following compound exposure. Rats were divided into three groups. One group was not subjected to handling procedures and served as controls for both handled groups. Animals in the other two groups were weighed, subjected to restraint in Broome restrainers, and administered water via oral gavage daily for 1 or 4 days with tail vein blood collections at 1, 2, 4, and 8 hours postdose on days 1 and 4. Significantly altered genes were identified in livers of animals following 1 or 4 days of handling when compared to the unhandled animals. Gene changes in animals handled for 4 days were similar to those handled for 1 day, suggesting a lack of habituation. The altered genes were primarily immune function related genes. These findings, along with a correlating increase in corticosterone levels suggest that common handling procedures may cause a minor immune system perturbance.

A multifactorial approach to hepatobiliary transporter assessment enables improved therapeutic compound development.

The bile salt export pump (BSEP) is expressed at the canalicular domain of hepatocytes, where it serves as the primary route of elimination for monovalent bile acids (BAs) into the bile canaliculi. The most compelling evidence linking dysfunction in BA transport with liver injury in humans is found with carriers of mutations that render BSEP nonfunctional. Based on mounting evidence, there appears to be a strong association between drug-induced BSEP interference and liver injury in humans; however, causality has not been established. For this reason, drug-induced BSEP interference is best considered a susceptibility factor for liver injury as other host- or drug-related properties may contribute to the development of hepatotoxicity. To better understand the association between BSEP interference and liver injury in humans, over 600 marketed or withdrawn drugs were evaluated in BSEP expressing membrane vesicles. The example of a compound that failed during phase 1 human trials is also described, AMG 009. AMG 009 showed evidence of liver injury in humans that was not predicted by preclinical safety studies, and BSEP inhibition was implicated. For 109 of the drugs with some effect on in vitro BSEP function, clinical use, associations with hepatotoxicity, pharmacokinetic data, and other information were annotated. A steady state concentration (C(ss)) for each of these annotated drugs was estimated, and a ratio between this value and measured IC50 potency values were calculated in an attempt to relate exposure to in vitro potencies. When factoring for exposure, 95% of the annotated compounds with a C(ss)/BSEP IC50 ratio ≥ 0.1 were associated with some form of liver injury. We then investigated the relationship between clinical evidence of liver injury and effects to multidrug resistance-associated proteins (MRPs) believed to play a role in BA homeostasis. The effect of 600+ drugs on MRP2, MRP3, and MRP4 function was also evaluated in membrane vesicle assays. Drugs with a C(ss)/BSEP IC50 ratio ≥ 0.1 and a C(ss)/MRP IC50 ratio ≥ 0.1 had almost a 100% correlation with some evidence of liver injury in humans. These data suggest that integration of exposure data, and knowledge of an effect to not only BSEP but also one or more of the MRPs, is a useful tool for informing the potential for liver injury due to altered BA transport.

Raf inhibition causes extensive multiple tissue hyperplasia and urinary bladder neoplasia in the rat.

Seven novel and potent Raf small molecule kinase inhibitors (C1-7) were evaluated in seven-day oral repeat dose rat toxicity studies. All compounds tested induced hyperplasia in multiple tissues. Consistently affected was stratified squamous epithelium at a number of sites and transitional epithelium of urinary bladder and kidney. A seven-day time course study in rats showed morphologic evidence of epithelial proliferation in the nonglandular stomach within four to five hours after a single dose of C-1. Similar indications of cellular proliferation were observed in the urinary bladder by day 2 and in the heart, kidney, and liver by day 3. Transcriptional evidence of proliferation in the urinary bladder was detected within four to five hours after a single dose consistent with activation of the PI3K/AKT and ERK/MAPK pathways. In a twenty-eight-day rat toxicity study of C-1, hyperplasia was observed in the esophagus, nonglandular stomach, skin, urinary bladder, kidney, and heart. Hyperplasia of transitional epithelium of the urinary bladder was particularly severe and in one female rat was accompanied by the presence of a transitional cell carcinoma. These results suggest that these Raf inhibitors induce early transcriptional changes driving unchecked cell proliferation, resulting in marked tissue hyperplasia that can progress to carcinoma within a short time frame.

The evolution of bioinformatics in toxicology: advancing toxicogenomics.

As one reflects back through the past 50 years of scientific research, a significant accomplishment was the advance into the genomic era. Basic research scientists have uncovered the genetic code and the foundation of the most fundamental building blocks for the molecular activity that supports biological structure and function. Accompanying these structural and functional discoveries is the advance of techniques and technologies to probe molecular events, in time, across environmental and chemical exposures, within individuals, and across species. The field of toxicology has kept pace with advances in molecular study, and the past 50 years recognizes significant growth and explosive understanding of the impact of the compounds and environment to basic cellular and molecular machinery. The advancement of molecular techniques applied in a whole-genomic capacity to the study of toxicant effects, toxicogenomics, is no doubt a significant milestone for toxicological research. Toxicogenomics has also provided an avenue for advancing a joining of multidisciplinary sciences including engineering and informatics in traditional toxicological research. This review will cover the evolution of the field of toxicogenomics in the context of informatics integration its current promise, and limitations.

Interference with bile salt export pump function is a susceptibility factor for human liver injury in drug development.

The bile salt export pump (BSEP) is an efflux transporter, driving the elimination of endobiotic and xenobiotic substrates from hepatocytes into the bile. More specifically, it is responsible for the elimination of monovalent, conjugated bile salts, with little or no assistance from other apical transporters. Disruption of BSEP activity through genetic disorders is known to manifest in clinical liver injury such as progressive familial intrahepatic cholestasis type 2. Drug-induced disruption of BSEP is hypothesized to play a role in the development of liver injury for several marketed or withdrawn therapeutics. Unfortunately, preclinical animal models have been poor predictors of the liver injury associated with BSEP interference observed for humans, possibly because of interspecies differences in bile acid composition, differences in hepatobiliary transporter modulation or constitutive expression, as well as other mechanisms. Thus, a BSEP-mediated liver liability may go undetected until the later stages of drug development, such as during clinical trials or even postlicensing. In the absence of a relevant preclinical test system for BSEP-mediated liver injury, the toxicological relevance of available in vitro models to human health rely on the use of benchmark compounds with known clinical outcomes, such as marketed or withdrawn drugs. In this study, membrane vesicles harvested from BSEP-transfected insect cells were used to assess the activity of more than 200 benchmark compounds to thoroughly investigate the relationship between interference with BSEP function and liver injury. The data suggest a relatively strong association between the pharmacological interference with BSEP function and human hepatotoxicity. Although the most accurate translation of risk would incorporate pharmacological potency, pharmacokinetics, clearance mechanisms, tissue distribution, physicochemical properties, indication, and other drug attributes, the additional understanding of a compound's potency for BSEP interference should help to limit or avoid BSEP-related liver liabilities in humans that are not often detected by standard preclinical animal models.

Application of genomics for identification of systemic toxicity triggers associated with VEGF-R inhibitors.

The key to the discovery of new pharmaceuticals is to develop molecules that interact with the intended target and minimize interaction with unintended molecular targets, therefore minimizing toxicity. This is aided by the use of various in vitro selectivity assays that are used to select agents most potent for the desired target. Typically, molecules from similar chemical series, with similar in vitro potencies, are expected to yield comparable in vivo pharmacological and toxicological profiles, predictive of target effects. However, in this study, we investigated the in vivo effects of two analogue compounds that similarly inhibit several receptor tyrosine kinases such as vascular endothelial growth factor receptor 1 (VEGFR/Flt1), vascular endothelial growth factor 2 (VEGFR2/kinase domain receptor/Flk-1), vascular endothelial growth factor receptor 3 (VEGFR3/Flt4), platelet-derived growth factor receptor (PDGFR), and Kit receptors, which bear similar chemical structures, have comparable potencies, but differ markedly in their rodent toxicity profiles. Global gene expression data were used to generate hypotheses regarding the existence of toxicity triggers that would reflect the perturbation of signaling in multiple organs such as the liver, adrenal glands, and the pancreas in response to compound treatment. We concluded that differences in pharmacokinetic properties of the two analogues, such as volume of distribution, half-life, and organ concentrations, resulted in marked differences in the chemical burden on target organs and may have contributed to the vast differences in toxicity profiles observed with the two otherwise similar molecules. We propose including select toxicokinetic parameters such as V(ss), T(1/2), and T(max) as additional criteria that could be used to rank order compounds from the same pharmacological series to possibly minimize organ toxicity. Assessment of toxicokinetics is not an atypical activity on toxicology studies, even in early screening studies; however, these data may not always be used in decision making for selecting or eliminating one compound over another. Finally, we illustrate that in vivo gene expression profiles can serve as a complementary assessor of this activity and simultaneously help provide an assessment of on or off-target biological activity.

Current and future applications of toxicogenomics: Results summary of a survey from the HESI Genomics State of Science Subcommittee.

BACKGROUND: In spite of the application of toxicogenomic (TGx) data to the field of toxicology for the past 10 years, the broad implementation and full impact of TGx for chemical and drug evaluation to improve decision making within organizations and by policy makers has not been achieved. OBJECTIVES: The goal of the Health and Environmental Sciences Institute (HESI) Committee on the Application of Genomics to Mechanism-based Risk Assessment was to construct and summarize a multisector survey, addressing key issues and perspectives on the current and future practical uses and challenges of implementing TGx data to facilitate discussions for decision making within organizations and by policy makers. METHODS: An online survey to probe the current status and future challenges facing the field of TGx for drug and chemical evaluation in experimental and nonclinical models was taken by scientists and scientific decision/policy makers actively engaged in the field of TGx within industrial, academic, and regulatory sectors of the United States, Europe, and Japan. For this survey, TGx refers specifically to the analysis of gene expression responses to evaluate xenobiotic exposure in experimental and preclinical models. RESULTS: The survey results are summarized from questions covering broad areas including technology used, organizational capacity and resource allocation, experimental approaches, data storage and exchange, perceptions of benefits and hurdles, and future expectations. CONCLUSIONS: The survey findings provide valuable information on the current state of the science of TGx applications and identify key areas in which TGx will have an impact as well as the key hurdles in applying TGx data to address issues. The findings serve as a public resource to facilitate discussions on the focus of future TGx efforts to ensure that a maximal benefit can be obtained from toxicogenomic studies.

A novel statistical algorithm for gene expression analysis helps differentiate pregnane X receptor-dependent and independent mechanisms of toxicity.

Genome-wide gene expression profiling has become standard for assessing potential liabilities as well as for elucidating mechanisms of toxicity of drug candidates under development. Analysis of microarray data is often challenging due to the lack of a statistical model that is amenable to biological variation in a small number of samples. Here we present a novel non-parametric algorithm that requires minimal assumptions about the data distribution. Our method for determining differential expression consists of two steps: 1) We apply a nominal threshold on fold change and platform p-value to designate whether a gene is differentially expressed in each treated and control sample relative to the averaged control pool, and 2) We compared the number of samples satisfying criteria in step 1 between the treated and control groups to estimate the statistical significance based on a null distribution established by sample permutations. The method captures group effect without being too sensitive to anomalies as it allows tolerance for potential non-responders in the treatment group and outliers in the control group. Performance and results of this method were compared with the Significant Analysis of Microarrays (SAM) method. These two methods were applied to investigate hepatic transcriptional responses of wild-type (PXR(+/+)) and pregnane X receptor-knockout (PXR(-/-)) mice after 96 h exposure to CMP013, an inhibitor of ß-secretase (ß-site of amyloid precursor protein cleaving enzyme 1 or BACE1). Our results showed that CMP013 led to transcriptional changes in hallmark PXR-regulated genes and induced a cascade of gene expression changes that explained the hepatomegaly observed only in PXR(+/+) animals. Comparison of concordant expression changes between PXR(+/+) and PXR(-/-) mice also suggested a PXR-independent association between CMP013 and perturbations to cellular stress, lipid metabolism, and biliary transport.

Assessment of repeated microarray experiments using mixed tissue RNA reference samples.

Genome-scale gene expression technologies are increasingly being applied for biological research as a whole and toxicological screening in particular. In order to monitor data quality and process drift, we adopted the use of two rat-tissue mixtures (brain, liver, kidney, and testis) previously introduced as RNA reference samples. These samples were processed every time a microarray experiment was hybridized, thereby verifying the comparability of the resulting expression data for cross-study comparison. This study presents the analysis of 21 technical replicates of these two mixed-tissue samples using Affymetrix RAE230_2 GeneChip over a period of 12 months. The results show that detection sensitivity, measured by the number of present and absent sequences, is robust, and data correlation, indicated by scatter plots, varies little over time. Receiver operating characteristic (ROC) curves show the sensitivity and specificity of the current measurements are consistent with arrays previously classified as well performing. Overall, this paper shows that the inclusion of standard samples during microarray labeling and hybridization experiments is useful to benchmark the performance of microarray experiments over time and allows discovery of any process drift that, if it occurs, may confound the comparison of these datasets.

Insight into redox-regulated gene networks in vascular cells.

To understand the complex nature of the atherogenic response initiated by oxidative stress in vascular smooth muscle cells (vSMCs), computational prediction methodology was employed to define putative gene-gene and gene-environment interactions in vSMCs subjected to oxidative chemical stress. Computational relationships were derived from the global gene expression profiles of murine cells challenged with a chemical pro-oxidant to cause oxidative stress or cells treated with anti-oxidant prior to oxidative injury. Target clones were chosen based on their biological relevance within the context of the atherogenic response and included lysyl oxidase, matrix metalloproteinase 2, insulin like growth factor binding protein 5, and lymphocyte antigen 6c. Established biological relationships were derived computationally confirming the usefulness of the algorithm in uncovering novel biological relationships worthy of future investigation. Thus, the predictive algorithm can be a useful tool to advance the frontiers of biological discovery.

Use of a mixed tissue RNA design for performance assessments on multiple microarray formats.

The comparability and reliability of data generated using microarray technology would be enhanced by use of a common set of standards that allow accuracy, reproducibility and dynamic range assessments on multiple formats. We designed and tested a complex biological reagent for performance measurements on three commercial oligonucleotide array formats that differ in probe design and signal measurement methodology. The reagent is a set of two mixtures with different proportions of RNA for each of four rat tissues (brain, liver, kidney and testes). The design provides four known ratio measurements of >200 reference probes, which were chosen for their tissue-selectivity, dynamic range coverage and alignment to the same exemplar transcript sequence across all three platforms. The data generated from testing three biological replicates of the reagent at eight laboratories on three array formats provides a benchmark set for both laboratory and data processing performance assessments. Close agreement with target ratios adjusted for sample complexity was achieved on all platforms and low variance was observed among platforms, replicates and sites. The mixed tissue design produces a reagent with known gene expression changes within a complex sample and can serve as a paradigm for performance standards for microarrays that target other species.

Prediction of nephrotoxicant action and identification of candidate toxicity-related biomarkers.

A vast majority of pharmacological compounds and their metabolites are excreted via the urine, and within the complex structure of the kidney,the proximal tubules are a main target site of nephrotoxic compounds. We used the model nephrotoxicants mercuric chloride, 2-bromoethylamine hydrobromide, hexachlorobutadiene, mitomycin, amphotericin, and puromycin to elucidate time- and dose-dependent global gene expression changes associated with proximal tubular toxicity. Male Sprague-Dawley rats were dosed via intraperitoneal injection once daily for mercuric chloride and amphotericin (up to 7 doses), while a single dose was given for all other compounds. Animals were exposed to 2 different doses of these compounds and kidney tissues were collected on day 1, 3, and 7 postdosing. Gene expression profiles were generated from kidney RNA using 17K rat cDNA dual dye microarray and analyzed in conjunction with histopathology. Analysis of gene expression profiles showed that the profiles clustered based on similarities in the severity and type of pathology of individual animals. Further, the expression changes were indicative of tubular toxicity showing hallmarks of tubular degeneration/regeneration and necrosis. Use of gene expression data in predicting the type of nephrotoxicity was then tested with a support vector machine (SVM)-based approach. A SVM prediction module was trained using 120 profiles of total profiles divided into four classes based on the severity of pathology and clustering. Although mitomycin C and amphotericin B treatments did not cause toxicity, their expression profiles were included in the SVM prediction module to increase the sample size. Using this classifier, the SVM predicted the type of pathology of 28 test profiles with 100% selectivity and 82% sensitivity. These data indicate that valid predictions could be made based on gene expression changes from a small set of expression profiles. A set of potential biomarkers showing a time- and dose-response with respect to the progression of proximal tubular toxicity were identified. These include several transporters (Slc21a2, Slc15, Slc34a2), Kim 1, IGFbp-1, osteopontin, alpha-fibrinogen, and Gstalpha.

Integrating phenotypic and expression profiles to map arsenic-response networks.

BACKGROUND: Arsenic is a nonmutagenic carcinogen affecting millions of people. The cellular impact of this metalloid in Saccharomyces cerevisiae was determined by profiling global gene expression and sensitivity phenotypes. These data were then mapped to a metabolic network composed of all known biochemical reactions in yeast, as well as the yeast network of 20,985 protein-protein/protein-DNA interactions. RESULTS: While the expression data unveiled no significant nodes in the metabolic network, the regulatory network revealed several important nodes as centers of arsenic-induced activity. The highest-scoring proteins included Fhl1, Msn2, Msn4, Yap1, Cad1 (Yap2), Pre1, Hsf1 and Met31. Contrary to the gene-expression analyses, the phenotypic-profiling data mapped to the metabolic network. The two significant metabolic networks unveiled were shikimate, and serine, threonine and glutamate biosynthesis. We also carried out transcriptional profiling of specific deletion strains, confirming that the transcription factors Yap1, Arr1 (Yap8), and Rpn4 strongly mediate the cell's adaptation to arsenic-induced stress but that Cad1 has negligible impact. CONCLUSIONS: By integrating phenotypic and transcriptional profiling and mapping the data onto the metabolic and regulatory networks, we have shown that arsenic is likely to channel sulfur into glutathione for detoxification, leads to indirect oxidative stress by depleting glutathione pools, and alters protein turnover via arsenation of sulfhydryl groups on proteins. Furthermore, we show that phenotypically sensitive pathways are upstream of differentially expressed ones, indicating that transcriptional and phenotypic profiling implicate distinct, but related, pathways.

Exploration of low-dose estrogen effects: identification of No Observed Transcriptional Effect Level (NOTEL).

Identifying a minimal dose capable of eliciting a biological response is a fundamental issue in a number of scientific fields, including: drug development, signal transduction research, and environmental toxicology. Frequently, proliferation, viability, and other assays based on the cellular response to a treatment are used to assess the threshold dose for minimal activity. Here we propose a novel approach for identifying the effects of low dose treatments and pinpointing the threshold dose. Using microarrays, we examined the transcriptional response of a hormone responsive breast cancer cell line (MCF-7) stimulated with various concentrations of estrogen. Previous studies have focused on transcriptional responses to physiologically relevant concentrations of estrogen. However, relatively few studies have examined the transcriptional effects of concentrations below normal physiologic levels. These doses may not stimulate the expression of any genes or, alternatively, may regulate a different subset of genes that had not been previously characterized as estrogen responsive. We used gene expression profiling, coupled with a detailed analysis of replicates, to measure estrogen effects on many transcriptional targets and found that only physiologically relevant doses of estrogen (1 x 10(-10) M and higher) were capable of inducing a transcriptional response. This study demonstrates the utility of gene expression profiling as a means to identify concentrations that do not elicit a change in gene expression, or simply a No Observed Transcriptional Effect Level (NOTEL). The identification of a NOTEL for a given compound may be beneficial in several different scientific disciplines. For example, in the development of therapeutic drugs, a NOTEL could be used to identify doses of pharmaceutical compounds that are no longer effective at modulating the expression of biomarkers of efficacy.