Structure Guided Discovery of Novel Pan Metallo-ß-Lactamase Inhibitors with Improved Gram-Negative Bacterial Cell Penetration.

The use of ß-lactam (BL) and ß-lactamase inhibitor combination to overcome BL antibiotic resistance has been validated through clinically approved drug products. However, unmet medical needs still exist for the treatment of infections caused by Gram-negative (GN) bacteria expressing metallo-ß-lactamases. Previously, we reported our effort to discover pan inhibitors of three main families in this class: IMP, VIM, and NDM. Herein, we describe our work to improve the GN coverage spectrum in combination with imipenem and relebactam. This was achieved through structure- and property-based optimization to tackle the GN cell penetration and efflux challenges. A significant discovery was made that inhibition of both VIM alleles, VIM-1 and VIM-2, is essential for broad GN coverage, especially against VIM-producing P. aeruginosa. In addition, pharmacokinetics and nonclinical safety profiles were investigated for select compounds. Key findings from this drug discovery campaign laid the foundation for further lead optimization toward identification of preclinical candidates.

Nonclinical Cardiovascular Assessment of the Soluble Guanylate Cyclase Stimulator Vericiguat.

Vericiguat and its metabolite M-1 were assessed for proarrhythmic risk in nonclinical in vitro and in vivo studies. In vitro manual voltage-clamp recordings at room temperature determined the effect of vericiguat on human Ether-a-go-go Related Gene (hERG) K+ channels. Effects of vericiguat and M-1 on hERG K+, Nav1.5, hCav1.2, hKvLQT1/1minK, and hKv4.3 channels were investigated via automated voltage-clamp recordings at ambient temperature. Effects of vericiguat and M-1 on hERG K+ and Nav1.5 channels at pathophysiological conditions were explored via manual voltage-clamp recordings at physiologic temperature. Single oral doses of vericiguat (0.6, 2.0, and 6.0 mg/kg) were assessed for in vivo proarrhythmic risk via administration to conscious telemetered dogs; electrocardiogram (ECG) and hemodynamic parameters were monitored. ECG recordings were included in 4- and 39-week dog toxicity studies. In manual voltage-clamp recordings, vericiguat inhibited hERG K+-mediated tail currents in a concentration-dependent manner (20% threshold inhibitory concentration ∼1.9 µM). In automated voltage-clamp recordings, neither vericiguat nor M-1 were associated with biologically relevant inhibition (>20%) of hNav1.5, hCav1.2, hKvLQT1, and hKv4.3. No clinically relevant observations were made for hNav1.5 and hKvLQT1 under simulated pathophysiological conditions. Vericiguat was associated with expected mode-of-action-related dose-dependent changes in systolic arterial blood pressure (up to -20%) and heart rate (up to +53%). At maximum vericiguat dose, corrected QT (QTc) interval changes from baseline varied slightly (-6 to +1%) depending on correction formula. Toxicity studies confirmed absence of significant QTc interval changes. There was no evidence of an increased proarrhythmic risk from nonclinical studies with vericiguat or M-1. SIGNIFICANCE STATEMENT: There was no evidence of an increased proarrhythmic risk from in vitro and in vivo nonclinical studies with vericiguat or M-1. The integrated risk assessment of these nonclinical data combined with existing clinical data demonstrate administration of vericiguat 10 mg once daily in patients with heart failure with reduced ejection fraction is not associated with a proarrhythmic risk.

Early-Onset albuminuria and Associated Renal Pathology in Leucine-Rich Repeat Kinase 2 Knockout Rats.

Activating mutations of the leucine-rich repeat kinase 2 (LRRK2) gene are associated with Parkinson disease (PD), prompting development of LRRK2 inhibitors as potential treatment for PD. However, kidney safety concerns have surfaced from LRRK2 knockout (KO) mice and rats and from repeat-dose studies in rodents administered LRRK2 inhibitors. To support drug development of this therapeutic target, we conducted a study of 26 weeks' duration in 2-month-old wild-type and LRRK2 KO Long-Evans Hooded rats to systematically examine the performance of urinary safety biomarkers and to characterize the nature of the morphological changes in the kidneys by light microscopy and by ultrastructural evaluation. Our data reveal the time course of early-onset albuminuria at 3 and 4 months in LRRK2 KO female and male rats, respectively. The increases in urine albumin were not accompanied by concurrent increases in serum creatinine, blood urea nitrogen, or renal safety biomarkers such as kidney injury molecule 1 or clusterin, although morphological alterations in both glomerular and tubular structure were identified by light and transmission electron microscopy at 8 months of age. Diet optimization with controlled food intake attenuated the progression of albuminuria and associated renal changes.

