Dose-Response Analysis in Risk Assessment: Evaluation of Behavioral Specificity John R. Glowa Laboratory of Medicinal Chemistry, National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health, Bethesda, Maryland Several methods of quantitative risk assessment that have been described recently are particularly applicable to neurotoxic end points. These methods can be broadly divided into two types of approaches based on their treatment of dose-response data to estimate risks. Benchmark approaches estimate risks using variability in response to a fixed dose level in comparison with background control variability. Probabilistic approaches estimate risks using the variability in the dose to produce a small effect in the sample population. The current report seeks to extend the development of probabilistic approaches for neurotoxic end points. Because behavioral data are often used to assess therapeutic efficacy as well as toxicity (unwanted effects), this analysis focused on the relative risks of producing these effects with the same agent. The therapeutic potential of GBR 12909 was determined by its ability to decrease cocaine-maintained responding in monkeys. The effects of this agent were also assessed in the same monkeys using food-maintained responding to provide an indication of behavioral toxicity. GBR 12909 decreased both behaviors, with complete decreases on drug-seeking behavior occurring at doses that had minimal effects on food-maintained responding. The difference in the estimates of doses to decrease drug-seeking and food-maintained behavior suggested that specific therapeutic effects could be obtained in the absence of unwanted side effects for a definable proportion of the population. These results also suggest that multiple behavioral end points can be useful for identifying specific effects of chemicals for the purposes of risk assessment. - Environ Health Perspect 104(Suppl 21:391-396 (1996) Key words: risk assessment, drug abuse, rhesus monkeys, GBR 12909, schedule-controlled behavior, neurotoxicology Introduction Risk assessment is the attempt to estimate agent that produce adverse effects in the chance of obtaining an adverse effect of humans, risk assessments are typically based exposure to an agent. Although it is possible on effects obtained over a range of non- to predict virtually safe levels if sufficient toxic-to-toxic doses determined in animals. information is known about doses of an Such methods were developed initially to This paper was prepared as background for the Workshop on Risk Assessment Methodology for Neurobehavioral Toxicity convened by the Scientific Group on Methodologies for the Safety Evaluation of Chemicals (SGOMSEC) held 12-17 June 1994 in Rochester, New York. Manuscript received 1 February 1995; manuscript accepted 17 December 1995. This work was supported, in part, by a NIDA Interagency Agreement #RA-ND-94-24 (JR Glowa, principal investigator). Portions of the data were presented at the College for Problems of Drug Dependence, Toronto, Canada, 1993. All animal procedures conformed to the Guide for Care and Use of Laboratory Animals endorsed by the National Institutes of Health. The author gratefully acknowledges the counsel and support of PB Dews and RC MacPhail. Their continued discussions of risk assessment methodology, as well as comments on ear- lier versions of this manuscript, are greatly appreciated. Address correspondence to Dr. John R. Glowa, Chief, Behavioral Pharmacology Unit, LMC/NIDDK/NIH, Building 14D, Room 311, Bethesda, MD 20892. Telephone: (301) 496-5059. Fax: (301) 402-0378. E-mail: G7W@cu.nih.gov Abbreviations used: ADI, acceptable daily intake; RfD, reference dose; NOAEL, no observable adverse effect level; LOAEL, lowest observable adverse effect level; DA, dopamine; TLVs, threshold limit values; FRF, fixed- ratio 3G-response food presentation; FRC, fixed-ratio 30-response cocaine injection; ANOVA, analysis of vari- ance; LD50, median lethal dose; ED5, median effective dose. predict the chance of producing carcino- genic and mutagenic effects with exposure to ionizing radiation. Higher doses were often associated with an increased inci- dence of cancers, so the focus of risk assess- ment shifted toward the attempt to predict effects produced by very low doses. Two problems developed. First, the spontaneous rate of cancer production called into doubt the ability to measure an increase in risks associated with very low doses of an agent. Second, several low-dose extrapolation approaches were developed based on the assumption that any increment in dose would increase the probability of an effect. These issues have questioned the relevance of the low-dose extrapolation approach for other types of adverse effect. For example, few neurotoxic events occur spontaneously, although disease, age, or other processes may exacerbate their effects. Likewise, mechanisms such as cellular repair, plastic- ity, system redundancy, or tolerance may contribute to the lack of measurable effect following exposure to low doses of neuro- toxic agents. Alternatively, some agents have been simply found to lack effects at low doses. As a result, risk assessment methodol- ogy for noncancer end points shifted away from the estimation of effects of very low doses toward the attempt to estimate levels of an agent that might be considered safe. The earliest practice used to predict safe levels was the acceptable daily intake approach (ADI), which attempted to predict a dose (within an order of magni- tude of uncertainty) that could be tolerated over the lifetime without producing harm. However, while some argue that it was vir- tually impossible to determine an entirely safe dose, others pointed to an apparent inappropriateness of using the word "accept- able" in the context of a poison. This led to a revised terminology called the reference dose (RfD) approach. This approach first finds a dose level with no effect (i.e., no observable adverse effect level, or NOAEL) or a minimal effect (i.e., lowest observable adverse effect level, or LOAEL) in animals and then divides this dose by a series of uncertainty factors (invariably 10-fold each) to calculate a RfD for humans. The result- ing RfD could be several orders of magni- tude smaller than the original NOAEL or LOAEL (1). There have been a number of criticisms of the RfD approach (2-6), including a) its dependence on the actual dose spacing used, b) its failure to incorpo- rate the slope of the dose-response function, Environmental Health Perspectives * Vol 104, Supplement 2 - April 1996 391 J.R. GLOWA c) a lack of biological justification for uncertainty factors, d) its questionable applicability for some types of data, e) the fact that the RfD is not a risk figure but rather a point estimate of an exposure at which risk is considered negligible, and f) its perverse ability to reward biased risk assessors with less conservative risk figures for using small sample sizes. In response, a new trend in risk assessment that uses more quantitative approaches has emerged. These approaches are of interest because they circumvent many of the problems asso- ciated with the RfD and low-dose extrapola- tion approaches and thus provide true quantitative risk estimates. Given the recent increase in concerns over other end points such as neurotoxicity (7,8), and the equivo- cal performance of the older approaches to estimate risks (9,10), the time is ripe to scrutinize these quantitative methods. Quantitative approaches generally describe dose-effect functions in mathe- matical terms and specifically incorporate the variability in different parameters of this function to assess risks (6,11). One subfam- ily of these approaches can be characterized by its focus on variability in an effect at a fixed-exposure level (effect-tolerance or benchmark approaches). One advantage of this approach is that statistical techniques to fit functions and characterize variability are readily available. Another subfamily of approaches is characterized by a focus on variability in exposure levels that can be shown to produce a fixed effect (dose-toler- ance or probabilistic approaches). The advantage of the dose-tolerance approaches is that they characterize individual differ- ences, thus providing direct information as to the potential adequacy of the use of an uncertainty factor for individual variability (6). The current report extends the use of the probabilistic (dose-tolerance) approach to multiple behavioral end points in the same organisms to address the specificity of the effect of the agent. The agent used for these analyses is a drug that may have therapeutic merit. Though drugs are not commonly thought of as neurotoxicants, any pharmacologist or toxicologist will acknowledge that high doses of therapeutic agents will have adverse effects. It has been a common practice over the years to compare doses expected to pro- duce a therapeutic effect with those expected to produce a toxic effect. The cur- rent report departs from this approach by attempting to compare the risk of producing these effects. Another small departure from standard methods is the use of a different end point in the current report. Behavioral toxicology has been primarily concerned with the direct effects of toxic agents on learned or acquired behaviors. Most often, these behaviors are experimentally devel- oped in the laboratory, typically by train- ing animals to respond under schedules of food reinforcement. Stable rates and pat- terns of behavior provide a baseline from which a dose-effect function can be deter- mined. While almost all of these studies have employed behaviors that lead to the delivery of food as a reinforcer, some others have shown that the event used to main- tain behavior can determine different effects for some agents (12). It has been known for some time that drugs can be used to maintain responding (e.g., self- administration studies). Interestingly, some toxicants have been shown to serve as rein- forcers. For example, toluene has been shown to support responding in monkeys (13). Little is known of solvent abuse; however, growing knowledge suggests that, if left unchecked, drug abuse results in changes in behavior that are sufficient to meet the criteria of a psychiatric disorder (14). Neurotoxic changes are also seen with prolonged exposure to drugs of abuse (15). Thus, drug abuse may be the most prevalent form of behavioral toxicology observable in humans today. One question is whether a therapeutic agent could be developed that would decrease drug abuse but not alter other normal behavior. Previous studies have shown that oper- ant responding can be maintained by cocaine delivery (self-administration), which provides an animal model of drug abuse (16). Left unrestrained, cocaine self-admin- istration can occur to the point of produc- ing serious harm in animals (17). Such excessive behavior may result in neurotoxic effects in central mesolimbic dopamine (DA) pathways associated with reinforce- ment (18). In contrast to some forms of drug abuse, there is no currently accepted strategy to treat cocaine abuse. However, clues from successful treatment approaches suggest that targeting mechanisms associ- ated with the abuse potential of a drug may lead to a treatment drug. Recently, a specific neuronal site (the DA reuptake or transporter) has been identified as a recep- tor that is strongly associated with the rein- forcing effects of cocaine (19). Studies have shown that drugs that bind to this site need not exhibit reinforcing effects in humans (20) or in animals (21). These observations have led to the suggestion that long-acting DA reuptake inhibitors may attenuate cocaine-seeking behavior in a manner similar to the ability of methadone to decrease opiate-seeking behavior or nicotine patches to decrease smoking (22). However, there is no clear consensus as to the type of effect a medication should exhibit on cocaine self-administration in animals that would predict its efficacy to decrease cocaine-seeking behavior in humans. The current research is meant to address this issue from a therapeutic point of view; however, it is obvious that some concern should be directed at the potential toxicity of such an agent. Thus, the current studies are designed to compare the relative therapeutic actions of an agent being devel- oped to treat cocaine abuse with other unwanted effects. The current studies focus on one compound, GBR 12909, which has previously been shown to exhibit a number of effects similar to those of cocaine and has a high affinity for the DA reuptake site. In contrast, because GBR 12909 is long acting and has a chemical structure unlike cocaine, it may block the effects of cocaine. Interestingly, this drug has been reported to exhibit a nonstimulant profile of action in normal human volunteers following oral administration (23). Thus, this agent could serve as a prototype from which even more therapeutic and less toxic analogs could be developed. The primary focus of the cur- rent work was to separate those effects. For these purposes, the effects of GBR 12909 on food-maintained behavior were taken to represent a behaviorally toxic effect, and its effects on cocaine self- administration were used as an indication of its therapeutic potential. To accomplish the risk assessment analysis of these data, methods previously developed to assess the risks of exposure to agents using a single behavioral end point were used (24). This approach was origi- nally described by Dews (3,4). Data from individual subjects are first described by means of a mathematical function, and the dose resulting in a small but measurable effect is estimated from each function. Through replications in different subjects, a distribution of individual point estimates is obtained. The variability in these esti- mates is used to predict population toler- ances, from which the proportion of the population that would be expected to be affected at lower levels can be determined. This method has been used extensively to assess the risks associated with exposure to a variety of organic solvents (25-27). These estimates of a small decrement in neurobe- havioral functioning in relatively small Environmental Health Perspectives * Vol 104. Supplement 2 * April 1996 392 NEUROBEHAVIORAL RISK ASSESSMENTS FOR THERAPEUTICS sample sizes (n= 10-12) agreed well with established threshold limit values (TLVs). A subsequent report (28) compared estimates from small (n= 10) and large (n= 40) sample sizes to illustrate the effectiveness of these rel- atively small experiments, as well as to com- pare estimates obtained with this approach to those of others. This approach has also been compared with other approaches in terms of its use of the dose-effect data (6,10), applicability to cross-species extrapo- lation (29), and ability to characterize the neurotoxic effects of the pesticide carbaryl (6). The potential advantage of using two end points in the present studies is that the nonspecific (toxic) effects of the agent can be directly compared to its intended therapeutic effects, thereby producing a thzerapeutic ratio. Methods This report uses methods that were previ- ously published (30,31). Briefly, rhesus monkeys were surgically implanted with a subcutaneous port/catheter system (32) and trained to respond under a multiple fixed-ratio 30-response food presentation (FRF), fixed-ratio 30-response cocaine injection (FRC) schedule to establish a baseline where approximately equal rates of responding were maintained by food and cocaine. The acute effects of GBR 12909 were then determined on these perfor- mances using eight monkeys. Individual control rates of responding for both perfor- mances were obtained for several noninjec- tion control sessions and served as a baseline for each monkey. The effects of pretreat- ment with GBR 12909 were assessed on this baseline by comparing the mean over- all rates of responding during drug ses- sions with rates during control sessions. Two advantages of these data are that the effects of the agent were well characterized (doses were assessed 2-3 times per monkey) and the individual differences in effect occurring among the monkeys could serve to model individual differences in the human population. GBR 12909 (1-1.7 mg/kg) decreased responding maintained by cocaine almost completely without affecting food-main- tained responding. This type of effect showed that cocaine-maintained respond- ing could be decreased by low doses of this drug. With higher doses of GBR 12909, food-maintained responding decreased. This effect may have been associated with a drug-induced disruption of behavior or because some aspect of food-maintained responding provided a greater resistance to the drug effect than that seen with cocaine-maintained responding. A mean dose-effect function for each of these behav- ioral end points was then constructed using the individual data from the eight animals. While there was no overlap of these curves (i.e., GBR 12909 always decreased cocaine- maintained responding at lower doses than food-maintained responding), some animals were more sensitive to the effects of GBR 12909 than others. However, traditional statistical approaches confirmed that GBR 12909 had different effects on these two behaviors. For example, single factor repeated measures analysis of variance (ANOVA) indicated a significant main effect of maintaining event (F= 18.428; df = 1, 14; p< 0.0007), dose (p < 0.)000 1), and an interaction. Figure 1 shows the essential elements of the method of risk assessment applied to these data. First, the effect of GBR 12909 on each type of response was determined relative to each monkey's baseline rate for each response by dividing rates during drug sessions by rates during control sessions for each component. This represents a normal dose-effect function for a single animal, only two different end points could be described. Next, a straight line was fit to the linear portion (only those doses that pro- duced greater than a 20% decrement but less than an 80% decrement in responding) of each curve using linear regression. Effect 100 - 80 - 60 - 40 - 20 - I I I 2 -1 0 1 In dose GBR 12909, mg/kg Figure 1. Transformed dose-effect summaries of the rate-decreasing effects of GBR 12909 on responding maintained under multiple FR30 schedules of cocaine delivery (solid line) and food-pellet delivery (dotted line). Effects are the mean of double-determined obser- vations of the rate-decreasing effects (see text for a description of the restriction of the range) in eight monkeys, expressed as a percentage of the individual baseline rate of responding. Doses were transformed to the natural logarithm (In). Variability (1 z-score) in the dose to produce 90 and 10% effects, based on dif- ferences in the individual curves, is indicated by the horizontal bars. (percent of control) was plotted against the natural logarithm (ln) of the dose. The figure shows only the mean circumscribed descending limb, from 90 to 10% of con- trol. The variability (SD) in the 90 and 10% of control was determined from the individ- ual functions (horizontal lines). For exam- ple, the mean ln dose of GBR 12909 to produce a 10% decrement in cocaine-main- tained responding was -1, with a SD of about 0.5. As is apparent from Figure 1, the doses of GBR 12909 required to decrease food-maintained responding were high. In fact, for three animals complete decreases in food-maintained responding were not pro- duced by the highest dose of GBR 12909 studied. Since previous studies suggested that a higher dose (5.6 mg/kg) would com- pletely abolish responding, it was conserva- tively assumed that similar effects would be obtained in these animals and that it would be in the best interest of the animals to not test higher doses. The estimates of risk were then obtained by calculating the successive prob- abilities (i.e., p=0.1, p=0.01, p=0.001) of obtaining a 10% decrement in this group of animals using a z-score approach. In the past, only estimates from the low end of these distributions have been assessed because of the specific interest in the effects of low doses. In this study, the overlap of the effects of similar doses of GBR 12909 on both cocaine- and food-maintained responding could be directly compared. Transformation then allowed determina- tion of doses (in milligram per kilogram) expected to decrease responding by 10% in the specified proportion (i.e., 1 out of 10, 1 out of 100, etc.) of the population (24). Figure 1 characterizes the descending limbs of the dose-effect curve for each end point. (A mean curve was never used for the analysis; it simply represents the mean slope joining the mean 90 and 10% point esti- mates.) The mean parameters (a=intercept; b=slope) for the individual curves, the 90% effect (in In dose and in milligram per kilo- gram), and the SD of In dose for cocaine- maintained responding were FRC, mg/kg a b 90% (In) 90% Mean 18.103 -74.52 -1.0538 (SD) (0.4226) 0.37 The resulting risk figures for estimation of a 90% effect in successively smaller propor- tions (1/10 [p= 0.1], 1/100 [p= 0.0I], and 1/1000 [p= 0.001]) of the population based on cocaine-maintained responding were Environmental Health Perspectives * Vol 104, Supplement 2 - April 1 996 tm .s a 0 CL U) 9) -6 r 0 L) 0 a a) P CD CL 393 J.R. GLOWA In dose Dose, mg/kg The mean par Probabllity In addition, the effects of GBR 12909 on 0.1 0.01 0.001 cocaine-maintained responding may also -0.652 -2.321 -3.076 be similar to those of delivering free food 0.19 0.10 0.05 on behavior maintained by food delivery. ameters for food-maintained However, in contrast to the decrease in responding were FRF cocaine-maintained behavior produced by GBR 12909, the decrease in food-main- a b 90% (in) 90%, mg/kg tained behavior by this drug was consid- 170.45 -134.6 0.5579 2.02 ered adverse (i.e., an effect that might (0.6258) interfere with the normal behavioral func- tioning). The variability in the likelihood of producing each effect in individual ani- tulting risk figures for estimation of mals was then used to specify the probabil- effect in successively smaller propor- ity of seeing that effects at doses that would /10 [p= 0.1], 1/100 [p= 0.01], and be expected to affect certain proportions of D [p = 0.001 O]) of the population the population. in food-maintained responding were Further analysis of the current data was Probability used to explore potential therapeutic 0.1 0.01 0.001 indices for GBR 12909. Traditional thera- peutic ratios have been constructed using glkg 0.72 0.27 0.09 comparable levels of effect on different end points (e.g., LD50/ED50). The current approach could extend these methods to mparisons of the distributions of compare the extent to which similar levels o decrease each type of responding of effect on both behaviors overlapped. For y t-tests indicated (p <0.001) the example, the low end of the distribution of ral specificity of the effects of GBR effects for food-maintained responding (i.e., the two behaviors were affected suggests that 1 out of 10 individuals would rent doses). exhibit a 10% decrease in responding when lusions given 0.72 mg/kg GBR 12909. Likewise, 1 out of 100 individuals would be expected irrent report is not intended to to exhibit this effect when given 0.27 tively explore the relationships mg/kg. For cocaine-maintained responding, a different effects of an agent on dif- the high end of the distribution suggests ,nd points but rather to extend the that 9 out of 10 individuals would exhibit ment of a risk assessment technique at least a 10% decrease in drug-seeking L well-defined data set. The current behavior at 0.45 mg/kg. Thus, a dose could Is are of interest because they direcdy be conceived that diminished drug-seeking e sources of variability that have been in at least 90% of the population and with- d in other risk assessment processes. out producing adverse effects in 90 to 99% Lmple, the use of the sample variance of the population. nate risks provides data directly rele- In contrast, it may be of interest to ) one source of uncertainty in risk compare the overlap of a highly therapeutic lent (individual differences). In the dose (i.e., the ED90 for cocaine-maintained iproach this uncertainty factor is typi- responding) with a minimally toxic dose sumed to be an order of magnitude. (i.e., EDIO for food-maintained respond- irrent studies also suggested that, ing) in order to assess the probabilities of :here was a clear difference in the obtaining the most benefit at the least cost. Df GBR 12909 on these two types of Alternatively, it might be argued that only Dr, the variability in effects on each a single dose of GBR 12909 could be given ed relevant data from which risks at one time, focusing the emphasis of such oe assessed. analyses on distributions of two different ! current methods established dose- levels of effect (e.g., a dose of 1.7 mg/kg elationships for GBR 12909 on two decreased drug-maintained responding it behaviors. These dose-effect func- almost completely while having little effect ,uggested that each behavior was on food-maintained responding). Obviously, i differently by this drug, as expected. the extent to which there is no apparent 2909 decreased drug-seeking behav- overlap between the mean desired and ese effects may be similar effects to undesired effects suggests the drug may be of methadone to treat heroin or nico- relatively safe, although a 1/10 chance of an tches to treat tobacco abuse (28,29). adverse effect may warrant concern. The present data suggest that there is some overlap of effects with GBR 12909, although it is by no means complete. It may be constructive to compare the degree of overlap with that of other therapeutic agents (e.g., diazepam) when assessed on relevant behaviors (e.g., punished and nonpunished responding). In the present case, some of this overlap may have resulted from the treat- ment of the data. The slope of the function for food-maintained responding was steeper than that for cocaine-maintained respond- ing, suggesting that the artificial decreases imposed for monkeys that did not exhibit complete decreases in food-maintained responding may have inflated variability. One interesting feature of the current approach is the variability of an effect (and hence, risks estimates) that was almost con- stant over the range of the dose-effect function. For example, examination of the variability at different effect levels (i.e., 90, 50, 10%) on each curve demonstrated a remarkable similarity at different points for food-maintained responding (from 0.639 at 90% to 0.626 at 10%). Variability increased slightly with the increasing effect of GBR 12909 on cocaine-maintained responding (from 0.423 at 90% to 0.724 at 10%). These observations illustrate that dose-toler- ance and benchmark approaches character- ize variability in a completely different manner. Confidence intervals invariably flare at the high and low ends of the dose-effect function, even if the observed variability is greatest at the intermediate doses. One fea- ture of these hyperbolidike functions is that, with sufficient variability in the dose-effect data, this flare may preclude the determina- tion of a benchmark dose because the entire lower confidence interval may never inter- cept the 10% effect level. A feature of quantitative approaches that appears to be gaining acceptance is the use of logarithmic conversions of dose. While log-normal distributions have been used and accepted in pharmacology for many years (33) primarily because receptor interactions can span orders or magnitude, their use in risk assessment has received less attention. Previous studies have pointed to the advan- tage of preventing negative numbers in risk assessment (3,4) and to the validity of log- normal distributions of individual effects when sufficient numbers of subjects are used (27). Standard benchmark approaches may also convert doses to log doses in unex- posed algorithms (2). In the current data this conversion had an unusual effect on the upper end of the ED estimates. When the normally distributed ln point estimates were Environmental Health Perspectives * Vol 104, Supplement 2 * April 1996 Mean (SD) The res a 90% ( tions (1 1/100( based o In dose Dose, mj Cor doses t4 10% b behavio 12909 4 by diffe Conc The ct exhaus betweer ferent e develop using a method measurt assume( For exa to estin vant to assessm RID ap cally as The cu while t effects ( behavic provid could b Thc effect ri differen tions s affected GBR 1 ior. Th the use tine pai 394 NEUROBEHAVIORAL RISK ASSESSMENTS FOR THERAPEUTICS converted back to real numbers, the con- verted distribution dispersed at a logarith- mic rate of change. The result should be that observed effects rarely exceed those predicted at either end of the distribution, although in practice predicted effects at the low end seem to occur with more regular- ity than at the high end. This could dimin- ish some concern of the predicted overlap of these two effects described above. Another feature of the current method that deserves consideration is the use of linear functions to fit individual dose-effect data. There may be several issues including the likelihood that other types of functions may fit the data better. However, previous authors have noted that the type of the function used to fit the linear region usually has little effect on the risk estimate (24,34). In fact, traditional statistical techniques (ANOVA, etc.) may bias the slope of the function by incorporating excessive determi- nations of 0 and 100% effects. Eliminating excessive determination of doses without effect seems desirable, but it requires an arbitrary rule for data exclusion. The use of functions that can effectively use all the data to more accurately fix the linear portion (e.g., logistic curves) should be explored. Linear functions are preferred, of course, because they strongly suggest that a direct relationship exists between dose and effect. While it is obvious that agents may have many effects, the most parsimonious conclusion upon observation of a linear function is that the one that is being mea- sured is directly affected by the agent. Unfortunately, not all dose-effect data can be accurately characterized by a linear func- tion. For example, there is often an increase in responding with some agents when assessed on behavioral measures. These increases inevitably give way to decreases, creating an inverted U-shaped function. Wood and Cox (35) recently reported risk assessments using data in which low doses of toluene increased rates and higher doses decreased rates. They fit a linear function to the ascending limb (rate-increasing effects) alone. While this approach clearly avoids the issues of the interpretation of a complex function and uses the most sensitive effect of toluene, it questions the issue of whether any effect of a toxicant is adverse. In the current studies, effects of GBR 12909 on drug-seeking behavior were the most sensi- tive but were not considered toxic. The current report applies a quantita- tive approach to assess the risks of obtaining two different effects with the same agent. Drug-seeking behavior was targeted for selective elimination by doses of GBR 12909 that were not expected to produce toxic effects. It should be clear from these results that a decrease (or increase) in one behavior alone may not be the best indica- tion of behavioral toxicity. Likewise, selec- tive decreases in food-seeking behavior produced by other agents [e.g., appetite suppressants (36)] may be an intended effect. Literature is incomplete on the dif- ferent behavioral effects that agents can have on behaviors maintained by different events. In the absence of specific knowledge of the effect of an agent on behaviors main- tained by certain events, complete decreases in responding may warrant concern. When a single behavior is studied, it must be con- sidered that doses which decrease respond- ing completely might decrease all types of responding completely, a clear reason for concern. However, when one behavior is more affected than another, it could be because the susceptible behavior is less well maintained. This possibility points to the importance of reinforcement processes in the assessment of behavioral toxicology. In the current case, a decrease in the relatively well-maintained food-presentation behavior is clear reason for concern. The dose-tolerance method was easily adapted to demonstrate a selective effect of GBR 12909 on two different behaviors, as well as to estimate the risks of obtaining those effects in a population. This method should be considered positively over other approaches since it directly employs indi- vidual differences in sensitivity to an agent as a metric of risk. However, not all types of end points can be studied using single-sub- ject, repeated-exposure designs. Recently, we have developed a variant of the dose-tol- erance approach for use with groups-design single-exposure data (Bogdan et al., unpub- lished data). This approach uses an iterating line-fitting program to calculate all possible dose-effect combinations that could be cre- ated from individual effects when different groups of animals are treated with different doses (e.g., three dose groups with six ani- mals each would produce a 6 x 6 x 6 matrix yielding 216 possible dose-effect functions). This approach was very effective when applied to a data set addressing the effects of carbaryl on motor activity of rats (6). The risk figures it produced were less conserva- tive than those of the current (within-sub- ject) approach but were more conservative than those produced by the maximal likeli- hood estimation/confidence-interval approach. Future studies will continue to compare the results of these approaches to provide the practitioner a range of options for quantitative risk assessment. REFERENCES 1. Barnes G, Dourson M. Reference dose (RfD): description and use in health risk assessments. Reg Toxicol Pharmacol 8:471-486 (1988). 2. Crump KS. A new method for determining allowable daily intakes. Fundam Appl Toxicol 4:854-871 (1984). 3. Dews PB. On the assessment of risk. In: Developmental Behavioral Pharmacology (Krasnegor N, Gray J, Thompson T, eds). Hillsdale, NJ:Lawrence Erlbaum, 1986;53-65. 4. Dews PB. Some general problems in behavioral toxicology. In: Neurobehavioraf Toxicology (Annau Z, ed). Baltimore:The Johns Hopkins University Press, 1986;424-434. 5. Dourson ML. New approaches in the derivation of acceptable daily intake (ADI). Comments Toxicol 1:35-48 (1986). 6. Glowa JR, MacPhail RC. Dose-response models for quantitative risk assessment in neurotoxicology. In: Handbook of Neurotoxicology, Vol 3: Approaches and Methodologies (Chang L, Tilson H, eds). New York:Raven Press, 1995; 777-787. 7. McMillan DE. Risk assessment for neurobehavioral toxicity. Environ Health Perspect 76:155-161 (1987). 8. Raffaele KC, Rees DC. Neurotoxicology dose/response assess- ment for several cholinesterase inhibitors: use of uncertainty factors. Neurotoxicology 11:237-256 (1990). 9. Weiss B. Neurobehavioral toxicity as a basis for risk assessment. Trends Pharmacol Sci 91:59-62 (1988). 10. National Research Council. Environmental Neurotoxicology. Washington:National Academy Press (1992). 11. Cote I, Rees DC, Glowa JR. An introduction to the principles of risk assessment. In: The Vulnerable Brain and Environmental Risks, Vol 3: Toxins in Air and Water (Isaacson R, Jensen K, eds). New York:Plenum Press, 1994;183-207. 12. Barrett JE. Effects of alcohol, chlordiazepoxide, cocaine and pentobarbital on responding maintained under fixed-interval schedules of food or shock presentation. J Pharmacol Exp Ther 196:605-615 (1976). Environmental Health Perspectives * Vol 104, Supplement 2 - April 1996 395 J.R. GLOWA 13. Wood RW. Stimulus properties of inhaled substances: an update. In: Nervous System Toxicology (Mitchell CL, ed). New York:Raven Press 1982; 199-212. 14. American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders, 3rd ed rev. Washington:APA Press, 1987. 15. Nestler EJ. Molecular mechanisms of drug addiction. J Neurosci 12:2439-2450 (1992) . 16. Balster RL. Drug abuse potential evaluation in animals. Br J Addict 86:1549-1558 (1991). 17. Aigner T, Balster RL. Choice behavior in rhesus monkeys: cocaine versus food. Science 201:534-535 (1978). 18. Koob GF, Bloom FE. Cellular and molecular mechanisms of drug dependence. Science 242:715-723 (1988). 19. Ritz MC, Lamb RJ, Goldberg SR, Kuhar MJ. Cocaine recep- tors on dopamine transporters are related to self-administration of cocaine. Science 237:1219-1223 (1987). 20. Peck AW, Hamilton M. Psychopharmacology of bupropion in normal volunteers. J Clin Psychiat 44:202-205 (1983). 21. Bergman J, Madras BK, Johnson SE, Spealman RD. Effects of cocaine and related drugs in nonhuman primates: III. Self- administration by squirrel monkeys. J Pharmacol Exp Ther 251:150-155 (1989). 22. Rothman RB. High affinity dopamine reuptake blockers as potential cocaine antagonists: a strategy for drug development. Life Sci 46:PL17-21 (1990). 23. S0gaard U, Michalow J, Butler B, Lund Laursen A, Ingersen SH, Skrumsager BK, Rafaelsen 0. A tolerance study of single and multiple dosing of the selective dopamine uptake inhibitor GBR 12909 in healthy subjects. J Int Clin Psychopharmacol 5:237-251 (1990). 24. Glowa JR, DeWeese J, Natale ME, Holland JJ, Dews PB. Behavioral toxicology of volatile organic solvents: I. Methods: acute effects of to1uene. J Environ Pathol Toxicol Oncol 6:153-168 (1986). 25. Glowa JR. Behavioral toxicity of n-octane. In: Advances in Modern Environmental Toxicology, Vol 6 (Mehlman MA, ed). Princeton, NJ:Princeton Scientific Publishers, 1984;245-253. 26. Glowa JR, Dews PB. Behavioral toxicology of volatile organic solvents: IV. Comparisons of the behavioral effects of acetone, methyl ethyl ketone, ethyl acetate, carbon disulfide, and toluene of the responding of mice. J Am Coll Toxicol 6:461-469 (1987). 27. Glowa JR. Behavioral toxicology of volatile organic solvents: V. Comparisons of the behavioral and neuroendocrine effects among n-alkanes. J Am Coll Toxicol 10:639-646 (1992). 28. Glowa JR. A comparison of dose-effect based risk assessment methods. Neurosci Biobehav Rev 15:153-158 (1991). 29. Rees DC, Glowa JR. Extrapolations to humans for neurotoxi- cants. In: The Vulnerable Brain and Environmental Risks, Vol 3: Toxins in Air and Water (Isaacson R, Jensen K, eds). New York:Plenum Press, 1994;207-230. 30. Glowa JR, Wojnicki FHE, Matecka D, Bacher J, Mansbach RS, Balster RL, Rice KC. Effects of dopamine reuptake inhibitors on food and cocaine maintained responding: I. Dependence on unit dose of cocaine. Exp Clin Psychopharmacol 3:219-231 (1995). 31. Glowa JR, Wojnicki FH, Matecka D, Rice K, Rothman RB. Effects of dopamine reuptake inhibitors on food and cocaine maintained responding: II. Comparisons of effects of psychomo- tor stimulants and subchronic efects. Exp Clin Psychopharmacol 3:231-239 (1995). 32. Wojnicki FHE, Bacher JD, Glowa JR. The use of subcuta- neous vascular access port in rhesus monkeys. Lab Anim Med 44:491-494 (1994). 33. Gaddum JH. Lognornal distributions. Nature 156:463-466 (1945). 34. Crump KS, Allen BC, Faustman EM. The Use of the Benchmark Dose Approach in Health Risk Assessment. Final report. Washington:U.S. Environmental Protection Agency, 1992. 35. Wood RW, Cox CC. A repeated-measures approach to the detection of the minimal acute effects of toluene. Toxicology 6:221 (1986). 36. Corwin R, Schuster CR. Anorectic specificity as measured in a choice paradigm in rhesus monkeys. Pharmacol Biochem Behav45:131-141 (1993). 396 Environmental Health Perspectives - Vol 104, Supplement 2 - April 1996