A long-acting GDF15 analog causes robust, sustained weight loss and reduction of food intake in an obese nonhuman primate model.

Growth Differentiation Factor-15 (GDF15) is a circulating polypeptide linked to cellular stress and metabolic adaptation. GDF15's half-life is ~3 h and activates the glial cell line-derived neurotrophic factor family receptor alpha-like (GFRAL) receptor expressed in the area postrema. To characterize sustained GFRAL agonism on food intake (FI) and body weight (BW), we tested a half-life extended analog of GDF15 (Compound H [CpdH]) suitable for reduced dosing frequency in obese cynomolgus monkeys. Animals were chronically treated once weekly (q.w.) with CpdH or long-acting GLP-1 analog dulaglutide. Mechanism-based longitudinal exposure-response modeling characterized effects of CpdH and dulaglutide on FI and BW. The novel model accounts for both acute, exposure-dependent effects reducing FI and compensatory changes in energy expenditure (EE) and FI occurring over time with weight loss. CpdH had linear, dose-proportional pharmacokinetics (terminal half-life ~8 days) and treatment caused exposure-dependent reductions in FI and BW. The 1.6 mg/kg CpdH reduced mean FI by 57.5% at 1 week and sustained FI reductions of 31.5% from weeks 9-12, resulting in peak reduction in BW of 16 ± 5%. Dulaglutide had more modest effects on FI and peak BW loss was 3.8 ± 4.0%. Longitudinal modeling of both the FI and BW profiles suggested reductions in BW observed with both CpdH and dulaglutide were fully explained by exposure-dependent reductions in FI without increase in EE. Upon verification of the pharmacokinetic/pharmacodynamic relationship established in monkeys and humans for dulaglutide, we predicted that CpdH could reach double digit BW loss in humans. In summary, a long-acting GDF15 analog led to sustained reductions in FI in overweight monkeys and holds potential for effective clinical obesity pharmacotherapy.

Multivalent human antibody-centyrin fusion protein to prevent and treat Staphylococcus aureus infections.

Treating and preventing infections by antimicrobial-resistant bacterial pathogens is a worldwide problem. Pathogens such as Staphylococcus aureus produce an array of virulence determinants, making it difficult to identify single targets for the development of vaccines or monoclonal therapies. We described a human-derived anti-S. aureus monoclonal antibody (mAb)-centyrin fusion protein ("mAbtyrin") that simultaneously targets multiple bacterial adhesins, resists proteolysis by bacterial protease GluV8, avoids Fc engagement by S. aureus IgG-binding proteins SpA and Sbi, and neutralizes pore-forming leukocidins via fusion with anti-toxin centyrins, while maintaining Fc- and complement-mediated functions. Compared with the parental mAb, mAbtyrin protected human phagocytes and boosted phagocyte-mediated killing. The mAbtyrin also reduced pathology, reduced bacterial burden, and protected from different types of infections in preclinical animal models. Finally, mAbtyrin synergized with vancomycin, enhancing pathogen clearance in an animal model of bacteremia. Altogether, these data establish the potential of multivalent mAbs for treating and preventing S. aureus diseases.

The quest for balance between capturing data and model complexity: A quantitative clinical pharmacology approach applied to monoclonal antibodies.

The main objective of this tutorial is to provide the readers with a roadmap of how to establish increasingly complex target-mediated drug disposition (TMDD) models for monoclonal antibodies. To this end, we built mathematical models, each with a detailed visualization, starting from the basic TMDD model by Mager and Jusko to the well-established, physiologically based model by Li et al. in a step-wise fashion to highlight the relative importance of key physiological processes that impact mAb kinetics and system dynamics. As the models become more complex, the question of structural and parameter identifiability arises. To address this question, we work through a trastuzumab case example to guide the modeler's choice for model and parameter optimization in light of the context of use. We leave the readers of this tutorial with a brief summary of the advantages and limitations of each model expansion, as well as the model source codes for further self-guided exploration and hands-on analysis.

