Model-Informed Precision Dosing for Personalized Ustekinumab Treatment in Plaque Psoriasis.

BACKGROUND/OBJECTIVES: Implementing model-informed precision dosing (MIPD) strategies guided by population pharmacokinetic/pharmacodynamic (PK/PD) models could enhance the management of inflammatory diseases such as psoriasis. However, the extent of individual experimental data gathered during MIPD significantly influences the uncertainty in estimating individual PK/PD parameters, affecting clinical dose selection decisions. METHODS: This study proposes a methodology to individualize ustekinumab (UTK) dosing strategies for 23 Spanish patients with moderate to severe chronic plaque psoriasis., considering the uncertainty of individual parameters within a population PK/PD model. RESULTS: An indirect response model from previous research was used to describe the PK/PD relationship between UTK serum concentrations and the Psoriasis Area and Severity Index (PASI) score. A maximum inhibition drug effect (Imax) model was selected, and a first-order remission constant rate of psoriatic skin lesion (kout = 0.016 d-1) was estimated. CONCLUSIONS: The MIPD approach predicted that 35% and 26% of the patients would need an optimized and intensified dosage regimen, respectively, compared to the regimen typically used in clinical practice. This analysis demonstrated its utility as a tool for selecting personalized UTK dosing regimens in clinical practice in order to optimize the probability of achieving targeted clinical outcomes in patients with psoriasis.

Population Pharmacokinetic-Pharmacodynamic Analysis of a Reserpine-Induced Myalgia Model in Rats.

(1) Background: Fibromyalgia syndrome (FMS) is a chronic pain condition with widespread pain and multiple comorbidities, for which conventional therapies offer limited benefits. The reserpine-induced myalgia (RIM) model is an efficient animal model of FMS in rodents. This study aimed to develop a pharmacokinetic-pharmacodynamic (PK-PD) model of reserpine in rats, linking to its impact on monoamines (MAs). (2) Methods: Reserpine was administered daily for three consecutive days at dose levels of 0.1, 0.5, and 1 mg/kg. A total of 120 rats were included, and 120 PK and 828 PD observations were collected from 48 to 96 h after the first dose of reserpine. Non-linear mixed-effect data analysis was applied for structural PK-PD model definition, variability characterization, and covariate analysis. (3) Results: A one-compartment model best described reserpine in rats (V = 1.3 mL/kg and CL = 4.5 × 10-1 mL/h/kg). A precursor-pool PK-PD model (kin = 6.1 × 10-3 mg/h, kp = 8.6 × 10-4 h-1 and kout = 2.7 × 10-2 h-1) with a parallel transit chain (k0 = 1.9 × 10-1 h-1) characterized the longitudinal levels of MA in the prefrontal cortex, spinal cord, and amygdala in rats. Reserpine stimulates the degradation of MA from the pool compartment (Slope1 = 1.1 × 10-1 h) and the elimination of MA (Slope2 = 1.25 h) through the transit chain. Regarding the reference dose (1 mg/kg) of the RIM model, the administration of 4 mg/kg would lead to a mean reduction of 65% (Cmax), 80% (Cmin), and 70% (AUC) of MA across the brain regions tested. (4) Conclusions: Regional brain variations in neurotransmitter depletion were identified, particularly in the amygdala, offering insights for therapeutic strategies and biomarker identification in FMS research.

Impact of Time on Parameters for Assessing the Microstructure Equivalence of Topical Products: Diclofenac 1% Emulsion as a Case Study.

