Whole-Body Disposition and Physiologically Based Pharmacokinetic Modeling of Adeno-Associated Viruses and the Transgene Product.

To facilitate model-informed drug development (MIDD) of adeno-associated virus (AAV) therapy, here we have developed a physiologically based pharmacokinetic (PBPK) model for AAVs following preclinical investigation in mice. After 2E11 Vg/mouse dose of AAV8 and AAV9 encoding a monoclonal antibody (mAb) gene, whole-body disposition of both the vector and the transgene mAb was evaluated over 3 weeks. At steady-state, the following tissue-to-blood (T/B) concentration ratios were found for AAV8/9: ∼50 for liver; ∼10 for heart and muscle; ∼2 for brain, lung, kidney, adipose, and spleen; ≤1 for bone, skin, and pancreas. T/B values for mAb were compared with the antibody biodistribution coefficients, and five different clusters of organs were identified based on their transgene expression profile. All the biodistribution data were used to develop a novel AAV PBPK model that incorporates: (i) whole-body distribution of the vector; (ii) binding, internalization, and intracellular processing of the vector; (iii) transgene expression and secretion; and (iv) whole-body disposition of the secreted transgene product. The model was able to capture systemic and tissue PK of the vector and the transgene-produced mAb reasonably well. Pathway analysis of the PBPK model suggested that liver, muscle, and heart are the main contributors for the secreted transgene mAb. Unprecedented PK data and the novel PBPK model developed here provide the foundation for quantitative systems pharmacology (QSP) investigations of AAV-mediated gene therapies. The PBPK model can also serve as a quantitative tool for preclinical study design and preclinical-to-clinical translation of AAV-based gene therapies.

Pharmacokinetics of AAV9 Mediated Trastuzumab Expression in Rat Brain Following Systemic and Local Administration.

INTRODUCTION: Recombinant adeno-associated viruses(rAAVs) are an attractive tool to ensure long-term expression monoclonal antibody(mAb) in the central nervous system(CNS). It is still unclear whether systemic injection or local CNS administration of AAV9 is more beneficial for the exposure of the expressed mAb in the brain. Hence, we compared the biodistribution and transgene expression following AAV9-Trastuzumab administration through different routes. METHODS AND RESULT: In-house generated AAV9-Trastuzumab vectors were administered at 5E+11 Vgs/rat through intravenous(IV), intracerebroventricular(ICV), intra-cisterna magna(ICM) and intrastriatal(IST) routes. Vector and trastuzumab blood/plasma concentrations were assessed at different time points up to the terminal time point of 21 days. Different brain regions in addition to the spinal cord, cerebrospinal fluid(CSF) and interstitial fluid(ISF), were also analyzed at the terminal time point. Our results show that vector biodistribution and Trastuzumab expression in the brain could the ranked as follows: IST>ICM>ICV>IV. Rapid clearance of vector was observed after administration via the ICM and ICV routes. The ICV route produced similar expression levels across different brain regions, while the ICM route had better expression in the hindbrain and spinal cord region. The IST route had higher expression in the forebrain region compared to the hindbrain region. A sharp decline in trastuzumab plasma concentration was observed across all routes of administration due to anti-trastuzumab antibody response. CONCLUSION: In this study we have characterized vector biodistribution and transgene mAb expression after AAV9 vector administration through different routes in rats. IST and ICM represent the best administration routes to deliver antibody genes to the brain.

Assessment of AAV9 distribution and transduction in rats after administration through Intrastriatal, Intracisterna magna and Lumbar Intrathecal routes.