Lead Optimization to Advance Protease-Activated Receptor-1 Antagonists in Early Discovery.

Vorapaxar is an approved drug for the reduction of thrombotic cardiovascular events in patients with a history of myocardial infarction or with peripheral arterial disease. Subsequent to the discovery of Vorapaxar, medicinal chemistry efforts were continued to identify structurally differentiated leads. Toward this goal, extensive structure-activity relationship studies using a C-ring-truncated version of Vorapaxar culminated in the discovery of three leads, represented as 13, 14, and 23. Among these leads, compound 14 possessed favorable pharmacokinetic properties and an off-target profile, which supported additional profiling in an exploratory rat toxicology study.

Guiding Chemically Synthesized Peptide Drug Lead Optimization by Derisking Mast Cell Degranulation-Related Toxicities of a NaV1.7 Peptide Inhibitor.

Studies have shown that some peptides and small molecules can induce non IgE-mediated anaphylactoid reactions through mast cell activation. Upon activation, mast cells degranulate and release vasoactive and proinflammatory mediators, from cytoplasmic granules into the extracellular environment which can induce a cascade of severe adverse reactions. This study describes a lead optimization strategy to select NaV1.7 inhibitor peptides that minimize acute mast cell degranulation (MCD) toxicities. Various in vitro, in vivo, and PKPD models were used to screen candidates and guide peptide chemical modifications to mitigate this risk. Anesthetized rats dosed with peptides demonstrated treatment-related decreases in blood pressure and increases in plasma histamine concentrations which were reversible with a mast cell stabilizer, supporting the MCD mechanism. In vitro testing in rat mast cells with NaV1.7 peptides demonstrated a concentration-dependent increase in histamine. Pharmacodynamic modeling facilitated establishing an in vitro to in vivo correlation for histamine as a biomarker for blood pressure decline via the MCD mechanism. These models enabled assessment of structure-activity relationship (SAR) to identify substructures that contribute to peptide-mediated MCD. Peptides with hydrophobic and cationic characteristics were determined to have an elevated risk for MCD, which could be reduced or avoided by incorporating anionic residues into the protoxin II scaffold. Our analyses support that in vitro MCD assessment in combination with PKPD modeling can guide SAR to improve peptide lead optimization and ensure an acceptable early in vivo tolerability profile with reduced resources, cycle time, and animal use.

Development of ProTx-II Analogues as Highly Selective Peptide Blockers of Nav1.7 for the Treatment of Pain.

Inhibitor cystine knot peptides, derived from venom, have evolved to block ion channel function but are often toxic when dosed at pharmacologically relevant levels in vivo. The article describes the design of analogues of ProTx-II that safely display systemic in vivo blocking of Nav1.7, resulting in a latency of response to thermal stimuli in rodents. The new designs achieve a better in vivo profile by improving ion channel selectivity and limiting the ability of the peptides to cause mast cell degranulation. The design rationale, structural modeling, in vitro profiles, and rat tail flick outcomes are disclosed and discussed.

A Series of Novel, Highly Potent, and Orally Bioavailable Next-Generation Tricyclic Peptide PCSK9 Inhibitors.

Proprotein convertase subtilisin-like/kexin type 9 (PCSK9) is a key regulator of plasma LDL-cholesterol (LDL-C) and a clinically validated target for the treatment of hypercholesterolemia and coronary artery disease. Starting from second-generation lead structures such as 2, we were able to refine these structures to obtain extremely potent bi- and tricyclic PCSK9 inhibitor peptides. Optimized molecules such as 44 demonstrated sufficient oral bioavailability to maintain therapeutic levels in rats and cynomolgus monkeys after dosing with an enabled formulation. We demonstrated target engagement and LDL lowering in cynomolgus monkeys essentially identical to those observed with the clinically approved, parenterally dosed antibodies. These molecules represent the first report of highly potent and orally bioavailable macrocyclic peptide PCSK9 inhibitors with overall profiles favorable for potential development as once-daily oral lipid-lowering agents. In this manuscript, we detail the design criteria and multiparameter optimization of this novel series of PCSK9 inhibitors.

