# Characterization of Toxicokinetics and Toxicodynamics with Linear Systems Theory: Application to Lead-Associated Cognitive Decline.

We present a theoretical approach to analysis of toxicokinetics and toxicodynamics using linear systems theory. In our approach, we define two impulse response functions that characterize the kinetic behavior of an environmental agent in the body and the dynamic time-course behavior of its effect on the body. This approach provides a formalism for understanding the relation among exposure, dose, and cumulative biologically effective dose and for understanding the implications of an effect time-course on cross-sectional and longitudinal data analyses. We use lead-associated cognitive decline as a specific example where the approach may be applied. Key word: bone lead, exposure assessment, linear systems model, toxicodynamics, toxicokinetics. Environ Health Perspect 109:361-368 (2001). [Online 16 March 2001]http://ehpnet1.niehs.nih.gov/docs/2001/109p361-3681inks/abstract.html

Environmental and occupational epidemiologic research involves the identification of relations between past exposures to putative toxicants and subsequent adverse health effects in individuals within study populations. Such relations are often hard to fully characterize because of difficulties in accurately quantifying exposure, dose, and effect. Meaningful quantification is particularly difficult because the exposure--effect relation arises from a multistage process, often referred to as the toxicologic paradigm. As a result, the use of biomarkers in molecular epidemiology research has gained widespread attention (1).

In using biomarkers in environmental epidemiology, it is critical to consider the actual process or parameter a given biomarker reflects. For example, the concepts of exposure, internal dose, and biologically effective dose are often blurred in practice, and the loss of these distinctions can influence interpretation of data. Measures of exposure or internal dose are often assumed to be surrogates of the biologically effective dose; their use implies a set of assumptions that is not usually fully articulated or considered. These assumptions involve considerations of the toxicokinetics of the agent and "exposure" biomarker, which are influenced by varying exposure intensity and duration, the residence time of the active form of the agent at the sensitive target, saturation effects, and release from body stores. Similarly, measures of health effects depend on the time course of response to a given exposure (the toxicodynamics of the agent and "response" biomarker) and are influenced by varying response magnitude and duration, dose-dependent repair mechanisms, and multiagent synergistic effects.

In practice, simplifying assumptions are usually made to directly relate exposure to biologically effective dose and to relate this surrogate of dose to health effect. These simplifying assumptions generally overlook the potentially important influence of bidirectional transfer of the agent from one anatomic or physiologic compartment to another. Less frequently, more sophisticated approaches to relating exposure to internal dose and internal dose to biologically effective dose at the sensitive target have been used. In general, these approaches are based on multicompartmental pharmacokinetic modeling (2-4). Even with sophisticated approaches to characterizing the toxicokinetics of a given agent--effect paradigm, simplified models of response are usually used in epidemiologic investigations (e.g., the assumption of an irreversible, static response to a given exposure). To help understand these issues, we present a conceptual framework based on linear systems theory and its application to the analysis of lead-associated neurocognitive decline. We specifically consider issues of residence time of the agent at the sensitive target, later release of the agent from body stores (with corresponding re-residence at the sensitive target), and the time course of response.

Theory

A fundamental assumption in dose--response studies of toxic agents is that the active form of the agent at the sensitive target site causes the effect (4). Thus, to characterize the relation between exposure to an environmental agent and subsequent development of an adverse health effect, two types of processes must be modeled: the toxicokinetics that describe the relation between environmental exposure and ultimate cumulative dose at the sensitive target, and the toxicodynamics that describe the relation between this cumulative dose at the sensitive target and the adverse effects. In practice, it is frequently assumed that the total adverse effect is proportional to the area under the curve ([AUC.sub.T]) of the time-concentration relation for the active form of the agent at the target site, (4):

[1] [AUC.sub.T] = [integral of] [C.sub.T] (t)dt

Thus, it is highly desirable in environmental epidemiologic investigations to be able to estimate [AUC.sub.T] and to be able to characterize the relation between [AUC.sub.T] and the observed effect at any given measurement time.

Toxicokinetics: Measures of [AUC.sub.T]

A common surrogate for [AUC.sub.T] is cumulative exposure, E, defined as the integral exposure to a certain time-dependent environmental concentration of a toxic agent, [C.sub.E](t):

[2] E = [integral of] [C.sub.E](t)dt

Another common surrogate for [AUC.sub.T] is the cumulative dose, D, which is frequently based on the assumption of a linear relation between exposure and dose:

[3] D(t) = kE(t) D = k [integral of] [C.sub.E](t)dt

The use of either cumulative exposure or cumulative dose as a surrogate for [AUC.sub.T] implies the assumption of a linear relation between [AUC.sub.T] and either E or D, with y-intercept equal to zero:

[4] [AUC.sub.T] =[k.sub.E]E = [k.sub.D]D

In other words, Equation 4 implies that the toxicokinetics are strictly linear; that is, the transfer of a toxic agent from the environment to the sensitive target (including bioactivation, if relevant) follows a linear dependence. Of perhaps greater importance, Equation 4 implies that the residence time of the agent at the sensitive target is relatively brief, and that, if the agent is concentrated in body stores, it is not subsequently released and made bioavailable to the sensitive target. This conclusion follows directly from characterization of the transfer of an environmental agent to the body [i.e., from exposure, [C.sub.E](t), to internal dose, [C.sub.D](t)] and from the initial biodistribution in the body to the sensitive target [i.e., biologically effective dose, [C.sub.T](t)] with linear systems theory (5), as discussed below. If these assumptions are incorrect, [C.sub.T](t) will be significantly different in shape than either [C.sub.E](t) or [C.sub.D](t) [i.e., [C.sub.T](t) will be spread out compared with either [C.sub.E](t) or [C.sub.D](t)], even in the case of a truly linear relation between exposure and internal dose, and the use of Equation 4 will lead to an underestimation of the actual [AUC.sub.T]. Thus, the common use of the cumulative exposure index [i.e., integrated exposure over time (Equation 2)] (4) may underestimate the [AUC.sub.T].

