New directions in incidence-dose modeling

Acute Oral, Subacute, and Developmental Toxicity Profiling of Naphthalene 2-Yl, 2-Chloro, 5-Nitrobenzoate: Assessment Based on Stress Response, Toxicity, and Adverse Outcome Pathways

The U.S. National Research Council (NRC) introduced new approaches to report toxicity studies. The NRC vision is to explore the toxicity pathways leading to the adverse effects in intact organisms by the exposure of the chemicals. This study examines the toxicity profiling of the naphthalene-2-yl 2-chloro-5-dinitrobenzoate (SF5) by adopting the vision of NRC that moves from traditional animal studies to the cellular pathways. Acute, subacute, and developmental toxicity studies were assayed according to the Organization for Economic Cooperation and Development (OECD) guidelines. The stress response pathway, toxicity pathway, and adverse effects outcome parameters were analyzed by using their standard protocols. The results showed that the acute toxicity study increases the liver enzyme levels. In a subacute toxicity study, alkaline phosphatase (ALP) levels were raised in both male and female animals. SF5 significantly increases the normal sperm count in the male animals corresponding to a decrease in the abnormality count. Developmental toxicity showed the normal skeletal and morphological parameters, except little hydrocephalus was observed in developmental toxicity. Doses of 20 mg/kg in males and 4 mg/kg in females showed decreased glutathione (GSH) levels in the kidney and liver. MDA levels were also increased in the kidney and liver. However, histopathological studies did not show any cellular change in these organs. No statistical difference was observed in histamine levels, testosterone, nuclear factor erythroid two-related factor-2 (Nrf2), and nuclear factor-kappa B (NF-κB), which showed no initiation of the stress response, toxicity, and adverse effect pathways. Immunomodulation was observed at low doses in subacute toxicity studies. It was concluded that SF5 did not produce abrupt and high-toxicity levels in organs and biochemical parameters. So, it is safe for further studies.

Multi-omics phenotyping of the gut-liver axis reveals metabolic perturbations from a low-dose pesticide mixture in rats

Health effects of pesticides are not always accurately detected using the current battery of regulatory toxicity tests. We compared standard histopathology and serum biochemistry measures and multi-omics analyses in a subchronic toxicity test of a mixture of six pesticides frequently detected in foodstuffs (azoxystrobin, boscalid, chlorpyrifos, glyphosate, imidacloprid and thiabendazole) in Sprague-Dawley rats. Analysis of water and feed consumption, body weight, histopathology and serum biochemistry showed little effect. Contrastingly, serum and caecum metabolomics revealed that nicotinamide and tryptophan metabolism were affected, which suggested activation of an oxidative stress response. This was not reflected by gut microbial community composition changes evaluated by shotgun metagenomics. Transcriptomics of the liver showed that 257 genes had their expression changed. Gene functions affected included the regulation of response to steroid hormones and the activation of stress response pathways. Genome-wide DNA methylation analysis of the same liver samples showed that 4,255 CpG sites were differentially methylated. Overall, we demonstrated that in-depth molecular profiling in laboratory animals exposed to low concentrations of pesticides allows the detection of metabolic perturbations that would remain undetected by standard regulatory biochemical measures and which could thus improve the predictability of health risks from exposure to chemical pollutants.

Stem cells in natural product and medicinal plant drug discovery-An overview of new screening approaches

Natural products remain a rich source of new drugs, and the search for bioactive molecules from nature continues to play an important role in the development of new medicines. Also, there is increasing use of herbal medicines for the treatment of a plethora of diseases, and demands for more scientific evidence for their efficacy and safety remains a huge challenge. The propensity of stem cells to differentiate into almost every cell type not only holds promise for the delivery of cell-based therapies for currently incurable diseases or a useful tool in studying cell physiology and pathophysiology. Increasingly, stem cells are becoming an important tool in preclinical drug screening and toxicity testing. In this review, we examine the scientific advances made towards the use of pluripotent stem cells as a model for the screening of plant-based medicines. The combination of well-established in vitro electrophysiological and a plethora of toxicogenomic technologies, together with the optimisation of culture methods of herbal plants and pluripotent stem cells can be explored to establish the basis for efficacy, and tissue/organ-based toxicities of many currently used medicinal plants whose efficacies and toxicities remain unknown.

Embracing Systems Toxicology at Single-Cell Resolution

