A Randomized, Double-Blind, Multicenter, Noninferiority Study Comparing Paliperidone Palmitate 6-Month Versus the 3-Month Long-Acting Injectable in Patients With Schizophrenia.

BACKGROUND: This double-blind (DB), randomized, parallel-group study was designed to evaluate efficacy and safety of paliperidone palmitate 6-month (PP6M) formulation relative to paliperidone palmitate 3-month (PP3M) formulation in patients with schizophrenia. METHODS: Following screening, patients entered an open-label (OL) maintenance phase and received 1 injection cycle of paliperidone palmitate 1-month (PP1M; 100 or 150 mg eq.) or PP3M (350 or 525 mg eq.). Clinically stable patients were randomized (2:1) to receive PP6M (700 or 1000 mg eq., gluteal injections) or PP3M (350 or 525 mg eq.) in a 12-month DB phase; 2 doses of PP6M (corresponding to doses of PP1M and PP3M) were chosen. RESULTS: Overall, 1036 patients were screened, 838 entered the OL phase, and 702 (mean age: 40.8 years) were randomized (PP6M: 478; PP3M: 224); 618 (88.0%) patients completed the DB phase (PP6M: 416 [87.0%]; PP3M: 202 [90.2%]). Relapse rates were PP6M, 7.5% (n = 36) and PP3M, 4.9% (n = 11). The Kaplan-Meier estimate of the difference (95% CI) between treatment groups (PP6M - PP3M) in the percentages of patients who remained relapse free was -2.9% (-6.8%, 1.1%), thus meeting noninferiority criteria (95% CI lower bound is larger than the pre-specified noninferiority margin of -10%). Secondary efficacy endpoints corroborated the primary analysis. Incidences of treatment-emergent adverse events were similar between PP6M (62.1%) and PP3M (58.5%). No new safety concerns emerged. CONCLUSIONS: The efficacy of a twice-yearly dosing regimen of PP6M was noninferior to that of PP3M in preventing relapse in patients with schizophrenia adequately treated with PP1M or PP3M. TRIAL REGISTRATION: Clinical Trials.gov identifier: NCT03345342.

The impact of CYP3A5*3 polymorphism on sirolimus pharmacokinetics: insights from predictions with a physiologically-based pharmacokinetic model.

AIMS: Sirolimus is an mTOR inhibitor metabolized by CYP3A4 and CYP3A5. Reported effects of CYP3A5 polymorphisms on sirolimus pharmacokinetics (PK) have shown unexplained discrepancies across studies. We quantitatively assessed the effect of CYP3A5*3 status on sirolimus PK by in vitro assessment and simulation using a physiologically-based PK (PBPK) model. In addition, we explored designs for an adequately powered pharmacogenetic association study. METHOD: In vitro metabolism studies were conducted to confirm individual CYP contribution to sirolimus metabolism. PK profiles were simulated in CYP3A5 expressers and non-expressers with a PBPK model. The pre-dose concentration predictions were used as the outcome parameter to estimate the required sample size for a pharmacogenetic association study. RESULTS: Sirolimus metabolism was inhibited by over 90% by ketoconazole, a CYP3A specific inhibitor. The PBPK model developed based on CL(int) of recombinant CYP3A4, CYP3A5 and CYP2C8 predicted a small CYP3A5*3 effect on simulated sirolimus PK profiles. A subsequent power analysis based on these findings indicated that at least 80 subjects in an enrichment design, 40 CYP3A5 expressers and 40 non-expressers, would be required to detect a significant difference in the predicted trough concentrations at 1 month of therapy (P < 0.05, 80% power). CONCLUSIONS: This study suggests that CYP3A5 contribution to sirolimus metabolism is much smaller than that of CYP3A4. Observed discrepancies across studies could be explained as the result of inadequate sample size. PBPK model simulations allowed mechanism-based evaluation of the effects of CYP3A5 genotype on sirolimus PK and provided preliminary data for the design of a future prospective study.

Population pharmacokinetic-pharmacodynamic modeling and dosing simulation of propofol maintenance anesthesia in severely obese adolescents.

