Carrie L. Moland

QSAR Models Using a Large Diverse Set of Estrogens

Endocrine disruptors (EDs) have a variety of adverse effects in humans and animals. About 58,000 chemicals, most having little safety data, must be tested in a group of tiered assays. As assays will take years, it is important to develop rapid methods to help in priority setting. For application to large data sets, we have developed an integrated system that contains sequential four phases to predict the ability of chemicals to bind to the estrogen receptor (ER), a prevalent mechanism for estrogenic EDs. Here we report the results of evaluating two types of QSAR models for inclusion in phase III to quantitatively predict chemical binding to the ER. Our data set for the relative binding affinities (RBAs) to the ER consists of 130 chemicals covering a wide range of structural diversity and a 6 orders of magnitude spread of RBAs. CoMFA and HQSAR models were constructed and compared for performance. The CoMFA model had a r2 = 0.91 and a q2LOO = 0.66. HQSAR showed reduced performance compared to CoMFA with r2 = 0.76 and q2LOO = 0.59. A number of parameters were examined to improve the CoMFA model. Of these, a phenol indicator increased the q2LOO to 0.71. When up to 50% of the chemicals were left out in the leave-N-out cross-validation, the q2 remained significant. Finally, the models were tested by using two test sets; the q2pred for these were 0.71 and 0.62, a significant result which demonstrates the utility of the CoMFA model for predicting the RBAs of chemicals not included in the training set. If used in conjunction with phases I and II, which reduced the size of the data set dramatically by eliminating most inactive chemicals, the current CoMFA model (phase III) can be used to predict the RBA of chemicals with sufficient accuracy and to provide quantitative information for priority setting.

A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages

The rat has been used extensively as a model for evaluating chemical toxicities and for understanding drug mechanisms. However, its transcriptome across multiple organs, or developmental stages, has not yet been reported. Here we show, as part of the SEQC consortium efforts, a comprehensive rat transcriptomic BodyMap created by performing RNA-Seq on 320 samples from 11 organs of both sexes of juvenile, adolescent, adult and aged Fischer 344 rats. We catalogue the expression profiles of 40,064 genes, 65,167 transcripts, 31,909 alternatively spliced transcript variants and 2,367 non-coding genes/non-coding RNAs (ncRNAs) annotated in AceView. We find that organ-enriched, differentially expressed genes reflect the known organ-specific biological activities. A large number of transcripts show organ-specific, age-dependent or sex-specific differential expression patterns. We create a web-based, open-access rat BodyMap database of expression profiles with crosslinks to other widely used databases, anticipating that it will serve as a primary resource for biomedical research using the rat model.

Underlying mitochondrial dysfunction triggers flutamide-induced oxidative liver injury in a mouse model of idiosyncratic drug toxicity.

Flutamide, a widely used nonsteroidal anti-androgen, but not its bioisostere bicalutamide, has been associated with idiosyncratic drug-induced liver injury. Although the susceptibility factors are unknown, mitochondrial injury has emerged as a putative hazard of flutamide. To explore the role of mitochondrial sensitization in flutamide hepatotoxicity, we determined the effects of superimposed drug stress in a murine model of underlying mitochondrial abnormalities. Male wild-type or heterozygous Sod2(+/-) mice were injected intraperitoneously with flutamide (0, 30 or 100 mg/kg/day) for 28 days. A kinetic pilot study revealed that flutamide (100 mg/kg/day) caused approximately 10-fold greater exposure than the reported therapeutic mean plasma levels. Mutant (5/10), but not wild-type, mice in the high-dose group exhibited small foci of hepatocellular necrosis and an increased number of apoptotic hepatocytes. Hepatic GSSG/GSH, protein carbonyl levels, and serum lactate levels were significantly increased, suggesting oxidant stress and mitochondrial dysfunction. Measurement of mitochondrial superoxide in cultured hepatocytes demonstrated that mitochondria were a significant source of flutamide-enhanced oxidant stress. Indeed, mitochondria isolated from flutamide-treated Sod2(+/-) mice exhibited decreased aconitase activity as compared to vehicle controls. A transcriptomics analysis using MitoChips revealed that flutamide-treated Sod2(+/-) mice exhibited a selective decrease in the expression of all complexes I and III subunits encoded by mitochondrial DNA. In contrast, Sod2(+/-) mice receiving bicalutamide (50 mg/kg/day) did not reveal any hepatic changes. These results are compatible with our concept that flutamide targets hepatic mitochondria and exerts oxidant stress that can lead to overt hepatic injury in the presence of an underlying mitochondrial abnormality.

Transgenic expression of proximal tubule peroxisome proliferator–activated receptor-α in mice confers protection during acute kidney injury