As systems biology expands its multi-omic spectrum to increasing resolutions, distinguishing cells based on single-cell profiles becomes feasible. Unlike traditional bulk assays that average cellular responses and blur the distinct identities of responsive cells, single-cell technologies enable sensitive detection of small cellular changes and precise identification of those cells perturbed by toxicants. Among the suite of omic technologies that continue to expand and become affordable, single-cell RNA sequencing (scRNA-seq) is at the cutting edge and leading the way to transform systems toxicology. Single-cell systems toxicology can provide a wealth of information to elucidate cell-specific alterations and response trajectories, detect points-of-departure, map and develop dynamical models of toxicity pathways.

From the exposome to mechanistic understanding of chemical-induced adverse effects

The exposome encompasses an individual's exposure to exogenous chemicals, as well as endogenous chemicals that are produced or altered in response to external stressors. While the exposome concept has been established for human health, its principles can be extended to include broader ecological issues. The assessment of exposure is tightly interlinked with hazard assessment. Here, we explore if mechanistic understanding of the causal links between exposure and adverse effects on human health and the environment can be improved by integrating the exposome approach with the adverse outcome pathway (AOP) concept that structures and organizes the sequence of biological events from an initial molecular interaction of a chemical with a biological target to an adverse outcome. Complementing exposome research with the AOP concept may facilitate a mechanistic understanding of stress-induced adverse effects, examine the relative contributions from various components of the exposome, determine the primary risk drivers in complex mixtures, and promote an integrative assessment of chemical risks for both human and environmental health.

A tiered, integrated biological and chemical monitoring framework for contaminants of emerging concern in aquatic ecosystems

The chemical-specific risk-based paradigm that informs monitoring and assessment of environmental contaminants does not apply well to the many thousands of new chemicals that are being introduced into ambient receiving waters. We propose a tiered framework that incorporates bioanalytical screening tools and diagnostic nontargeted chemical analysis to more effectively monitor for contaminants of emerging concern (CECs). The framework is based on a comprehensive battery of in vitro bioassays to first screen for a broad spectrum of CECs and nontargeted analytical methods to identify bioactive contaminants missed by the currently favored targeted analyses. Water quality managers in California have embraced this strategy with plans to further develop and test this framework in regional and statewide pilot studies on waterbodies that receive discharge from municipal wastewater treatment plants and stormwater runoff. In addition to directly informing decisions, the data obtained using this framework can be used to construct and validate models that better predict CEC occurrence and toxicity. The adaptive interplay among screening results, diagnostic assessment and predictive modeling will allow managers to make decisions based on the most current and relevant information, instead of extrapolating from parameters with questionable linkage to CEC impacts. Integr Environ Assess Manag 2016;12:540-547. © 2015 SETAC.

Systems Biology and Biomarkers of Early Effects for Occupational Exposure Limit Setting

In a recent National Research Council document, new strategies for risk assessment were described to enable more accurate and quicker assessments. This report suggested that evaluating individual responses through increased use of bio-monitoring could improve dose-response estimations. Identification of specific biomarkers may be useful for diagnostics or risk prediction as they have the potential to improve exposure assessments. This paper discusses systems biology, biomarkers of effect, and computational toxicology approaches and their relevance to the occupational exposure limit setting process. The systems biology approach evaluates the integration of biological processes and how disruption of these processes by chemicals or other hazards affects disease outcomes. This type of approach could provide information used in delineating the mode of action of the response or toxicity, and may be useful to define the low adverse and no adverse effect levels. Biomarkers of effect are changes measured in biological systems and are considered to be preclinical in nature. Advances in computational methods and experimental -omics methods that allow the simultaneous measurement of families of macromolecules such as DNA, RNA, and proteins in a single analysis have made these systems approaches feasible for broad application. The utility of the information for risk assessments from -omics approaches has shown promise and can provide information on mode of action and dose-response relationships. As these techniques evolve, estimation of internal dose and response biomarkers will be a critical test of these new technologies for application in risk assessment strategies. While proof of concept studies have been conducted that provide evidence of their value, challenges with standardization and harmonization still need to be overcome before these methods are used routinely.

A novel model for neuroendocrine toxicology: neurobehavioral effects of BPA exposure in a prosocial species, the prairie vole (Microtus ochrogaster)