BACKGROUND: Optimal dosing of propofol to maintain appropriate anesthetic depth is challenging in severely obese (SO) adolescents. We previously reported that total body weight (TBW) is predictive of propofol clearance. This study was aimed at characterizing pharmacokinetics (PK) and pharmacodynamics (PD) of propofol in SO adolescents, using bispectral index (BIS), and toward developing PK/PD model-based dosing guidelines. METHODS: A prospective PK/PD study was conducted in 26 SO children and adolescents aged 9-18 years (body mass index 31-69 kg·m(-2)), undergoing surgery with intravenous propofol anesthesia clinically titrated by providers blinded to BIS. BIS data and propofol infusion schemes were recorded. Venous blood samples collected during and after propofol infusion were assayed for propofol concentrations. A propofol PK/PD model was developed using NONMEM and model-based simulations were performed to determine propofol dosing regimens targeting BIS of 50 ± 10. RESULTS: A three-compartment PK model linked to a sigmoidal inhibitory Emax PD model by a first-order rate constant, adequately described the propofol concentration (n = 375) and BIS (n = 3334) data. TBW was the most predictive covariate for propofol clearance [CL (l·min(-1) ) = 1.65 × (TBW/70)(0.75)]. An effect-site propofol concentration of 3.19 µg·ml(-1) was estimated for half-maximal effect, with no identifiable predictive covariates. The proposed maintenance dosing regimen targeted to a BIS of 50 ± 10, based on our PK/PD model, was able to predict desired propofol concentrations and BIS in a representative obese teen when used in conjunction with accepted PK/PD models for children/obese adults (PK:Eleveld/PD: Cortinez), further supporting evidence for the dosing based on TBW. CONCLUSION: This is the first study to describe the PK/PD of propofol in SO adolescents. The proposed maintenance dosing regimen for propofol uses TBW in an allometric function as a dosing scalar, with an exponent of 0.75. Our results suggest no relevant effect of obesity on the propofol concentration-BIS relationship.

Model-Informed Clinical Development of Once-Every-6-Month Injection of Paliperidone Palmitate in Patients with Schizophrenia: A Pharmacometric Bridging Approach (Part I).