The demonstration of bioequivalence proposed in the European Medicines Agency's (EMA's) draft guideline for topical products with the same qualitative and quantitative composition requires the confirmation of the internal structure equivalence. The impact of the shelf-life on the parameters proposed for internal structure comparison has not been studied. The objectives of this work were: (1) to quantify the effect of the time since manufacturing on the mean value and variability of the parameters proposed by the EMA to characterize the internal structure and performance of topical formulations of a complex topical formulation, and (2) to evaluate the impact of these changes on the assessment of the microstructure equivalence. A total of 5 batches of a topical emulgel containing 1% diclofenac diethylamine were evaluated 5, 14, and 23 months after manufacture. The zero-shear viscosity (η0), viscosity at 100 s-1 (η100), yield stress (σ0), elastic (G') and viscous (G″) moduli, internal phase droplet size and in vitro release of the active ingredient were characterized. While no change in variability over time was detected, the mean value of all the parameters changed, especially the droplet size and in vitro release. Thus, combining data from batches of different manufacturing dates may compromise the determination of bioequivalence. The results confirm that to assess the microstructural similarity of complex formulations (such as emulgel), the 90% confidence interval limit for the mean difference in rheological and in vitro release parameters should be 20% and 25%, respectively.

Bioequivalence risk assessment of oral formulations containing racemic ibuprofen through a chiral physiologically based pharmacokinetic model of ibuprofen enantiomers.

The characterization of the time course of ibuprofen enantiomers can be useful in the selection of the most sensitive analyte in bioequivalence studies. Physiologically based pharmacokinetic (PBPK) modelling and simulation represents the most efficient methodology to virtually assess bioequivalence outcomes. In this work, we aim to develop and verify a PBPK model for ibuprofen enantiomers administered as a racemic mixture with different immediate release dosage forms to anticipate bioequivalence outcomes based on different particle size distributions. A PBPK model incorporating stereoselectivity and non-linearity in plasma protein binding and metabolism as well as R-to-S unidirectional inversion has been developed in Simcyp®. A dataset composed of 11 Phase I clinical trials with 54 scenarios (27 per enantiomer) and 14,452 observations (7129 for R-ibuprofen and 7323 for S-ibuprofen) was used. Prediction errors for AUC0-t and Cmax for both enantiomers fell within the 0.8-1.25 range in 50/54 (93 %) and 42/54 (78 %) of scenarios, respectively. Outstanding model performance, with 10/10 (100 %) of Cmax and 9/10 (90 %) of AUC0-t within the 0.9-1.1 range, was demonstrated for oral suspensions, which strongly supported its use for bioequivalence risk assessment. The deterministic bioequivalence risk assessment has revealed R-ibuprofen as the most sensitive analyte to detect differences in particle size distribution for oral suspensions containing 400 mg of racemic ibuprofen, suggesting that achiral bioanalytical methods would increase type II error and declare non-bioequivalence for formulations that are bioequivalent for the eutomer.

Towards precision medicine of long-acting aripiprazole through population pharmacokinetic modelling.

Population pharmacokinetic (popPK) models constitute a valuable tool for characterizing the pharmacokinetic properties of once-monthly long-acting injectable aripiprazole (LAI aripiprazole) and quantifying the sources of variability in drug exposure. Our aim is to develop a popPK model of both aripiprazole and its metabolite dehydro-aripiprazole in patients treated with LAI aripiprazole, and to personalize the dosing regimen of aripiprazole across different sub-groups of patients. This is a prospective study investigating the pharmacokinetics of LAI aripiprazole. A total of 93 patients were included, 21 for model development and 71 for external model evaluation. A one-compartment model with linear absorption and elimination adequately described both aripiprazole and dehydro-aripiprazole concentrations. The weight of the patients has been shown to be the factor that most influences the absorption. However, the metabolizing phenotype for CYP2D6 and the concomitant treatment with strong inhibitors of this cytochrome have been shown to be the covariates that most influence total drug exposure. This is the first popPK model developed for LAI aripiprazole that includes aripiprazole and its main active metabolite, dehydroaripiprazole. It provides a personalized dosage recommendation that maximizes the probability of achieving optimal therapeutic concentrations and minimizes the difficulties associated with trial-and-error therapeutic strategies carried out in clinical practice.

Integrating Mechanistic and Toxicokinetic Information in Predictive Models of Cholestasis.

Drug development involves the thorough assessment of the candidate's safety and efficacy. In silico toxicology (IST) methods can contribute to the assessment, complementing in vitro and in vivo experimental methods, since they have many advantages in terms of cost and time. Also, they are less demanding concerning the requirements of product and experimental animals. One of these methods, Quantitative Structure-Activity Relationships (QSAR), has been proven successful in predicting simple toxicity end points but has more difficulties in predicting end points involving more complex phenomena. We hypothesize that QSAR models can produce better predictions of these end points by combining multiple QSAR models describing simpler biological phenomena and incorporating pharmacokinetic (PK) information, using quantitative in vitro to in vivo extrapolation (QIVIVE) models. In this study, we applied our methodology to the prediction of cholestasis and compared it with direct QSAR models. Our results show a clear increase in sensitivity. The predictive quality of the models was further assessed to mimic realistic conditions where the query compounds show low similarity with the training series. Again, our methodology shows clear advantages over direct QSAR models in these situations. We conclude that the proposed methodology could improve existing methodologies and could be suitable for being applied to other toxicity end points.

New insights into the role of VKORC1 polymorphisms for optimal warfarin dose selection in Caribbean Hispanic patients through an external validation of a population PK/PD model.

Warfarin, an oral anticoagulant, has been used for decades to prevent thromboembolic events. The complex interplay between CYP2C9 and VKORC1 genotypes on warfarin PK and PD properties is not fully understood in special sub-groups of patients. This study aimed to externally validate a population pharmacokinetic/pharmacodynamic (PK/PD) model for the effect of warfarin on international normalized ratio (INR) and to evaluate optimal dosing strategies based on the selected covariates in Caribbean Hispanic patients. INR, and CYP2C9 and VKORC1 genotypes from 138 patients were used to develop a population PK/PD model in NONMEM. The structural definition of a previously published PD model for INR was implemented. A numerical evaluation of the parameter-covariate relationship was performed. Simulations were conducted to determine optimal dosing strategies for each genotype combinations, focusing on achieving therapeutic INR levels. Findings revealed elevated IC50 for G/G, G/A, and A/A VKORC1 haplotypes (11.76, 10.49, and 9.22 mg/L, respectively), in this population compared to previous reports. The model-guided dosing analysis recommended daily warfarin doses of 3-5 mg for most genotypes to maintain desired INR levels, although subjects with combination of CYP2C9 and VKORC1 genotypes * 2/* 2-, * 2/* 3- and * 2/* 5-A/A would require only 1 mg daily. This research underscores the potential of population PK/PD modeling to inform personalized warfarin dosing in populations typically underrepresented in clinical studies, potentially leading to improved treatment outcomes and patient safety. By integrating genetic factors and clinical data, this approach could pave the way for more effective and tailored anticoagulation therapy in diverse patient groups.

Mechanistic characterization of oscillatory patterns in unperturbed tumor growth dynamics: The interplay between cancer cells and components of tumor microenvironment.

Mathematical modeling of unperturbed and perturbed tumor growth dynamics (TGD) in preclinical experiments provides an opportunity to establish translational frameworks. The most commonly used unperturbed tumor growth models (i.e. linear, exponential, Gompertz and Simeoni) describe a monotonic increase and although they capture the mean trend of the data reasonably well, systematic model misspecifications can be identified. This represents an opportunity to investigate possible underlying mechanisms controlling tumor growth dynamics through a mathematical framework. The overall goal of this work is to develop a data-driven semi-mechanistic model describing non-monotonic tumor growth in untreated mice. For this purpose, longitudinal tumor volume profiles from different tumor types and cell lines were pooled together and analyzed using the population approach. After characterizing the oscillatory patterns (oscillator half-periods between 8-11 days) and confirming that they were systematically observed across the different preclinical experiments available (p<10-9), a tumor growth model was built including the interplay between resources (i.e. oxygen or nutrients), angiogenesis and cancer cells. The new structure, in addition to improving the model diagnostic compared to the previously used tumor growth models (i.e. AIC reduction of 71.48 and absence of autocorrelation in the residuals (p>0.05)), allows the evaluation of the different oncologic treatments in a mechanistic way. Drug effects can potentially, be included in relevant processes taking place during tumor growth. In brief, the new model, in addition to describing non-monotonic tumor growth and the interaction between biological factors of the tumor microenvironment, can be used to explore different drug scenarios in monotherapy or combination during preclinical drug development.

Usage of model combination in computational toxicology.

New Approach Methodologies (NAMs) have ushered in a new era in the field of toxicology, aiming to replace animal testing. However, despite these advancements, they are not exempt from the inherent complexities associated with the study's endpoint. In this review, we have identified three major groups of complexities: mechanistic, chemical space, and methodological. The mechanistic complexity arises from interconnected biological processes within a network that are challenging to model in a single step. In the second group, chemical space complexity exhibits significant dissimilarity between compounds in the training and test series. The third group encompasses algorithmic and molecular descriptor limitations and typical class imbalance problems. To address these complexities, this work provides a guide to the usage of a combination of predictive Quantitative Structure-Activity Relationship (QSAR) models, known as metamodels. This combination of low-level models (LLMs) enables a more precise approach to the problem by focusing on different sub-mechanisms or sub-processes. For mechanistic complexity, multiple Molecular Initiating Events (MIEs) or levels of information are combined to form a mechanistic-based metamodel. Regarding the complexity arising from chemical space, two types of approaches were reviewed to construct a fragment-based chemical space metamodel: those with and without structure sharing. Metamodels with structure sharing utilize unsupervised strategies to identify data patterns and build low-level models for each cluster, which are then combined. For situations without structure sharing due to pharmaceutical industry intellectual property, the use of prediction sharing, and federated learning approaches have been reviewed. Lastly, to tackle methodological complexity, various algorithms are combined to overcome their limitations, diverse descriptors are employed to enhance problem definition and balanced dataset combinations are used to address class imbalance issues (methodological-based metamodels). Remarkably, metamodels consistently outperformed classical QSAR models across all cases, highlighting the importance of alternatives to classical QSAR models when faced with such complexities.

Application of population physiologically based pharmacokinetic modelling to optimize target expression and clearance mechanisms of therapeutic monoclonal antibodies.

AIMS: To use population physiologically based pharmacokinetic (PopPBPK) modelling to optimize target expression, kinetics and clearance of HER1/2 directed therapeutic monoclonal antibodies (mAbs). Thus, to propose a general workflow of PopPBPK modelling and its application in clinical pharmacology. METHODS: Full PBPK model of pertuzumab (PTZ) was developed in patient population using Simcyp V21R1 incorporating mechanistic targeted-mediated drug disposition process by fitting known clinical PK and sparse receptor proteomics data to optimize target expression and kinetics of HER2 receptor. Trastuzumab (TTZ) PBPK modelling was used to validate the optimized HER2 target. Additionally, the simulator was also used to develop a full PBPK model for the HER1-directed mAb cetuximab (CTX) to assess the underlying targeted-mediated drug disposition-independent elimination mechanisms. RESULTS: HER2 final parameterisation coming from the PBPK modelling of PTZ was successfully cross validated through PBPK modelling of TTZ with average fold error (AFE), absolute AFE and percent prediction error values for area under the concentration-time curve (AUC) and maximum plasma concentration (Cmax ) of 1.13, 1.16 and 16, and 1.01, 1.07 and 7, respectively. CTX PBPK model performance was validated after the incorporation of an additional systemic clearance of 0.033 L/h as AFE and absolute AFE showed an acceptable predictive power of AUC and Cmax with percent prediction error of 13% for AUC and 10% for Cmax . CONCLUSIONS: Optimisation of both system and drug related parameters were performed through PBPK modelling to improve model performance of therapeutic mAbs (PTZ, TTZ and CTX). General workflow was proposed to develop and apply PopPBPK to support clinical development of mAbs targeting same receptor.

Alternative Pharmacokinetic Metrics in Single-Dose Studies to Ensure Bioequivalence of Prolonged-Release Products at Steady State-A Case Study.

(1) Background: this article investigates which PK metrics in a single-dose study (concentration at the end of posology interval, Cτ, partial areas under the curve, pAUCs, or half-value duration, HVD) are more sensitive and less variable for predicting the failure of a prolonged-release product at steady-state that was the bioequivalent for Cmax, AUC0-t and AUC0-inf, in the single-dose study; (2) Methods: a cross-over study was performed in 36 subjects receiving desvenlafaxine 100 mg prolonged-release tablets. Conventional (Cmax, AUC0-t and AUC0-inf) and additional (Cτ, pAUCs and HVD) PK metrics were considered after single-dose conditions. Predicted PK metrics at steady state (AUC0-τ, Cmax,ss, and Cτ,ss) were derived using a population PK model approach; (3) Results: the existing differences in the shape of the concentration-time curves precluded to show equivalence for Cτ,ss in the simulated study at steady state. This failure to show equivalence at steady state was predicted by Cτ, pAUCs and HVD in the single-dose study. Cτ was the most sensitive metric for detecting the different shape, with a lower intra-subject variability than HVD; (4) Conclusions: conventional PK metrics for single-dose studies (Cmax, AUC0-t and AUC0-inf) are not enough to guarantee bioequivalence at steady state for prolonged-release products.

Regulatory Requirements for the Development of Second-Entry Semisolid Topical Products in the European Union.

The development of second-entry topical products is hampered by several factors. The excipient composition should be similar to the reference product because excipients may also contribute to efficacy. Conventional pharmacokinetic bioequivalence studies were not considered acceptable because drug concentrations are measured downstream after the site of action. There was no agreed methodology to characterize the microstructure of semisolids, and waivers of therapeutic equivalence studies with clinical endpoints were not possible. Only the vasoconstrictor assay for corticosteroids was accepted as a surrogate. This paper describes the implementation of the European Union's stepwise approach for locally acting products to cutaneous products, discusses the equivalence requirements of the EMA Draft Guideline on the Quality and Equivalence of Topical Products, and compares them with the US Food and Drug Administration recommendations. Step 1 includes the possibility of waivers for simple formulations based on in vitro data only (Q1 + Q2 + Q3 + IVRT). Step 2 includes step 1 requirements plus a kinetic study (TS/IVPT/PKBE) to compare the local availability of complex formulations. Step 3 refers to clinical studies with pharmacodynamic/clinical endpoints. As excipients may affect the local tolerability and efficacy of the products, the similarity of excipient composition is required in all steps, except where clinical endpoints are compared.

Population Pharmacokinetic/Pharmacodynamic Modelling of Daptomycin for Schedule Optimization in Patients with Renal Impairment.

The aims of this study are (i) to develop a population pharmacokinetic/pharmacodynamic model of daptomycin in patients with normal and impaired renal function, and (ii) to establish the optimal dose recommendation of daptomycin in clinical practice. Several structural PK models including linear and non-linear binding kinetics were evaluated. Monte Carlo simulations were conducted with a fixed combination of creatinine clearance (30-90 mL/min/1.73 m2) and body weight (50-100 kg). The final dataset included 46 patients and 157 daptomycin observations. A two-compartment model with first-order peripheral distribution and elimination kinetics assuming non-linear protein-binding kinetics was selected. The bactericidal effect for Gram+ strains with MIC ≤ 0.5 mg/L could be achieved with 5-12 mg/kg daily daptomycin based on body weight and renal function. The administration of 10-17 mg/kg q48 h daptomycin allows to achieve bactericidal effect for Gram+ strains with MIC ≤ 1 mg/L. Four PK samples were selected as the optimal sampling strategy for an accurate AUC estimation. A quantitative framework has served to characterize the non-linear binding kinetics of daptomycin in patients with normal and impaired renal function. The impact of different dosing regimens on the efficacy and safety outcomes of daptomycin treatment based on the unbound exposure of daptomycin and individual patient characteristics has been evaluated.

A physiologically based pharmacokinetic model for open acid and lactone forms of atorvastatin and metabolites to assess the drug-gene interaction with SLCO1B1 polymorphisms.

Atorvastatin is the most prescribed 3-hydroxy-3-methylglutaryl coenzyme A reductase inhibitor used to lower cardiovascular risk and constitutes one of the best-selling drugs world-wide. Several physiologically based pharmacokinetic (PBPK) models have been developed to assess its non-straightforward pharmacokinetics (PK) as well as that of its metabolites and have been only applied to assess drug-drug interactions (DDI). Here we present a full PBPK model for atorvastatin and its metabolites able to predict within a 2-fold error their PK after the administration of a solid oral dosage form containing the calcium salt of atorvastatin in single and multiple dosing schedules at 20, 40, and 80 mg and 10 mg dose levels, respectively. Internal validation with data from Phase 1 clinical trials as well as external validation in predicting clinically relevant DDIs consolidated model structure and parameterization. The model has been used to quantitatively assess the drug-gene interaction (DGI) between SLCO1B1 polymorphisms and atorvastatin exposure and revealed that patients with a reduced activity in hepatic uptake of atorvastatin are at increased risk of suffering muscle discomfort because of a 30% lower clearance (p < 0.01), leading to a 40% and 33% higher (p < 0.05) atorvastatin AUC and Cmax, respectively. These findings could explain the reported hazard ratio of 1.4 (95% CI: 1.1-1.7, p = 0.02) for suffering statin-induced myopathies and the treatment discontinuation among these patients (odds ratio 1.67, p = 0.0001) observed in the context of routine clinical care.

A Review of Population Pharmacokinetic Analyses of Linezolid.

In recent years, many studies on population pharmacokinetics of linezolid have been conducted. This comprehensive review aimed to summarize population pharmacokinetic models of linezolid, by focusing on dosage optimization to maximize the probability of attaining a certain pharmacokinetic-pharmacodynamic parameter in special populations. We searched the PubMed and EMBASE databases for population pharmacokinetic analyses of linezolid using a parametric non-linear mixed-effect approach, including both observational and prospective trials. Of the 32 studies, 26 were performed in adults, four in children, and one in both adults and children. High between-subject variability was determined in the majority of the models, which was in line with the variability of linezolid concentrations previously detected in observational studies. Some studies found that patients with renal impairment, hepatic failure, advanced age, or low body weight had higher exposure and adverse reactions rates. In contrast, lower concentrations and therapeutic failure were associated with obese patients, young patients, and patients who had undergone renal replacement techniques. In critically ill patients, the inter-individual and intra-individual variability was even greater, suggesting that this population is at an even higher risk of underexposure and overexposure. Therapeutic drug monitoring may be warranted in a large proportion of patients given that the Monte Carlo simulations demonstrated that the one-size-fits-all labeled dosing of 600 mg every 12 h could lead to toxicity or therapeutic failure for high values of the minimum inhibitory concentration of the target pathogen. Further research on covariates, including renal function, hepatic function, and drug-drug interactions related to P-glycoprotein could help to explain variability and improve linezolid dosing regimens.

Pharmacometric characterization of entero-hepatic circulation processes of orally administered formulations of amiodarone under complex binding kinetics.

AIMS: The aims of this work are (i) to characterize the absorption properties of orally administered formulations at different dose levels, and (ii) to evaluate the impact of entero-hepatic circulation on the time-course of amiodarone (AM) in rats in order to optimize the development of new oral (OR) formulations. METHODS: Intravenous (IV) formulation consisted on a solution of a commercial injectable of AM chlorhydrate. OR formulations included the IV commercial formulation (Trangorex®) (Solution I), an aqueous supramicellar solution of AM chlorhydrate with Polysorbate at 5% (Solution II) and a suspension from Trangorex® tablets (Tablet). Data from 96 male Wistar rats, including 985 AM observations, were analyzed using NONMEM v7.4. RESULTS: The population pharmacokinetic (PK) model assumes linear absorption processes, showing ka of AM from Solution II (Polysorbate 80, 5%) and Solution I increased by 2.5- and 1.62-fold compared to Tablet formulation. OR bioavailability of AM from Tablet, Solution I and Solution II was 37%, 40%, and 50%, respectively. The structural model of AM disposition was adapted from a previously population PK model and expanded by incorporating entero-hepatic reabsorption (EHR) processes, which estimated a 12.3% biliary excretion of AM and complete re-absorption from lumen. CONCLUSIONS: The current population PK model of AM demonstrated the absorption rate enhancement when AM is formulated with supramicellar concentrations of Polysorbate 80. The study design allowed to characterize the EHR of AM and its contribution in the overall AM disposition.

Developing Clinically Relevant Dissolution Specifications (CRDSs) for Oral Drug Products: Virtual Webinar Series.

A webinar series that was organised by the Academy of Pharmaceutical Sciences Biopharmaceutics focus group in 2021 focused on the challenges of developing clinically relevant dissolution specifications (CRDSs) for oral drug products. Industrial scientists, together with regulatory and academic scientists, came together through a series of six webinars, to discuss progress in the field, emerging trends, and areas for continued collaboration and harmonisation. Each webinar also hosted a Q&A session where participants could discuss the shared topic and information. Although it was clear from the presentations and Q&A sessions that we continue to make progress in the field of CRDSs and the utility/success of PBBM, there is also a need to continue the momentum and dialogue between the industry and regulators. Five key areas were identified which require further discussion and harmonisation.

Impact of Pharmacokinetic and Pharmacodynamic Properties of Monoclonal Antibodies in the Management of Psoriasis.

The treatment of psoriasis has been revolutionized by the emergence of biological therapies. Monoclonal antibodies (mAb) generally have complex pharmacokinetic (PK) properties with nonlinear distribution and elimination. In recent years, several population pharmacokinetic/pharmacodynamic (PK/PD) models capable of describing different types of mAb have been published. This study aims to summarize the findings of a literature search about population PK/PD modeling and therapeutic drug monitoring (TDM) of mAb in psoriasis. A total of 22 articles corresponding to population PK/PD models of tumor necrosis factor (TNF)-α inhibitors (adalimumab and golimumab), interleukin (IL)-23 inhibitors (guselkumab, tildrakizumab, and risankizumab), IL-23/IL-12 inhibitor (ustekinumab), and IL-17 inhibitors (secukinumab, ixekizumab, and brodalumab) were collected. A summary of the clinical trials conducted so far in psoriasis was included, together with the current structural population PK and PD models. The most significant and clinical covariates were body weight (BW) and the presence of immunogenicity on clearance (CL). The lack of consensus on PK/PD relationships has prevented establishing an adequate dosage and, therefore, accentuates the need for TDM in psoriasis.