Modeling Sex Differences in Pharmacokinetics, Pharmacodynamics, and Disease Progression Effects of Naproxen in Rats with Collagen-Induced Arthritis

The Role of Pharmacometrics in Advancing the Therapies for Autoimmune Diseases

Autoimmune diseases (AIDs) are a group of disorders in which the immune system attacks the body's own tissues, leading to chronic inflammation and organ damage. These diseases are difficult to treat due to variability in drug PK among individuals, patient responses to treatment, and the side effects of long-term immunosuppressive therapies. In recent years, pharmacometrics has emerged as a critical tool in drug discovery and development (DDD) and precision medicine. The aim of this review is to explore the diverse roles that pharmacometrics has played in addressing the challenges associated with DDD and personalized therapies in the treatment of AIDs. Methods: This review synthesizes research from the past two decades on pharmacometric methodologies, including Physiologically Based Pharmacokinetic (PBPK) modeling, Pharmacokinetic/Pharmacodynamic (PK/PD) modeling, disease progression (DisP) modeling, population modeling, model-based meta-analysis (MBMA), and Quantitative Systems Pharmacology (QSP). The incorporation of artificial intelligence (AI) and machine learning (ML) into pharmacometrics is also discussed. Results: Pharmacometrics has demonstrated significant potential in optimizing dosing regimens, improving drug safety, and predicting patient-specific responses in AIDs. PBPK and PK/PD models have been instrumental in personalizing treatments, while DisP and QSP models provide insights into disease evolution and pathophysiological mechanisms in AIDs. AI/ML implementation has further enhanced the precision of these models. Conclusions: Pharmacometrics plays a crucial role in bridging pre-clinical findings and clinical applications, driving more personalized and effective treatments for AIDs. Its integration into DDD and translational science, in combination with AI and ML algorithms, holds promise for advancing therapeutic strategies and improving autoimmune patients' outcomes.

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.

Indigenous Nigeria medicinal herbal remedies: A potential source for therapeutic against rheumatoid arthritis

Rheumatoid arthritis (RA) is a debilitating disease associated with locomotion impairment, and conventional therapeutic drugs are not optimal for managing RA. There is an avalanche of medications used for the management of RA. Still, studies have shown that they are associated with severe side effects, including hepatotoxicity, retinopathy, and cardiotoxicity disorders of the central nervous system (CNS), skin, blood, and infections. Complementary and alternative medicine (CAM) is currently gaining attention as a novel panacea for managing debilitating diseases, such as RA. Nigerian folk herbal remedies are replete with a plethora of curative medicine, albeit unvalidated scientifically but with seemingly miraculous provenance. Studies of the identification of bioactive compounds present in these botanicals using advanced spectral analytical techniques have enhanced our understanding of the role of Nigerian herbal remedies in the treatment and management of RA. Interestingly, experimental studies abound that the bioactive compounds present in the extracts of plant botanicals protected animals from the development of RA in different experimental models and reduced the toxicity associated with conventional therapeutics. Validated mechanisms of RA amelioration in human and animal models include suppression of the expression of NF-κB, IL-1β, TNF-α, IL-6, IL-8, IL-17, IL-23, chemokines, TGF-β, RANKL, RANK, iNOS, arginase, COX-2, VEGFA, VEGFR, NFATC1, and TRAP in the synoviocytes. Decreased ROS, NO, MDA, carbonyl groups, and PGE₂ in the synovial fluid increased the expression of PPARα/γ; antioxidant and anti-inflammatory molecules also improve RA etiology. In this mini-review, we discuss the global burden of RA, the novel role of plant-based botanicals as potential therapeutics against signaling pathways in RA. Also addressed is the possible repurposing/reprofiling of plant botanicals to increase their therapeutic index among RA patients that patronize traditional healers in Nigeria with a global projection.

Teaching computational systems biology with an eye on quantitative systems pharmacology at the undergraduate level: Why do it, who would take it, and what should we teach?