Diversity in Clinical Pharmacology: A Call to Action.

Prioritization of diversity, equity, and inclusion in all facets of our work is long overdue for the clinical pharmacology community. Increasing diversity in clinical research will deepen our understanding of nuanced patient populations and help improve all patient outcomes. Fostering an inclusive and diverse workforce will lead to broader perspectives that can better inform critical decisions and create work environments where everyone can thrive. In this call to action, we invite you to join us.

How resilient were we in 2021? Results of a LinkedIn Survey including biomedical and pharmaceutical professionals using the Benatti Resiliency Model.

Enhancing resiliency should elevate innovation and efficiency in biomedical research and development (R&D); however, compared with other professions, data on practice of resilience is lacking. Using the Benatti Resiliency Model (5 anchors: Well-Being, Self-Awareness, Brand, Connection, and Innovation), we surveyed professionals, including those in biomedical and pharmaceutical R&D. A structured LinkedIn questionnaire (March 16-May 23, 2021), surveyed each model anchor using five categories. One hundred fifty-eight participants (~6% student/trainee, 18%, 27%, and 49% in 1-5, 5-15 or >15 years post-terminal degree) took the survey (90 in biomedical and pharmaceutical R&D). Over 50% chose "always"/"often" across questions, except external influence or engagement. The question with one of the lowest "always" scores (~15%) was "I get feedback on my influence and impact in my career" in Brand, highlighting areas for leadership development and coaching. In the anchor of Well-being, nutrition and stress management also received some lowest "always" scores (~15% for both). Connection and Innovation scores trended slightly higher in biomedical and pharmaceutical R&D. No students/trainees chose "always" in Brand, indicating evolution of brand maturity over time. Self- and survey-assessed resiliency scores were associated (rs  = 0.37, p < 0.0001). Our survey yielded actionable insights on Resilience, including "best practices" through an open-ended question for one thing most useful to boost resilience in the survey and is the first application of the Benatti Model for crowdsourced research.

Changing the drug development and therapeutic paradigm with biologic drug combinations and bispecifics: How to choose between these two approaches?

Biologics are increasingly being co-developed in combination or as novel constructs like bispecific antibodies (BsAbs) with the goal of targeting multiple, non-redundant mechanisms of action. Rational design of combinations and dual-targeting approaches that consider disease complexities have the potential to improve efficacy and safety, to increase duration of clinical benefit, and to minimize clinical resistance mechanisms. Here we summarize examples of BsAbs and biologic combinations that have been approved by health authorities and present drug development considerations when deciding between these two strategies. These include an understanding of target biology, nonclinical safety risks, dose optimization strategies, the regulatory framework, pharmacokinetic, immunogenicity, and bioanalytical assay considerations. The disease biology, target dynamics, and pharmacology objectives were identified as important factors in early drug development to decide between a BsAb versus a combination. Nonclinical safety assessment and dose optimization strategies can also pose challenges for BsAb versus combinations. High unmet medical needs and lack of treatment options are often the common denominators for deciding to develop a BsAb or a combination. Future development of biologic triple combinations and BsAbs combinations with other biologics will further increase drug development complexities and hold promise for more effective treatment options for patients.

Design, synthesis and preclinical evaluation of bio-conjugated amylinomimetic peptides as long-acting amylin receptor agonists.