Our previous studies suggest that peroxisome proliferator-activated receptor-alpha (PPARalpha) plays a critical role in regulating fatty acid beta-oxidation in kidney tissue and this directly correlated with preservation of kidney morphology and function during acute kidney injury. To further study this, we generated transgenic mice expressing PPARalpha in the proximal tubule under the control of the promoter of KAP2 (kidney androgen-regulated protein 2). Segment-specific upregulation of PPARalpha expression by testosterone treatment of female transgenic mice improved kidney function during cisplatin or ischemia-reperfusion-induced acute kidney injury. Ischemia-reperfusion injury or treatment with cisplatin in wild-type mice caused inhibition of fatty-acid oxidation, reduction of mitochondrial genes of oxidative phosphorylation, mitochondrial DNA, fatty-acid metabolism, and the tricarboxylic acid cycle. Similar injury in testosterone-treated transgenic mice resulted in amelioration of these effects. Similarly, there were increases in the levels of 4-hydroxy-2-hexenal-derived lipid peroxidation products in wild-type mice, which were also reduced in the transgenic mice. Similarly, necrosis of the S3 segment was reduced in the two injury models in transgenic mice compared to wild type. Our results suggest proximal tubule PPARalpha activity serves as a metabolic sensor. Its increased expression without the use of an exogenous PPARalpha ligand in the transgenic mice is sufficient to protect kidney function and morphology, and to prevent abnormalities in lipid metabolism associated with acute kidney injury.

The Liver Toxicity Biomarker Study: Phase I Design and Preliminary Results

Drug-induced liver injury (DILI) is the primary adverse event that results in withdrawal of drugs from the market and a frequent reason for the failure of drug candidates in development. The Liver Toxicity Biomarker Study (LTBS) is an innovative approach to investigate DILI because it compares molecular events produced in vivo by compound pairs that (a) are similar in structure and mechanism of action, (b) are associated with few or no signs of liver toxicity in preclinical studies, and (c) show marked differences in hepatotoxic potential. The LTBS is a collaborative preclinical research effort in molecular systems toxicology between the National Center for Toxicological Research and BG Medicine, Inc., and is supported by seven pharmaceutical companies and three technology providers. In phase I of the LTBS, entacapone and tolcapone were studied in rats to provide results and information that will form the foundation for the design and implementation of phase II. Molecular analysis of the rat liver and plasma samples combined with statistical analyses of the resulting datasets yielded marker analytes, illustrating the value of the broad-spectrum, molecular systems analysis approach to studying pharmacological or toxicological effects.

Elimination of laboratory ozone leads to a dramatic improvement in the reproducibility of microarray gene expression measurements

BackgroundEnvironmental ozone can rapidly degrade cyanine 5 (Cy5), a fluorescent dye commonly used in microarray gene expression studies. Cyanine 3 (Cy3) is much less affected by atmospheric ozone. Degradation of the Cy5 signal relative to the Cy3 signal in 2-color microarrays will adversely reduce the Cy5/Cy3 ratio resulting in unreliable microarray data.ResultsOzone in central Arkansas typically ranges between ~22 ppb to ~46 ppb and can be as high as 60–100 ppb depending upon season, meteorological conditions, and time of day. These levels of ozone are common in many areas of the country during the summer. A carbon filter was installed in the laboratory air handling system to reduce ozone levels in the microarray laboratory. In addition, the airflow was balanced to prevent non-filtered air from entering the laboratory. These modifications reduced the ozone within the microarray laboratory to ~2–4 ppb. Data presented here document reductions in Cy5 signal on both in-house produced microarrays and commercial microarrays as a result of exposure to unfiltered air. Comparisons of identically hybridized microarrays exposed to either carbon-filtered or unfiltered air demonstrated the protective effect of carbon-filtration on microarray data as indicated by Cy5 and Cy3 intensities. LOWESS normalization of the data was not able to completely overcome the effect of ozone-induced reduction of Cy5 signal. Experiments were also conducted to examine the effects of high humidity on microarray quality. Modest, but significant, increases in Cy5 and Cy3 signal intensities were observed after 2 or 4 hours at 98–99% humidity compared to 42% humidity.ConclusionSimple installation of carbon filters in the laboratory air handling system resulted in low and consistent ozone levels. This allowed the accurate determination of gene expression by microarray using Cy5 and Cy3 fluorescent dyes.

Development of doxorubicin-induced chronic cardiotoxicity in the B6C3F1 mouse model

Serum levels of cardiac troponins serve as biomarkers of myocardial injury. However, troponins are released into the serum only after damage to cardiac tissue has occurred. Here, we report development of a mouse model of doxorubicin (DOX)-induced chronic cardiotoxicity to aid in the identification of predictive biomarkers of early events of cardiac tissue injury. Male B6C3F(1) mice were administered intravenous DOX at 3mg/kg body weight, or an equivalent volume of saline, once a week for 4, 6, 8, 10, 12, and 14weeks, resulting in cumulative DOX doses of 12, 18, 24, 30, 36, and 42mg/kg, respectively. Mice were sacrificed a week following the last dose. A significant reduction in body weight gain was observed in mice following exposure to a weekly DOX dose for 1week and longer compared to saline-treated controls. DOX treatment also resulted in declines in red blood cell count, hemoglobin level, and hematocrit compared to saline-treated controls after the 2nd weekly dose until the 8th and 9th doses, followed by a modest recovery. All DOX-treated mice had significant elevations in cardiac troponin T concentrations in plasma compared to saline-treated controls, indicating cardiac tissue injury. Also, a dose-related increase in the severity of cardiac lesions was seen in mice exposed to 24mg/kg DOX and higher cumulative doses. Mice treated with cumulative DOX doses of 30mg/kg and higher showed a significant decline in heart rate, suggesting drug-induced cardiac dysfunction. Altogether, these findings demonstrate the development of DOX-induced chronic cardiotoxicity in B6C3F(1) mice.