Challenges in obtaining efficient transduction of brain and spinal cord following systemic AAV delivery have led to alternative administration routes being used in clinical trials that directly infuse the virus into the CNS. However, data comparing different direct AAV injections into the brain remain limited making it difficult to choose optimal routes. Here we tested both AAV9-egfp and AAV9-fLuc delivery via intrastriatal (IST), intracisterna magna (ICM) and lumbar intrathecal (LIT) routes in adult rats and assessed vector distribution and transduction in brain, spinal cord and peripheral tissues. We find that IST infusion leads to robust transgene expression in the striatum, thalamus and cortex with lower peripheral tissue transduction and anti-AAV9 capsid titers compared to ICM or LIT. ICM delivery provided strong GFP and luciferase expression across more brain regions than the other routes and similar expression in the spinal cord to LIT injections, which itself largely failed to transduce the rat brain. Our data highlight the strengths and weaknesses of each direct CNS delivery route which will help with future clinical targeting.

Hallmarks of neurodegenerative disease: A systems pharmacology perspective.

Age-related central neurodegenerative diseases, such as Alzheimer's and Parkinson's disease, are a rising public health concern and have been plagued by repeated drug development failures. The complex nature and poor mechanistic understanding of the etiology of neurodegenerative diseases has hindered the discovery and development of effective disease-modifying therapeutics. Quantitative systems pharmacology models of neurodegeneration diseases may be useful tools to enhance the understanding of pharmacological intervention strategies and to reduce drug attrition rates. Due to the similarities in pathophysiological mechanisms across neurodegenerative diseases, especially at the cellular and molecular levels, we envision the possibility of structural components that are conserved across models of neurodegenerative diseases. Conserved structural submodels can be viewed as building blocks that are pieced together alongside unique disease components to construct quantitative systems pharmacology (QSP) models of neurodegenerative diseases. Model parameterization would likely be different between the different types of neurodegenerative diseases as well as individual patients. Formulating our mechanistic understanding of neurodegenerative pathophysiology as a mathematical model could aid in the identification and prioritization of drug targets and combinatorial treatment strategies, evaluate the role of patient characteristics on disease progression and therapeutic response, and serve as a central repository of knowledge. Here, we provide a background on neurodegenerative diseases, highlight hallmarks of neurodegeneration, and summarize previous QSP models of neurodegenerative diseases.

Current progress and limitations of AAV mediated delivery of protein therapeutic genes and the importance of developing quantitative pharmacokinetic/pharmacodynamic (PK/PD) models.

While protein therapeutics are one of the most successful class of drug molecules, they are expensive and not suited for treating chronic disorders that require long-term dosing. Adeno-associated virus (AAV) mediated in vivo gene therapy represents a viable alternative, which can deliver the genes of protein therapeutics to produce long-term expression of proteins in target tissues. Ongoing clinical trials and recent regulatory approvals demonstrate great interest in these therapeutics, however, there is a lack of understanding regarding their cellular disposition, whole-body disposition, dose-exposure relationship, exposure-response relationship, and how product quality and immunogenicity affects these important properties. In addition, there is a lack of quantitative studies to support the development of pharmacokinetic-pharmacodynamic models, which can support the discovery, development, and clinical translation of this delivery system. In this review, we have provided a state-of-the-art overview of current progress and limitations related to AAV mediated delivery of protein therapeutic genes, along with our perspective on the steps that need to be taken to improve clinical translation of this therapeutic modality.

From the prion-like propagation hypothesis to therapeutic strategies of anti-tau immunotherapy.