Pramlintide is an equipotent amylin analogue that reduces food intake and body weight in obese subjects and has been clinically approved as an adjunctive therapy for the treatment of adult diabetic patients. However, due to its extremely short half-life in vivo, a regimen of multiple daily administrations is required for achieving clinical effectiveness. Herein is described the development of prototypical long-acting pramlintide bioconjugates, in which pramlintide's disulfide-linked macrocycle was replaced by a cyclic thioether motif. This modification enabled stable chemical conjugation to a half-life extending antibody. In contrast to pramlintide (t1/2 < 0.75 h), bioconjugates 35 and 38 have terminal half-lives of ∼2 days in mice and attain significant exposure levels that are maintained up to 7 days. Single dose subcutaneous administration of 35 in lean mice, given 18-20 h prior to oral acetaminophen (AAP) administration, significantly reduced gastric emptying (as determined by plasma AAP levels). In a separate study, similar administration of 35 in fasted lean mice effected a reduction in food intake for up to 48 h. These data are consistent with durable amylinomimetic responses and provide the basis for further development of such long-acting amylinomimetic conjugates for the potential treatment of obesity and associated pathologies.

Minimal Physiologically Based Pharmacokinetic-Pharmacodynamic (mPBPK-PD) Model of N-Acetylgalactosamine-Conjugated Small Interfering RNA Disposition and Gene Silencing in Preclinical Species and Humans.

Conjugation of small interfering RNA (siRNA) to tris N-acetylgalactosamine [(GalNAc)3] can enable highly selective, potent, and durable knockdown of targeted proteins in the liver. However, potential knowledge gaps between in vitro experiments, preclinical species, and clinical scenarios remain. A minimal physiologically based pharmacokinetic-pharmacodynamic model for GalNAc-conjugated siRNA (GalNAc-siRNA) was developed using published data for fitusiran (ALN-AT3), an investigational compound targeting liver antithrombin (AT), to delineate putative determinants governing the whole-body-to-cellular pharmacokinetic (PK) and pharmacodynamic (PD) properties of GalNAc-siRNA and facilitate preclinical-to-clinical translation. The model mathematically linked relevant mechanisms: 1) hepatic biodistribution, 2) tris-GalNAc binding to asialoglycoprotein receptors (ASGPRs) on hepatocytes, 3) ASGPR endocytosis and recycling, 4) endosomal transport and escape of siRNA, 5) cytoplasmic RNA-induced silencing complex (RISC) loading, 6) degradation of target mRNA by bound RISC, and 7) knockdown of protein. Physiologic values for 36 out of 48 model parameters were obtained from the literature. Kinetic parameters governing (GalNAc)3-ASGPR binding and internalization were derived from published studies of uptake in hepatocytes. The proposed model well characterized reported pharmacokinetics, RISC dynamics, and knockdown of AT mRNA and protein by ALN-AT3 in mice. The model bridged multiple PK-PD data sets in preclinical species (mice, rat, monkey) and successfully captured reported plasma pharmacokinetics and AT knockdown in a phase I ascending-dose study. Estimates of in vivo potency were similar (∼2-fold) across species. Subcutaneous absorption and serum AT degradation rate constants scaled across species by body weight with allometric exponents of -0.29 and -0.22. The proposed mechanistic modeling framework characterizes the unique PK-PD properties of GalNAc-siRNA. SIGNIFICANCE STATEMENT: Tris N-acetylgalactosamine (GalNAc)3-conjugated small interfering RNA (siRNA) therapeutics enable liver-targeted gene therapy and precision medicine. Using a translational and systems-based minimal physiologically based pharmacokinetic-pharmacodynamic (mPBPK-PD) modeling approach, putative determinants influencing GalNAc-conjugated siRNA (GalNAc-siRNA) functionality in three preclinical species and humans were investigated. The developed model successfully integrated and characterized relevant published in vitro-derived biomeasures, mechanistic PK-PD profiles in animals, and observed clinical PK-PD responses for an investigational GalNAc-siRNA (fitusiran). This modeling effort delineates the disposition and liver-targeted pharmacodynamics of GalNAc-siRNA.

Model-Based Cellular Kinetic Analysis of Chimeric Antigen Receptor-T Cells in Humans.

