Minimal brain PBPK model to support the preclinical and clinical development of antibody therapeutics for CNS diseases

A minimal PBPK model to accelerate preclinical development of drugs against tuberculosis

Introduction: Understanding drug exposure at disease target sites is pivotal to profiling new drug candidates in terms of tolerability and efficacy. Such quantification is particularly tedious for anti-tuberculosis (TB) compounds as the heterogeneous pulmonary microenvironment due to the infection may alter lung permeability and affect drug disposition. Murine models have been a longstanding support in TB research so far and are here used as human surrogates to unveil the distribution of several anti-TB compounds at the site-of-action via a novel and centralized PBPK design framework. Methods: As an intermediate approach between data-driven pharmacokinetic (PK) models and whole-body physiologically based (PB) PK models, we propose a parsimonious framework for PK investigation (minimal PBPK approach) that retains key physiological processes involved in TB disease, while reducing computational costs and prior knowledge requirements. By lumping together pulmonary TB-unessential organs, our minimal PBPK model counts 9 equations compared to the 36 of published full models, accelerating the simulation more than 3-folds in Matlab 2022b. Results: The model has been successfully tested and validated against 11 anti-TB compounds-rifampicin, rifapentine, pyrazinamide, ethambutol, isoniazid, moxifloxacin, delamanid, pretomanid, bedaquiline, OPC-167832, GSK2556286 - showing robust predictability power in recapitulating PK dynamics in mice. Structural inspections on the proposed design have ensured global identifiability and listed free fraction in plasma and blood-to-plasma ratio as top sensitive parameters for PK metrics. The platform-oriented implementation allows fast comparison of the compounds in terms of exposure and target attainment. Discrepancies in plasma and lung levels for the latest BPaMZ and HPMZ regimens have been analyzed in terms of their impact on preclinical experiment design and on PK/PD indices. Conclusion: The framework we developed requires limited drug- and species-specific information to reconstruct accurate PK dynamics, delivering a unified viewpoint on anti-TB drug distribution at the site-of-action and a flexible fit-for-purpose tool to accelerate model-informed drug design pipelines and facilitate translation into the clinic.

Local Identifiability Analysis, Parameter Subset Selection and Verification for a Minimal Brain PBPK Model

Physiologically-based pharmacokinetic (PBPK) modeling is important for studying drug delivery in the central nervous system, including determining antibody exposure, predicting chemical concentrations at target locations, and ensuring accurate dosages. The complexity of PBPK models, involving many variables and parameters, requires a consideration of parameter identifiability; i.e., which parameters can be uniquely determined from data for a specified set of concentrations. We introduce the use of a local sensitivity-based parameter subset selection algorithm in the context of a minimal PBPK (mPBPK) model of the brain for antibody therapeutics. This algorithm is augmented by verification techniques, based on response distributions and energy statistics, to provide a systematic and robust technique to determine identifiable parameter subsets in a PBPK model across a specified time domain of interest. The accuracy of our approach is evaluated for three key concentrations in the mPBPK model for plasma, brain interstitial fluid and brain cerebrospinal fluid. The determination of accurate identifiable parameter subsets is important for model reduction and uncertainty quantification for PBPK models.

A Minimal PBPK Model for Plasma and Cerebrospinal Fluid Pharmacokinetics of Trastuzumab after Intracerebroventricular Administration in Patients with HER2-Positive Brain Metastatic Localizations

BACKGROUND

Dosing regimens of trastuzumab administered by intracerebroventricular (icv) route to patients with HER2-positive brain localizations remain empirical. The objectives of this study were to describe pharmacokinetics (PK) of trastuzumab in human plasma and cerebrospinal fluid (CSF) after simultaneous icv and intravenous (iv) administration using a minimal physiologically-based pharmacokinetic model (mPBPK) and to perform simulations of alternative dosing regimens to achieve therapeutic concentrations in CSF.

METHODS

Plasma and CSF PK data were collected in two patients with HER2-positive brain localizations. A mPBPK model for mAbs consisting of four compartments (tight and leaky tissues, plasma and lymph) was enriched by an additional compartment for ventricular CSF. The comparison between observed and model-predicted concentrations was evaluated using prediction error (PE).