The term "propagon" is used to define proteins that may transmit misfolding in vitro, in tissues or in organisms. Among propagons, misfolded tau is thought to be involved in the pathogenic mechanisms of various "tauopathies" that include Alzheimer's disease, progressive supranuclear palsy, and argyrophilic grain disease. Here, we review the available data in the literature and point out how the prion-like tau propagation has been extended from Alzheimer's disease to tauopathies. First, in Alzheimer's disease, the progression of tau aggregation follows stereotypical anatomical stages which may be considered as spreading. The mechanisms of the propagation are now subject to intensive and controversial research. It has been shown that tau may be secreted in the interstitial fluid in an active manner as reflected by high and constant concentration of extracellular tau during Alzheimer's pathology. Animal and cell models have been devised to mimic tau seeding and propagation, and despite their limitations, they have further supported to the prion-like propagation hypothesis. Finally, such new ways of thinking have led to different therapeutic strategies in anti-tau immunotherapy among tauopathies and have stimulated new clinical trials. However, it appears that the prion-like propagation hypothesis mainly relies on data obtained in Alzheimer's disease. From this review, it appears that further studies are needed (1) to characterize extracellular tau species, (2) to find the right pathological tau species to target, (3) to follow in vivo tau pathology by brain imaging and biomarkers and (4) to interpret current clinical trial results aimed at reducing the progression of these pathologies. Such inputs will be essential to have a comprehensive view of these promising therapeutic strategies in tauopathies.

A translational platform PBPK model for antibody disposition in the brain.

In this manuscript, we have presented the development of a novel platform physiologically-based pharmacokinetic (PBPK) model to characterize brain disposition of mAbs in the mouse, rat, monkey and human. The model accounts for known anatomy and physiology of the brain, including the presence of distinct blood-brain barrier and blood-cerebrospinal fluid (CSF) barrier. CSF and interstitial fluid turnover, and FcRn mediated transport of mAbs are accounted for. The model was first used to characterize published and in-house pharmacokinetic (PK) data on the disposition of mAbs in rat brain, including the data on PK of mAb in different regions of brain determined using microdialysis. Majority of model parameters were fixed based on literature reported values, and only 3 parameters were estimated using rat data. The rat PBPK model was translated to mouse, monkey, and human, simply by changing the values of physiological parameters corresponding to each species. The translated PBPK models were validated by a priori predicting brain PK of mAbs in all three species, and comparing predicted exposures with observed data. The platform PBPK model was able to a priori predict all the validation PK profiles reasonably well (within threefold), without estimating any parameters. As such, the platform PBPK model presented here provides an unprecedented quantitative tool for prediction of mAb PK at the site-of-action in the brain, and preclinical-to-clinical translation of mAbs being developed against central nervous system (CNS) disorders. The proposed model can be further expanded to account for target engagement, disease pathophysiology, and novel mechanisms, to support discovery and development of novel CNS targeting mAbs.

Computer-assembled cross-species/cross-modalities two-pore physiologically based pharmacokinetic model for biologics in mice and rats.

Two-pore physiologically-based pharmacokinetic (PBPK) models can be expected to describe the tissue distribution and elimination kinetics of soluble proteins, endogenous or dosed, as function of their size. In this work, we amalgamated our previous two-pore PBPK model for an inert domain antibody (dAb) in mice with the cross-species platform PBPK model for monoclonal antibodies described in literature into a unified two-pore platform that describes protein modalities of different sizes and includes neonatal Fc receptor (FcRn) mediated recycling. This unified PBPK model was parametrized for organ-specific lymph flow rates and the endosomal recycling rate constant using an extended tissue distribution time-course dataset that included an inert dAb, albumin and IgG in rats and mice. The model was evaluated by comparing the ab initio predictions for the tissue distribution and elimination properties of albumin-binding dAbs (AlbudAbsTM) in mice and rats with the experimental observations. Due to the large number of molecular species and reactions involved in large-scale PBPK models, we have also developed and deployed a MatlabTM script for automating the assembly of SimBiologyTM-based two-pore biologics PBPK models which drastically cuts the time and effort required for model building.

Efficacy of Epratuzumab, an Anti-CD22 Monoclonal IgG Antibody, in Systemic Lupus Erythematosus Patients With Associated Sjögren's Syndrome: Post Hoc Analyses From the EMBODY Trials.

