Validity conditions of approximations for a target-mediated drug disposition model: A novel first-order approximation and its comparison to other approximations.

Target-mediated drug disposition (TMDD) is a phenomenon characterized by a drug's high-affinity binding to a target molecule, which significantly influences its pharmacokinetic profile within an organism. The comprehensive TMDD model delineates this interaction, yet it may become overly complex and computationally demanding in the absence of specific concentration data for the target or its complexes. Consequently, simplified TMDD models employing quasi-steady state approximations (QSSAs) have been introduced; however, the precise conditions under which these models yield accurate results require further elucidation. Here, we establish the validity of three simplified TMDD models: the Michaelis-Menten model reduced with the standard QSSA (mTMDD), the QSS model reduced with the total QSSA (qTMDD), and a first-order approximation of the total QSSA (pTMDD). Specifically, we find that mTMDD is applicable only when initial drug concentrations substantially exceed total target concentrations, while qTMDD can be used for all drug concentrations. Notably, pTMDD offers a simpler and faster alternative to qTMDD, with broader applicability than mTMDD. These findings are confirmed with antibody-drug conjugate real-world data. Our findings provide a framework for selecting appropriate simplified TMDD models while ensuring accuracy, potentially enhancing drug development and facilitating safer, more personalized treatments.

Physicochemical characterization and phase II metabolic profiling of echinochrome A, a bioactive constituent from sea urchin, and its physiologically based pharmacokinetic modeling in rats and humans.

Echinochrome A, a natural naphthoquinone pigment found in sea urchins, is increasingly being investigated for its nutritional and therapeutic value associated with antioxidant, anticancer, antiviral, antidiabetic, and cardioprotective activities. Although several studies have demonstrated the biological effects and therapeutic potential of echinochrome A, little is known regarding its biopharmaceutical behaviors. Here, we aimed to investigate the physicochemical properties and metabolic profiles of echinochrome A and establish a physiologically-based pharmacokinetic (PBPK) model as a useful tool to support its clinical applications. We found that the lipophilicity, color variability, ultraviolet/visible spectrometry, and stability of echinochrome A were markedly affected by pH conditions. Moreover, metabolic and pharmacokinetic profiling studies demonstrated that echinochrome A is eliminated primarily by hepatic metabolism and that four possible metabolites, i.e., two glucuronidated and two methylated conjugates, are formed in rat and human liver preparations. A whole-body PBPK model incorporating the newly identified hepatic phase II metabolic process was constructed and optimized with respect to chemical-specific parameters. Furthermore, model simulations suggested that echinochrome A could exhibit linear disposition profiles without systemic and local tissue accumulation in clinical settings. Our proposed PBPK model of echinochrome A could be a valuable tool for predicting drug interactions in previously unexplored scenarios and for optimizing dosage regimens and drug formulations.

Artificial neural network models: implementation of functional near-infrared spectroscopy-based spontaneous lie detection in an interactive scenario.

Deception is an inevitable occurrence in daily life. Various methods have been used to understand the mechanisms underlying brain deception. Moreover, numerous efforts have been undertaken to detect deception and truth-telling. Functional near-infrared spectroscopy (fNIRS) has great potential for neurological applications compared with other state-of-the-art methods. Therefore, an fNIRS-based spontaneous lie detection model was used in the present study. We interviewed 10 healthy subjects to identify deception using the fNIRS system. A card game frequently referred to as a bluff or cheat was introduced. This game was selected because its rules are ideal for testing our hypotheses. The optical probe of the fNIRS was placed on the subject's forehead, and we acquired optical density signals, which were then converted into oxy-hemoglobin and deoxy-hemoglobin signals using the Modified Beer-Lambert law. The oxy-hemoglobin signal was preprocessed to eliminate noise. In this study, we proposed three artificial neural networks inspired by deep learning models, including AlexNet, ResNet, and GoogleNet, to classify deception and truth-telling. The proposed models achieved accuracies of 88.5%, 88.0%, and 90.0%, respectively. These proposed models were compared with other classification models, including k-nearest neighbor, linear support vector machines (SVM), quadratic SVM, cubic SVM, simple decision trees, and complex decision trees. These comparisons showed that the proposed models performed better than the other state-of-the-art methods.

A Robust Computer-Aided Automated Brain Tumor Diagnosis Approach Using PSO-ReliefF Optimized Gaussian and Non-Linear Feature Space.