Impacts on brain and behavior have been reported in laboratory rodents after developmental exposure to bisphenol A (BPA), raising concerns about possible human effects. Epidemiological data suggest links between prenatal BPA exposure and altered affective behaviors in children, but potential mechanisms are unclear. Disruption of mesolimbic oxytocin (OT)/vasopressin (AVP) pathways have been proposed, but supporting evidence is minimal. To address these data gaps, we employed a novel animal model for neuroendocrine toxicology: the prairie vole (Microtus ochrogaster), which are more prosocial than lab rats or mice. Male and female prairie vole pups were orally exposed to 5-μg/kg body weight (bw)/d, 50-μg/kg bw/d, or 50-mg/kg bw/d BPA or vehicle over postnatal days 8-14. Subjects were tested as juveniles in open field and novel social tests and for partner preference as adults. Brains were then collected and assessed for immunoreactive (ir) tyrosine hydroxylase (TH) (a dopamine marker) neurons in the principal bed nucleus of the stria terminalis (pBNST) and TH-ir, OT-ir, and AVP-ir neurons in the paraventricular nucleus of the hypothalamus (PVN). Female open field activity indicated hyperactivity at the lowest dose and anxiety at the highest dose. Effects on social interactions were also observed, and partner preference formation was mildly inhibited at all dose levels. BPA masculinized principal bed nucleus of the stria terminalis TH-ir neuron numbers in females. Additionally, 50-mg/kg bw BPA-exposed females had more AVP-ir neurons in the anterior PVN and fewer OT-ir neurons in the posterior PVN. At the 2 lowest doses, BPA eliminated sex differences in PVN TH-ir neuron numbers and reversed this sex difference at the highest dose. Minimal behavioral effects were observed in BPA-exposed males. These data support the hypothesis that BPA alters affective behaviors, potentially via disruption of OT/AVP pathways.

High throughput nano-liter RT-qPCR to classify soil contamination using a soil arthropod

Background

To incorporate genomics data into environmental assessments a mechanistic perspective of interactions between chemicals and induced biological processes needs to be developed. Since chemical compounds with structural similarity often induce comparable biological responses in exposed animals, gene expression signatures can serve as a starting point for the assessment of chemicals and their toxicity, but only when relevant and stable gene panels are available. To design such a panel, we isolated differentially expressed gene fragments from the soil arthropod Folsomia candida, a species often used for ecotoxicological testing. Animals were exposed to two chemically distinct compounds, being a metal (cadmium) and a polycyclic aromatic hydrocarbon (phenanthrene). We investigated the affected molecular responses resulting from either treatment and developed and validated 44 qPCR assays for their responses using a high throughput nano-liter RT-qPCR platform for the analysis of the samples.

Results

Suppressive subtractive hybridization (SSH) was used to retrieve stress-related gene fragments. SSH libraries revealed pathways involved in mitochondrial dysfunction and protein degradation for cadmium and biotransformation for phenanthrene to be overrepresented. Amongst a small cluster of SSH-derived cadmium responsive markers were an inflammatory response protein and an endo-glucanase. Conversely, cytochrome P450 family 6 or 9 was specifically induced by phenanthrene. Differential expressions of these candidate biomarkers were also highly significant in the independently generated test sample set. Toxicity levels in different training samples were not reflected by any of the markers' intensity of expressions. Though, a model based on partial least squares differential analysis (PLS-DA) (with RMSEPs between 9 and 22% and R²s between 0.82 and 0.97) using gene expressions of 25 important qPCR assays correctly predicted the nature of exposures of test samples.

Conclusions

For the application of molecular bio-indication in environmental assessments, multivariate analyses obviously have an added value over univariate methods. Our results suggest that compound discrimination can be achieved by PLS-DA, based on a hard classification of the within-class rankings of samples from a test set. This study clearly shows that the use of high throughput RT-qPCR could be a valuable tool in ecotoxicology combining high throughput with analytical sensitivity.

Predicting chemical impacts on vertebrate endocrine systems