RESULTS

The developed mPBPK model described plasma and CSF trastuzumab concentrations reasonably well with mean PE for plasma and CSF data of 41.8% [interquartile range, IQR = -9.48; 40.6] and 18.3% [-36.7; 60.6], respectively, for patient 1 and 11.4% [-10.8; 28.7] and 22.5% [-27.7; 77.9], respectively, for patient 2. Trastuzumab showed fast clearance from CSF to plasma with Cmin,ss of 0.56 and 0.85 mg/L for 100 and 150 mg q1wk, respectively. Repeated dosing of 100 and 150 mg q3day resulted in Cmin,ss of 10.3 and 15.4 mg/L, respectively. Trastuzumab CSF target concentrations are achieved rapidly and maintained above 60 mg/L from 7 days after a continuous perfusion at 1.0 mg/h.

CONCLUSION

Continuous icv infusion of trastuzumab at 1.0 mg/h could be an alternative dosing regimen to rapidly achieve intraventricular CSF therapeutic concentrations.

Cerebrospinal Fluid Pharmacokinetics of Nicardipine Following Intrathecal Administration in Subarachnoid Hemorrhage Patients

Subarachnoid hemorrhage (SAH) is a devastating type of stroke, leading to high mortality and morbidity rates. Cerebral vasospasm and delayed cerebral ischemia (DCI) are common complications following SAH and contribute significantly to the poor outcomes observed in these patients. Intrathecal (IT) nicardipine delivered via an existing external ventricular drain has been shown to be correlated with reduced DCI and improved patient outcomes. The current study aims to characterize population pharmacokinetic (popPK) properties of intermittent IT nicardipine. Following informed consent, serial cerebrospinal fluid (CSF) samples were obtained from 16 SAH patients (50.4 ± 9.3 years old; 12 females) treated with IT nicardipine every 6 hours (n=8) or every 8 hours (n=8), which were subject to high-performance liquid chromatography for measurement of its CSF concentration. Our popPK analysis showed that the CSF PK of IT nicardipine in the cohort was adequately described by a two-compartment model with a lag time, with reliable parameter estimates (relative standard error < 50%). The intracranial pressure influenced both the total clearance and the central volume. Calculated PK parameters were similar between q6h and q8h dosing regimens. Despite a small cohort of SAH patients, we successfully developed a popPK model to describe the nicardipine disposition kinetics in the CSF following IT administration. These findings may help inform future clinical trials designed to examine the optimal dosing of IT nicardipine.

Development of a physiologically based toxicokinetic model for lead in pregnant women: The role of bone tissue in the maternal and fetal internal exposure

Epidemiological studies have shown associations between prenatal exposure to lead (Pb) and neurodevelopmental effects in young children. Prenatal exposure is generally characterized by measuring the concentration in the umbilical cord at delivery or in the maternal blood during pregnancy. To assess internal Pb exposure during prenatal life, we developed a pregnancy physiologically based pharmacokinetic (p-PBPK) model that to simulates Pb levels in blood and target tissues in the fetus, especially during critical periods for brain development. An existing Pb PBPK model was adapted to pregnant women and fetuses. Using data from literature, both the additional maternal bone remodeling, that causes Pb release into the blood, and the Pb placental transfers were estimated by Bayesian inference. Additional maternal bone remodeling was estimated to start at 21.6 weeks. Placental transfers were estimated between 4.6 and 283 L.day-1 at delivery with high interindividual variability. Once calibrated, the p-PBPK model was used to simulate fetal exposure to Pb. Internal fetal exposure greatly varies over the pregnancy with two peaks of Pb levels in blood and brain at the end of the 1st and 3rd trimesters. Sensitivity analysis shows that the fetal blood lead levels are affected by the maternal burden of bone Pb via maternal bone remodeling and by fetal bone formation at different pregnancy stages. Coupling the p-PBPK model with an effect model such as an adverse outcome pathway could help to predict the effects on children's neurodevelopment.

Utility of Minimal Physiologically Based Pharmacokinetic Models for Assessing Fractional Distribution, Oral Absorption, and Series-Compartment Models of Hepatic Clearance