Chimeric antigen receptor (CAR)-T cell therapy has achieved considerable success in treating B-cell hematologic malignancies. However, the challenges of extending CAR-T therapy to other tumor types, particularly solid tumors, remain appreciable. There are substantial variabilities in CAR-T cellular kinetics across CAR-designs, CAR-T products, dosing regimens, patient responses, disease types, tumor burdens, and lymphodepletion conditions. As a "living drug," CAR-T cellular kinetics typically exhibit four distinct phases: distribution, expansion, contraction, and persistence. The cellular kinetics of CAR-T may correlate with patient responses, but which factors determine CAR-T cellular kinetics remain poorly defined. Herein, we developed a cellular kinetic model to retrospectively characterize CAR-T kinetics in 217 patients from 7 trials and compared CAR-T kinetics across response status, patient populations, and tumor types. Based on our analysis results, CAR-T cells exhibited a significantly higher cell proliferation rate and capacity but a lower contraction rate in patients who responded to treatment. CAR-T cells proliferate to a higher degree in hematologic malignancies than in solid tumors. Within the assessed dose ranges (107 -109 cells), CAR-T doses were weakly correlated with CAR-T cellular kinetics and patient response status. In conclusion, the developed CAR-T cellular kinetic model adequately characterized the multiphasic CAR-T cellular kinetics and supported systematic evaluations of the potential influencing factors, which can have significant implications for the development of more effective CAR-T therapies.

A minimal physiologically based pharmacokinetic model to characterize colon TNF suppression and treatment effects of an anti-TNF monoclonal antibody in a mouse inflammatory bowel disease model.

Biotherapeutic drugs against tumor necrosis factor (TNF) are effective treatments for moderate to severe inflammatory bowel disease (IBD). Here, we evaluated CNTO 5048, an antimurine TNF surrogate monoclonal antibody (mAb), in a CD45RBhigh adoptive T cell transfer mouse colitis model, which allows examination of the early immunological events associated with gut inflammation and the therapeutic effects. The study was designed to quantitatively understand the effects of IBD on CNTO 5048 disposition, the ability of CNTO 5048 to neutralize pathogenic TNF at the colon under disease conditions, and the impact of dosing regimen on CNTO 5048 treatment effect. CNTO 5048 and TNF concentrations in both mice serum and colon homogenate were also measured. Free TNF concentrations in colon, but not in serum, were shown to correlate well with the colon pharmacodynamic readout, such as the summed histopathology score and neutrophil score. A minimal physiologically based pharmacokinetic (mPBPK) model was developed to characterize CNTO 5048 PK and disposition, as well as colon soluble TNF target engagement (TE). The mPBPK/TE model reasonably captured the observed data and provided a quantitative understanding of an anti-TNF mAb on its colon TNF suppression and therapeutic effect in a physiologically relevant IBD animal model. These results also provided insights into the potential benefits of using induction doses for the treatment of IBD patients.

Characterization of concurrent target suppression by JNJ-61178104, a bispecific antibody against human tumor necrosis factor and interleukin-17A.

Tumor necrosis factor (TNF) and interleukin (IL)-17A are pleiotropic cytokines implicated in the pathogenesis of several autoimmune diseases including rheumatoid arthritis (RA) and psoriatic arthritis (PsA). JNJ-61178104 is a novel human anti-TNF and anti-IL-17A monovalent, bispecific antibody that binds to both human TNF and human IL-17A with high affinities and blocks the binding of TNF and IL-17A to their receptors in vitro. JNJ-61178104 also potently neutralizes TNF and IL-17A-mediated downstream effects in multiple cell-based assays. In vivo, treatment with JNJ-61178104 resulted in dose-dependent inhibition of cellular influx in a human IL-17A/TNF-induced murine lung neutrophilia model and the inhibitory effects of JNJ-61178104 were more potent than the treatment with bivalent parental anti-TNF or anti-IL-17A antibodies. JNJ-61178104 was shown to engage its targets, TNF and IL-17A, in systemic circulation measured as drug/target complex formation in normal cynomolgus monkeys (cyno). Surprisingly, quantitative target engagement assessment suggested lower apparent in vivo target-binding affinities for JNJ-61178104 compared to its bivalent parental antibodies, despite their similar in vitro target-binding affinities. The target engagement profiles of JNJ-61178104 in humans were in general agreement with the predicted profiles based on cyno data, suggesting similar differences in the apparent in vivo target-binding affinities. These findings show that in vivo target engagement of monovalent bispecific antibody does not necessarily recapitulate that of the molar-equivalent dose of its bivalent parental antibody. Our results also offer valuable insights into the understanding of the pharmacokinetics/pharmacodynamics and target engagement of other bispecific biologics against dimeric and/or trimeric soluble targets in vivo.