Computational systems biology (CSB) is a field that emerged primarily as the product of research activities. As such, it grew in several directions in a distributed and uncoordinated manner making the area appealing and fascinating. The idea of not having to follow a specific path but instead creating one fueled innovation. As the field matured, several interdisciplinary graduate programs emerged attempting to educate future generations of computational systems biologists. These educational initiatives coordinated the dissemination of information across student populations that had already decided to specialize in this field. However, we are now entering an era where CSB, having established itself as a valuable research discipline, is attempting the next major step: Entering undergraduate curricula. As interesting as this endeavor may sound, it has several difficulties, mainly because the field is not uniformly defined. In this manuscript, we argue that this diversity is a significant advantage and that several incarnations of an undergraduate-level CSB biology course could, and should, be developed tailored to programmatic needs. In this manuscript, we share our experiences creating a course as part of a Biomedical Engineering program.

Sex-Related Differences in Drugs with Anti-Inflammatory Properties

There is increasing evidence of sex differences in the action of anti-inflammatory drugs, with women being at significantly higher risk of adverse effects. Nevertheless, clinicians' awareness of the implications of these sex differences on dosing and adverse event monitoring in routine practice is still in need of improvement. We reviewed the literature evaluating sex differences in terms of pharmacokinetics and pharmacodynamics of anti-inflammatory drugs. The anti-thrombotic activity of selective and non-selective COX-inhibitors tends to be stronger in men than women. Side effect profiles differ with regards to gastro-intestinal, renal and hepatic complications. Glucocorticosteroids were found to be more effective in men; women were more sensitive to corticosteroids when their oestradiol levels were high, a finding important for women taking hormonal contraception. TNF-alpha inhibitors have a longer half-life in men, leading to stronger immunosuppression and this a higher incidence of infections as side effects. Although research on sex differences in the effectiveness and safety of drugs is increasing, findings are often anecdotal and controversial. There is no systematic sex-differentiated reporting from clinical trials, and women are often under-represented. As personalized medicine is gaining in importance, sex, and gender aspects need to become integral parts of future research and policy making.

Transitioning from Basic toward Systems Pharmacodynamic Models: Lessons from Corticosteroids

Technology in bioanalysis, -omics, and computation have evolved over the past half century to allow for comprehensive assessments of the molecular to whole body pharmacology of diverse corticosteroids. Such studies have advanced pharmacokinetic and pharmacodynamic (PK/PD) concepts and models that often generalize across various classes of drugs. These models encompass the "pillars" of pharmacology, namely PK and target drug exposure, the mass-law interactions of drugs with receptors/targets, and the consequent turnover and homeostatic control of genes, biomarkers, physiologic responses, and disease symptoms. Pharmacokinetic methodology utilizes noncompartmental, compartmental, reversible, physiologic [full physiologically based pharmacokinetic (PBPK) and minimal PBPK], and target-mediated drug disposition models using a growing array of pharmacometric considerations and software. Basic PK/PD models have emerged (simple direct, biophase, slow receptor binding, indirect response, irreversible, turnover with inactivation, and transduction models) that place emphasis on parsimony, are mechanistic in nature, and serve as highly useful "top-down" methods of quantitating the actions of diverse drugs. These are often components of more complex quantitative systems pharmacology (QSP) models that explain the array of responses to various drugs, including corticosteroids. Progressively deeper mechanistic appreciation of PBPK, drug-target interactions, and systems physiology from the molecular (genomic, proteomic, metabolomic) to cellular to whole body levels provides the foundation for enhanced PK/PD to comprehensive QSP models. Our research based on cell, animal, clinical, and theoretical studies with corticosteroids have provided ideas and quantitative methods that have broadly advanced the fields of PK/PD and QSP modeling and illustrates the transition toward a global, systems understanding of actions of diverse drugs. SIGNIFICANCE STATEMENT: Over the past half century, pharmacokinetics (PK) and pharmacokinetics/pharmacodynamics (PK/PD) have evolved to provide an array of mechanism-based models that help quantitate the disposition and actions of most drugs. We describe how many basic PK and PK/PD model components were identified and often applied to the diverse properties of corticosteroids (CS). The CS have complications in disposition and a wide array of simple receptor-to complex gene-mediated actions in multiple organs. Continued assessments of such complexities have offered opportunities to develop models ranging from simple PK to enhanced PK/PD to quantitative systems pharmacology (QSP) that help explain therapeutic and adverse CS effects. Concurrent development of state-of-the-art PK, PK/PD, and QSP models are described alongside experimental studies that revealed diverse CS actions.

ATLAS mPBPK: A MATLAB-Based Tool for Modeling and Simulation of Minimal Physiologically-Based Pharmacokinetic Models