Minimal physiologically based pharmacokinetic (mPBPK) models are physiologically relevant, require less information than full PBPK models, and offer flexibility in pharmacokinetics (PK). The well-stirred hepatic model (WSM) is commonly used in PBPK, whereas the more plausible dispersion model (DM) poses computational complexities. The series-compartment model (SCM) mimics the DM but is easier to operate. This work implements the SCM and mPBPK models for assessing fractional tissue distribution, oral absorption, and hepatic clearance using literature-reported blood and liver concentration-time data in rats for compounds mainly cleared by the liver. Further handled were various complexities, including nonlinear hepatic binding and metabolism, differing absorption kinetics, and sites of administration. The SCM containing one to five (n) liver subcompartments yields similar fittings and provides comparable estimates for hepatic extraction ratio (ER), prehepatic availability (Fg ), and first-order absorption rate constants (ka ). However, they produce decreased intrinsic clearances (CLint ) and liver-to-plasma partition coefficients (Kph ) with increasing n as expected. Model simulations demonstrated changes in intravenous and oral PK profiles with alterations in Kph and ka and with hepatic metabolic zonation. The permeability (PAMPA P) of the various compounds well explained the fitted fractional distribution (fd ) parameters. The SCM and mPBPK models offer advantages in distinguishing systemic, extrahepatic, and hepatic clearances. The SCM allows for incorporation of liver zonation and is useful in assessing changes in internal concentration gradients potentially masked by similar blood PK profiles. Improved assessment of intraorgan drug concentrations may offer insights into active moieties driving metabolism, biliary excretion, pharmacodynamics, and hepatic toxicity. SIGNIFICANCE STATEMENT: The minimal physiologically based pharmacokinetic model and the series-compartment model are useful in assessing oral absorption and hepatic clearance. They add flexibility in accounting for various drug- or system-specific complexities, including fractional distribution, nonlinear binding and saturable hepatic metabolism, and hepatic zonation. These models can offer improved insights into the intraorgan concentrations that reflect physiologically active moieties often driving disposition, pharmacodynamics, and toxicity.

Analysis of clinical failure of anti-tau and anti-synuclein antibodies in neurodegeneration using a quantitative systems pharmacology model

Misfolded proteins in Alzheimer's disease and Parkinson's disease follow a well-defined connectomics-based spatial progression. Several anti-tau and anti-alpha synuclein (aSyn) antibodies have failed to provide clinical benefit in clinical trials despite substantial target engagement in the experimentally accessible cerebrospinal fluid (CSF). The proposed mechanism of action is reducing neuronal uptake of oligomeric protein from the synaptic cleft. We built a quantitative systems pharmacology (QSP) model to quantitatively simulate intrasynaptic secretion, diffusion and antibody capture in the synaptic cleft, postsynaptic membrane binding and internalization of monomeric and oligomeric tau and aSyn proteins. Integration with a physiologically based pharmacokinetic (PBPK) model allowed us to simulate clinical trials of anti-tau antibodies gosuranemab, tilavonemab, semorinemab, and anti-aSyn antibodies cinpanemab and prasineuzumab. Maximal target engagement for monomeric tau was simulated as 45% (semorinemab) to 99% (gosuranemab) in CSF, 30% to 99% in ISF but only 1% to 3% in the synaptic cleft, leading to a reduction of less than 1% in uptake of oligomeric tau. Simulations for prasineuzumab and cinpanemab suggest target engagement of free monomeric aSyn of only 6-8% in CSF, 4-6% and 1-2% in the ISF and synaptic cleft, while maximal target engagement of aggregated aSyn was predicted to reach 99% and 80% in the synaptic cleft with similar effects on neuronal uptake. The study generates optimal values of selectivity, sensitivity and PK profiles for antibodies. The study identifies a gradient of decreasing target engagement from CSF to the synaptic cleft as a key driver of efficacy, quantitatively identifies various improvements for drug design and emphasizes the need for QSP modelling to support the development of tau and aSyn antibodies.

Minimal physiologically-based hybrid model of pharmacokinetics in pregnant women: Application to antenatal corticosteroids