Amino acids are sensitive glucagon receptor-specific biomarkers for glucagon-like peptide-1 receptor/glucagon receptor dual agonists.

AIM: The aim of this study was to evaluate amino acids as glucagon receptor (GCGR)-specific biomarkers in rodents and cynomolgus monkeys in the presence of agonism of both glucagon-like peptide-1 receptor (GLP1R) and GCGR with a variety of dual agonist compounds. MATERIALS AND METHODS: Primary hepatocytes, rodents (normal, diet-induced obese and GLP1R knockout) and cynomolgus monkeys were treated with insulin (hepatocytes only), glucagon (hepatocytes and cynomolgus monkeys), the GLP1R agonist, dulaglutide, or a variety of dual agonists with varying GCGR potencies. RESULTS: A long-acting dual agonist, Compound 2, significantly decreased amino acids in both wild-type and GLP1R knockout mice in the absence of changes in food intake, body weight, glucose or insulin, and increased expression of hepatic amino acid transporters. Dulaglutide, or a variant of Compound 2 lacking GCGR agonism, had no effect on amino acids. A third variant with ~31-fold less GCGR potency than Compound 2 significantly decreased amino acids, albeit to a significantly lesser extent than Compound 2. Dulaglutide (with saline infusion) had no effect on amino acids, but an infusion of glucagon dose-dependently decreased amino acids on the background of GLP1R engagement (dulaglutide) in cynomolgus monkeys, as did Compound 2. CONCLUSIONS: These results show that amino acids are sensitive and translatable GCGR-specific biomarkers.

Toward Progress in Quantitative Translational Medicine: A Call to Action.

Quantitative translational medicine (QTM) is envisioned as a multifaceted discipline that will galvanize the path from idea to medicine through quantitative translation across the discovery, development, regulatory, and utilization spectrum. Here, we summarize results of an American Society for Clinical Pharmacology and Therapeutics (ASCPT) survey on barriers relevant to the advancement of QTM and propose opportunities for its deployment. Importantly, we offer a call to action to break down these barriers through patient-centered stewardship, effective communication, cross-sector collaboration, and a modernized educational curriculum.

A Long-Acting PYY3-36 Analog Mediates Robust Anorectic Efficacy with Minimal Emesis in Nonhuman Primates.

The gut hormone PYY3-36 reduces food intake in humans and exhibits at least additive efficacy in combination with GLP-1. However, the utility of PYY analogs as anti-obesity agents has been severely limited by emesis and rapid proteolysis, a profile similarly observed with native PYY3-36 in obese rhesus macaques. Here, we found that antibody conjugation of a cyclized PYY3-36 analog achieved high NPY2R selectivity, unprecedented in vivo stability, and gradual infusion-like exposure. These properties permitted profound reduction of food intake when administered to macaques for 23 days without a single emetic event in any animal. Co-administration with the GLP-1 receptor agonist liraglutide for an additional 5 days further reduced food intake with only one animal experiencing a single bout of emesis. This antibody-conjugated PYY analog therefore may enable the long-sought potential of GLP-1/PYY-based combination treatment to achieve robust, well-tolerated weight reduction in obese patients.

Rapid Purification of Human Bispecific Antibodies via Selective Modulation of Protein A Binding.