OBJECTIVE: EMBODY 1 (ClinicalTrials.gov identifier: NCT01262365) and EMBODY 2 (ClinicalTrials.gov identifier: NCT01261793) investigated the efficacy and safety of epratuzumab, a CD22-targeted humanized monoclonal IgG antibody, in patients with systemic lupus erythematosus (SLE). The studies showed no significant difference from placebo in primary or secondary clinical outcome measures but did demonstrate B cell-specific immunologic activity. The aim of this post hoc analysis was to determine whether epratuzumab had a different clinical efficacy profile in SLE patients with versus those without an associated diagnosis of Sjögren's syndrome (SS). METHODS: The efficacy and safety of epratuzumab were compared between 2 patient subpopulations randomized in EMBODY 1 and 2: SLE patients with and those without a diagnosis of associated SS. British Isles Lupus Assessment Group (BILAG) total score, BILAG-based Combined Lupus Assessment (BICLA) clinical response to treatment, biologic markers (including B cells, IgG, IgM, and IgA), and safety were assessed. RESULTS: A total of 1,584 patients were randomized in the EMBODY 1 and EMBODY 2 trials; 113 patients were anti-SSA positive and had a diagnosis of associated SS, and 1,375 patients (86.8%) had no diagnosis of associated SS (918 patients were randomized to receive epratuzumab and 457 to receive placebo). For patients with associated SS, but not those without associated SS, a higher proportion of patients receiving epratuzumab achieved a BICLA response and a reduction from baseline in BILAG total score. B cell reduction was faster in patients with associated SS. The sensitivity of B cells to epratuzumab as measured by the mean concentration producing 50% of the maximum B cell count depletion was lower for patients with associated SS (9.5 µg/ml) versus the total EMBODY population (87.1 µg/ml). No difference in the frequency of adverse events in those receiving placebo was reported. CONCLUSION: Patients with SLE and associated SS treated with epratuzumab showed improvement in SLE disease activity, which was associated with bioactivity, such as decreases in B cell number and IgM level.

Development of a physiologically based pharmacokinetic model for a domain antibody in mice using the two-pore theory.

Domain antibodies (dAbs) are the smallest antigen-binding fragments of immunoglobulins. To date, there is limited insight into the pharmacokinetics of dAbs, especially their distribution into tissues and elimination. The objective of this work was to develop a physiologically-based pharmacokinetic model to investigate the biodisposition of a non-specific dAb construct in mice. Following a single IV administration of 10 mg/kg dummy dAb protein to twenty four female mice, frequent blood samples were collected and whole body lateral sections were analyzed by quantitative whole-body autoradiography. The model is based on the two-pore hypothesis of extravasation where organ-specific isogravimetric flow rates (Jorg,ISO) and permeability-surface area products (PSorg) are expressed as linear functions of the lymph flow rate (Jorg) and the kidney compartment is modified to account for glomerular filtration of dAb. As a result, only Jorg, glomerular filtration coefficient and the combined volume of Bowman's capsule, proximal and distal renal tubules and loop of Henle were optimized by fitting simultaneously all blood and organ data to the model. Our model captures the pharmacokinetic profiles of dAb in blood and all organs and shows that extravasation into interstitial space is a predominantly diffusion-driven process. The parameter values were estimated with good precision (%RMSE ≈ 30) and low cross-correlation (R(2) < 0.2). We developed a flexible model with a limited parameter number that may be applied to other biotherapeutics after adapting for size-related effects on extravasation and renal elimination processes.

Safety, pharmacokinetic, and functional effects of the nogo-a monoclonal antibody in amyotrophic lateral sclerosis: a randomized, first-in-human clinical trial.