Brain tumors are among the deadliest diseases in the modern world. This study proposes an optimized machine-learning approach for the detection and identification of the type of brain tumor (glioma, meningioma, or pituitary tumor) in brain images recorded using magnetic resonance imaging (MRI). The Gaussian features of the image are extracted using speed-up robust features (SURF), whereas its non-linear features are obtained using KAZE, owing to their high performance against rotation, scaling, and noise problems. To retrieve local-level information, all brain MRI images are segmented into an 8 × 8 pixel grid. To enhance the accuracy and reduce the computational time, the variance-based k-means clustering and PSO-ReliefF algorithms are employed to eliminate the redundant features of the brain MRI images. Finally, the performance of the proposed hybrid optimized feature vector is evaluated using various machine learning classifiers. An accuracy of 96.30% is obtained with 169 features using a support vector machine (SVM). Furthermore, the computational time is also reduced to 1 min compared to the non-optimized features used for training of the SVM. The findings are also compared with previous research, demonstrating that the suggested approach might assist physicians and doctors in the timely detection of brain tumors.

Fractional transit compartment model for describing drug delayed response to tumors using Mittag-Leffler distribution on age-structured PKPD model.

The response of a cell population is often delayed relative to drug injection, and individual cells in a population of cells have a specific age distribution. The application of transit compartment models (TCMs) is a common approach for describing this delay. In this paper, we propose a TCM in which damaged cells caused by a drug are given by a single fractional derivative equation. This model describes the delay as a single equation composed of fractional and ordinary derivatives, instead of a system of ODEs expressed in multiple compartments, applicable to the use of the PK concentration in the model. This model tunes the number of compartments in the existing model and expresses the delay in detail by estimating an appropriate fractional order. We perform model robustness, sensitivity analysis, and change of parameters based on the amount of data. Additionally, we resolve the difficulty in parameter estimation and model simulation using a semigroup property, consisting of a system with a mixture of fractional and ordinary derivatives. This model provides an alternative way to express the delays by estimating an appropriate fractional order without determining the pre-specified number of compartments.

Extended transit compartment model to describe tumor delay using Coxian distribution.

The measured response of cell population is often delayed relative to drug injection, and individuals in a population have a specific age distribution. Common approaches for describing the delay are to apply transit compartment models (TCMs). This model reflects that all damaged cells caused by drugs suffer transition processes, resulting in death. In this study, we present an extended TCM using Coxian distribution, one of the phase-type distributions. The cell population attacked by a drug is described via age-structured models. The mortality rate of the damaged cells is expressed by a convolution of drug rate and age density. Then applying to Erlang and Coxian distribution, we derive Erlang TCM, representing the existing model, and Coxian TCMs, reflecting sudden death at all ages. From published data of drug and tumor, delays are compared after parameter estimations in both models. We investigate the dynamical changes according to the number of the compartments. Model robustness and equilibrium analysis are also performed for model validation. Coxian TCM is an extended model considering a realistic case and captures more diverse delays.

Development and application of a physiologically based pharmacokinetic model for entrectinib in rats and scale-up to humans: Route-dependent gut wall metabolism.

Entrectinib (Rozlytrek®) is an oral antineoplastic agent approved by the U.S. Food and Drug Administration in 2019 for the treatment of c-ros oncogene 1 (ROS1)-positive non-small cell lung cancer and neurotrophic tyrosine receptor kinase (NTRK) fusion-positive solid tumors. Although there have been a few studies on the pharmacokinetics of entrectinib, the relative contributions of several kinetic factors determining the oral bioavailability and systemic exposure of entrectinib are still worthy of investigation. Experimental data on the intestinal absorption and disposition of entrectinib in rats were acquired from studies on in vitro protein binding/tissue S9 metabolism, in situ intestinal perfusion, and in vivo dose-escalation/hepatic extraction. Using these datasets, an in-house whole-body physiologically based pharmacokinetic (PBPK) model incorporating the QGut model concepts and segregated blood flow in the gut was constructed and optimized with respect to drug-specific parameters. The established rat PBPK model was further extrapolated to humans through relevant physiological scale-up and parameter optimization processes. The optimized rat and human PBPK models adequately captured the impact of route-dependent gut metabolism on the systemic exposure to entrectinib and closely mirrored various preclinical and clinical observations. Our proposed PBPK model could be useful in optimizing dosage regimens and predicting drug interaction potential in various clinical conditions, after partial modification and validation.