Age and sex dependent changes in liver gene expression during the life cycle of the rat.

BackgroundAge- and sex-related susceptibility to adverse drug reactions and disease is a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of hepatic genes expressed at various life cycle stages will impact susceptibility to adverse drug reactions. Understanding the basal liver gene expression patterns is a necessary first step in addressing this hypothesis and will inform our assessments of adverse drug reactions as the liver plays a central role in drug metabolism and biotransformation. Untreated male and female F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 52, 78, and 104 weeks of age. Liver tissues were collected for histology and gene expression analysis. Whole-genome rat microarrays were used to query global expression profiles.ResultsAn initial list of differentially expressed genes was selected using criteria based upon p-value (p < 0.05) and fold-change (+/- 1.5). Three dimensional principal component analyses revealed differences between males and females beginning at 2 weeks with more divergent profiles beginning at 5 weeks. The greatest sex-differences were observed between 8 and 52 weeks before converging again at 104 weeks. K-means clustering identified groups of genes that displayed age-related patterns of expression. Various adult aging-related clusters represented gene pathways related to xenobiotic metabolism, DNA damage repair, and oxidative stress.ConclusionsThese results suggest an underlying role for genes in specific clusters in potentiating age- and sex-related differences in susceptibility to adverse health effects. Furthermore, such a comprehensive picture of life cycle changes in gene expression deepens our understanding and informs the utility of liver gene expression biomarkers.

Improvement in the Reproducibility and Accuracy of DNA Microarray Quantification by Optimizing Hybridization Conditions

BackgroundDNA microarrays, which have been increasingly used to monitor mRNA transcripts at a global level, can provide detailed insight into cellular processes involved in response to drugs and toxins. This is leading to new understandings of signaling networks that operate in the cell, and the molecular basis of diseases. Custom printed oligonucleotide arrays have proven to be an effective way to facilitate the applications of DNA microarray technology. A successful microarray experiment, however, involves many steps: well-designed oligonucleotide probes, printing, RNA extraction and labeling, hybridization, and imaging. Optimization is essential to generate reliable microarray data.ResultsHybridization and washing steps are crucial for a successful microarray experiment. By following the hybridization and washing conditions recommended by an oligonucleotide provider, it was found that the expression ratios were compressed greater than expected and data analysis revealed a high degree of non-specific binding. A series of experiments was conducted using rat mixed tissue RNA reference material (MTRRM) and other RNA samples to optimize the hybridization and washing conditions. The optimized hybridization and washing conditions greatly reduced the non-specific binding and improved the accuracy of spot intensity measurements.ConclusionThe results from the optimized hybridization and washing conditions greatly improved the reproducibility and accuracy of expression ratios. These experiments also suggested the importance of probe designs using better bioinformatics approaches and the need for common reference RNA samples for platform performance evaluation in order to fulfill the potential of DNA microarray technology.

Early biomarkers of doxorubicin-induced heart injury in a mouse model

Cardiac troponins, which are used as myocardial injury markers, are released in plasma only after tissue damage has occurred. Therefore, there is a need for identification of biomarkers of earlier events in cardiac injury to limit the extent of damage. To accomplish this, expression profiling of 1179 unique microRNAs (miRNAs) was performed in a chronic cardiotoxicity mouse model developed in our laboratory. Male B6C3F1 mice were injected intravenously with 3mg/kg doxorubicin (DOX; an anti-cancer drug), or saline once a week for 2, 3, 4, 6, and 8weeks, resulting in cumulative DOX doses of 6, 9, 12, 18, and 24mg/kg, respectively. Mice were euthanized a week after the last dose. Cardiac injury was evidenced in mice exposed to 18mg/kg and higher cumulative DOX dose whereas examination of hearts by light microscopy revealed cardiac lesions at 24mg/kg DOX. Also, 24 miRNAs were differentially expressed in mouse hearts, with the expression of 1, 1, 2, 8, and 21 miRNAs altered at 6, 9, 12, 18, and 24mg/kg DOX, respectively. A pro-apoptotic miR-34a was the only miRNA that was up-regulated at all cumulative DOX doses and showed a significant dose-related response. Up-regulation of miR-34a at 6mg/kg DOX may suggest apoptosis as an early molecular change in the hearts of DOX-treated mice. At 12mg/kg DOX, up-regulation of miR-34a was associated with down-regulation of hypertrophy-related miR-150; changes observed before cardiac injury. These findings may lead to the development of biomarkers of earlier events in DOX-induced cardiotoxicity that occur before the release of cardiac troponins.