In vitro and in vivo butadiene (BD) metabolism data from laboratory animals were integrated into a rodent physiologically based pharmacokinetic (PBPK) model with flow- and diffusion-limited compartments. The resulting model describes experimental data from multiple sources under scenarios such as closed chamber inhalation and nose-only flow-through inhalation exposures. Incorporation of diurnal glutathione (GSH) variation allows accurate simulation of GSH changes observed in air control nose-only exposures and BD exposures. An isolated tissue model based on rate parameters determined in vitro predicts the decrease in epoxide concentrations in intact animals during the time lag between exsanguination and tissue removal for tissues capable of epoxide biotransformation, providing a better indication of in vivo dosimetry. Further refinements of the model are required relative to model predictions of an important BD metabolite, diepoxybutane.