UNLABELLED: The neurite outgrowth inhibitor, Nogo-A, has been shown to be overexpressed in skeletal muscle in amyotrophic lateral sclerosis (ALS); it is both a potential biomarker and therapeutic target. We performed a double-blind, two-part, dose-escalation study, in subjects with ALS, assessing safety, pharmacokinetics (PK) and functional effects of ozanezumab, a humanized monoclonal antibody against Nogo-A. In Part 1, 40 subjects were randomized (3∶1) to receive single dose intravenous ozanezumab (0.01, 0.1, 1, 5, or 15 mg/kg) or placebo. In Part 2, 36 subjects were randomized (3∶1) to receive two repeat doses of intravenous ozanezumab (0.5, 2.5, or 15 mg/kg) or placebo, approximately 4 weeks apart. The primary endpoints were safety and tolerability (adverse events [AEs], vital signs, electrocardiogram (ECG), and clinical laboratory tests). Secondary endpoints included PK, immunogenicity, functional endpoints (clinical and electrophysiological), and biomarker parameters. Overall, ozanezumab treatment (0.01-15 mg/kg) was well tolerated. The overall incidence of AEs in the repeat dose 2.5 mg/kg and 15 mg/kg ozanezumab groups was higher than in the repeat dose placebo group and repeat dose 0.5 mg/kg ozanezumab group. The majority were considered not related to study drug by the investigators. Six serious AEs were reported in three subjects receiving ozanezumab; none were considered related to study drug. No study drug-related patterns were identified for ECG, laboratory, or vital signs parameters. One subject (repeat dose 15 mg/kg ozanezumab) showed a weak, positive anti-ozanezumab-antibody result. PK results were generally consistent with monoclonal antibody treatments. No apparent treatment effects were observed for functional endpoints or muscle biomarkers. Immunohistochemical staining showed dose-dependent co-localization of ozanezumab with Nogo-A in skeletal muscle. In conclusion, single and repeat dose ozanezumab treatment was well tolerated and demonstrated co-localization at the site of action. These findings support future studies with ozanezumab in ALS. TRIAL REGISTRATION: ClinicalTrials.gov NCT00875446 GSK-ClinicalStudyRegister.com GSK ID 111330.

A randomized, controlled trial of gabapentin enacarbil in subjects with neuropathic pain associated with diabetic peripheral neuropathy.

BACKGROUND: Gabapentin enacarbil (GEn), a transported prodrug of gabapentin, provides sustained, dose-proportional gabapentin exposure. The purpose of this study was to investigate the dose response of GEn to select the optimal dose(s) for clinical use in subsequent diabetic peripheral neuropathy (DPN) trials. METHODS: This was a multicenter, randomized, double-blind, double-dummy, parallel group, placebo-controlled trial with a study duration of approximately 20 weeks (Clinicaltrials.gov database, Identifier ! NCT00643760). Pregabalin (PGB) (Lyrica(®) ; Pfizer Inc.) was used as an active control to provide assay sensitivity of the trial. A total of 421 adult subjects with DPN were randomized in a ratio of 2:1:1:1:2 to receive oral GEn 3,600 mg/day, GEn 2,400 mg/day, GEn 1,200 mg/day, PGB 300 mg/day, or matching placebo, respectively. The primary efficacy endpoint was change from baseline to end of maintenance treatment with respect to the mean 24-hour average pain intensity score based on an 11-point Pain Intensity Numerical Rating Scale (PI-NRS). Safety and tolerability assessments included treatment-emergent adverse events (TEAEs), laboratory evaluations, vital signs, electrocardiograms (ECG), neurological examination, and pedal edema. RESULTS: The adjusted mean difference vs. placebo at the end of maintenance treatment with respect to the mean 24-hour average PI-NRS pain intensity score for GEn 1,200 mg (-0.35; [95% CI: -1.02, 0.31]; P = 0.295), GEn 2,400 mg (-0.02; [95% CI: -0.71, 0.66]; P = 0.946), and GEn 3,600 mg (-0.55; [95% CI: -1.10, 0.01]; P = 0.105) was not statistically significant. The active control, PGB (300 mg/day), did not differentiate from placebo. CONCLUSION: Overall, none of the GEn treatment groups differentiated from placebo. Analyses of the secondary endpoints showed comparable results across treatment groups. However, the majority of the endpoints, including all of the pain endpoints, showed the largest numerical treatment difference was between GEn 3,600 mg and placebo. The active control, PGB (300 mg/day), did not differentiate from placebo.

