httk (version 1.8)

solve_pbtk: Solve_PBTK

Description

This function solves for the amounts or concentrations in uM of a chemical in different tissues as functions of time based on the dose and dosing frequency.

Usage

solve_pbtk(chem.name = NULL, chem.cas = NULL, times=NULL, parameters=NULL,
           days=10, tsteps = 4, daily.dose=1,dose = NULL,doses.per.day=NULL,
           initial.values=NULL,plots=F,suppress.messages=F,species="Human", 
           iv.dose=F,output.units='uM',method="lsoda",rtol=1e-8,atol=1e-12,
           default.to.human=F,recalc.blood2plasma=F,recalc.clearance=F,
           dosing.matrix=NULL,adjusted.Funbound.plasma=T,regression=T, 
           restrictive.clearance = T,...)

Arguments

chem.name

Either the chemical name, CAS number, or the parameters must be specified.

chem.cas

Either the chemical name, CAS number, or the parameters must be specified.

times

Optional time sequence for specified number of days. Dosing sequence begins at the beginning of times.

parameters

Chemical parameters from parameterize_pbtk function, overrides chem.name and chem.cas.

days

Length of the simulation.

tsteps

The number time steps per hour.

daily.dose

Total daily dose, mg/kg BW.

dose

Amount of a single dose, mg/kg BW. Overwrites daily.dose.

doses.per.day

Number of doses per day.

initial.values

Vector containing the initial concentrations or amounts of the chemical in specified tissues with units corresponding to output.units. Defaults are zero.

plots

Plots all outputs if true.

suppress.messages

Whether or not the output message is suppressed.

species

Species desired (either "Rat", "Rabbit", "Dog", "Mouse", or default "Human").

iv.dose

Simulates a single i.v. dose if true.

output.units

Desired units (either "mg/L", "mg", "umol", or default "uM").

method

Method used by integrator (deSolve).

rtol

Argument passed to integrator (deSolve).

atol

Argument passed to integrator (deSolve).

default.to.human

Substitutes missing animal values with human values if true (hepatic intrinsic clearance or fraction of unbound plasma).

recalc.blood2plasma

Recalculates the ratio of the amount of chemical in the blood to plasma using the input parameters, calculated with hematocrit, Funbound.plasma, and Krbc2pu.

recalc.clearance

Recalculates the the hepatic clearance (Clmetabolism) with new million.cells.per.gliver parameter.

dosing.matrix

Vector of dosing times or a matrix consisting of two columns or rows named "dose" and "time" containing the time and amount, in mg/kg BW, of each dose.

adjusted.Funbound.plasma

Uses adjusted Funbound.plasma when set to TRUE along with partition coefficients calculated with this value.

regression

Whether or not to use the regressions in calculating partition coefficients.

restrictive.clearance

Protein binding not taken into account (set to 1) in liver clearance if FALSE.

...

Additional arguments passed to the integrator.

Value

A matrix of class deSolve with a column for time(in days), each compartment, the area under the curve, and plasma concentration and a row for each time point.

Details

Note that the model parameters have units of hours while the model output is in days.

Default NULL value for doses.per.day solves for a single dose.

The compartments used in this model are the gutlumen, gut, liver, kidneys, veins, arteries, lungs, and the rest of the body.

The extra compartments include the amounts or concentrations metabolized by the liver and excreted by the kidneys through the tubules.

AUC is the area under the curve of the plasma concentration.

Examples

# NOT RUN {
solve_pbtk(chem.name='Bisphenol-A',dose=.5,days=1,doses.per.day=2,tsteps=2)
out <- solve_pbtk(chem.name='bisphenola',dose=0,output.units='mg', 
                  plots=TRUE,initial.values=c(Agut=200))
# }
# NOT RUN {
parameters <- parameterize_pbtk(chem.name = "triclosan", species = "rat")
parameters["Funbound.plasma"] <- 0.1
out <- solve_pbtk(parameters=parameters)

library("ggplot2")
out <- solve_pbtk(chem.name = "Bisphenol A", days = 50, doses.per.day = 3)
plot.data <- as.data.frame(out)
css <- calc_analytic_css(chem.name = "Bisphenol A")
c.vs.t <- ggplot(plot.data,aes(time, Cplasma)) + geom_line() +
geom_hline(yintercept = css) + ylab("Plasma Concentration (uM)") +
xlab("Day") + theme(axis.text = element_text(size = 16), axis.title =
element_text(size = 16), plot.title = element_text(size = 17)) +
ggtitle("Bisphenol A")
print(c.vs.t)
# }