Modeling Corticosteroid Pharmacogenomics and Proteomics in Rat Liver

Mechanistic Pharmacokinetics and Pharmacodynamics of GalNAc-siRNA: Translational Model Involving Competitive Receptor-Mediated Disposition and RISC-Dependent Gene Silencing Applied to Givosiran

Triantennary N-acetyl-D galactosamine (GalNAc)3-conjugated small interfering RNA (siRNA) have majorly advanced the development of RNA-based therapeutics. Chemically stabilized GalNAc-siRNAs exhibit extensive albeit capacity-limited (nonlinear) distribution into hepatocytes with additional complexities in intracellular liver disposition and pharmacology. A mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) model of GalNAc-siRNA was developed to i) quantitate ASGPR-mediated disposition and downstream RNA-induced silencing complex (RISC)-dependent pharmacology following intravenous (IV) and subcutaneous (SC) dosing, ii) assess the kinetics of formed active metabolite, iii) leverage, as an example, published experimental data for givosiran, and iv) demonstrate PK translation across two preclinical species (rat and monkey) with subsequent prediction of human plasma PK. The structural model is based on competition between parent and formed active metabolite for occupancy and uptake via ASGPR into hepatocytes, intracellular sequestration and degradation, and downstream engagement of RNA-induced silencing complex (RISC) governing target mRNA degradation. The model jointly and accurately captured available concentration-time profiles of givosiran and/or AS(N-1)3' givosiran in rat and/or monkey plasma, liver, and/or kidney following givosiran administered both IV and SC. RISC-dependent gene silencing of ALAS1 mRNA was well-characterized. The model estimated an in vivo affinity (KD) value of 27.7 nM for GalNAc-ASGPR and weight-based allometric exponents of -0.27 and -0.24 for SC absorption and intracellular (endolysosomal) degradation rate constants. The model well-predicted reported givosiran plasma PK profiles in humans. PK simulations revealed net-shifts in liver-to-kidney distribution ratios with increasing IV and SC dose. Importantly, decreases in the relative liver uptake efficiency were demonstrated following IV and, to a lesser extent, following SC dosing explained by differential ASGPR occupancy profiles over time.

Acute serum protein and cytokine response of single dose of prednisone in adult volunteers

Pharmacological glucocorticoids are the most prescribed anti-inflammatory medications, and are chemical variants of cortisol, the circadian and stress hormone. Both endogenous and pharmacological glucocorticoids bind the glucocorticoid receptor (NR3C1) with high affinity, and both then bind downstream gene promoter elements (GRE) to drive positive gene transcription of many proteins. Glucocorticoid/GR complexes also bind distinct negative gene promoter elements (nGRE) to inhibit expression of genes involved in NF-κB innate immunity signaling. We sought to define the acute response of a single dose of prednisone (0.2 mg/kg) in young adult volunteers, with blood samples taken at baseline, 2, 3, 4 and 6 h post-oral dose. To control for circadian morning cortisol hitting the same molecular pathways, a day of blood draws was done without oral prednisone (same time of day), one day prior to drug day. Serum samples were processed for steroid hormone profiles (mass spectrometry; 9 steroidal hormones), proteomics (SOMAscan aptamer panels, 1,305 proteins), and inflammatory markers (Meso Scale Discovery; 10 pro-inflammatory cytokines). The pharmacological effect of the prednisone dose was shown by significant declines of adrenal steroids by 3 h after dosing. IL-10 showed drug-related increase to 4 hrs, then decrease to 6 hrs. IL-8 showed drug-related decrease in serum by 4 h, consistent with direct negative action of GR/ligand on IL-8 gene promoter. Proteomics data showed beta-2 microglobulin, TNFSF15, TSH, CST3, NBL1 to show time-related decreases with prednisone, while CXCL13 showed increases, although these require validation. In summary, a single low dose of prednisone leads to broad suppression of the adrenal axis within 3 h, and down-regulation of inflammatory serum proteins by 6 h.

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.

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

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

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.

Modeling Pathway Dynamics of the Skeletal Muscle Response to Intravenous Methylprednisolone (MPL) Administration in Rats: Dosing and Tissue Effects

A model-based approach for the assessment of pathway dynamics is explored to characterize metabolic and signaling pathway activity changes characteristic of the dosing-dependent differences in response to methylprednisolone in muscle. To consistently compare dosing-induced changes we extend the principles of pharmacokinetics and pharmacodynamics and introduce a novel representation of pathway-level dynamic models of activity regulation. We hypothesize the emergence of dosing-dependent regulatory interactions is critical to understanding the mechanistic implications of MPL dosing in muscle. Our results indicate that key pathways, including amino acid and lipid metabolism, signal transduction, endocrine regulation, regulation of cellular functions including growth, death, motility, transport, protein degradation, and catabolism are dependent on dosing, exhibiting diverse dynamics depending on whether the drug is administered acutely of continuously. Therefore, the dynamics of drug presentation offer the possibility for the emergence of dosing-dependent models of regulation. Finally, we compared acute and chronic MPL response in muscle with liver. The comparison revealed systematic response differences between the two tissues, notably that muscle appears more prone to adapt to MPL.