Minimal physiologically-based pharmacokinetic (mPBPK) models are an alternative to full physiologically-based pharmacokinetic (PBPK) models as they offer reduced complexity while maintaining the physiological interpretation of key model components. Full PBPK models have been developed for pregnancy, but a mPBPK model eases the ability to perform a "top-down" meta-analysis melding all available pharmacokinetic (PK) data in the mother and fetus. Our hybrid mPBPK model consists of mPBPK models for the mother and fetus with connection by the placenta. This model was applied to describe the rich PK data of antenatal corticosteroid betamethasone (BET) jointly with the limited data for dexamethasone (DEX) in the mother and fetus. Physiologic model parameters were obtained from the literature while drug-dependent parameters were estimated by the simultaneous fitting of all available data for DEX and BET. Maternal clearances of DEX and BET confirmed the literature values, and the expected fetal-to-maternal plasma ratios ranged from 0.3 to 0.4 for both drugs. Simulations of maternal plasma concentrations for the dosing regimens of BET and DEX recommended by the World Health Organization based on our findings revealed up to 60% lower exposures than found in nonpregnant women and offers a means of devising alternative dosing regimens. Our hybrid mPBPK model and meta-analysis approach could facilitate assessment of other classes of drugs indicated for the treatment of pregnant women.

A combined physiologically-based pharmacokinetic and quantitative systems pharmacology model for modeling amyloid aggregation in Alzheimer's disease

Antibody-mediated removal of aggregated β-amyloid (Aβ) is the current, most clinically advanced potential disease-modifying treatment approach for Alzheimer's disease. We describe a quantitative systems pharmacology (QSP) approach of the dynamics of Aβ monomers, oligomers, protofibrils, and plaque using a detailed microscopic model of Aβ₄₀ and Aβ₄₂ aggregation and clearance of aggregated Aβ by activated microglia cells, which is enhanced by the interaction of antibody-bound Aβ. The model allows for the prediction of Aβ positron emission tomography (PET) imaging load as measured by a standardized uptake value ratio. A physiology-based pharmacokinetic model is seamlessly integrated to describe target exposure of monoclonal antibodies and simulate dynamics of cerebrospinal fluid (CSF) and plasma biomarkers, including CSF Aβ₄₂ and plasma Aβ₄₂ /Aβ₄₀ ratio biomarkers. Apolipoprotein E genotype is implemented as a difference in microglia clearance. By incorporating antibody-bound, plaque-mediated macrophage activation in the perivascular compartment, the model also predicts the incidence of amyloid-related imaging abnormalities with edema (ARIA-E). The QSP platform is calibrated with pharmacological and clinical information on aducanumab, bapineuzumab, crenezumab, gantenerumab, lecanemab, and solanezumab, predicting adequately the change in PET imaging measured amyloid load and the changes in the plasma Aβ₄₂ /Aβ₄₀ ratio while slightly overestimating the change in CSF Aβ₄₂ . ARIA-E is well predicted for all antibodies except bapineuzumab. This QSP model could support the clinical trial design of different amyloid-modulating interventions, define optimal titration and maintenance schedules, and provide a first step to understand the variability of biomarker response in clinical practice.

Application of physiologically-based pharmacokinetic models for therapeutic proteins and other novel modalities

The past two decades have seen diversification of drug development pipelines and approvals from traditional small molecule therapies to alternative modalities including monoclonal antibodies, engineered proteins, antibody drug conjugates (ADCs), oligonucleotides and gene therapies. At the same time, physiologically-based pharmacokinetic (PBPK) models for small molecules have seen increased industry and regulatory acceptance.This review focusses on the current status of the application of PBPK models to these newer modalities and give a perspective on the successes, challenges and future directions of this field.There is greatest experience in the development of PBPK models for therapeutic proteins, and PBPK models for ADCs benefit from prior experience for both therapeutic proteins and small molecules. For other modalities, the application of PBPK models is in its infancy.Challenges are discussed and a common theme is lack of availability of physiological and experimental data to characterise systems and drug parameters to enable a priori prediction of pharmacokinetics. Furthermore, sufficient clinical data are required to build confidence in developed models.The PBPK modelling approach provides a quantitative framework for integrating knowledge and data from multiple sources and can be built on as more data becomes available.