BACKGROUND AND OBJECTIVE: A model-informed drug development (MIDD) approach was implemented for paliperidone palmitate (PP) 6-month (PP6M) clinical development, using pharmacokinetics and pharmacokinetic/pharmacodynamic model-based simulations. METHODS: PP6M pharmacokinetics were simulated by extending the PP 3-month (PP3M) pharmacokinetic model to account for increased injection volume, and hence dose. Contribution of the MIDD approach to the design of the pivotal PP6M phase-3 study (PP6M/PP3M noninferiority study, NCT03345342) investigating schizophrenia relapse rates was twofold: (1) PP6M dose selection, and (2) hypothesis generation that lower trough concentrations (Ctrough) associated with PP6M, relative to PP3M, were not associated with lower efficacy, which was to be evaluated in the phase-3 study. Moreover, accompanied by an intense sampling scheme to adequately characterize paliperidone pharmacokinetics and to elucidate the potential relationship between concentration and safety/efficacy, the bridging strategy eliminated the need for additional phase-1/phase-2 clinical studies. RESULTS: Using a MIDD bridging strategy, PP6M doses were selected that, compared with PP3M, were expected to have a similar range of exposures and a noninferior relapse rate and safety profile. Clinical data from PP6M/PP3M noninferiority study confirmed that PP6M, compared with PP3M, had a similar range of exposures (T'jollyn et al. in Eur J Drug Metab Pharmacokinet 2024), as well as a noninferior relapse rate and safety profile (this manuscript). CONCLUSIONS: Consistency of the MIDD approach with observed clinical outcomes confirmed the hypothesis that lower Ctrough did not lead to increased relapse rates at the doses administered. Although higher paliperidone peak concentrations are achieved with corresponding doses of PP6M relative to PP3M in the phase-3 clinical study, types and incidences of treatment-related adverse events were comparable between PP6M and PP3M groups and no new safety concerns emerged for PP6M (Najarian et al. in Int J Neuropsychopharmacol 25(3):238-251, 2022).

Model-Informed Clinical Development of 6-Monthly Injection of Paliperidone Palmitate in Patients with Schizophrenia: Dosing Strategies Guided by Population Pharmacokinetic Modeling and Simulation (Part II).

BACKGROUND AND OBJECTIVE: Paliperidone palmitate 6-month (PP6M) intramuscular (IM) injection is the longest-acting treatment available for patients with schizophrenia. A population pharmacokinetic (popPK) modeling and simulation approach was deployed to inform dosing strategies for PP6M. METHODS: The extensive analysis database included 15,932 paliperidone samples from 700 patients receiving gluteal paliperidone palmitate 3-month (PP3M) or PP6M injections in the double-blind phase of a phase-3 noninferiority study (NCT03345342). Exposure parameters for paliperidone appeared to increase dose-proportionally within each dosing schedule (PP3M/PP6M). The range of paliperidone exposures after IM administration of PP6M overlaps with that of corresponding doses of oral paliperidone extended release, PP 1-month (PP1M), and PP3M. Model-based simulations were performed to investigate paliperidone exposures in different PP6M dosing scenarios and relevant subpopulations. RESULTS: A dosing window of ≤ 2 weeks earlier and ≤ 3 weeks later than the target 6-month interval for maintenance treatment with PP6M dosing maintains paliperidone exposures at levels that are not expected to substantially impact its safety and efficacy. For missed-dose scenarios, tailored re-initiation regimens are proposed that should be applied before resuming PP6M maintenance dosing. Regarding subpopulations, PP6M 700 mg eq. is the highest dose recommended in mild renal-impairment patients; the paliperidone pharmacokinetics after PP6M administration is not affected by sex, body mass index, or age in a clinically meaningful way. CONCLUSION: Paliperidone concentration-time profiles after PP6M and PP3M dosing were adequately described by the popPK model. Model-based simulation results provide guidance for clinicians on initiating PP6M therapy, transitioning between paliperidone formulations, the dosing windows to use for maintenance dosing, and managing missed PP6M doses.

Pharmacokinetic Characteristics of Long-Acting Injectable Antipsychotics for Schizophrenia: An Overview.

The availability of long-acting injectable (LAI) antipsychotics for the treatment of schizophrenia provides clinicians with options that deliver continuous drug exposure and may improve adherence compared with daily oral antipsychotics. However, all LAI antipsychotics have unique formulations and pharmacokinetic characteristics that have implications for medication selection, administration interval, and injection site. This review outlines key differences in drug formulations and pharmacokinetics among LAI antipsychotics. A systematic search of the PubMed database was conducted to identify physical and formulation properties and pharmacokinetic data of commercially available LAI antipsychotics, including flupentixol decanoate, fluphenazine decanoate, haloperidol decanoate, zuclopenthixol decanoate, aripiprazole monohydrate, aripiprazole lauroxil, olanzapine pamoate, paliperidone palmitate, risperidone microspheres, and risperidone polymeric microspheres. Additional information was obtained from package inserts and product monographs. Relevant data on drug properties, administration details, pharmacokinetic parameters, and oral dose equivalencies of LAI antipsychotics are summarized. Based on our analysis, formulation characteristics (e.g., vehicle medium) and administration characteristics (e.g., injection site) can affect rate of absorption and adverse effects and may factor into whether oral supplementation or an additional injection is needed. Dose adjustments may be necessary based on potential drug-drug interactions, and approximate dose equivalence with oral formulations can help inform titration when switching from oral to LAI formulations. Clinicians administering LAI antipsychotics should consider these formulation and pharmacokinetic factors to maximize clinical impact and to adjust to an individual patient's needs and treatment goals.

Semi-mechanistic modelling platform to assess cardiac contractility and haemodynamics in preclinical cardiovascular safety profiling of new molecular entities.

BACKGROUND AND PURPOSE: Cardiovascular safety is one of the most frequent causes of safety-related attrition both preclinically and clinically. Preclinical cardiovascular safety is routinely assessed using dog telemetry monitoring key cardiovascular functions. The present research was to develop a semi-mechanistic modelling platform to simultaneously assess changes in contractility (dPdtmax ), heart rate (HR) and mean arterial pressure (MAP) in preclinical studies. EXPERIMENTAL APPROACH: Data from dPdtmax , HR, preload (left ventricular end-diastolic pressure [LVEDP]) and MAP were available from dog telemetry studies after dosing with atenolol (n = 27), salbutamol (n = 5), L-NG -nitroarginine methyl ester (L-NAME; n = 4), milrinone (n = 4), verapamil (n = 12), dofetilide (n = 8), flecainide (n = 4) and AZ001 (n = 14). Literature model for rat CV function was used for the structural population pharmacodynamic model development. LVEDP was evaluated as covariate to account for the effect of preload on dPdtmax . KEY RESULTS: The model was able to describe drug-induced changes in dPdtmax , HR and MAP for all drugs included in the developed framework adequately, by incorporating appropriate drug effects on dPdtmax , HR and/or total peripheral resistance. Consistent with the Starling's law, incorporation of LVEDP as a covariate on dPdtmax to correct for the preload effect was found to be statistically significant. CONCLUSIONS AND IMPLICATIONS: The contractility and haemodynamics semi-mechanistic modelling platform accounts for diurnal variation, drug-induced changes and inter-animal variation. It can be used to hypothesize and evaluate pharmacological effects and provide a holistic cardiovascular safety profile for new drugs.

Integrating cell-cycle progression, drug penetration and energy metabolism to identify improved cancer therapeutic strategies.

The effectiveness of chemotherapeutic drugs in tumors is reduced by multiple effects including drug diffusion and variable susceptibility of local cell populations. We hypothesized that quantifying the interactions between drugs and tumor microenvironments could be used to identify more effective anti-cancer strategies. To test this hypothesis we created a mathematical model that integrated intracellular metabolism, nutrient and drug diffusion, cell-cycle progression, cellular drug effects, and drug pharmacokinetics. To our knowledge, this is the first model that combines these elements and has coupled them to experimentally derived parameters. Drug cytotoxicity was assumed to be cell-cycle phase specific, and progression through the cell cycle was assumed to be dependent on ATP generation. The model consisted of a coupled set of nonlinear partial differential, ordinary differential and algebraic equations with an outer free boundary, which was solved using orthogonal collocation on a moving grid of finite elements. Model simulations showed the existence of an optimum drug diffusion coefficient: a low diffusivity prevents effective penetration before the drug is cleared from the blood and a high diffusivity limits drug retention. This result suggests that increasing the molecular weight of the anti-cancer drug paclitaxel from 854 to approximately 20,000 by nanoparticle conjugation would improve its efficacy. The simulations also showed that fast growing tumors are less responsive to therapy than are slower tumors with more quiescent cells, demonstrating the competing effects of regrowth and cytotoxicity. The therapeutic implications of the simulation results are that (1) monolayer cultures are inadequate for accurately determining therapeutic effects in vitro, (2) decreasing the diffusivity of paclitaxel could increase its efficacy, and (3) measuring the proliferation fraction in tumors could enhance the prediction of therapeutic efficacy.

CYP2D6 pharmacogenetic and oxycodone pharmacokinetic association study in pediatric surgical patients.

AIM: Oxycodone is partly metabolized to the active metabolite oxymorphone by hepatic CYP2D6 in the liver. Significant genetic variability in CYP2D6 activity affects oxymorphone formation. This study aimed to associate CYP2D6 genotype and oxycodone's metabolism. METHODS: 30 children were administered oral oxycodone postoperatively. Plasma levels of oxycodone and oxymorphone, and CYP2D6 genotype were analyzed. CYP2D6 genotype and oxycodone metabolism phenotype were determined based on CYP2D6 total activity score (TAS) and metabolism phenotype: poor metabolizer (PM), intermediate metabolizer (IM), extensive metabolizer (EM) or ultrarapid metabolizer (UM). RESULTS: Compared with PM/IM subjects, significantly greater oxymorphone exposure was seen in EM subjects (p = 0.02 for Cmax, p = 0.016 for AUC0-6 and p = 0.026 for AUC0-24). Similarly, higher TAS value was found to be associated with greater oxymorphone exposure. Higher conversion of oxycodone to oxymorphone was observed in EM subjects compared with PM/IM subjects (p = 0.0007 for Cmax, p = 0.001 for AUC0-6 and p = 0.004 for AUC0-24). CONCLUSION: CYP2D6 phenotypes explain metabolism of oxycodone in children, and oxymorphone exposure is higher in CYP2D6 EM phenotype. Further studies are needed to predict the occurrence of adverse event and tailor oxycodone dose for a specific CYP2D6 phenotype.

Fatty acid amide hydrolase-morphine interaction influences ventilatory response to hypercapnia and postoperative opioid outcomes in children.

AIM: Fatty acid amide hydrolase (FAAH) degrades anandamide, an endogenous cannabinoid. We hypothesized that FAAH variants will predict risk of morphine-related adverse outcomes due to opioid-endocannabinoid interactions. PATIENTS & METHODS: In 101 postsurgical adolescents receiving morphine analgesia, we prospectively studied ventilatory response to 5% CO2 (HCVR), respiratory depression (RD) and vomiting. Blood was collected for genotyping and morphine pharmacokinetics. RESULTS: We found significant FAAH-morphine interaction for missense (rs324420) and several regulatory variants, with HCVR (p < 0.0001) and vomiting (p = 0.0339). HCVR was more depressed in patients who developed RD compared with those who did not (p = 0.0034), thus FAAH-HCVR association predicts risk of impending RD from morphine use. CONCLUSION: FAAH genotypes predict risk for morphine-related adverse outcomes.

Propofol Pharmacokinetics and Estimation of Fetal Propofol Exposure during Mid-Gestational Fetal Surgery: A Maternal-Fetal Sheep Model.

BACKGROUND: Measuring fetal drug concentrations is extremely difficult in humans. We conducted a study in pregnant sheep to simultaneously describe maternal and fetal concentrations of propofol, a common intravenous anesthetic agent used in humans. Compared to inhalational anesthesia, propofol supplemented anesthesia lowered the dose of desflurane required to provide adequate uterine relaxation during open fetal surgery. This resulted in better intraoperative fetal cardiac outcome. This study describes maternal and fetal propofol pharmacokinetics (PK) using a chronically instrumented maternal-fetal sheep model. METHODS: Fetal and maternal blood samples were simultaneously collected from eight mid-gestational pregnant ewes during general anesthesia with propofol, remifentanil and desflurane. Nonlinear mixed-effects modeling was performed by using NONMEM software. Total body weight, gestational age and hemodynamic parameters were tested in the covariate analysis. The final model was validated by bootstrapping and visual predictive check. RESULTS: A total of 160 propofol samples were collected. A 2-compartment maternal PK model with a third fetal compartment appropriately described the data. Mean population parameter estimates for maternal propofol clearance and central volume of distribution were 4.17 L/min and 37.7 L, respectively, in a typical ewe with a median heart rate of 135 beats/min. Increase in maternal heart rate significantly correlated with increase in propofol clearance. The estimated population maternal-fetal inter-compartment clearance was 0.0138 L/min and the volume of distribution of propofol in the fetus was 0.144 L. Fetal propofol clearance was found to be almost negligible compared to maternal clearance and could not be robustly estimated. CONCLUSIONS: For the first time, a maternal-fetal PK model of propofol in pregnant ewes was successfully developed. This study narrows the gap in our knowledge in maternal-fetal PK model in human. Our study confirms that maternal heart rate has an important influence on the pharmacokinetics of propofol during pregnancy. Much lower propofol concentration in the fetus compared to maternal concentrations explain limited placental transfer in in-vivo paired model, and less direct fetal cardiac depression we observed earlier with propofol supplemented inhalational anesthesia compared to higher dose inhalational anesthesia in humans and sheep.

Rapid uptake of glucose and lactate, and not hypoxia, induces apoptosis in three-dimensional tumor tissue culture.

The spatial arrangement of cellular metabolism in tumor tissue critically affects the treatment of cancer. However, little is known about how diffusion and cellular uptake relate to intracellular metabolism and cell death in three dimensions. To quantify these mechanisms, fluorescent microscopy and multicellular tumor cylindroids were used to measure pH and oxygen profiles, and quantify the distribution of viable, apoptotic and necrotic cells. Spheroid dissociation, enzymatic analysis, and mass spectrometry were used to measure concentration profiles of glucose, lactate and glutamine. A mathematical model was used to integrate these measurements and calculate metabolic rate parameters. It was found that large cylindroids, >500 µm in diameter, contained apoptotic and necrotic cells, whereas small cylindroids contained apoptotic but not necrotic cells. The center of cylindroids was found to be acidic but not hypoxic. From the edge to the center, concentrations of glucose, lactate and glutamine decreased rapidly. Throughout the cell masses lactate was consumed and not produced. These measurements indicate that apoptosis was the primary mechanism of cell death; acidity was not caused by lactic acid; and cell death was caused by depletion of carbon sources and not hypoxia. The mathematical model showed that the transporter enzymes for glucose and lactate were not saturated; oxygen uptake was limited by intracellular metabolism; and oxygen uptake was not limited by membrane-transport or diffusion. Unsaturated transmembrane uptake may be the cause of both proliferative and apoptotic regimes in cancer. These results suggest that transporter enzymes are excellent targets for treating well oxygenated tumors.

ABCC3 and OCT1 genotypes influence pharmacokinetics of morphine in children.

AIM: Large interindividual variability in morphine pharmacokinetics could contribute to variability in morphine analgesia and adverse events. METHODS: Influence of weight, genetic polymorphisms, race and sex on morphine clearance and metabolite formation from 220 children undergoing outpatient adenotonsillectomy was studied. A nonlinear mixed effects model was developed in NONMEM to describe morphine and morphine glucuronide pharmacokinetics. RESULTS: Children with ABCC3 -211C>T polymorphism C/C genotype had significantly higher levels of morphine-6-glucuronide and morphine-3-glucuronide formation (∼40%) than C/T+T/T genotypes (p < 0.05). In this extended cohort similar to our earlier report, OCT1 homozygous genotypes (n = 13, OCT1*2-*5/*2-*5) had lower morphine clearance (14%; p = 0.06), and in addition complementing lower metabolite formation (∼39%) was observed. ABCB1 3435C>T TT genotype children had lower levels of morphine-3-glucuronide formation though no effect was observed on morphine and morphine-6-glucuronide pharmacokinetics. CONCLUSION: Our data suggest that besides bodyweight, OCT1 and ABCC3 genotypes play a significant role in the pharmacokinetics of intravenous morphine and its metabolites in children.

OCT1 genetic variants influence the pharmacokinetics of morphine in children.

AIM: Large interindividual variability in morphine disposition could contribute to unpredictable variability in morphine analgesia and adverse events. Caucasian children have more adverse effects and slower morphine clearance than African-American children. To study variations in intravenous morphine pharmacokinetics in children, we examined the influence of genetic polymorphisms in OCT1. METHODS: In 146 children undergoing adenotonsillectomy, 146 concentration-time profiles (2-4 measurements per patient) were available. Population pharmacokinetic analysis characterized the profiles in NONMEM(®) and tested OCT1 variants as covariates. RESULTS: Allometrically scaled post hoc Bayesian morphine clearance in homozygotes of loss-of-function OCT1 variants (n = 9, OCT1*2-*5/*2-*5) was significantly lower (20%) than in wild-type (n = 85, OCT1*1/*1) and heterozygotes (n = 52, OCT1*1/*2-*5; p < 0.05). CONCLUSION: Besides bodyweight, OCT1 genotypes play a significant role in intravenous morphine pharmacokinetics. Relatively high allelic frequencies of defective OCT1 variants among Caucasians may explain their lower morphine clearance and possibly higher frequencies of adverse events compared with African-American children. Original submitted 21 December 2012; Revision submitted 7 May 2013.

Micrometer-scale oxygen delivery rearranges cells and prevents necrosis in tumor tissue in vitro.

Oxygen availability plays a critical role in cancer progression and is correlated with poor prognosis. Despite this connection, the independent effects of oxygen gradients on tumor tissues have not been measured. To address this, we developed an oxygen delivery device that uses microelectrodes to generate oxygen directly underneath three-dimensional tumor cylindroids composed of colon carcinoma cells. The extent of cell death was measured using fluorescence staining. Supplying oxygen for 60 h eliminated the necrotic region typically found in the center of cylindroids despite the continued presence of other nutrient gradients. A mathematical model of cylindroid growth showed that the rate of cell death was more sensitive to oxygen than the growth rate. After oxygenation, a ring of dead cells was observed at the outside edge of cylindroids, and dead cells were observed moving outward from cylindroid centers. This movement suggests that dead cells were pushed by viable cells migrating in response to oxygen gradients, a mechanism that may connect transient oxygen gradients to metastasis formation. These measurements show that oxygen gradients are a primary factor governing cell viability and rearrange cells in tumors.

Incorporating energy metabolism into a growth model of multicellular tumor spheroids.

Diffusion limitations in tumors create regions that are deficient in essential nutrients and contain a large number of quiescent and dying cells. Chemotherapeutic compounds are not effective against quiescent cells and therefore have reduced efficacy against tumors with extensive quiescence. We have formulated a mathematical model that predicts the extent and location of quiescence in multicellular spheroids. Multicellular spheroids are in vitro models of in vivo tumor growth that have proven to be useful experimental systems for studying radiation therapy, drug penetration, and novel chemotherapeutic strategies. Our model incorporates a realistic description of primary energy metabolism within reaction-diffusion equations to predict local glucose, oxygen, and lactate concentrations and an overall spheroid growth rate. The model development is based on the assumption that local cellular growth and death rates are determined by local ATP production generated by intracellular energy metabolism. Dynamic simulation and parametric sensitivity studies are used to evaluate model behavior, including the spatial distribution of proliferating, quiescent, and dead cells for different cellular characteristics. Using this model we have determined the critical cell survival parameters that have the greatest impact on overall spheroid physiology, and we have found that oxygen transport has a greater effect than glucose transport on the distribution of quiescent cells. By predicting the extent of quiescence based on individual cellular characteristic alone this model has the potential to predict therapeutic efficiency and can be used to design effective chemotherapeutic strategies.