Invention of MK-8262, a Cholesteryl Ester Transfer Protein (CETP) Inhibitor Backup to Anacetrapib with Best-in-Class Properties.

Cholesteryl ester transfer protein (CETP) represents one of the key regulators of the homeostasis of lipid particles, including high-density lipoprotein (HDL) and low-density lipoprotein (LDL) particles. Epidemiological evidence correlates increased HDL and decreased LDL to coronary heart disease (CHD) risk reduction. This relationship is consistent with a clinical outcomes trial of a CETP inhibitor (anacetrapib) combined with standard of care (statin), which led to a 9% additional risk reduction compared to standard of care alone. We discuss here the discovery of MK-8262, a CETP inhibitor with the potential for being the best-in-class molecule. Novel in vitro and in vivo paradigms were integrated to drug discovery to guide optimization informed by a critical understanding of key clinical adverse effect profiles. We present preclinical and clinical evidence of MK-8262 safety and efficacy by means of HDL increase and LDL reduction as biomarkers for reduced CHD risk.

Series of Novel and Highly Potent Cyclic Peptide PCSK9 Inhibitors Derived from an mRNA Display Screen and Optimized via Structure-Based Design.

Proprotein convertase subtilisin-like/kexin type 9 (PCSK9) is a key regulator of plasma LDL-cholesterol (LDL-C) and a clinically validated target for the treatment of hypercholesterolemia and coronary artery disease. In this paper, we describe a series of novel cyclic peptides derived from an mRNA display screen which inhibit the protein-protein interaction between PCSK9 and LDLR. Using a structure-based drug design approach, we were able to modify our original screening lead 2 to optimize the potency and metabolic stability and minimize the molecular weight to provide novel bicyclic next-generation PCSK9 inhibitor peptides such as 78. These next-generation peptides serve as a critical foundation for continued exploration of potential oral, once-a-day PCSK9 therapeutics for the treatment of cardiovascular disease.

A T-cell-dependent antibody response study using a murine surrogate anti-PD-1 monoclonal antibody as an alternative to a non-human primate model.

The programmed cell death 1 (PD-1) pathway represents a major immune checkpoint which may be engaged by cells in a tumor microenvironment to overcome active T-cell immune surveillance. Pembrolizumab (Keytruda®) is a potent and highly selective humanized monoclonal antibody (mAb) of the IgG4/κ isotype designed to directly block the interaction between PD-1 and its ligands, PD-L1 and PD-L2. The current work was focused on developing a mouse T-Dependent Antibody Response (TDAR) model using a murinized rat anti-mouse PD-1 antibody (muDX400; a rodent surrogate for pembrolizumab) to evaluate the potential impact of treatment with a PD-1 inhibitor on immune responses to an antigen challenge (e.g. HBsAg in Hepatitis B vaccine). Despite the lower binding affinity and T1/2 compared to pembrolizumab, ligand blocking data indicated muDX400 had appropriate pharmacological activity and demonstrated efficacy in mouse tumor models, thus was suitable for pharmacodynamic and vaccination studies in mice. In a vaccination study in which mice were concomitantly administered muDX400 and the Hepatitis B vaccine, muDX400 was well-tolerated and did not result in any immune-mediated adverse effects. The treatment with muDX400 was associated with a shift in the ratio between naive and memory cells in both CD4+ and CD8+ T-lymphocytes in the spleen but did not affect anti-HBsAg antibody response profile. The mouse TDAR model using the Hepatitis B vaccine and the surrogate anti-PD1 monoclonal antibody was a useful tool in the evaluation of the potential immune-mediated effects of pembrolizumab following vaccination and appears to be a suitable alternative for the nonhuman primate TDAR models utilized for other checkpoint inhibitors.

Translational Safety Biomarkers of Kidney Injury.

Acute kidney injury continues to be a common problem and there continues to be a medical need for sensitive translational biomarkers for clinical monitoring. The past decade has yielded unprecedented progress in fundamental research into novel kidney biomarker evaluation and the mechanistic understanding of kidney injury; as such, these novel biomarkers increasingly are being used in preclinical drug development and in early clinical trials of drug candidates on a case-by-case basis, as well as in medical and veterinary practice. With the recent successful clinical qualification of a subset of novel accessible biomarker candidates for use in early phase clinical trials, continued clinical evaluation may enable expanded regulatory qualification for more generalized clinical use. This review provides a comprehensive overview about the discovery and development of kidney safety biomarkers with a focus on current progress in nonclinical research, progress toward translation to the clinic, and perspectives on future opportunities.