Physiologically based pharmacokinetic (PBPK) modeling of everolimus (RAD001) in rats involving non-linear tissue uptake.

Everolimus is a novel macrolide immunosuppressant developed for the prophylaxis of allogeneic renal or cardiac transplant rejection. Treatments with immunosuppressants are often associated with organ toxicity that is linked to high organ exposure. Therefore, gaining insight into the pharmacokinetics of everolimus in various organs is highly desirable especially those organs of therapeutic interest or those that pose safety concerns. The aim of this work was to characterize the disposition kinetics of everolimus in rats by physiologically based pharmacokinetic (PBPK) modeling. Blood and tissue samples were collected from male Wistar rats over 24 hr following intravenous (iv) bolus and iv infusion of 1 mg/kg and 10 mg/kg/2 hr of everolimus. Further blood samples were collected between 1 and 170 hr from a third group of rats, which received iv infusion of 1 mg/kg/2 hr of everolimus. Drug concentrations in blood and tissues were determined by a liquid chromatography reverse dilution method. Distribution of everolimus between blood fractions was determined in vitro at 37 degrees C. The results of the study demonstrated that everolimus exhibited moderate non-linear binding to red blood cells. Also, the tissue-to-blood concentration ratio decreased in all tissues as blood concentration increased. A PBPK model involving non-linear tissue binding was able to successfully describe the observed data in blood and all the organs investigated. The highest binding potential was observed in thymus, lungs, and spleen with the greatest tissue affinity observed in thymus, skin, and muscle as compared to other tissues. Everolimus exhibited a high clearance rate that was limited to the hepatic blood flow (47.2 ml/min/kg). The PBPK model was also able to predict the venous blood concentration reasonably well following oral administration. The oral bioavailability value, as estimated with the PBPK, was 12% and was similar to the value obtained by non-compartmental analysis. In conclusion, A PBPK model has been developed that successfully predicts the time course of everolimus in blood and a variety of organs. This model takes into account the non- linear binding of everolimus to red blood cells and tissues. This model may be used to predict everolimus concentration-time course in organs from other species including humans.

Physiologically based pharmacokinetic modeling of FTY720 (2-amino-2[2-(-4-octylphenyl)ethyl]propane-1,3-diol hydrochloride) in rats after oral and intravenous doses.

FTY720 (2-amino-2[2-(-4-octylphenyl)ethyl]propane-1,3-diol hydrochloride) is a new sphingosine-1-phosphate receptor agonist being developed for multiple sclerosis and prevention of solid organ transplant rejection. A physiologically based pharmacokinetic model was developed to predict the concentration of FTY720 in various organs of the body. Single oral and intravenous doses of FTY720 were administered to male Wistar rats, with blood and tissue sampling over 360 h analyzed by liquid chromatography/tandem mass spectrometry. A well stirred model (perfusion rate-limited) described FTY720 kinetics in heart, lungs, spleen, muscle, kidneys, bone, and liver, with a permeability rate-limited model being required for brain, thymus, and lymph nodes. Tissue-to-blood partition coefficients (RT) ranged from 4.69 (muscle) to 41.4 (lungs). In lymph nodes and spleen, major sites for FTY720-induced changes in sequestration of lymphocytes, RT values were 22.9 and 34.7, respectively. Permeability-surface area products for brain, thymus, and lymph nodes were 39.3, 122, and 176 ml/min. Intrinsic hepatic clearance was 23,145 l/h/kg for the free drug in blood (f(ub) 0.000333); systemic clearance was 0.748 l/h/kg and terminal half-life was 23.4 h. The fraction orally absorbed was 71%. The model characterized well FTY720 disposition for this extensive dosing and tissue collection study in the rat. On scaling the model to dogs and humans, good agreement was found between the actual and predicted blood concentration-time profiles. More importantly, brain concentrations in dogs were well predicted from those of the rat. In absolute terms, the predictions were slightly lower than observed values, just under a 1.5-fold deviation, but the model accurately predicted the terminal elimination of FTY720 from the brain.