Methods to rapidly generate high quality bispecific antibodies (BsAb) having normal half-lives are critical for therapeutic programs. Here, we identify 3 mutations (T307P, L309Q, and Q311R or "TLQ") in the Fc region of human IgG1 which disrupt interaction with protein A while enhancing interaction with FcRn. The mutations are shown to incrementally alter the pH at which a mAb elutes from protein A affinity resin. A BsAb comprised of a TLQ mutant and a wild-type IgG1 can be efficiently separated from contaminating parental mAbs by differential protein A elution starting from either a) purified parental mAbs, b) in-supernatant crossed parental mAbs, or c) co-transfected mAbs. We show that the Q311R mutation confers enhanced FcRn interaction in vitro, and Abs harboring either the Q311R or TLQ mutations have serum half-lives as long as wild-type human IgG1. The mutant Abs have normal thermal stability and Fcγ receptor interactions. Together, the results lead to a method for high-throughput generation of BsAbs suitable for in vivo studies.

Modulation of protein A binding allows single-step purification of mouse bispecific antibodies that retain FcRn binding.

The increased number of bispecific antibodies (BsAb) under therapeutic development has resulted in a need for mouse surrogate BsAbs. Here, we describe a one-step method for generating highly pure mouse BsAbs suitable for in vitro and in vivo studies. We identify two mutations in the mouse IgG2a and IgG2b Fc region: one that eliminates protein A binding and one that enhances protein A binding by 8-fold. We show that BsAbs harboring these mutations can be purified from the residual parental monoclonal antibodies in one step using protein A affinity chromatography. The structural basis for the effects of these mutations was analyzed by X-ray crystallography. While the mutation that disrupted protein A binding also inhibited FcRn interaction, a bispecific mutant in which one subunit retained the ability to bind protein A could still interact with FcRn. Pharmacokinetic analysis of the serum half-lives of the mutants showed that the mutant BsAb had a serum half-life comparable to a wild-type Ab. The results describe a rapid method for generating panels of mouse BsAbs that could be used in mouse studies.

Cross-arm binding efficiency of an EGFR x c-Met bispecific antibody.

Multispecific proteins, such as bispecific antibodies (BsAbs), that bind to two different ligands are becoming increasingly important therapeutic agents. Such BsAbs can exhibit markedly increased target binding and target residence time when both pharmacophores bind simultaneously to their targets. The cross-arm binding efficiency (χ) describes an increase in apparent affinity when a BsAb binds to the second target or receptor (R2) following its binding to the first target or receptor (R1) on the same cell. χ is an intrinsic characteristic of a BsAb mostly related to the binding epitopes on R1 and R2. χ can have significant impacts on the binding to R2 for BsAbs targeting two receptors on the same cell. JNJ-61186372, a BsAb that targets epidermal growth factor receptor (EGFR) and c-Met, was used as the model compound for establishing a method to characterize χ. The χ for JNJ-61186372 was successfully determined via fitting of in vitro cell binding data to a ligand binding model that incorporated χ. The model-derived χ value was used to predict the binding of JNJ-61186372 to individual EGFR and c-Met receptors on tumor cell lines, and the results agreed well with the observed IC50 for EGFR and c-Met phosphorylation inhibition by JNJ-61186372. Consistent with the model, JNJ-61186372 was shown to be more effective than the combination therapy of anti-EGFR and anti-c-Met monovalent antibodies at the same dose level in a mouse xenograft model. Our results showed that χ is an important characteristic of BsAbs, and should be considered for rationale design of BsAbs targeting two membrane bound targets on the same cell.

Evaluation of changes in oral drug absorption in preterm and term neonates for Biopharmaceutics Classification System (BCS) class I and II compounds.