PBPK-PD modeling for the preclinical development and clinical translation of tau antibodies for Alzheimer’s disease

Disrupted tau proteostasis and transneuronal spread is a pathological hallmark of Alzheimer’s disease. Neurodegenerative diseases remain an unmet medical need and novel disease modifying therapeutics are paramount. Our objective was to develop a mechanistic mathematical model to enhance our understanding of tau antibody pharmacokinetics and pharmacodynamics in animals and humans. A physiologically-based pharmacokinetic-pharmacodynamic (PBPK-PD) modeling approach was employed to support the preclinical development and clinical translation of therapeutic antibodies targeting tau for the treatment of Alzheimer’s disease. The pharmacokinetics of a tau antibody was evaluated in rat and non-human primate microdialysis studies. Model validation for humans was performed using publicly available clinical data for gosuranemab. In-silico analyses were performed to predict tau engagement in human brain for a range of tau antibody affinities and various dosing regimens. PBPK-PD modeling enabled a quantitative understanding for the relationship between dose, affinity, and target engagement, which supported lead candidate optimization and predictions of clinically efficacious dosing regimens.

Investigating brain uptake of a non-targeting monoclonal antibody after intravenous and intracerebroventricular administration

Monoclonal antibodies play an important role in the treatment of various diseases. However, the development of these drugs against neurological disorders where the drug target is located in the brain is challenging and requires a good understanding of the local drug concentration in the brain. In this original research, we investigated the systemic and local pharmacokinetics in the brain of healthy rats after either intravenous (IV) or intracerebroventricular (ICV) administration of EGFRvIII-T-Cell bispecific (TCB), a bispecific monoclonal antibody. We established an experimental protocol that allows serial sampling in serum, cerebrospinal fluid (CSF) and interstitial fluid (ISF) of the prefrontal cortex in freely moving rats. For detection of drug concentration in ISF, a push-pull microdialysis technique with large pore membranes was applied. Brain uptake into CSF and ISF was characterized and quantified with a reduced brain physiologically-based pharmacokinetic model. The model allowed us to interpret the pharmacokinetic processes of brain uptake after different routes of administration. The proposed model capturing the pharmacokinetics in serum, CSF and ISF of the prefrontal cortex suggests a barrier function between the CSF and ISF that impedes free antibody transfer. This finding suggests that ICV administration may not be better suited to reach higher local drug exposure as compared to IV administration. The model enabled us to quantify the relative contribution of the blood-brain barrier (BBB) and Blood-CSF-Barrier to the uptake into the interstitial fluid of the brain. In addition, we compared the brain uptake of three monoclonal antibodies after IV dosing. In summary, the presented approach can be applied to profile compounds based on their relative uptake in the brain and provides quantitative insights into which pathways are contributing to the net exposure in the brain.

Revisiting Cerebrospinal Fluid Flow Direction and Rate in Physiologically Based Pharmacokinetic Model

The bidirectional pulsatile movement of cerebrospinal fluid (CSF), instead of the traditionally believed unidirectional and constant CSF circulation, has been demonstrated. In the present study, the structure and parameters of the CSF compartments were revisited in our comprehensive and validated central nervous system (CNS)-specific, physiologically based pharmacokinetic (PBPK) model of healthy rats (LeiCNS-PK3.0). The bidirectional and site-dependent CSF movement was incorporated into LeiCNS-PK3.0 to create the new LeiCNS-PK“3.1” model. The physiological CSF movement rates in healthy rats that are unavailable from the literature were estimated by fitting the PK data of sucrose, a CSF flow marker, after intra-CSF administration. The capability of LeiCNS-PK3.1 to describe the PK profiles of other molecules was compared with that of the original LeiCNS-PK3.0 model. LeiCNS-PK3.1 demonstrated superior description of the CSF PK profiles of a range of small molecules after intra-CSF administration over LeiCNS-PK3.0. LeiCNS-PK3.1 also retained the same level of predictability of CSF PK profiles in cisterna magna after intravenous administration. These results support the theory of bidirectional and site-dependent CSF movement across the entire CSF space over unidirectional and constant CSF circulation in healthy rats, pointing out the need to revisit the structures and parameters of CSF compartments in CNS-PBPK models.