Animals have evolved diverse protective mechanisms for responding to toxic chemicals of both natural and anthropogenic origin. From a governmental regulatory perspective, these protective responses complicate efforts to establish acceptable levels of chemical exposure. To explore this issue, we considered vertebrate endocrine systems as potential targets for environmental contaminants. Using the hypothalamic-pituitary-thyroid (HPT), hypothalamic-pituitary-gonad (HPG), and hypothalamic-pituitary-adrenal (HPA) axes as case examples, we identified features of these systems that allow them to accommodate and recover from chemical insults. In doing so, a distinction was made between effects on adults and those on developing organisms. This distinction was required because endocrine system disruption in early life stages may alter development of organs and organ systems, resulting in permanent changes in phenotypic expression later in life. Risk assessments of chemicals that impact highly regulated systems must consider the dynamics of these systems in relation to complex environmental exposures. A largely unanswered question is whether successful accommodation to a toxic insult exerts a fitness cost on individual animals, resulting in adverse consequences for populations. Mechanistically based mathematical models of endocrine systems provide a means for better understanding accommodation and recovery. In the short term, these models can be used to design experiments and interpret study findings. Over the long term, a set of validated models could be used to extrapolate limited in vitro and in vivo testing data to a broader range of untested chemicals, species, and exposure scenarios. With appropriate modification, Tier 2 assays developed in support of the U.S. Environmental Protection Agency's Endocrine Disruptor Screening Program could be used to assess the potential for accommodation and recovery and inform the development of mechanistically based models.

Application of toxicogenomics in hepatic systems toxicology for risk assessment: acetaminophen as a case study

Hepatic systems toxicology is the integrative analysis of toxicogenomic technologies, e.g., transcriptomics, proteomics, and metabolomics, in combination with traditional toxicology measures to improve the understanding of mechanisms of hepatotoxic action. Hepatic toxicology studies that have employed toxicogenomic technologies to date have already provided a proof of principle for the value of hepatic systems toxicology in hazard identification. In the present review, acetaminophen is used as a model compound to discuss the application of toxicogenomics in hepatic systems toxicology for its potential role in the risk assessment process, to progress from hazard identification towards hazard characterization. The toxicogenomics-based parallelogram is used to identify current achievements and limitations of acetaminophen toxicogenomic in vivo and in vitro studies for in vitro-to-in vivo and interspecies comparisons, with the ultimate aim to extrapolate animal studies to humans in vivo. This article provides a model for comparison of more species and more in vitro models enhancing the robustness of common toxicogenomic responses and their relevance to human risk assessment. To progress to quantitative dose-response analysis needed for hazard characterization, in hepatic systems toxicology studies, generation of toxicogenomic data of multiple doses/concentrations and time points is required. Newly developed bioinformatics tools for quantitative analysis of toxicogenomic data can aid in the elucidation of dose-responsive effects. The challenge herein is to assess which toxicogenomic responses are relevant for induction of the apical effect and whether perturbations are sufficient for the induction of downstream events, eventually causing toxicity.

Modeling Dose-response at Low Dose: A Systems Biology Approach for Ionization Radiation

For ionization radiation (IR) induced cancer, a linear non-threshold (LNT) model at very low doses is the default used by a number of national and international organizations and in regulatory law. This default denies any positive benefit from any level of exposure. However, experimental observations and theoretical biology have found that both linear and J-shaped IR dose-response curves can exist at those very low doses. We develop low dose J-shaped dose-response, based on systems biology, and thus justify its use regarding exposure to IR. This approach incorporates detailed, molecular and cellular descriptions of biological/toxicological mechanisms to develop a dose-response model through a set of nonlinear, differential equations describing the signaling pathways and biochemical mechanisms of cell cycle checkpoint, apoptosis, and tumor incidence due to IR. This approach yields a J-shaped dose response curve while showing where LNT behaviors are likely to occur. The results confirm the hypothesis of the J-shaped dose response curve: the main reason is that, at low-doses of IR, cells stimulate protective systems through a longer cell arrest time per unit of IR dose. We suggest that the policy implications of this approach are an increasingly correct way to deal with precautionary measures in public health.

Toxicity testing in the 21st century: a vision and a strategy