For the sake of our discussion, a system is considered linear if the relation between the input [e.g., [C.sub.E](t)] and the output [e.g., [C.sub.T](t)] has the properties of additivity and scaling. A system has the property of additivity if the sum of the outputs from two independent inputs equals the unified output from the sum of the inputs. That is, if the function f that describes a system is linear, then:

f[[x.sub.1](t) + [x.sub.2](t)] = f[[x.sub.1](t)] + f[[x.sub.2](t)].

In our case, this implies that the sum of the two actual time--concentration curves at the sensitive target resulting from two separate exposures would be the same as the single time--concentration curve resulting from the sum of the exposures. That is, we assume that a single, complex time-varying exposure may be conceptualized as a series of intensity-scaled instantaneous exposures, and that the actual observed time--concentration curves resulting from the single, complex time-varying exposure can be modeled as the sum of individual time--concentration curves from the series of intensity-scaled instantaneous exposures.

A system has the property of scaling if the output is scalable by the input:

f[ax(t)] = af[x(t)].

In our case, this implies that the time--concentration curve at the sensitive target multiplicatively scales with a multiplicative change in exposure. For example, if the exposure intensity doubles, we assume that the time--concentration curve will double (i.e., the concentration will double, and the curve will retain the same shape).

Additivity and scaling are usually combined into the principle of superposition. In general, biological systems follow this super-position principle, up to the point of saturation or other mass effects at very high concentrations.

An equally important consideration for a linear system is time invariance. In a time-invariant system, a time shift or delay of the input produces a corresponding shift in the output, without any other change. In our case, this implies that the toxicokinetics do not change over time (i.e., the same exposure would always produce the same time--concentration curve). This assumption is frequently made in assessments of environmental concentrations of toxicants (6), but has not always explicitly been applied to toxicokinetic behavior. In practice, this assumption implies that ongoing exposure does not influence the toxicokinetics per se. Such is the case when the adverse effects (i.e., the toxicodynamics) do not influence physiologic processes that govern the toxicokinetics. For example, if lead does not influence blood flow or organ extraction from blood, chronic exposure does not change the biodistribution of an additional dose relative to the first dose. As a contrasting example, initial exposure to a respiratory irritant may influence airway caliber or mucociliary clearance, thus changing the toxicokinetics of subsequent exposures. For the example given below (lead-associated neurocognitive decline), we assume that the system is time invariant.

Given a linear, time-invariant system, the output [e.g., [C.sub.T](t)] for any given arbitrary input [e.g., [C.sub.E](t)] can be directly predicted from knowledge of the input and the system's impulse response function (IRF). This IRF characterizes the system's response (i.e., the output) to an infinitely short duration input (mathematically equivalent to a delta or Dirac function in time, denoted as [Delta]) (5):

[5] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

From a toxicokinetic perspective, [IRF.sub.TK] is the [C.sub.T](t) curve that would be observed from a single, infinitely short duration exposure. Because [IRF.sub.TK] is a time-varying curve and is conceptually derived from a single pulse exposure, its use implies that the ongoing toxicokinetics of even a single molecule of a toxicant could be captured by this formalism. For example, [IRF.sub.TK] could represent the behavior of a single lead molecule that enters the brain (and produces an effect), is cleared and stored in bone, is then released back to the blood, and reenters the brain (and produces a second effect).

Assuming that the toxicokinetics can be described by a linear, time-invariant system, the observed target site concentration--time curve from an arbitrary exposure time course is given by the mathematical convolution (denoted by [cross product]) of the actual exposure time course with the target site [IRF.sub.TK](5):

[6] [C.sub.T](t) = [k.sub.E][C.sub.E](t) [cross product] [IRF.sub.TK](t),

where [C.sub.E](t) is the exposure time--concentration curve, and [k.sub.E] is a constant relating units of exposure to units of biologically effective dose. Equation 6 may be recast as the standard convolution integral:

[7] [C.sub.T](t) = [k.sub.E] [integral of] [C.sub.E]([Tau])[IRF.sub.TK](t-[Tau])d[Tau],

where [Tau] is a dummy variable of integration. Of major importance, if either of the two terms in the convolution [e.g., [C.sub.E](t) or [IRF.sub.TK] Equations 6 and 7] is a [Delta] function, then that term drops out of the convolution, and the relation reduces to a straight equivalency. That is, if the residence time of a toxicant at the sensitive target is very short, and if there is no subsequent bioavailability due to release from body stores, then [IRF.sub.TK] is essentially a [Delta] function, and the assumptions implicit in Equation 4 are valid. On the other hand, if the residence time is significant, then [IRF.sub.TK] is not a [Delta] function, but a curve with some spread in time, and Equation 7 must be used (with Equation 1) instead of Equations 2 and 4 to determine [AUC.sub.T]. In a similar fashion, if there is significant subsequent release from body stores, [IRF.sub.TK] will be multi-peaked or have an initial peak (representing the initial transfer from the environment), followed by non-zero values over time (representing the release from body stores), and Equation 7 must again be used to compute [AUC.sub.T] via Equation 1.

In essence, one extreme scenario for the IRF and Equations 6 and 7 is that [IRF.sub.TK] is a [Delta] function, implying that the residence time is essentially zero (i.e., the active agent rapidly transits through the sensitive target) and that no release from body stores occurs. The opposite scenario from a toxicokinetic perspective is that the active agent permanently resides in the sensitive target; that is, [IRF.sub.TK] is a constant with time (either because the agent never clears from the sensitive target or because ongoing significant biorelease from body stores constantly replenishes that amount of agent cleared from the target). In general, Equations 1 and 7 can be combined, with [AUC.sub.T] expressed as a double integral, with integration limits from time zero to the current observation time, T:

[8] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

If [IRF.sub.TK] is a [Delta] function, the inner integral drops out, and Equation 8 reduces to that relation implied by Equations 3 and 4:

[9] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