Overview on the application of transcription profiling using selected nephrotoxicants for toxicology assessment.

Microarrays allow for the simultaneous measurement of changes in the levels of thousands of messenger RNAs within a single experiment. As such, the potential for the application of transcription profiling to preclinical safety assessment and mechanism-based risk assessment is profound. However, several practical and technical challenges remain. Among these are nomenclature issues, platform-specific data formats, and the lack of uniform analysis methods and tools. Experiments were designed to address biological, technical, and methodological variability, to evaluate different approaches to data analysis, and to understand the application of the technology to other profiling methodologies and to mechanism-based risk assessment. These goals were addressed using experimental information derived from analysis of the biological response to three mechanistically distinct nephrotoxins: cisplatin, gentamicin, and puromycin aminonucleoside. In spite of the technical challenges, the transcription profiling data yielded mechanistically and topographically valuable information. The analyses detailed in the articles from the Nephrotoxicity Working Group of the International Life Sciences Institute Health and Environmental Sciences Institute suggest at least equal sensitivity of microarray technology compared to traditional end points. Additionally, microarray analysis of these prototypical nephrotoxicants provided an opportunity for the development of candidate bridging biomarkers of nephrotoxicity. The potential future extension of these applications for risk assessment is also discussed.

Identification of platform-independent gene expression markers of cisplatin nephrotoxicity.

Within the International Life Sciences Institute Committee on Genomics, a working group was formed to focus on the application of microarray technology to preclinical assessments of drug-induced nephrotoxicity. As part of this effort, Sprague-Dawley rats were treated with the nephrotoxicant cisplatin at doses of 0.3-5 mg/kg over a 4- to 144-hr time course. RNA prepared from these animals was run on a variety of microarray formats at multiple sites. A set of 93 differentially expressed genes associated with cisplatin-induced renal injury was identified on the National Institute of Environmental Health Sciences (NIEHS) custom cDNA microarray platform using quadruplicate measurements of pooled animal RNA. The reproducibility of this profile of statistically significant gene changes on other platforms, in pooled and individual animal replicate samples, and in an independent study was investigated. A good correlation in response between platforms was found among the 48 genes in the NIEHS data set that could be matched to probes on the Affymetrix RGU34A array by UniGene identifier or sequence alignment. Similar results were obtained with genes that could be linked between the NIEHS and Incyte or PHASE-1 arrays. The degree of renal damage induced by cisplatin in individual animals was commensurate with the number of differentially expressed genes in this data set. These results suggest that gene profiles linked to specific types of tissue injury or mechanisms of toxicity and identified in well-performed replicated microarray experiments may be extrapolatable across platform technologies, laboratories, and in-life studies.

Identification of putative gene based markers of renal toxicity.

This study, designed and conducted as part of the International Life Sciences Institute working group on the Application of Genomics and Proteomics, examined the changes in the expression profile of genes associated with the administration of three different nephrotoxicants--cisplatin, gentamicin, and puromycin--to assess the usefulness of microarrays in the understanding of mechanism(s) of nephrotoxicity. Male Sprague-Dawley rats were treated with daily doses of puromycin (5-20 mg/kg/day for 21 days), gentamicin (2-240 mg/kg/day for 7 days), or a single dose of cisplatin (0.1-5 mg/kg). Groups of rats were sacrificed at various times after administration of these compounds for standard clinical chemistry, urine analysis, and histological evaluation of the kidney. RNA was extracted from the kidney for microarray analysis. Principal component analysis and gene expression-based clustering of compound effects confirmed sample separation based on dose, time, and degree of renal toxicity. In addition, analysis of the profile components revealed some novel changes in the expression of genes that appeared to be associated with injury in specific portions of the nephron and reflected the mechanism of action of these various nephrotoxicants. For example, although puromycin is thought to specifically promote injury of the podocytes in the glomerulus, the changes in gene expression after chronic exposure of this compound suggested a pattern similar to the known proximal tubular nephrotoxicants cisplatin and gentamicin; this prediction was confirmed histologically. We conclude that renal gene expression profiling coupled with analysis of classical end points affords promising opportunities to reveal potential new mechanistic markers of renal toxicity.

ATM-dependent and -independent gene expression changes in response to oxidative stress, gamma irradiation, and UV irradiation.