With the release of the landmark report Toxicity Testing in the 21st Century: A Vision and a Strategy, the U.S. National Academy of Sciences, in 2007, precipitated a major change in the way toxicity testing is conducted. It envisions increased efficiency in toxicity testing and decreased animal usage by transitioning from current expensive and lengthy in vivo testing with qualitative endpoints to in vitro toxicity pathway assays on human cells or cell lines using robotic high-throughput screening with mechanistic quantitative parameters. Risk assessment in the exposed human population would focus on avoiding significant perturbations in these toxicity pathways. Computational systems biology models would be implemented to determine the dose-response models of perturbations of pathway function. Extrapolation of in vitro results to in vivo human blood and tissue concentrations would be based on pharmacokinetic models for the given exposure condition. This practice would enhance human relevance of test results, and would cover several test agents, compared to traditional toxicological testing strategies. As all the tools that are necessary to implement the vision are currently available or in an advanced stage of development, the key prerequisites to achieving this paradigm shift are a commitment to change in the scientific community, which could be facilitated by a broad discussion of the vision, and obtaining necessary resources to enhance current knowledge of pathway perturbations and pathway assays in humans and to implement computational systems biology models. Implementation of these strategies would result in a new toxicity testing paradigm firmly based on human biology.

A bistable switch underlying B-cell differentiation and its disruption by the environmental contaminant 2,3,7,8-tetrachlorodibenzo-p-dioxin

The differentiation of B cells into antibody-secreting plasma cells upon antigen stimulation, a crucial step in the humoral immune response, is disrupted by 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). Several key regulatory proteins in the B-cell transcriptional network have been identified, with two coupled mutually repressive feedback loops among the three transcription factors B-cell lymphoma 6 (Bcl-6), B lymphocyte-induced maturation protein 1(Blimp-1), and paired box 5 (Pax5) forming the core of the network. However, the precise mechanisms underlying B-cell differentiation and its disruption by TCDD are not fully understood. Here we show with a computational systems biology model that coupling of the two feedback loops at the Blimp-1 node, through parallel inhibition of Blimp-1 gene activation by Bcl-6 and repression of Blimp-1 gene deactivation by Pax5, can generate a bistable switch capable of directing B cells to differentiate into plasma cells. We also use bifurcation analysis to propose that TCDD may suppress the B-cell to plasma cell differentiation process by raising the threshold dose of antigens such as lipopolysaccharide required to trigger the bistable switch. Our model further predicts that high doses of TCDD may render the switch reversible, thus causing plasma cells to lose immune function and dedifferentiate to a B cell-like state. The immunotoxic implications of these predictions are twofold. First, TCDD and related compounds would disrupt the initiation of the humoral immune response by reducing the proportion of B cells that respond to antigen and differentiate into antibody-secreting plasma cells. Second, TCDD may also disrupt the maintenance of the immune response by depleting the pool of available plasma cells through dedifferentiation.

Incorporating information from the U.S. Department of Energy low-dose program into regulatory decision-making: three policy integration challenges

The U.S. Department of Energy (U.S. DOE) sponsors a research program aimed at gaining a better understanding of how low-dose radiation affects cellular functioning and progression toward disease. There have been calls to incorporate into regulatory decision-making the scientific information that this program has produced. After a discussion of the evolution of radiation protection law and the weight-of-evidence approach that agencies employ, this paper offers some preliminary thoughts about how to approach this complex and important policy question. Three implementation challenges are identified and discussed. The first implementation challenge involves explaining low-dose effects in a systems biology model. The second challenge arises when issues of population susceptibility are juxtaposed against molecular and mechanistic studies, such as those that make up much of the U.S. DOE low-dose program. The third challenge concerns integrating the results of radiation epidemiology, especially epidemiologic studies among cohorts that are exposed to low dose and low-dose rate radiation, with the results of U.S. DOE low-dose studies.

On the incorporation of chemical-specific information in risk assessment

This paper describes the evolution of chemical risk assessment from its early dependence on generic default approaches to the current situation in which mechanistic and biokinetic data are routinely incorporated to support a more chemical-specific approach. Two methodologies that have played an important role in this evolution are described: mode-of-action evaluation and physiologically based biokinetic (PBBK) modelling. When used together, these techniques greatly increase the opportunity for the incorporation of biokinetic and mechanistic data in risk assessment. The resulting risk assessment approaches are more appropriately tailored to the specific chemical and are more likely to provide an accurate assessment of the potential hazards associated with human exposures. The appropriate application of PBBK models in risk assessment demands well-formulated statements about the chemical mode of action. It is this requirement for an explicit, mechanistic hypothesis that gives biologically motivated models their power, but at the same time serves as the greatest impediment to the acceptance of a chemical-specific risk assessment approach by regulators. The chief impediment to the regulatory acceptance and application of PBBK models in risk assessment is concern about uncertainties associated with their use. To some extent such concerns can be addressed by the development of generally accepted approaches for model evaluation and quantitative uncertainty analysis. In order to assure the protection of public health while limiting the economic and social consequences of over-regulation, greater dialogue between researchers and regulators is crucially needed to foster an increased use of emerging scientific information and innovative methods in chemical risk assessments.