If [IRF.sub.TK] a constant function, it drops out of the inner integral in Equation 8, but the double integration remains:

[10] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

[Note that [IRF.sub.TK] is causal, which means that it results from an event in real time; thus, [IRF.sub.TK] (t [is less than] 0) = 0. This assumption facilitates the simplification of Equation 8 to either Equation 9 or Equation 10.]

Equations 9 and 10 define opposite scenarios. If the actual [IRF.sub.TK] is closer to a [Delta] function than a constant, the use of Equation 9 in environmental epidemiologic investigations should provide biologically effective dose estimates that are better predictors of health outcomes than would the use of Equation 10. If the actual [IRF.sub.TK] is significantly spread, either because of significant residence time or significant release over time from body stores, then the use of Equation 10 should provide better biologically effective dose estimates, and thus stronger and more consistent associations with the health outcomes under study.

Toxicodynamics: Time-Dependent Measures of Response

In the previous section, we used linear systems theory to characterize the time course of the active form of the agent at the sensitive target, [C.sub.T](t). We now want to characterize the time course of the health outcome or response, R(t). As with [C.sub.T](t), R(t) can be conceptualized as a linear system characterized by an IRF. Since we designated the IRF for the toxicokinetic relation [IRF.sub.TK], a similar IRF, designated [IRF.sub.TD], can be used to characterize the toxicodynamics:

[11] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

From a toxicodynamic perspective, [IRF.sub.TD] is the R(t) curve that would be observed from a single, infinitely short duration [C.sub.T](t) curve [i.e., in the case where both [C.sub.E](t) and [IRF.sub.TK] are (5 functions]. The shape of [IRF.sub.TD] directly indicates reversible versus irreversible (persistent) versus progressive effects from a single exposure, as shown in Figure 1. A reversible effect would yield an [IRF.sub.TD] that goes from zero at the point of exposure to some response value and back to zero. An irreversible effect would yield an [IRF.sub.TD] that goes from zero at the point of exposure to some response value that persists, independent of time. A progressive (increasing) effect (from a previous exposure) would yield an [IRF.sub.TD] that continuously increases with time, starting from zero at the point of exposure.

[ILLUSTRATION OMITTED]

The utility of Equation 11 is that it identifies issues in the design of an epidemiologic study in which response is measured either cross-sectionally at some time after exposure ceases, or longitudinally. For example, a cross-sectional design uses a surrogate of cumulative exposure and a single measure of response at some later time. Equation 11 tells us that such a design is only applicable if [IRF.sub.TK] is a [Delta] function and if [IRF.sub.TD] is a step function (i.e., the effect is persistent rather than either reversible or progressive). In a longitudinal study, it is imperative to consider possible shapes for [IRF.sub.TD], because the observed R(t) implies underlying [IRF.sub.TK] and [IRF.sub.TD] functions that are usually not directly obtainable. For example, the observation of a progressive increase in effect with time could be the result of a) the ongoing presence of active agent at the sensitive target, due to either long residence time or ongoing release from body stores, either of which would cause [AUC.sub.T] to continue to increase with time; b) progressive response from past exposures; or c) a combination of the two. Equation 11 suggests that the observed R(t) could be the result of a significantly non-[Delta] function IRF for either the toxicokinetic or toxicodynamic portion of the toxicologic paradigm. For example, from a purely mathematical point of view, the same R(t) would be observed if either [IRF.sub.TK] or [IRF.sub.TD] were a [Delta] function and the other were a step function.

Application to Longitudinal Data Analysis

Bandeen-Roche et al. (7) have previously described a general model for analysis of data from prospective observational studies with multiple outcome measures over time. Their model includes a family of exposure summaries whose mathematical formalism is a convolution integral similar to Equation 7, although they did not approach the overall model's conceptualization with linear systems. In their data analysis model, outcome or response is a function of exposure history plus a random error. With our notation:

[12] R(t) = F[[C.sub.T](t)] + [Epsilon](t)

In their approach, F[.] represents a complex function representing the effects of the family of exposure summaries. They use this construct to create a generalized linear model that separates cross-sectional, longitudinal, historic, and regression-to-the-mean effects. Here, "cross-sectional," "longitudinal," and "regression-to-the-mean" have their convention. al epidemiologic definitions, and "historic" addresses the influence of previous exposures on the toxicodynamics resulting from subsequent exposures. If we use separate [Beta] coefficients to designate cross-sectional ([[Beta].sub.cs]), longitudinal ([[Beta].sub.l]), historic ([[Beta].sub.h]), and regression-to-the-mean ([[Beta].sub.rm]) terms in the model, then initial (i.e., cross-sectional) and subsequent (i.e., longitudinal) relations can be defined (with our notation):

[13] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

It is common in data and statistical analyses of longitudinal data to use the difference between successive measurements. In such a case, the model above reduces to:

[14] [Delta]R = [R.sub.2] - [R.sub.1] = [Beta]* + [[Beta].sub.1] ([AUC.sub.2] - [AUC.sub.1]) + [[Beta].sub.h][AUC.sub.1] + [B.sub.rm][R.sub.1] + [Epsilon]*

This equation is important in generalizing the use of the linear systems approach we propose in longitudinal data analysis. Such analyses go directly to the heart of the toxicokinetic and toxicodynamic implications of the IRF model we have invoked.

Empirical Validation: Cognitive Effects of Lead

Lead is widely recognized as a significant neurotoxicant, and the development of biomarkers of lead exposure has been vigorously pursued. X-ray fluorescence (XRF) measurement of lead in bone has been adopted as the method of choice to assess cumulative exposure (8,9) because lead in blood has a clearance half-time of 30 days, whereas lead in bone has a clearance half-time of 15-30 years. XRF measures of bone lead highly correlate with the integral of the time-course of blood lead concentration [also called the "cumulative blood lead index" (10,11)]. Independent data suggest an association between cumulative blood lead level and cumulative brain uptake (12-13). Because multiple blood-lead level measurements as a function of time, which are necessary to compute the cumulative index, are usually not available in epidemiologic studies and only rarely in occupational studies, single XRF measures of lead in bone are taken to represent "cumulative exposure" (via Equation 2) or "cumulative dose" (via Equation 3).