Ataxia telangiectasia (AT) is an autosomal recessive disorder characterized by progressive cerebellar degeneration, immunodeficiencies, telangiectasias, sensitivity to ionizing radiation, and high predisposition for malignancies. The ataxia telangiectasia mutated (ATM) gene encodes a protein (ATM) with serine/threonine kinase activity. DNA-double strand breaks are known to increase its kinase activity. While cells from individuals with AT are attenuated in their G(1)-, S- and G(2)-phase cell cycle checkpoint functions in response to gamma irradiation and oxidative stress, their response to UV irradiation appears to be equivalent to that of wild-type cells. In this study, we investigated changes in gene expression in response to gamma irradiation, oxidative stress, and UV irradiation, focusing on the dependence on ATM. Doses for all three treatments were selected that resulted in roughly an equivalent induction of a G(1) checkpoint response and inhibition of progression through S phase. To investigate gene expression changes, logarithmically growing wild-type and AT dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 micro M t-butyl-hydroperoxide), or UV radiation (7.5 J/m(2)), and RNA was harvested 6 h after treatment. Gene expression analysis was performed using the NIEHS Human ToxChip 2.0 with approximately 1900 cDNA clones representing known genes and ESTs. All three treatments resulted in distinct patterns of gene expression changes, as shown previously. ATM-dependent and ATM-independent components were detected within these patterns, as were novel indications of involvement of ATM in regulation of transcription factors such as SP1, AP1 and MTF1.

Identification of distinct and common gene expression changes after oxidative stress and gamma and ultraviolet radiation.

The human genome is exposed to many different kinds of DNA-damaging agents. While most damage is detected and repaired through complex damage recognition and repair machineries, some damage has the potential to escape these mechanisms. Unrepaired DNA damage can give rise to alterations and mutations in the genome in an individual cell, which can result in malignant transformation, especially when critical genes are deregulated. In this study, we investigated gene expression changes in response to oxidative stress, gamma (gamma) radiation, and ultraviolet (UV) radiation and their potential implications in cancer development. Doses were selected for each of the three treatments, based on their ability to cause a similar G(1) checkpoint induction and slow down in early S-phase progression, as reflected by a comparable reduction in cyclin E-associated kinase activity of at least 75% in logarithmically growing human dermal diploid fibroblasts. To investigate gene expression changes, logarithmically growing dermal diploid fibroblasts were exposed to either gamma radiation (5 Gy), oxidative stress (75 microM of tert-butyl hydroperoxide (t-butyl-OOH)), or UV radiation (UVC) (7.5 J/m(2)) and RNA was harvested 6 h after treatment. Gene expression was analyzed using the NIEHS Human ToxChip 2.0 with approximately 1901 cDNA clones representing known genes and expressed sequence tags (ESTs). We were able to identify common and distinct responses in dermal diploid fibroblasts to the three different stimuli used. Within our analysis, gene expression profiles in response to gamma radiation and oxidative stress appeared to be more similar than profiles expressed after UV radiation. Interestingly, equivalent cyclin E-associated kinase activity reduction with all the three treatments was associated with greater transcriptional changes after UV radiation than after gamma radiation and oxidative stress. While samples treated with UV radiation displayed modulations of their mitogen activated protein kinase (MAPK) pathway, gamma radiation had its major influence on cell-cycle progression in S-phase and mitosis. In addition, cell cultures from different individuals displayed significant differences in their gene expression responses to DNA damage.

Genomic interrogation of mechanism(s) underlying cellular responses to toxicants.

Assessment of the impact of xenobiotic exposure on human health and disease progression is complex. Knowledge of mode(s) of action, including mechanism(s) contributing to toxicity and disease progression, is valuable for evaluating compounds. Toxicogenomics, the subdiscipline which merges genomics with toxicology, holds the promise to contributing significantly toward the goal of elucidating mechanism(s) by studying genome-wide effects of xenobiotics. Global gene expression profiling, revolutionized by microarray technology and a crucial aspect of a toxicogenomic study, allows measuring transcriptional modulation of thousands of genes following exposure to a xenobiotic. We use our results from previous studies on compounds representing two different classes of xenobiotics (barbiturate and peroxisome proliferator) to discuss the application of computational approaches for analyzing microarray data to elucidate mechanism(s) underlying cellular responses to toxicants. In particular, our laboratory demonstrated that chemical-specific patterns of gene expression can be revealed using cDNA microarrays. Transcript profiling provides discrimination between classes of toxicants, as well as, genome-wide insight into mechanism(s) of toxicity and disease progression. Ultimately, the expectation is that novel approaches for predicting xenobiotic toxicity in humans will emerge from such information.