A Tumor-Immune Interaction Model for Synergistic Combinations of Anti PD-L1 and Ionizing Irradiation Treatment.

Combination therapy with immune checkpoint blockade and ionizing irradiation therapy (IR) generates a synergistic effect to inhibit tumor growth better than either therapy does alone. We modeled the tumor-immune interactions occurring during combined IT and IR based on the published data from Deng et al. The mathematical model considered programmed cell death protein 1 and programmed death ligand 1, to quantify data fitting and global sensitivity of critical parameters. Fitting of data from control, IR and IT samples was conducted to verify the synergistic effect of a combination therapy consisting of IR and IT. Our approach using the model showed that an increase in the expression level of PD-1 and PD-L1 was proportional to tumor growth before therapy, but not after initiating therapy. The high expression level of PD-L1 in T cells may inhibit IT efficacy. After combination therapy begins, the tumor size was also influenced by the ratio of PD-1 to PD-L1. These results highlight that the ratio of PD-1 to PD-L1 in T cells could be considered in combination therapy.

An immune therapy model for effective treatment on inflammatory bowel disease.

Inflammatory bowel disease (IBD) is a disease that causes inflammation throughout the digestive tract. Repeated inflammation and frequent relapses cause intestinal damage and expose the patient to a higher risk. In this work, we proposed an immune therapy model for effective treatment strategy through mathematical modeling for patients with IBD. We evaluated the ability of the patient's immune system to recover during treatment. For this, we defined the interval of healthy individual, and examined the frequency of compartments such as T cells and cytokines considered in the model maintain the normal state. Based on the fact that each patient has a unique immune system, we have shown at the same drug works differently, depending on the individual immune system characteristics for every patient. It is known that IBD is related to an imbalance between pro- and anti- inflammatory cytokines as the cause of the disease. So the ratios of pro- to anti- inflammatory cytokines are used as an indicator of patient's condition and inflammation status in various diseases. We compared the ratios of pro- to anti- inflammatory cytokine according to patient's individual immune system and drugs. Since the effects of biological drugs are highly dependent on the patient's own immune system, it is essential to define the immune system status before selecting and using a biological drug.

Evaluation of COVID-19 epidemic outbreak caused by temporal contact-increase in South Korea.

OBJECTIVES: On March 15, 2020, 61.3% of the confirmed cases of COVID-19 infection in South Korea are associated with the worship service that was organized on February 9 in the Shincheonji Church of Jesus in Daegu. We aim to evaluate the effects of mass infection in South Korea and assess the preventive control intervention. METHOD: Using openly available data of daily cumulative confirmed cases and deaths, the basic and effective reproduction numbers was estimated using a modified susceptible-exposed-infected-recovered-type epidemic model. RESULTS: The basic reproduction number was estimated to be R0=1.77. The effective reproduction number increased approximately 20 times after the mass infections from the 31 st patient, which was confirmed on February 9 in the Shincheonji Church of Jesus, Daegu. However, the effective reproduction number decreased to less than unity after February 28 owing to the implementation of high-level preventive control interventions in South Korea, coupled with voluntary prevention actions by citizens. CONCLUSION: Preventive action and control intervention were successfully established in South Korea.

A mathematical model for assessing the effectiveness of controlling relapse in Plasmodium vivax malaria endemic in the Republic of Korea.

Malaria has persisted as an endemic near the Demilitarized Zone in the Republic of Korea since the re-emergence of Plasmodium vivax malaria in 1993. The number of patients affected by malaria has increased recently despite many controls tools, one of the reasons behind which is the relapse of malaria via liver hypnozoites. Tafenoquine, a new drug approved by the United States Food and Drug Administration in 2018, is expected to reduce the rate of relapse of malaria hypnozoites and thereby decrease the prevalence of malaria among the population. In this work, we have developed a new transmission model for Plasmodium vivax that takes into account a more realistic intrinsic distribution from existing literature to quantify the current values of relapse parameters and to evaluate the effectiveness of the anti-relapse therapy. The model is especially suitable for estimating parameters near the Demilitarized Zone in Korea, in which the disease follows a distinguishable seasonality. Results were shown that radical cure could significantly reduce the prevalence level of malaria. However, eradication would still take a long time (over 10 years) even if the high-level treatment were to persist. In addition, considering that the vector's behavior is manipulated by the malaria parasite, relapse repression through vector control at the current level may result in a negative effect in containing the disease. We conclude that the use of effective drugs should be considered together with the increased level of the vector control to reduce malaria prevalence.