Minimal physiologically-based pharmacokinetic (mPBPK) models are frequently used to model plasma pharmacokinetic (PK) data and utilize and yield physiologically relevant parameters. Compared with classical compartment and whole-body physiologically-based pharmacokinetic modeling approaches, mPBPK models maintain a structure of intermediate physiological complexity that can be adequately informed by plasma PK data. In this tutorial, we present a MATLAB-based tool for the modeling and simulation of mPBPK models (ATLAS mPBPK) of small and large molecules. This tool enables the users to perform the following: (i) PK data visualization, (ii) simulation, (iii) parameter optimization, and (iv) local sensitivity analysis of mPBPK models in a simple and efficient manner. In addition to the theoretical background and implementation of the different tool functionalities, this tutorial includes simulation and sensitivity analysis showcases of small and large molecules with and without target-mediated drug disposition.

Modeling Sex Differences in Anti-inflammatory Effects of Dexamethasone in Arthritic Rats

PURPOSE

Collagen-induced arthritic (CIA) rats are used commonly for preclinical pharmacologic research into rheumatoid arthritis (RA). Dexamethasone (DEX), a potent corticosteroid (CS), remains an important component in combination therapy for RA. Although sex differences in RA and CS pharmacokinetics/pharmacodynamics (PK/PD) have been documented in humans, there has been no such comprehensive evaluation of sex differences in CIA rats.

METHODS

Paw size measurements were obtained for males and females from four groups of animals: healthy controls, non-drug treated arthritic animals, and both 0.225 and 2.25 mg/kg DEX-treated arthritic animals. A turnover model for disease progression, minimal PBPK model for drug concentrations, and inhibitory indirect response model were applied using population PK/PD modeling.

RESULTS

The clearances of DEX were 43% greater in males, but other PK parameters were similar. The temporal profiles of paw swelling exhibited earlier progression, peak edema times, and disease remission in females. DEX suppressed paw edema well in both males and females with similar capacity (I_max) values (=1.0), but DEX potency was less in females with higher IC₅₀ values (0.101 versus 0.015 ng/mL).

CONCLUSIONS

The pharmacology of DEX was well characterized in CIA rats. This study addresses knowledge gaps about sex differences and can be a guide for more mechanistic assessment of sex, drug, and disease differences in RA.

Pharmacological characterisation of CR6086, a potent prostaglandin E2 receptor 4 antagonist, as a new potential disease-modifying anti-rheumatic drug

Background

Prostaglandin E₂ (PGE₂) acts via its EP4 receptor as a cytokine amplifier (e.g., interleukin [IL]-6) and induces the differentiation and expansion of inflammatory T-helper (Th) lymphocytes. These mechanisms play a key role in the onset and progression of rheumatoid arthritis (RA). We present the pharmacological characterisation of CR6086, a novel EP4 receptor antagonist, and provide evidence for its potential as a disease-modifying anti-rheumatic drug (DMARD).

Methods

CR6086 affinity and pharmacodynamics were studied in EP4-expressing HEK293 cells by radioligand binding and cyclic adenosine monophosphate (cAMP) production, respectively. In immune cells, IL-6 and vascular endothelial growth factor (VEGF) expression were analysed by RT-PCR, and IL-23 and IL-17 release were measured by enzyme-linked immunosorbent assay (ELISA). In collagen-induced arthritis (CIA) models, rats or mice were immunised with bovine collagen type II. Drugs were administered orally (etanercept and methotrexate intraperitoneally) starting at disease onset. Arthritis progression was evaluated by oedema, clinical score and histopathology. Anti-collagen II immunoglobulin G antibodies were measured by ELISA.

Results

CR6086 showed selectivity and high affinity for the human EP4 receptor (K_i = 16.6 nM) and functioned as a pure antagonist (half-maximal inhibitory concentration, 22 nM) on PGE₂-stimulated cAMP production. In models of human immune cells in culture, CR6086 reduced key cytokine players of RA (IL-6 and VEGF expression in macrophages, IL-23 release from dendritic cells, IL-17 release from Th17 cells). In the CIA model of RA in rats and mice, CR6086 significantly improved all features of arthritis: severity, histology, inflammation and pain. In rats, CR6086 was better than the selective cyclooxygenase-2 inhibitor rofecoxib and at least as effective as the Janus kinase inhibitor tofacitinib. In mice, CR6086 and the biologic DMARD etanercept were highly effective, whereas the non-steroidal anti-inflammatory drug naproxen was ineffective. Importantly, in a study of CR6086/methotrexate, combined treatment greatly improved the effect of a fully immunosuppressive dose of methotrexate.