Biomedical implications of information processing in chemical systems: Non-classical approach to photochemistry of coordination compounds

Analogies between photoactive nitric oxide generators and various electronic devices: logic gates and operational amplifiers are presented. These analogies have important biological consequences: application of control parameters allows for better targeting and control of nitric oxide drugs. The same methodology may be applied in the future for other therapeutic strategies and at the same time helps to understand natural regulatory and signaling processes in biological systems.

Research toward the development of a biologically based dose response assessment for inorganic arsenic carcinogenicity: a progress report

Cancer risk assessments for inorganic arsenic have been based on human epidemiological data, assuming a linear dose response below the range of observation of tumors. Part of the reason for the continued use of the linear approach in arsenic risk assessments is the lack of an adequate biologically based dose response (BBDR) model that could provide a quantitative basis for an alternative nonlinear approach. This paper describes elements of an ongoing collaborative research effort between the CIIT Centers for Health Research, the U.S. Environmental Protection Agency, ENVIRON International, and EPRI to develop BBDR modeling approaches that could be used to inform a nonlinear cancer dose response assessment for inorganic arsenic. These efforts are focused on: (1) the refinement of physiologically based pharmacokinetic (PBPK) models of the kinetics of inorganic arsenic and its metabolites in the mouse and human; (2) the investigation of mathematical solutions for multi-stage cancer models involving multiple pathways of cell transformation; (3) the review and evaluation of the literature on the dose response for the genomic effects of arsenic; and (4) the collection of data on the dose response for genomic changes in the urinary bladder (a human target tissue for arsenic carcinogenesis) associated with in vivo drinking water exposures in the mouse as well as in vitro exposures of both mouse and human cells. An approach is proposed for conducting a biologically based margin of exposure risk assessment for inorganic arsenic using the in vitro dose response for the expression of genes associated with the obligatory precursor events for arsenic tumorigenesis.

Dose-Incidence Modeling: Consequences of Linking Quantal Measures of Response to Depletion of Critical Tissue Targets

In developing mechanistic PK-PD models, incidence of toxic responses in a population has to be described in relation to measures of biologically effective dose (BED). We have developed a simple dose-incidence model that links incidence with BED for compounds that cause toxicity by depleting critical cellular target molecules. The BED in this model was the proportion of target molecule adducted by the dose of toxic compound. Our modeling approach first estimated the proportion depleted for each dose and then calculated the tolerance distribution for toxicity in relation to either administered dose or log of administered dose. We first examined cases where the mean of the tolerance distribution for toxicity occurred when a significant proportion of target had been adducted (i.e., more than half). When a normal distribution was assumed to exist for the relationship of incidence and BED, the tolerance distribution based on administered dose for these cases becomes asymmetrical and logarithmic transformations of the administered dose axis lead to a more symmetrical distribution. These linked PK-PD models for tissue reactivity, consistent with conclusions from other work for receptor binding models (Lutz et al., 2005), indicate that log normal distributions with administered dose may arise from normal distributions for BED and nonlinear kinetics between BED and administered dose. These conclusions are important for developing biologically based dose response (BBDR) models that link incidences of toxicity or other biological responses to measures of BED.

Dose-response modeling in reproductive toxicology in the systems biology era

Systems biology approaches for modeling cellular signaling networks affected by chemical exposures should soon produce integrated methodologies capable of predicting dose-response relationships for developmental toxicants and for other toxic responses. This paper outlines an emerging strategy for systems biology approaches in dose-response modeling. Genome-wide functional screens, bioinformatic tools, and network mapping technologies together can provide directed graph representations of the cellular signaling networks. The graphical representations can be converted into mathematical models that permit predicting the shapes of dose-response curves for altered cell signaling by test compounds during development. Systems biology approaches require interdisciplinary teams with expertise in reproduction, cell biology, signal transduction, mathematical/biomedical modeling, and risk assessment. In addition to outlining a systems approach for dose-response research, this paper discusses initial stages of application of this strategy to examine inhibition of steroidogenesis in testes by phthalate esters.