Recent advances in physiologically based pharmacokinetic and pharmacodynamic models for anticancer nanomedicines.

Nanoparticles (NPs) have distinct pharmacokinetic (PK) properties and can potentially improve the absorption, distribution, metabolism, and elimination (ADME) of small-molecule drugs loaded therein. Owing to the unwanted toxicities of anticancer agents in healthy organs and tissues, their precise delivery to the tumor is an essential requirement. There have been numerous advancements in the development of nanomedicines for cancer therapy. Physiologically based PK (PBPK) models serve as excellent tools for describing and predicting the ADME properties and the efficacy and toxicity of drugs, in combination with pharmacodynamic (PD) models. The recent preliminary application of these modeling approaches to NPs demonstrated their potential benefits in research and development processes relevant to the ADME and pharmacodynamics of NPs and nanomedicines. Here, we comprehensively review the pharmacokinetics of NPs, the developed PBPK models for anticancer NPs, and the developed PD model for anticancer agents.

Modeling to capture bystander-killing effect by released payload in target positive tumor cells.

BACKGROUND: Antibody-drug conjugates (ADCs) are intended to bind to specific positive target antigens and eradicate only tumor cells from an intracellular released payload through the lysosomal protease. Payloads, such as MMAE, have the capacity to kill adjacent antigen-negative (Ag-) tumor cells, which is called the bystander-killing effect, as well as directly kill antigen-positive (Ag+) tumor cells. We propose that a dose-response curve should be independently considered to account for target antigen-positive/negative tumor cells. METHODS: A model was developed to account for the payload in Ag+/Ag- cells and the associated parameters were applied. A tumor growth inhibition (TGI) effect was explored based on an ordinary differential equation (ODE) after substituting the payload concentration in Ag+/Ag- cells into an Emax model, which accounts for the dose-response curve. To observe the bystander-killing effects based on the amount of Ag+/Ag- cells, the Emax model is used independently. TGI models based on ODE are unsuitable for describing the initial delay through a tumor-drug interaction. This was solved using an age-structured model based on the stochastic process. RESULTS: ß∈(0,1] is a fraction parameter that determines the proportion of cells that consist of Ag+/Ag- cells. The payload concentration decreases when the ratio of efflux to influx increases. The bystander-killing effect differs with varying amounts of Ag+ cells. The larger ß is, the less bystander-killing effect. The decrease of the bystander-killing effect becomes stronger as Ag+ cells become larger than the Ag- cells. Overall, the ratio of efflux to influx, the amount of released payload, and the proportion of Ag+ cells determine the efficacy of the ADC. The tumor inhibition delay through a payload-tumor interaction, which goes through several stages, may be solved using an age-structured model. CONCLUSIONS: The bystander-killing effect, one of the most important topics of ADCs, has been explored in several studies without the use of modeling. We propose that the bystander-killing effect can be captured through a mathematical model when considering the Ag+ and Ag- cells. In addition, the TGI model based on the age-structure can capture the initial delay through a drug interaction as well as the bystander-killing effect.

Mathematical modeling of antibody drug conjugates with the target and tubulin dynamics to predict AUC.

Antibody drug conjugates (ADCs)are one of the most recently developed chemotherapeutics to treat some types of tumor cells. They consist of monoclonal antibodies (mAbs), linkers, and potent cytotoxic drugs. Unlike common chemotherapies, ADCs combine selectively with a target at the surface of the tumor cell, and a potent cytotoxic drug (payload) effectively prevents microtubule polymerization. In this work, we construct an ADC model that considers both the target of antibodies and the receptor (tubulin) of the cytotoxic payloads. The model is simulated with brentuximab vedotin, one of ADCs, and used to investigate the pharmacokinetic (PK) characteristics of ADCs in vivo. It also predicts area under the curve (AUC) of ADCs and the payloads by identifying the half-life. The results show that dynamical behaviors fairly coincide with the observed data and half-life and capture AUC. Thus, the model can be used for estimating some parameters, fitting experimental observations, predicting AUC, and exploring various dynamical behaviors of the target and the receptor.