On the prediction of the human response: a recycled mechanistic pharmacokinetic/pharmacodynamic approach.

Although it is routine to predict the blood or plasma pharmacokinetics of compounds for man based upon preclinical studies, the real value of such predictions only comes when linked to drug effects. In the first example, the immunomodulator, FTY720, the first sphingosine-1-phosphate receptor agonist, stimulates the sequestration of lymphocytes into lymph nodes thus removing cells from blood circulation. A prior physiology-based pharmacokinetic model fitted the concentration-time course of FTY720 in rats. This was connected to an indirect response model of the lymphocyte system to characterise the cell trafficking effects. The IC(50) of FTY720 was different in the rat compared with the monkey; man was assumed to be similar to the monkey. The systemic lymphocyte half-lives were also different between species. To make predictions of the pharmacodynamic behaviour for man, two elements are required, i) systemic exposure, in this case from an upscaled physiology based model, and ii) an estimate of lymphocyte turnover in man, gained from the literature from other drug treatments. Predictions compared well with clinical results. The second example is the monoclonal antibody Xolair, designed to bind immunoglobulin E for atopic diseases. A mechanism based two-site binding model described the kinetics of both Xolair and endogenous IgE. This model has been reused for other monoclonal antibodies designed to bind fluid-phase ligands. Sensitivity analysis shows that if differences across species in the kinetics of the endogenous system are not accounted for, then pharmacokinetic/pharmacodynamic models may give misleading predictions of the time course and extent of the response.

Pharmacokinetics and cell trafficking dynamics of 2-amino-2-[2-(4-octylphenyl)ethyl]propane-1,3-diol hydrochloride (FTY720) in cynomolgus monkeys after single oral and intravenous doses.

The pharmacokinetics and cell trafficking dynamics of 2-amino-2-[2-(4-octylphenyl)ethyl]propane-1,3-diol hydrochloride (FTY720), a novel immunosuppressive agent, were examined in cynomolgus monkeys (three males and three females). After single doses of 0.1 mg/kg p.o. or i.v. bolus and 1 mg/kg p.o. were administered to the animals, the concentrations of FTY720, and the numbers of lymphocytes, CD20+CD2-B cells, and CD2+CD20-T cells in blood were measured over 23 days. A linear three-compartment model characterized the time course of FTY720 concentrations with a terminal half-life of about 31 h, clearance of about 0.53 l/h/kg, and bioavailability of about 38%. The dynamic responses were not area under the curve (or dose) proportional for either males or females. An indirect response model with a distribution pool captured the cell trafficking data for all doses for each cell type, where initial blood counts (R(0)) were about 7650, 2100, and 5250 cells/microl; maximum fractional inhibition (I(max)) about 0.88, 0.85, and 0.91; influx (k(in)) about 6014, 1312, and 5662 cells/microl/h; efflux (k(out)) about 0.798, 0.555, and 1.08 h(-1); intercompartmental k(cp) about 0.134, 0.192, and 0.082 h(-1); and intercompartmental k(pc) rate constants about 3.9 x 10(-4), and 0.016 and 8.9 x 10(-6) h(-1) for lymphocytes, B cells, and T cells, respectively. The inhibition concentration IC(50) was about 0.48 microg/l for all cells, which was remarkably low. The apparent distribution volumes of peripheral pool (V(p)) were markedly larger than blood volume (V(b)) for all cells. The I(max) for cell trafficking was achieved at doses smaller than that producing graft protection, indicating stronger central than peripheral effects of this drug. The profound cell trafficking effects of FTY720 can be readily captured and interpreted with an extended indirect response model.