AIMS: Evidence suggests that the rate of oral drug absorption changes during early childhood. Yet, respective clinical implications are currently unclear, particularly for preterm neonates. The objective of this study was to evaluate changes in oral drug absorption after birth for different Biopharmaceutics Classification System (BCS) class I and II compounds to better understand respective implications for paediatric pharmacotherapy. METHODS: Two paradigm compounds were selected for BCS class I (paracetamol (acetaminophen) and theophylline) and II (indomethacin and ibuprofen), respectively, based on the availability of clinical literature data following intravenous and oral dosing. A comparative population pharmacokinetic analysis was performed in a step-wise manner in NONMEM® 7.2 to characterize and predict changes in oral drug absorption after birth for paracetamol, theophylline and indomethacin. RESULTS: A one compartment model with an age-dependent maturation function for oral drug absorption was found appropriate to characterize the pharmacokinetics of paracetamol. Our findings indicate that the rate at which a drug is absorbed from the GI tract reaches adult levels within about 1 week after birth. The maturation function for paracetamol was found applicable to theophylline and indomethacin once solubility limitations were overcome via drug formulation. The influence of excipients on solubility and, hence, oral bioavailability was confirmed for ibuprofen, a second BCS class II compound. CONCLUSIONS: The findings of our study suggest that the processes underlying changes in oral drug absorption after birth are drug-independent and that the maturation function identified for paracetamol may be generally applicable to other BCS class I and II compounds for characterizing drug absorption in preterm as well as term neonates.

Utility of free and total target measurements as target engagement and efficacy biomarkers in biotherapeutic development--opportunities and challenges.

For biotherapeutics directed against soluble targets, most often monoclonal antibodies (mAbs), their therapeutic efficacy theoretically is driven by the magnitude and duration of free target suppression. However, for soluble targets of rapid turnover and low abundance, it can be technically challenging to directly measure the lowering of free target following treatment with biologics. The opportunities, challenges, and practical approaches to assess free and bound soluble targets and the utility of free and bound target measurements as biomarkers for target engagement and efficacy are covered in this review. In particular, case examples are presented to illustrate the interplay between drug and free/bound target, and how an integrated bioanalytical and pharmacokinetic/target engagement/pharmacodynamic (PK/TE/PD) modeling approach can be used to assess the target engagement for biologics directed against soluble targets with rapid turnover. Important caveats of the modeling approach in the absence of free target measurements are also discussed.

Development of a population pharmacokinetic model characterizing the tissue distribution of azithromycin in healthy subjects.

Recent clinical trials indicate that the use of azithromycin is associated with the emergence of macrolide resistance. The objective of our study was to simultaneously characterize free target site concentrations and correlate them with the MIC90s of clinically relevant pathogens. Azithromycin (500 mg once daily [QD]) was administered orally to 6 healthy male volunteers for 3 days. The free concentrations in the interstitial space fluid (ISF) of muscle and subcutaneous fat tissue as well as the total concentrations in plasma and polymorphonuclear leukocytes (PMLs) were determined on days 1, 3, 5, and 10. All concentrations were modeled simultaneously in NONMEM 7.2 using a tissue distribution model that accounts for nonlinear protein binding and ionization state at physiological pH. The model performance and parameter estimates were evaluated via goodness-of-fit plots and nonparametric bootstrap analysis. The model we developed described the concentrations at all sampling sites reasonably well and showed that the overall pharmacokinetics of azithromycin is driven by the release of the drug from acidic cell/tissue compartments. The model-predicted unionized azithromycin (AZM) concentrations in the cytosol of PMLs (6.0 ± 1.2 ng/ml) were comparable to the measured ISF concentrations in the muscle (8.7 ± 2.9 ng/ml) and subcutis (4.1 ± 2.4 ng/ml) on day 10, whereas the total PML concentrations were >1,000-fold higher (14,217 ± 2,810 ng/ml). The total plasma and free ISF concentrations were insufficient to exceed the MIC90s of the skin pathogens at all times. Our results indicate that the slow release of azithromycin from low pH tissue/cell compartments is responsible for the long terminal half-life of the drug and thus the extended period of time during which free concentrations reside at subinhibitory concentrations.