Coordination of altered DNA repair and damage pathways in arsenite-exposed keratinocytes.

Human exposure to arsenic, a ubiquitous and toxic environmental pollutant, is associated with an increased incidence of skin cancer. However, the mechanism(s) associated with AsIII-mediated toxicity and carcinogenesis at low levels of exposure remains elusive. Aberrations in cell proliferation, oxidative damage, and DNA-repair fidelity have been implicated in sodium arsenite (AsIII)-mediated carcinogenicity and toxicity, but these events have been examined in isolation in the majority of biological models of arsenic exposure. We hypothesized that the simultaneous interaction of these effects may be important in arsenic-mediated neoplasia in the skin. To evaluate this, normal human epidermal keratinocytes (NHEK) were exposed to nontoxic doses (0.005-5 micro M) of AsIII and monitored for several physiological endpoints at the times when cells were harvested for gene expression measurements (1-24 h). Two-fluor cDNA microarray analyses indicated that AsIII treatment decreased the expression of genes associated with DNA repair (e.g., p53 and Damage-specific DNA-binding protein 2) and increased the expression of genes indicative of the cellular response to oxidative stress (e.g., Superoxide dismutase 1, NAD(P)H quinone oxidoreductase, and Serine/threonine kinase 25). AsIII also modulated the expression of certain transcripts associated with increased cell proliferation (e.g., Cyclin G1, Protein kinase C delta), oncogenes, and genes associated with cellular transformation (e.g., Gro-1 and V-yes). These observations correlated with measurements of cell proliferation and mitotic measurements as AsIII treatment resulted in a dose-dependent increase in cellular mitoses at 24 h and an increase in cell proliferation at 48 h of exposure. Data in this manuscript demonstrates that AsIII exposure simultaneously modulates DNA repair, cell proliferation, and redox-related gene expression in nontransformed, normal NHEK. It is anticipated that data in this report will serve as a foundation for furthering our knowledge of AsIII-regulated gene expression in skin and other tissues and contribute to a better understanding of arsenic toxicity and carcinogenesis.

Oxidative stress and its role in skin disease.

Skin is a major target of oxidative stress due to reactive oxygen species (ROS) that originate in the environment and in the skin itself. ROS are generated during normal metabolism, are an integral part of normal cellular function, and are usually of little harm because of intracellular mechanisms that reduce their damaging effects. Antioxidants attenuate the damaging effects of ROS and can impair and/or reverse many of the events that contribute to epidermal toxicity and disease. However, increased or prolonged free radical action can overwhelm ROS defense mechanisms, contributing to the development of cutaneous diseases and disorders. Although ROS play a role in diseases such as skin cancer, their biological targets and pathogenic mode of action are still not fully understood. In addition, strategies useful in the therapeutic management of ROS action in human skin are still lacking. This review is intended to give investigators an introduction to ROS, antioxidants, two skin disorders influenced by ROS action (skin cancer and psoriasis), and relevant model systems used to study ROS action.

Methapyrilene toxicity: anchorage of pathologic observations to gene expression alterations.

Methapyrilene (MP) exposure of animals can result in an array of adverse pathological responses including hepatotoxicity. This study investigates gene expression and histopathological alterations in response to MP treatment in order to 1) utilize computational approaches to classify samples derived from livers of MP treated rats based on severity of toxicity incurred in the corresponding tissue, 2) to phenotypically anchor gene expression pattems, and 3) to gain insight into mechanism(s) of methapyrilene hepatotoxicity. Large-scale differential gene expression levels associated with the exposure of male Sprague-Dawley rats to the rodent hepatic carcinogen MP for 1, 3, or 7 days after daily dosage with 10 or 100 mg/kg/day were monitored. Hierarchical clustering and principal component analysis were successful in classifying samples in agreement with microscopic observations and revealed low-dose effects that were not observed histopathologically. Data from cDNA microarray analysis corroborated observed histopathological alterations such as hepatocellular necrosis, bile duct hyperplasia, microvesicular vacuolization, and portal inflammation observed in the livers of MP exposed rats and provided insight into the role of specific genes in the studied toxicological processes.