Hu et al. (14) considered two paradigms for the interpretation of skeletal lead, as measured by XRF: bone lead as an indicator of cumulative lead exposure, and bone lead as a source of body lead burden that can be mobilized into the circulation. The first paradigm considers bone lead as a surrogate marker for cumulative dose to sensitive targets, whereas the second considers bone lead as an important endogenous source of further exposure. With either paradigm, the dose--response relation could be linear or nonlinear and involve a threshold or not (14). In either case, XRF measures of bone lead could be predictive of a given health outcome (such as cognitive decline) as long as a strictly linear relation exists between the XRF measure and [AUC.sub.T]. As discussed above, if either the residence time (assuming the first paradigm) or release from bone stores (assuming the second paradigm) is significant, then the XRF measure, although clearly better than a blood measure, would not correlate as highly with health outcome as a truer [AUC.sub.T] metric. For example, a comparison of the strength of association with health outcomes between [AUC.sub.T] estimated from Equation 9 versus Equation 10 would provide some evidence for which of the two scenarios (short residence time and no release versus long residence time and/or subsequent release) is more likely. This distinction is particularly useful in making hypotheses about whether an effect is likely to be transient or persistent or progressive.

Methods

To empirically investigate these scenarios, we reanalyzed data from a longitudinal study of 535 former organolead manufacturing workers, for whom we have already reported results (15-19). Informed consent was obtained before a subject was enrolled in the study. As reported in more detail elsewhere (19), a battery of 19 cognitive tests was obtained annually (Table 1). The results of this battery were compared with blood and bone lead measurements. Current tibial lead was measured via XRF and used to estimate the peak tibial lead value at the time of cessation of occupational exposure. To do so, the clearance of lead from bone was modeled with a mono-exponential function, as has previously been demonstrated to fit longitudinal bone lead data (20); the clearance half-time was assumed to be 27 years (21). The [AUC.sub.T] was estimated in three ways: a) by assuming that the current tibial lead was proportional to [AUC.sub.T], b) by assuming that the peak tibial lead was proportional to [AUC.sub.T], and c) by forming and integrating the estimated tibial lead time course (which we designate AUC'). To do the latter, we back-extrapolated from the current tibial lead value to the time at which exposure ended (i.e., the same process used to get the peak tibial lead value), and then assumed a straight line between that peak value and a value of zero at the start of occupational exposure some years earlier, as shown in Figure 2. (We knew the date of start of occupational exposure and cessation of exposure for each subject.)

[ILLUSTRATION OMITTED]

Table 1. Generalized estimating equation linear modeling results identifying predictors of annual change in neurobehavioral test scores in 535 former organolead manufacturing workers comparing four different dose metrics, 1994-1998. Neurobehavioral measure(a) (used in separate regression models of Current PbB(b) Current TL(c) change) [Beta] (SE [Beta]) [Beta] (SE [Beta]) Block design (Wechsler 0.161 (0.135) -0.058 (0.160) Adult Intelligence Scale) Digit symbol (Wechsler -0.187 (0.126) -0.133 (0.135) Adult Intelligence Scale, revised) Symbol digit -0.012 (0.102) -0.099 (0.099) Serial digit learning -0.160 (0.132) -0.020 (0.153) Rey complex figure, copy -0.001 (0.097) -0.030 (0.091) Rey complex figure, -0.174 (0.088)(**) -0.077 (0.090) delayed recall Rey auditory verbal 0.049 (0.167) -0.255 (0.175) learning test, immediate recall, 5 trials Rey auditory verbal -0.055 (0.055) -0.128 (0.054)(**) learning test, delayed recall Rey auditory verbal -0.009 (0.053) 0.035 (0.062) learning test, recognition Trails A -0.225 (0.252) -0.586 (0.273)(**) Trails B -0.790 (0.629) 0.285 (0.740) Finger tapping, dominant -0.196 (0.124) -0.260 (0.155)(*) hand Finger tapping, -0.137 (0.099) -0.224 (0.128)(*) nondominant hand Pegboard, dominant hand 0.045 (0.087) -0.133 (0.092) Pegboard, nondominant hand -0.093 (0.087) -0.254 (0.100)(**) Pegboard, both hands 0.053 (0.074) -0.048 (0.095) Pegboard assembly 0.461 (0.235)(**) -0.034 (0.329) Stroop (C form - A form) -0.740 (0.389)(*) -0.676 (0.544) Choice reaction time -0.393 (2.541) -0.298 (2.703) average SIGNS of [Beta] 14/19 negative 20/22 negative coefficients Statistical significance -- -- 2 < 0.05 (1+) 3 < 0.05 1 < 0.01 2 < 0.10 Neurobehavioral measure(a) (used in separate regression models of Peak TL(d) AUC-lead(e) change) [Beta] (SE [Beta]) [Beta] (SE [Beta]) Block design (Wechsler -0.223 (0.165) -0.213 (0.173) Adult Intelligence Scale) Digit symbol (Wechsler -0.038 (0.140) 0.058 (0.139) Adult Intelligence Scale, revised) Symbol digit -0.206 (0.107)(**) -0.225 (0.108)(**) Serial digit learning -0.104 (0.161) -0.149 (0.165) Rey complex figure, copy -0.138 (0.094) -0.170 (0.096) Rey complex figure, -0.174 (0.096)(*) -0.196 (0.097)(**) delayed recall Rey auditory verbal -0.571 (0.193)(***) -0.671 (0.204)(***) learning test, immediate recall, 5 trials Rey auditory verbal -0.149 (0.058)(**) -0.134 (0.058)(**) learning test, delayed recall Rey auditory verbal -0.028 (0.072) -0.054 (0.078) learning test, recognition Trails A -0.503 (0.308) -0.451 (0.309) Trails B -0.061 (0.851) -0.212 (0.881) Finger tapping, dominant -0.170 (0.153) -0.056 (0.150) hand Finger tapping, -0.169 (0.135) -0.094 (0.137) nondominant hand Pegboard, dominant hand -0.138 (0.098) -0.122 (0.098) Pegboard, nondominant hand -0.250 (0.101)(**) -0.182 (0.102)(*) Pegboard, both hands -0.094 (0.087) -0.117 (0.086) Pegboard assembly -0.320 (0.320) -0.577 (0.301)(*) Stroop (C form - A form) -1.122 (0.606)(*) -1.366 (0.675)(**) Choice reaction time -2.121 (2.776) -2.968 (2.989) average SIGNS of [Beta] 22/22 negative 21/22 negative coefficients Statistical significance 1 < 0.01 1 < 0.01 3 < 0.05 4 < 0.05 2 < 0.10 2 < 0.10 (a) Adjusted for age, education, visit number, testing technician, and baseline score; tests were standardized for direction so that a negative coefficient indicates worsening performance with increasing blood or tibial lead. Beta coefficients are standardized so that they can be directly compared for each neurobehavioral test. The units of each [Beta] coefficient indicate change in neurobehavioral test score per SD unit increase in the lead biomarker. (b) Current PbB = current blood lead level. (c) Current TL = current tibial lead. (d) Peak TL = peak tibial lead, estimated from current tibial lead and years since last exposure, using an estimated half-time of lead in tibia of 27 years (see "Methods"). (e) AUC-lead = area under the curve of estimated tibial lead levels versus time (see "Methods"). (***) p-value s 0.01; (**) p-value < 0.05; (*) p-value < 0.10.

If the XRF measurement can be conceptualized via Equation 2 (i.e., the XRF value at a given time T represents the integral exposure from time zero to that time T), then the [AUC.sub.T] estimated by the third approach above (AUC') represents

[15] [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Equation 15 is thus identical in form to Equation 10 and implies permanent residence of lead at the sensitive target (central nervous system receptors in the brain, in our case). This is also seen by comparison of Equation 15 with Equation 8; the only way that Equation 8 (which is the general form) could be identical to Equation 15 is if [IRF.sub.TK] in Equation 8 is a constant.

In other words, we have four estimates of cumulative biologically effective dose or [AUC.sub.T]:

1. Current blood-lead level, which reflects both past exposures with a 30-day clearance half-time (i.e., it is a poor index of cumulative exposure) and the reintroduction of lead into the circulation via release from bone stores. This metric assumes a short residence time in brain and is mainly reflective of present release from bone and other stores because our population is no longer occupationally exposed.

2. Current tibial lead level, which reflects cumulative exposure (with a 27-year clearance half-time) and the magnitude of bone stores potentially available for release. This metric assumes a short residence time in brain and allows for the possibility of some release from body stores; it better reflects cumulative exposure than blood lead because of the much slower clearance.

3. Peak tibial lead level, which reflects cumulative exposure (and is corrected for clearance) and the magnitude of bone stores potentially available for release. This metric also assumes a short residence time in brain and allows for the possibility of some release from body stores.

4. AUC, which reflects cumulative biologically effective dose in the case of permanent residence (either because the initial residence time is long and/or there is significant ongoing release from bone stores).

Comparisons of the association of these four estimates of [AUC.sub.T] with each of the 19 cognitive tests were based on linear regression using generalized estimating equations methodology. Beta coefficients for the four estimates of [AUC.sub.T], for each of the 19 cognitive tests, were obtained and assessed for statistical significance. The cognitive test outcomes were z-transformed before modeling so that the [Beta] coefficients could be directly compared. The linear regression models controlled for age, education, visit number, testing technician, and baseline score on each test.

Two types of linear regression models were used. The first ([R.sub.1] in Equation 13) focused on cross-sectional data from the baseline measurements. Here, the baseline score on each of the 19 cognitive tests were the dependent variables in the linear regression models. The second focused on longitudinal data, starting with the baseline data, and adding three subsequent measurements, each 1 year apart, for each of the 19 cognitive tests. Here, the dependent variables were the annual change scores (via Equation 14), in practice defined as ([R.sub.i] - [R.sub.0])/[Delta] time rather than [R.sub.i] - [R.sub.i-1], and generalized estimating equations methods were used to examine associations of the lead measures with change in test scores over time.

Results

Cross-sectional analyses. We previously reported the associations of current and peak tibial lead levels with 19 neurobehavioral test scores (16); we now report associations for blood-lead level and AUC'. Taking a p-value [is less than] 0.65 as significant, current blood lead level was significantly associated with 4/19 tests, current tibial lead level with 9/19, peak tibial lead level with 11/19, and AUC' with 14/19. All significant [Beta] coefficients indicated that increasing lead levels were associated with lower neurobehavioral test scores.

Longitudinal analyses. Table 1 shows the results of the linear regression analyses using the baseline values for blood lead, current tibial lead, peak tibial lead, and AUC'. Taking a p-value [is less than] 0.05 as significant, current blood lead level was associated with 2/19 neurobehavioral test change scores (although one of the two significant [Beta] coefficients was in the opposite direction than expected), current tibial lead level with 3/19, peak tibial lead level with 4/19, and AUC' with 5/19. Of interest, current blood lead level had the highest [Beta] for 3/19 associations, current tibial lead level for 4/19 associations, peak tibial lead level for 3/19 associations, and AUC' for 9/19 associations. There was only one instance where AUC' did not have a significant [Beta] coefficient when one of the other estimates did; in this case, current tibial lead level produced the only [Beta] coefficient that achieved statistical significance.

Because the cognitive test outcomes were standardized, both the direction and magnitude of the associations ([Beta] coefficients) of each [AUC.sub.T] lead measure with each neurobehavioral test change score could be directly compared. Current blood-lead level had only two significant [Beta] coefficients, one of which was positive (i.e., in the opposite direction than expected). Use of current blood level would thus lead to the conclusion of no association between lead and cognitive decline. In contrast, progressing from current tibial lead level to peak tibial level to AUC' increased the number of significant [Beta] coefficients, and AUC' produced the largest [Beta] coefficient in more of the 19 tests than any of the other 3 measures.

Even though the [Beta] coefficients changed depending on which of the four estimates of AUC was used, some association between the four estimates was present. Pearson's correlations were significant (p [is less than] 0.01) between a) blood and current tibial lead (r = 0.44), peak tibial lead (r = 0.26), and AUC' (r = 0.18); b) current tibial lead and peak tibial lead (r = 0.86), and AUC' (r = 0.70); and c) peak tibial lead and AUC' (r = 0.94).

A critical issue in the interpretation of the longitudinal test outcomes is whether the change in cognitive performance over time could be completely (and thus solely) explained by the increase in AUC' with time (since the limit of integration for any AUC measure progressively increases with time). In other words, is any progressive cognitive decline simply the result of progressive cumulative dose? Accordingly, we also evaluated a generalized model in which [Delta]-AUC' ([Delta]-AUC' = AUC' end-of-interval -- AUC' baseline) was used instead of baseline AUC'. With this model, only 2/19 tests had significant beta coefficients, and the change in cognitive test outcomes could not be explained by the change in AUC alone. This suggests that a progressive model for [IRF.sub.TD] needs to be considered.

Discussion

In epidemiology, assessing the association between exposure to a putative toxicant and subsequent health outcome implies the existence of an underlying biologically based dose--response relation. The goal of exposure or internal dose assessment is thus to find an index that best represents the cumulative biologically effective dose of the active form of the toxicant at the sensitive target. In practice, this index should be proportional to the integral of the time course of the concentration of active agent at the sensitive target (i.e., Equation 1), designated [AUC.sub.T]. The use of cumulative exposure or internal dose as a surrogate for this time--concentration integral will only prove useful if the toxicokinetics are approximately linear over the concentration range expected, and if the effects are approximately cumulative (4,6). This is simply a restatement of Haber's rule: tissue damage should be related to the product of the mean exposure intensity and time (6). Thus, a major issue in the development and use of biomarkers is the degree to which these assumptions hold. In other words, it is important to understand whether the toxicokinetics and toxicodynamics are linear and time-invariant mathematically, and whether the toxicodynamics represent a reversible, persistent, or progressive process.

In this report, we present a conceptual framework based on linear systems theory as an aid to identifying and considering these issues. With respect to toxicokinetics, we have attempted to relate different surrogates of cumulative biologically effective dose and to identify the conditions under which certain assumptions are implicitly invoked. We have introduced the common linear systems concept of an impulse response function, [IRF.sub.TK], to describe the toxicokinetics following an infinitely short duration exposure. Use of linear systems theory and this concept allows us to define a general relation between the exposure time-course and the time-concentration curve of the active form of the agent at the sensitive target, whose integral, [AUC.sub.T], is likely best correlated with response. This general relation (Equation 8) could be transformed into two more specific relations, one representing the case where the residence time in the sensitive target is infinitely short and no biorelease occurs (Equation 9), and another representing the case where the effective residence time is infinitely long (either because the agent never dears from the sensitive target or because ongoing significant biorelease from body stores constantly replenishes that amount of agent cleared from the target; Equation 10).

We also sought to use linear systems theory to conceptualize the time-course of response, which was particularly important for the longitudinal data presented here, in which we studied the association of four estimates of [AUC.sub.T] with cognitive decline in an occupationally exposed cohort of 535 workers. We found that progressively more [Beta] coefficients were statistically significant as we moved from current blood lead level to current tibial lead level to peak tibial lead level to AUC' in both cross-sectional and longitudinal analyses.

Conceptualizing [AUC.sub.T] in this way and generating multiple estimates of [AUC.sub.T] is helpful in several ways: it clarifies the distinctions between exposure, internal dose, and biologically effective dose; it guides the development of different estimates of [AUC.sub.T]: and the results of a comparison of the association of these different estimates of [AUC.sub.T] with health outcomes provide indirect evidence of the underlying biological phenomena. For example, for those health outcomes in which a measure of recent or current dose has the highest association, the health outcome is likely an acute, reversible process, whereas for those health outcomes in which a measure of cumulative dose has the highest association, the health outcome is likely a persistent or progressive accumulative process. For those health outcomes in which a measure that includes consideration of residence time has the highest association, the agent's toxicokinetics likely include significant residence in the sensitive target and/or significant release from body stores.

Significant successes and utility have been reported with the use of pharmacokinetic modeling in predicting the toxicokinetics of environmental agents such as lead (4,22). Our approach differs from these previous mathematical efforts in that it is "model free" (i.e., it does not assume a certain "topology" or relationship among a series of anatomic or physiologic compartments, as do classical pharmacokinetic modeling approaches). Rather, our approach makes use of linear systems theory to describe the IRF of the system nonparametrically. This nonparametric IRF exists independent of compartment-based descriptions. In other words, there are no predetermined parameters whose presence, number, and character are fixed by an a priori hypothetical model.

Having stated this, we emphasize that we are not asserting that our approach is intrinsically better than multicompartment modeling; rather, we view the two approaches as highly complementary. For example, both may lead to a useful prediction of [AUC.sub.T], but they require different assumptions and independent data. The choice of approach depends on prior knowledge and on the types of information and relations desired. If compartments and relations can be identified and appropriate quantitative rate constants determined, then the multi-compartment toxicokinetic approach yields accurate predictions (4,22). On the other hand, when less is known about the potential compartments, and particularly when little is known about rate constant values, the conceptualization of the toxicokinetics via a linear system with an IRF may prove useful. Such is the case when an empirical time course in an organ or structure of interest is already known (23). In this regard, the empirical time course need not come from a [Delta] input function because deconvolution analysis can be used to obtain [IRF.sub.TK] from the combination of any arbitrary but known input function and the empirical time course (5). Once [IRF.sub.TK] is obtained in this way, the time-course for any other arbitrary but known input function can be predicted via Equation 6.

We also emphasize that our conceptual framework covers both toxicokinetics and toxicodynamics, whereas pharmacokinetic modeling only addresses predictions of toxicokinetics. In this regard, we highlight the complementary nature of pharmacokinetic modeling and linear systems analysis by suggesting that pharmacokinetic modeling can be used, when available, to predict [IRF.sub.TK] for subsequent use in linear systems analysis of the toxicodynamics. In essence, either pharmacokinetic modeling or our approach could be used to predict the needed kinetic time-course if the right data are in hand; it is important that the types of data are strikingly different. In practice, we used the toxicokinetic portion of our conceptual framework to create four different estimates of [AUC.sub.T]; we do not claim that any of these estimates is more accurate than one obtained from pharmacokinetic modeling. Rather, the conceptual framework provides a different understanding of the meaning of each of these estimates than one derived from pharmacokinetic modeling, and this understanding helps in the interpretation of the actual data.

The conceptual framework as presented, embodied in specific equations, and used so far, requires that the toxicokinetics and toxicodynamics be linear and time invariant. With respect to toxicokinetics, when the active form of the agent is produced by metabolism, nonlinear effects, especially at high doses, are expected as the processes saturate. Nonlinear kinetics could also arise from changes in individual uptake or susceptibility with time, from synergistic or antagonistic effects related to concurrent mixed exposures to other toxicants, from allergic responses, from changes produced by the initial exposure to that agent (e.g., upregulation of cytochrome P450), or from dose-rate effects (6). With respect to toxicodynamics, for stochastic processes like carcinogenesis, dose-response relations may be linear or nonlinear. In the case of a direct genotoxic carcinogen, a linear or linear-quadratic relation between [AUC.sub.T] and response is expected (4). For nonstochastic processes, a linear relation with a threshold is commonly observed. As with the toxicokinetics, up- or downregulation of receptors or tolerance effects can introduce nonlinearities. To the extent to which strict linearity is not present, the power of the conceptual framework and actual approach decreases; the degree of linearity may be different for the toxicokinetic and toxicodynamic portions of the analysis in a given application.

In contrast to the requirement for linearity, time invariance is not strictly required in our approach. For the sake of simplicity of presentation and implementation, we have invoked the assumption of time invariance; this assumption is what leads to the specific convolution integral given in Equation 7, and repeated below:

[17] [C.sub.T](t) = [k.sub.E] [integral of] [C.sub.E]([Tau])[IRF.sub.TK](t-[Tau])d[Tau]

If [IRF.sub.TK] is time varying, the convolution integral given in Equation 7 must be modified, as follows:

[7'] [C.sub.T](t) = [k.sub.E] [integral of] [C.sub.E]([Tau])[IRF.sub.TK](t; [Tau])d[Tau]

Equation 7' differs in a subtle way from Equation 7: [IRF.sub.TK] is now a function of both t and [Tau], not just t-[Tau]. An analogous situation holds for [IRF.sub.TD]. This implies the need to obtain a family of IRFs (as a function of t) rather than just a single IRF. In practice, if the IRF changes only slowly, then it becomes possible to treat segments of time as being time invariant and to use a single IRF during that time period.

In essence, our approach has permitted initial comparisons of different assumptions about the residence time of lead in brain, release from bone stores, and the persistence or progression of lead-associated neurobehavioral effects. Using both cross-sectional and longitudinal data analysis, we have been able to show, in preliminary form, that a) a measure of cumulative lead dose (AUC') that implies either long residence time of lead in brain or significant ongoing release of lead from body stores is the best predictor of both test scores at cross-section and test score declines over time; b) the change in this AUC metric over time is a poor predictor of longitudinal test score change; and c) the observed longitudinal change in test scores is consistent with a model of progressive neurobehavioral effect. Given our knowledge about the clearance half-time of lead in brain and the current blood and tibial lead levels in the former workers, we believe that the observed annual test score declines are likely due to a combination of newly induced effects from lead released from bony stores and, more significantly, progressive effects from past exposures to lead. In any event, the model requires that, at a minimum, effects persist for many years past the exposure that triggered the effect. These data thus support the hypothesis that this is not a transient neurochemical effect, that, by necessity, would depend on the continued presence of lead in brain to sustain the effect, but rather a persistent structural change (which may have been initially triggered by neurochemical events).

In general, elimination of the toxicant from the sensitive target site probably represents the rate-limiting step between exposure and response (6). For example, Rappaport (6) has alluded to the physiologic damping that "resulted from accumulation of lead over several months owing to the slow rate of elimination and distribution of this metal from the blood." If this elimination is slow, the accumulated burden is large relative to the amount of toxicant received (6). In such a case, knowledge of [AUC.sub.T] itself (compared with either integrated exposure or internal dose) will better predict response because even large short-term fluctuations in exposure will not directly provide important information on biologically effective dose.

Ultimately, we are interested in the development, validation, and application of biomarkers whose behavior we understand vis-a-vis the toxicologic paradigm. In this regard, the use of the IRF has already been described in noninvasive imaging (24) and may provide a means to obtain impulse response functions for the toxicokinetics and toxicodynamics of environmental agents of interest, such as lead, which has been radio-labeled (25).

REFERENCES AND NOTES

(1.) Links JM, Kensler TW, Groopman JD. Biomarkers and mechanistic approaches in molecular epidemiology. Annu Rev Public Health 16:83-103 (1995).

(2.) Smith TJ. Exposure-dose relationships. In: Exposure Assessment for Epidemiology and Hazard Control (Rappaport SM, Smith TJ, eds). Boca Raton, FL:Lewis Publishers, 1990;97-114.

(3.) Sato A, Endoh K, Johanson G. The use of models to investigate the toxicokinetic behavior of organic solvents. In: Exposure Assessment for Epidemiology and Hazard Control (Rappaport SM, Smith TJ, eds). Boca Raton, FL:Lewis Publishers, 1990;131-152.

(4.) Smith TJ. Pharmacokinetic models in the development of exposure indicators in epidemiology. Ann Occup Hyg 35:543-560 (1991).

(5.) Oppenheim AV, Schafer RW. Discrete-Time Signal Processing. Englewood Cliffs, NJ:Prentice Hall, 1989.

(6.) Rappaport SM. Assessment of long-term exposures to toxic substances in air. Ann Occup Hyg 35:61-121 (1991)

(7.) Bandeen-Roche K, Hall CB, Stewart WI, Zeger SL. Modelling disease progression in terms of exposure history. Stat Med 18:2899-2916 (1999).

(8.) Hu H, Milder FL, Burger DE. X-ray fluorescence: issues surrounding the application of a new tool for measuring burden of lead. Environ Res 49:295-317 (1989).

(9.) Todd AC, Chettle DR. In vivo X-ray fluorescence of lead in bone: review and current issues. Environ Health Perspect 102:172-177 (1994).

(10.) Hu H, Pepper L, Goldman R. Effect of repeated occupational exposure to lead, cessation of exposure, and chelation on levels of lead in bone. Am J Ind Med 20:723-735 (1991).

(11.) Fleming DEB, Chettle DR, Wetmur JG, Desnick Rd, Robin JP, Boulay D, Richard NS, Gordon CL, Webber CE. Effect of the d-aminolevulinate dehydratase polymorphism on the accumulation of lead in bone and blood in lead smelter workers. Environ Res 77:49-61 (1998).

(12.) Somervaille LJ, Chettle DR, Scott MC, Tennant DR, McKiernan M J, Skilbeck A, Trethowan WN. In vivo tibia lead measurements as an index of cumulative exposure in occupationally exposed subjects. Br J Ind Med 45:174-181 (1988).

(13.) Cremin JD, Luck ML, Laughlin NK, Smith DR. Efficacy of succimer chelation for reducing brain lead in a primate model of human lead exposure. Toxicol Appl Pharmacol 161:283-293 (1999).

(14.) Hu H, Rabinowitz M, Smith D. Bone lead as a biological marker in epidemiologic studies of chronic toxicity: conceptual paradigms. Environ Health Perspect 106:1-8 (1998).

(15.) Schwartz BS, Stewart WF, Todd AC, Links JM. Predictors of dimercaptosuccinic acid chelatable lead and tibial lead in former organolead manufacturing workers. Occup Environ Med 56:22-29 (1999).

(16.) Stewart WF, Schwartz BS, Simon D, Bolla KI, Todd AC, Links d. Neurobehavioral function and tibial and chelatable lead levels in 543 former organolead workers. Neurology 52:1610-1617 (1999).

(17.) Schwartz BS, Stewart WF, Todd AC, Simon D, Links JM. Different associations of blood lead, DMSA-chelateble lead, and tibial lead levels with blood pressure in 543 former organolead manufacturing workers. Arch Environ Health 55:85-92 (2000).

(18.) Schwartz BS, Stewart WF, Kelsey KT, Simon D, Park S, Links JM, Todd AC. Associations of tibial lead levels with Bsml polymorphisms in the vitamin D receptor in former organolead manufacturing workers. Environ Health Perspect 108:199-211 (2000).

(19.) Schwartz BS, Stewart WF, Bolla KI, Simon D, Bandeen-Roche K, Gordon B, Links JM, Todd AC. Past adult lead exposure is associated with longitudinal decline in neurobehavioral function in 535 former organolead manufacturing workers. Neurology 55:1144-1150 (2000).

(20.) Nilsson U, Attewell R, Christoffersson JO, Schutz A, Ahlgren L, Skerfving S, Mattsson S. Kinetics of lead in bone and blood after end of occupational exposure. Pharmacol Toxicol 68:477-484 (1991).

(21.) Gerhardsson L, Attewell R, Chettle DR, Englyst V, Lundstrom NG, Nordberg GF, Nyhlin H, Scott MC, Todd AC. In vivo measurements of lead in bone in long-term exposed lead smelter workers. Arch Environ Health 48:147-150 (1993).

(22.) O'Flaherty EJ, Inskip MJ, Franklin CA, Durbin PW, Manton WI, Baccanale CL. Evaluation and modification of a physiologically based model of lead kinetics using data from a sequential isotope study in cynomolgus monkeys. Toxicol Appl Pharmacol 149:1-16 (1998).

(23.) Widzowski DV, Cory-Slechta DA. Homogeneity of regional brain lead concentrations. Neurotoxicology 15:295-308 (1994).

(24.) Szabo Z, Vosberg H, Sondhaus CA, Feinendegen LE. Model identification and estimation of organ-function parameters using radioactive tracers and the impulse-response function. Eur J Nucl Med 11:265-274 (1985).

(25.) Lever SZ, Scheffel U. Regional distribution of [sup.203]Pb[Cl.sub.2] in the mouse after intravenous injection. Neurotoxicology 19:197-207 (1998).

Jonathan M. Links,(1) Brian S. Schwartz,(1,2) David Simon,(2) Karen Bandeen-Roche,(3) and Walter F. Stewart(2)

Departments of (1) Environmental Health Sciences, (2) Epidemiology, and (3) Biostatistics, The Johns Hopkins School of Hygiene and Public Health, Baltimore, Maryland, USA

Address correspondence to J.M. Links, Department of Environmental Health Sciences, Johns Hopkins School of Hygiene and Public Health, 615 North Wolfe Street, Baltimore, MD 21205-2179 USA. Telephone: (410) 955-9622. Fax: (410) 955-6222. E-mail: jlinks@jhsph.edu

We thank J. Prince and A. Todd for helpful discussions and the reviewers for important clarifying comments.

This research was supported in part by NIH grants R01 AG10785 from the National Institute on Aging and P30 ES03819 from the National Institute of Environmental Health Sciences.

Received 22 September 2000; accepted 25 October 2000.

## Reader Opinion