Conclusions

CR6086 is a novel, potent EP4 antagonist showing favourable immunomodulatory properties, striking DMARD effects in rodents, and anti-inflammatory activity targeted to immune-mediated inflammatory diseases and distinct from the general effects of cyclooxygenase inhibitors. These results support the clinical development of CR6086, both as a stand-alone DMARD and as a combination therapy with methotrexate. The proof-of-concept trial in patients with RA is ongoing.

Electronic supplementary material

The online version of this article (10.1186/s13075-018-1537-8) contains supplementary material, which is available to authorized users.

Modeling Combined Immunosuppressive and Anti-inflammatory Effects of Dexamethasone and Naproxen in Rats Predicts the Steroid-Sparing Potential of Naproxen

Dexamethasone (DEX), a widely prescribed corticosteroid, has long been the cornerstone of the treatment of inflammation and immunologic dysfunctions in rheumatoid arthritis. Corticosteroids are frequently used in combination with other antirheumatic agents such as nonsteroidal anti-inflammatory drugs (NSAIDs) and disease-modifying antirheumatic drugs to mitigate disease symptoms and minimize unwanted effects. We explored the steroid dose-sparing potential of the NSAID naproxen (NPX) with in vitro and in vivo studies. The single and joint suppressive effects of DEX and NPX on the in vitro mitogen-induced proliferation of T lymphocytes in blood and their anti-inflammatory actions on paw edema were investigated in female and male Lewis rats with collagen-induced arthritis (CIA). As expected, DEX was far more potent than NPX in these systems. Mathematical models incorporating an interaction term ψ were applied to quantitatively assess the nature and intensity of pharmacodynamic interactions between DEX and NPX. Modest synergistic effects of the two drugs were found in suppressing the mitogenic response of T lymphocytes. A pharmacokinetic/pharmacodynamic/disease progression model integrating dual drug inhibition quantitatively described the pharmacokinetics, time-course of single and joint anti-inflammatory effects (paw edema), and sex differences in CIA rats, and indicated additive effects of DEX and NPX. Further model simulations demonstrated the promising steroid-sparing potential of NPX in CIA rats, with the beneficial effects of the combination therapy more likely in males than females.

Effect of Disease-Related Changes in Plasma Albumin on the Pharmacokinetics of Naproxen in Male and Female Arthritic Rats

Naproxen (NPX) is used in the treatment of rheumatoid arthritis (RA) for alleviation of pain and inflammation. In view of the extensive albumin binding of NPX, this study investigates whether chronic inflammation and sex influence the physiologic albumin concentrations, plasma protein binding, and pharmacokinetics (PK) of NPX. The PK of NPX was evaluated in a rat model of RA [collagen-induced arthritis (CIA) in Lewis rats] and in healthy controls. These PK studies included 1) NPX in female and male CIA rats that received 10, 25, or 50 mg/kg NPX i.p.; and 2) NPX in healthy female and male rats after i.p. dosing of NPX at 50 mg/kg. Plasma albumin concentrations were quantified by enzyme-linked immunosorbent assay, and protein binding was assessed using ultrafiltration. The NPX concentrations in plasma and filtrates were determined by liquid chromatography-tandem mass spectrometry (LC-MS/MS). Plasma concentration-time data of NPX were first assessed by noncompartmental analysis (NCA). Nonlinear PK as indicated by dose-dependent NCA clearances and distribution volumes was observed. A two-compartment model with a first-order absorption process incorporating nonlinear protein binding in plasma and tissues jointly described the PK data of all groups. Saturable albumin binding accounts for the nonlinearity of NPX PK in all rats as well as part of the PK differences in arthritic rats. The CIA rats exhibited reduced albumin concentrations, reduced overall protein binding, and reduced clearances of unbound NPX, consistent with expectations during inflammation. The net effect of chronic inflammation was an elevation of the Cmax and area under the plasma concentration-time curve (AUC) of unbound drug.