Modeling receptor/gene-mediated effects of corticosteroids on hepatic tyrosine aminotransferase dynamics in rats: dual regulation by endogenous and exogenous corticosteroids

The Physiological and Pharmacological Significance of the Circadian Timing of the HPA Axis: A Mathematical Modeling Approach

As a potent endogenous regulator of homeostasis, the circadian time-keeping system synchronizes internal physiology to periodic changes in the external environment to enhance survival. Adapting endogenous rhythms to the external time is accomplished hierarchically with the central pacemaker located in the suprachiasmatic nucleus (SCN) signaling the hypothalamus-pituitary-adrenal (HPA) axis to release hormones, notably cortisol, which help maintain the body's circadian rhythm. Given the essential role of HPA-releasing hormones in regulating physiological functions, including immune response, cell cycle, and energy metabolism, their daily variation is critical for the proper function of the circadian timing system. In this review, we focus on cortisol and key fundamental properties of the HPA axis and highlight their importance in controlling circadian dynamics. We demonstrate how systems-driven, mathematical modeling of the HPA axis complements experimental findings, enhances our understanding of complex physiological systems, helps predict potential mechanisms of action, and elucidates the consequences of circadian disruption. Finally, we outline the implications of circadian regulation in the context of personalized chronotherapy. Focusing on the chrono-pharmacology of synthetic glucocorticoids, we review the challenges and opportunities associated with moving toward personalized therapies that capitalize on circadian rhythms.

Circadian rhythms: influence on physiology, pharmacology, and therapeutic interventions

Circadian rhythms are ubiquitous phenomena that recur daily in a self-sustaining, entrainable, and oscillatory manner, and orchestrate a wide range of molecular, physiological, and behavioral processes. Circadian clocks are comprised of a hierarchical network of central and peripheral clocks that generate, sustain, and synchronize the circadian rhythms. The functioning of the peripheral clock is regulated by signals from autonomic innervation (from the central clock), endocrine networks, feeding, and other external cues. The critical role played by circadian rhythms in maintaining both systemic and tissue-level homeostasis is well established, and disruption of the rhythm has direct consequence for human health, disorders, and diseases. Circadian oscillations in both pharmacokinetics and pharmacodynamic processes are known to affect efficacy and toxicity of several therapeutic agents. A variety of modeling approaches ranging from empirical to more complex systems modeling approaches have been applied to characterize circadian biology and its influence on drug actions, optimize time of dosing, and identify opportunities for pharmacological modulation of the clock mechanisms and their downstream effects. In this review, we summarize current understanding of circadian rhythms and its influence on physiology, pharmacology, and therapeutic interventions, and discuss the role of chronopharmacometrics in gaining new insights into circadian rhythms and its applications in chronopharmacology.

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.

Seeking Nonspecific Binding: Assessing the Reliability of Tissue Dilutions for Calculating Fraction Unbound

It has become commonplace (270+ article citations to date) to measure the fraction unbound (FrUn) of drugs in tissue homogenates and diluted plasma and then use a Correction Factor Equation (CFE) to extrapolate to the undiluted state. The CFE is based on assumptions of nonspecific binding with experimental use of very low drug concentrations. There are several possible determinants of apparent nonspecific binding as measured by methods such as equilibrium dialysis: true macromolecule binding and lipid partitioning along with receptor, enzyme, and transporter interactions. Theoretical calculations based on nonlinear protein binding indicate that the CFE will be most reliable to obtain FrUn when added drug concentration is small, binding constants are weak, protein concentrations are relatively high, and tissue dilution is minimal. When lipid partitioning is the sole factor determining apparent tissue binding, the CFE should be perfectly accurate. Use of very low drug concentrations, however, makes it more likely that specific binding to receptors and other targets may occur, and thus FrUn may reflect some binding to such components. Inclusion of trapped blood can clearly cause minor to marked discrepancies from purely tissue binding alone, which can be corrected. Furthermore, assessment of the occurrence of ionization/pH shifts, drug instability, and tissue metabolism may be necessary. Caution is needed in the use and interpretation of results from tissue dilution studies and other assessments of nonspecific binding, particularly for very strongly bound drugs with very small FrUn values and in tissues with metabolic enzymes, receptors, and trapped blood. SIGNIFICANCE STATEMENT: The use of tissue, plasma, and cell preparations to help obtain fraction unbound and tissue-to-plasma partition coefficients in pharmacokinetics has grown commonplace, especially for brain. This report examines theoretical, physiological, and experimental issues that need consideration before trusting such measurements and calculations.

Modeling Corticosteroid Pharmacokinetics and Pharmacodynamics, Part III: Estrous Cycle and Estrogen Receptor-Dependent Antagonism of Glucocorticoid-Induced Leucine Zipper (GILZ) Enhancement by Corticosteroids

Our previous report examined the pharmacokinetics (PK) of methylprednisolone (MPL) and adrenal suppression after a 50 mg/kg IM bolus in male and female rats, and we described in detail the development of a minimal physiologically based pharmacokinetic/pharmacodynamic (mPBPK/PD) model. In continuation of such assessments, we investigated sex differences in genomic MPL responses (PD). Message expression of the glucocorticoid-induced leucine zipper (GILZ) was chosen as a multitissue biomarker of glucocorticoid receptor (GR)-mediated drug response. Potential time-dependent interplay between sex hormone and glucocorticoid signaling in vivo was assessed by comparing the enhancement of GILZ by MPL in the uterus [high estrogen receptor (ER) density] and in liver (lower ER density) from male and female rats dosed within the proestrus (high estradiol/progesterone) and estrus (low estradiol/progesterone) phases of the rodent estrous cycle. An expanded-systems PD model of MPL considering circadian rhythms, multireceptor (ER and GR) control, and estrous variations delineated the determinants controlling receptor/gene-mediated steroid responses. Hepatic GILZ response was ∼3-fold greater in females, regardless of estrous stage, compared with males, driven predominantly by increased MPL exposure in females and a negligible influence of estrogen interaction. In contrast, GILZ response in the uterus during proestrus in females was 60% of that observed in estrus-phased females, despite no PK or receptor differences, providing in vivo support to the hypothesis of estrogen-mediated antagonism of glucocorticoid signaling. The developed model offers a mechanistic platform to assess the determinants of sex and tissue specificity in corticosteroid actions and, in turn, reveals a unique PD drug-hormone interaction occurring in vivo. SIGNIFICANCE STATEMENT: Mechanisms relating to sex-based pharmacodynamic variability in genomic responses to corticosteroids have been unclear. Using combined experimental and systems pharmacology modeling approaches, sex differences in both pharmacokinetic and pharmacodynamic mechanisms controlling the enhancement of a sensitive corticosteroid-regulated biomarker, the glucocorticoid-induced leucine zipper (GILZ), were clarified in vivo. The multiscale minimal physiologically based pharmacokinetics/pharmacodynamic model successfully captured the experimental observations and quantitatively discerned the roles of the rodent estrous cycle (hormonal variation) and tissue specificity in mediating the antagonistic coregulation of GILZ gene synthesis. These findings collectively support the hypothesis that estrogens antagonize pharmacodynamic signaling of genomic corticosteroid actions in vivo in a time- and estrogen receptor-dependent manner.

Modeling Corticosteroid Pharmacokinetics and Pharmacodynamics, Part II: Sex Differences in Methylprednisolone Pharmacokinetics and Corticosterone Suppression

Methylprednisolone (MPL), a corticosteroid of intermediate potency, remains an important immunomodulatory agent for autoimmune diseases. Although sex differences in corticosteroid pharmacokinetics/pharmacodynamics (PK/PD) have been documented in humans, comprehensive preclinical assessments of such differences have not been conducted. Limited in vitro evidence indicates possible sex differences in corticosteroid PK and PD. Therefore, it is hypothesized that comparative PK/PD assessments of MPL disposition and selected PD actions in both sexes will provide insights into factors controlling sex differences in steroid responses. This report focused on the plasma and tissue pharmacokinetics of MPL and its adrenal suppressive effects. Because time-dependent (estrous) regulation of sex hormones in females can influence drug responses, female rats were studied in the proestrus (high estradiol/progesterone) and estrus (low estradiol/progesterone) phases of the reproductive cycle. Cohorts of male and female rats were given a 50 mg/kg bolus dose of MPL intramuscularly. Plasma and liver concentrations of MPL as well as plasma corticosterone concentrations were assayed using high-performance liquid chromatography. An enhanced minimal physiologically-based PK/PD model was developed to characterize MPL kinetics and corticosterone dynamics. The clearance of MPL was ∼3-fold higher in males compared with females, regardless of estrous phase, likely attributable to sex-specific hepatic metabolism in males. Strong inhibitory effects on adrenal suppression were observed in all animals. These temporal steroid profiles in plasma and tissues will be used to drive receptor/gene-mediated PD effects of MPL in both sexes, as described in a companion article (Part III). SIGNIFICANCE STATEMENT: Sex is a relevant factor influencing the pharmacokinetics (PK) and pharmacodynamics (PD) of drugs. Few preclinical PK/PD studies, however, include sex as a variable. Sex differences in the PK and adrenal suppressive effects of the synthetic corticosteroid, methylprednisolone, were assessed in male and female rats as a function of the 4-day rodent reproductive cycle. Drug exposure was 3-fold higher in females, regardless of estrous stage, compared with males. An extended minimal physiologically-based PK/PD model utilizing in vitro and in vivo measurements was developed and applied. These studies provide a framework to account for sex-dependent variability in drug and endogenous agonist (corticosterone) exposures, serving as a prelude to more intricate assessments of sex-related variability in receptor/gene-mediated PD corticosteroid actions.

Modeling Corticosteroid Pharmacogenomics and Proteomics in Rat Liver

Corticosteroids (CS) regulate the expression of numerous genes at the mRNA and protein levels. The time course of CS pharmacogenomics and proteomics were examined in livers obtained from adrenalectomized rats given a 50-mg/kg bolus dose of methylprednisolone. Microarrays and mass spectrometry-based proteomics were employed to quantify hepatic transcript and protein dynamics. One-hundred, sixty-three differentially expressed mRNA and their corresponding proteins (163 genes) were clustered into two dominant groups. The temporal profiles of most proteins were delayed compared with their mRNA, attributable to synthesis delays and slower degradation kinetics. On the basis of our fifth-generation model of CS, mathematical models were developed to simultaneously describe the emergent time patterns for an array of steroid-responsive mRNA and proteins. The majority of genes showed time-dependent increases in mRNA and protein expression before returning to baseline. A model assuming direct, steroid-mediated stimulation of mRNA synthesis was applied. Some mRNAs and their proteins displayed down-regulation following CS. A model assuming receptor-mediated inhibition of mRNA synthesis was used. More complex patterns were observed for other genes (e.g., biphasic behaviors and opposite directionality in mRNA and protein). Models assuming either stimulation or inhibition of mRNA synthesis coupled with dual secondarily induced regulatory mechanisms affecting mRNA or protein turnover were derived. These findings indicate that CS-regulated gene expression manifested at the mRNA and protein levels are controlled via mechanisms affecting key turnover processes. Our quantitative models of CS pharmacogenomics were expanded from mRNA to proteins and provide extended hypotheses for understanding the direct, secondary, and downstream mechanisms of CS actions.

Receptor/gene/protein-mediated signaling connects methylprednisolone exposure to metabolic and immune-related pharmacodynamic actions in liver

A multiscale pharmacodynamic model was developed to characterize the receptor-mediated, transcriptomic, and proteomic determinants of corticosteroid (CS) effects on clinically relevant hepatic processes following a single dose of methylprednisolone (MPL) given to adrenalectomized (ADX) rats. The enhancement of tyrosine aminotransferase (TAT) mRNA, protein, and enzyme activity were simultaneously described. Mechanisms related to the effects of MPL on glucose homeostasis, including the regulation of CCAAT-enhancer binding protein-beta (C/EBPβ) and phosphoenolpyruvate carboxykinase (PEPCK) as well as insulin dynamics were evaluated. The MPL-induced suppression of circulating lymphocytes was modeled by coupling its effect on cell trafficking with pharmacogenomic effects on cell apoptosis via the hepatic (STAT3-regulated) acute phase response. Transcriptomic and proteomic time-course profiles measured in steroid-treated rat liver were utilized to model the dynamics of mechanistically relevant gene products, which were linked to associated systemic end-points. While time-courses of TAT mRNA, protein, and activity were well described by transcription-mediated changes, additional post-transcriptional processes were included to explain the lack of correlation between PEPCK mRNA and protein. The immune response model quantitatively discerned the relative roles of cell trafficking versus gene-mediated lymphocyte apoptosis by MPL. This systems pharmacodynamic model provides insights into the contributions of selected molecular events occurring in liver and explores mechanistic hypotheses for the multi-factorial control of clinically relevant pharmacodynamic outcomes.

On the analysis of complex biological supply chains: From Process Systems Engineering to Quantitative Systems Pharmacology

The use of models in biology has become particularly relevant as it enables investigators to develop a mechanistic framework for understanding the operating principles of living systems as well as in quantitatively predicting their response to both pathological perturbations and pharmacological interventions. This application has resulted in a synergistic convergence of systems biology and pharmacokinetic-pharmacodynamic modeling techniques that has led to the emergence of quantitative systems pharmacology (QSP). In this review, we discuss how the foundational principles of chemical process systems engineering inform the progressive development of more physiologically-based systems biology models.

Highlighting Vitamin D Receptor-Targeted Activities of 1α,25-Dihydroxyvitamin D3 in Mice via Physiologically Based Pharmacokinetic-Pharmacodynamic Modeling

We expanded our published physiologically based pharmacokinetic model (PBPK) on 1α,25-dihydroxyvitamin D₃ [1,25(OH)₂D₃], ligand of the vitamin D receptor (VDR), to appraise VDR-mediated pharmacodynamics in mice. Since 1,25(OH)₂D₃ kinetics was best described by a segregated-flow intestinal model (SFM) that described a low/partial intestinal (blood/plasma) flow to enterocytes, with feedback regulation of its synthesis (Cyp27b1) and degradation (Cyp24a1) enzymes, this PBPK(SFM) model was expanded to describe the VDR-mediated changes (altered/basal mRNA expression) of target genes/responses with the indirect response model. We examined data on 1) renal Trpv5 (transient receptor potential cation channel, subfamily V member 5) and Trpv6 and intestinal Trpv6 (calcium channels) for calcium absorption; 2) liver 3-hydroxy-3-methyl-glutaryl-coenzyme A reductase (Hmgcr) and cytochrome 7α-hydroxylase (Cyp7a1) for cholesterol synthesis and degradation, respectively; and 3) renal and brain Mdr1 (multidrug-resistance protein that encodes the P-glycoprotein) for digoxin disposition after repetitive intraperitoneal doses of 120 pmol 1,25(OH)₂D₃ Fitting, performed with modeling software, yielded reasonable prediction of a dominant role of intestinal Trpv6 in calcium absorption, circadian rhythm that is characterized by simple cosine models for Hmgcr and Cyp7a1 on liver cholesterol, and brain and renal Mdr1 on tissue efflux of digoxin. Fitted parameters on the E_max, EC₅₀, and turnover rate constants of VDR-target genes [zero-order production (k_in) and first-order degradation (k_out) rate constants] showed low coefficients of variation and acceptable median prediction errors (4.5%-40.6%). Sensitivity analyses showed that the E_max and EC₅₀ values are key parameters that could influence the pharmacodynamic responses. In conclusion, the PBPK(SFM)-pharmacodynamic model successfully characterized VDR gene activation and serves as a useful tool to predict the therapeutic effects of 1,25(OH)₂D₃.

Mechanistic Multi-Tissue Modeling of Glucocorticoid-Induced Leucine Zipper Regulation: Integrating Circadian Gene Expression with Receptor-Mediated Corticosteroid Pharmacodynamics

The glucocorticoid-induced leucine zipper (GILZ) is an important mediator of anti-inflammatory corticosteroid action. The pharmacokinetic/pharmacodynamic/pharmacogenomic effects of acute and chronic methylprednisolone (MPL) dosing on the tissue-specific dynamics of GILZ expression were examined in rats. A mechanism-based model was developed to investigate and integrate the role of MPL and circadian rhythms on the transcriptional enhancement of GILZ in multiple tissues. Animals received a single 50-mg/kg intramuscular bolus or a 7-day 0.3-mg/kg/h subcutaneous infusion of MPL and were euthanized at several time points. An additional group of rats were euthanized at several times and served as 24-hour light/dark (circadian) controls. Plasma MPL and corticosterone concentrations were measured by high-performance liquid chromatography. The expression of GILZ and glucocorticoid receptor (GR) mRNA was quantified in tissues using quantitative real-time reverse-transcription polymerase chain reaction. The pharmacokinetics of MPL were described using a two-compartment model. Mild-to-robust circadian oscillations in GR and GILZ mRNA expression were characterized in muscle, lung, and adipose tissues and modeled using Fourier harmonic functions. Acute MPL dosing caused significant down-regulation (40%-80%) in GR mRNA and enhancement of GILZ mRNA expression (500%-1080%) in the tissues examined. While GILZ returned to its rhythmic baseline following acute dosing, a new steady-state was observed upon enhancement by chronic dosing. The model captured the complex dynamics in all tissues for both dosing regimens. The model quantitatively integrates physiologic mechanisms, such as circadian processes and GR tolerance phenomena, which control the tissue-specific regulation of GILZ by corticosteroids. These studies characterize GILZ as a pharmacodynamic marker of corticosteroid actions in several tissues.

Functional proteomic analysis of corticosteroid pharmacodynamics in rat liver: Relationship to hepatic stress, signaling, energy regulation, and drug metabolism

Corticosteroids (CS) are anti-inflammatory agents that cause extensive pharmacogenomic and proteomic changes in multiple tissues. An understanding of the proteome-wide effects of CS in liver and its relationships to altered hepatic and systemic physiology remains incomplete. Here, we report the application of a functional pharmacoproteomic approach to gain integrated insight into the complex nature of CS responses in liver in vivo. An in-depth functional analysis was performed using rich pharmacodynamic (temporal-based) proteomic data measured over 66h in rat liver following a single dose of methylprednisolone (MPL). Data mining identified 451 differentially regulated proteins. These proteins were analyzed on the basis of temporal regulation, cellular localization, and literature-mined functional information. Of the 451 proteins, 378 were clustered into six functional groups based on major clinically-relevant effects of CS in liver. MPL-responsive proteins were highly localized in the mitochondria (20%) and cytosol (24%). Interestingly, several proteins were related to hepatic stress and signaling processes, which appear to be involved in secondary signaling cascades and in protecting the liver from CS-induced oxidative damage. Consistent with known adverse metabolic effects of CS, several rate-controlling enzymes involved in amino acid metabolism, gluconeogenesis, and fatty-acid metabolism were altered by MPL. In addition, proteins involved in the metabolism of endogenous compounds, xenobiotics, and therapeutic drugs including cytochrome P450 and Phase-II enzymes were differentially regulated. Proteins related to the inflammatory acute-phase response were up-regulated in response to MPL. Functionally-similar proteins showed large diversity in their temporal profiles, indicating complex mechanisms of regulation by CS.

SIGNIFICANCE

Clinical use of corticosteroid (CS) therapy is frequent and chronic. However, current knowledge on the proteome-level effects of CS in liver and other tissues is sparse. While transcriptomic regulation following methylprednisolone (MPL) dosing has been temporally examined in rat liver, proteomic assessments are needed to better characterize the tissue-specific functional aspects of MPL actions. This study describes a functional pharmacoproteomic analysis of dynamic changes in MPL-regulated proteins in liver and provides biological insight into how steroid-induced perturbations on a molecular level may relate to both adverse and therapeutic responses presented clinically.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Quantitative tissue-specific dynamics of in vivo GILZ mRNA expression and regulation by endogenous and exogenous glucocorticoids

Glucocorticoids (GC) are steroid hormones, which regulate metabolism and immune function. Synthetic GCs, or corticosteroids (CS), have appreciable clinical utility via their ability to suppress inflammation in immune-mediated diseases like asthma and rheumatoid arthritis. Recent work has provided insight to novel GC-induced genes that mediate their anti-inflammatory effects, including glucocorticoid-induced leucine zipper (GILZ). Since GILZ comprises an important part of GC action, its regulation by both drug and hormone will influence CS therapy. In addition, GILZ expression is often employed as a biomarker of GC action, which requires judicious selection of sampling time. Understanding the in vivo regulation of GILZ mRNA expression over time will provide insight into both the physiological regulation of GILZ by endogenous GC and the dynamics of its enhancement by CS. A highly quantitative qRT-PCR assay was developed for measuring GILZ mRNA expression in tissues obtained from normal and CS-treated rats. This assay was applied to measure GILZ mRNA expression in eight tissues; to determine its endogenous regulation over time; and to characterize its dynamics in adipose tissue, muscle, and liver following treatment with CS. We demonstrate that GILZ mRNA is expressed in several tissues. GILZ mRNA expression in adipose tissue displayed a robust circadian rhythm that was entrained with the circadian oscillation of endogenous corticosterone; and is strongly enhanced by acute and chronic dosing. Single dosing also enhanced GILZ mRNA in muscle and liver, but the dynamics varied. In conclusion, GILZ is widely expressed in the rat and highly regulated by endogenous and exogenous GCs.

Evidence for a glucocorticoid receptor beta splice variant in the rat and its physiological regulation in liver

Glucocorticoids are important regulators of metabolism and immune function. Synthetic glucocorticoids are extensively used for immunosuppression/anti-inflammatory therapy. Since the glucocorticoid receptor (GR) is central to most hormone effects; its in vivo regulation will influence hormone/drug action. An alternative splice variant, GRβ, is present in humans and may function as a dominant negative regulator of GR transcriptional activity. Recently, a similar splice variant was reported in mouse, although the mechanism of alternative splicing differs from that in humans. We present evidence that a splice variant of GR with an alternative C-terminus also occurs in the rat by a mechanism of intron inclusion. A highly quantitative qRT-PCR assay for the simultaneous measurement of both splice variants in a single sample was developed in order to accurately measure their regulation. We used this assay to assess the tissue specific expression of both mRNAs, and demonstrate that GRα is predominant in all tissues. In addition, the regulation of both GRα and GRβ mRNA by various physiological factors in rat liver was assessed. GRα showed a robust circadian rhythm, which was entrained with the circadian oscillation of the endogenous hormone. Time series experiments showed that both corticosteroids and LPS but not insulin dosing resulted in the transient down-regulation of GRα mRNA. LPS treatment also resulted in down-regulation of GRβ expression. A modest up-regulation in GRβ expression was observed only in animals having chronically elevated plasma insulin concentrations. However the expression of GRβ was significantly lower than that of GRα in all cases.

Mechanism-based pharmacodynamic modeling

Pharmacodynamic modeling is based on a quantitative integration of pharmacokinetics, pharmacological systems, and (patho-) physiological processes for understanding the intensity and time-course of drug effects on the body. Application of such models to the analysis of meaningful experimental data allows for the quantification and prediction of drug-system interactions for both therapeutic and adverse drug responses. In this chapter, commonly used mechanistic pharmacodynamic models are presented with respect to their important features, operable equations, and signature profiles. In addition, literature examples showcasing the utility of these models to adverse drug events are highlighted. Common model types that are covered include simple direct effects, biophase distribution, indirect effects, signal transduction, and irreversible effects.

Mechanism-based pharmacokinetic-pharmacodynamic modeling of the estrogen-like effect of ginsenoside Rb1 on neural 5-HT in ovariectomized mice

We sought to develop a mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) model to characterize the effects of ginsenoside Rb1 (Rb1) and estradiol (E(2)) on neural 5-hydroxytryptamine (5-HT) concentration in ovariectomized mice. PK data of Rb1 and E(2) were obtained in plasma and brain. Brain levels of 5-HT, tryptophan (TRP), 5-hydroxytryptophan (5-HTP), and 5-hydroxyindoleacetic acid (5-HIAA) were determined after a single intravenous injection of Rb1 (20mg/kg) and E(2) (0.2mg/kg) in ovariectomized mice. The activities of tryptophan hydroxylase (TPH), aromatic amino acid decarboxylase (AAAD), and monoamine oxidase (MAO) were also evaluated. Rb1 and E(2) elevated neural 5-HT levels via TPH activation and MAO inhibition, respectively. Effects were well described by the mechanism-based PK-PD model. The net effect of increased 5-HT induced by MAO inhibition is greater than TPH activation. The increased brain levels of 5-HT induced by Rb1 and E(2) were well described by the present PK-PD model, suggesting the use and further development of this mechanism-based model for the effects of ginsenoside on brain 5-HT levels.

Modeling autonomic regulation of cardiac function and heart rate variability in human endotoxemia

Heart rate variability (HRV), the quantification of beat-to-beat variability, has been studied as a potential prognostic marker in inflammatory diseases such as sepsis. HRV normally reflects significant levels of variability in homeostasis, which can be lost under stress. Much effort has been placed in interpreting HRV from the perspective of quantitatively understanding how stressors alter HRV dynamics, but the molecular and cellular mechanisms that give rise to both homeostatic HRV and changes in HRV have received less focus. Here, we develop a mathematical model of human endotoxemia that incorporates the oscillatory signals giving rise to HRV and their signal transduction to the heart. Connections between processes at the cellular, molecular, and neural levels are quantitatively linked to HRV. Rhythmic signals representing autonomic oscillations and circadian rhythms converge to modulate the pattern of heartbeats, and the effects of these oscillators are diminished in the acute endotoxemia response. Based on the semimechanistic model developed herein, homeostatic and acute stress responses of HRV are studied in terms of these oscillatory signals. Understanding the loss of HRV in endotoxemia serves as a step toward understanding changes in HRV observed clinically through translational applications of systems biology based on the relationship between biological processes and clinical outcomes.

Light-dark oscillations in the lung transcriptome: implications for lung homeostasis, repair, metabolism, disease, and drug action

Diurnal-nocturnal, or circadian-like, rhythms are 24-h variations in biological processes, evolved for the efficient functioning of living organisms. Such oscillations and their regulation in many peripheral tissues are still unclear. In this study, we used Affymetrix gene chips in a rich time-series experiment involving 54 animals killed at 18 time points within the 24-h cycle to examine light-dark cycle patterns of gene expression in rat lungs. Data mining identified 646 genes (represented by 1,006 probe sets) showing robust oscillations in expression in lung that were parsed into 8 distinct temporal clusters. Surprisingly, more than two-thirds of the probe sets showing cyclic expression peaked during the animal's light/inactive period. Six core clock genes and nine clock-related genes showed rhythmic oscillations in their expression in lung. Many of the genes that peaked during the inactive period included genes related to extracellular matrix, cytoskeleton, and protein processing and trafficking, which appear to be mainly involved in the repair and remodeling of the organ. Genes coding for growth factor ligands and their receptors, which play important roles in maintaining normal lung function, also showed rhythmic expression. In addition, genes involved in the metabolism and transport of endogenous compounds, xenobiotics, and therapeutic drugs, along with genes that are biomarkers or potential therapeutic targets for many lung diseases, also exhibited 24-h cyclic oscillations, suggesting an important role for such rhythms in regulating various aspects of the physiology and pathophysiology of lung.

Blood content of tyrosine is an index of glucocorticoid action on metabolism

Glucocorticoid hormones directly or indirectly control virtually all metabolic and physiological processes. Glucocorticoids are also shown to act on a multitude of genes, enzyme systems, and proinflammatory factors, but for these hormones there is no representative index of action on metabolism similar to glucose content in blood for insulin. The absence of such an index prevents the assessment of tissue provision with these hormones under various conditions and seems to be an essential cause of complications associated with the clinical use of glucocorticoid preparations. Considering specific features of tyrosine metabolism and data obtained experimentally and on a clinical model (adrenalectomy in rats and substitution therapy in endocrine disease), blood content of this amino acid seems promising as such an index. Based on comparing results of glucocorticoid treatment in patients with systemic lupus erythematosus with changes in their blood tyrosine contents, the pharmacological effect of glucocorticoid preparations is suggested to be mainly due to compensating a relative shortage of these hormones.

Circadian rhythms in gene expression: Relationship to physiology, disease, drug disposition and drug action

Circadian rhythms (24h cycles) are observed in virtually all aspects of mammalian function from expression of genes to complex physiological processes. The master clock is present in the suprachiasmatic nucleus (SCN) in the anterior part of the hypothalamus and controls peripheral clocks present in other parts of the body. Components of this core clock mechanism regulate the circadian rhythms in genome-wide mRNA expression, which in turn regulate various biological processes. Disruption of circadian rhythms can be either the cause or the effect of various disorders including metabolic syndrome, inflammatory diseases and cancer. Furthermore, circadian rhythms in gene expression regulate both the action and disposition of various drugs and affect therapeutic efficacy and toxicity based on dosing time. Understanding the regulation of circadian rhythms in gene expression plays an important role in both optimizing the dosing time for existing drugs and in the development of new therapeutics targeting the molecular clock.

Modeling corticosteroid effects in a rat model of rheumatoid arthritis I: mechanistic disease progression model for the time course of collagen-induced arthritis in Lewis rats

A mechanism-based model was developed to describe the time course of arthritis progression in the rat. Arthritis was induced in male Lewis rats with type II porcine collagen into the base of the tail. Disease progression was monitored by paw swelling, bone mineral density (BMD), body weights, plasma corticosterone (CST) concentrations, and tumor necrosis factor (TNF)-alpha, interleukin (IL)-1beta, IL-6, and glucocorticoid receptor (GR) mRNA expression in paw tissue. Bone mineral density was determined by PIXImus II dual energy X-ray densitometry. Plasma CST was assayed by high-performance liquid chromatography. Cytokine and GR mRNA were determined by quantitative real-time polymerase chain reaction. Disease progression models were constructed from transduction and indirect response models and applied using S-ADAPT software. A delay in the onset of increased paw TNF-alpha and IL-6 mRNA concentrations was successfully characterized by simple transduction. This rise was closely followed by an up-regulation of GR mRNA and CST concentrations. Paw swelling and body weight responses peaked approximately 21 days after induction, whereas bone mineral density changes were greatest at 23 days after induction. After peak response, the time course in IL-1beta, IL-6 mRNA, and paw edema slowly declined toward a disease steady state. Model parameters indicate TNF-alpha and IL-1beta mRNA most significantly induce paw edema, whereas IL-6 mRNA exerted the most influence on BMD. The model for bone mineral density captures rates of turnover of cancellous and cortical bone and the fraction of each in the different regions analyzed. This small systems model integrates and quantitates multiple factors contributing to arthritis in rats.

Modeling corticosteroid effects in a rat model of rheumatoid arthritis II: mechanistic pharmacodynamic model for dexamethasone effects in Lewis rats with collagen-induced arthritis

A mechanism-based model for pharmacodynamic effects of dexamethasone (DEX) was incorporated into our model for arthritis disease progression in the rat to aid in identification of the primary factors responsible for edema and bone loss. Collagen-induced arthritis was produced in male Lewis rats after injection of type II porcine collagen. DEX was given subcutaneously in single doses of 0.225 or 2.25 mg/kg or 7-day multiple doses of 0.045 or 0.225 mg/kg at 21 days postdisease induction. Effects on disease progression were measured by paw swelling, bone mineral density (BMD), body weights, plasma corticosterone (CST), and tumor necrosis factor (TNF)-alpha, interleukin (IL)-1beta, IL-6, and glucocorticoid receptor (GR) mRNA expression in paw tissue. Lumbar and femur BMD was determined by PIXImus II dual-energy X-ray absorptiometry. Plasma CST was assayed by high-performance liquid chromatography. Cytokine and GR mRNA were assayed by quantitative real-time polymerase chain reaction. Indirect response models, drug interaction models, transduction processes, and the fifth-generation model of corticosteroid dynamics were integrated and applied using S-ADAPT software to describe how dexamethasone binding to GR can regulate diverse processes. Cytokine mRNA, GR mRNA, plasma CST, and paw edema were suppressed after DEX administration. TNF-alpha mRNA expression and BMD seemed to increase immediately after dosing but were ultimately reduced. Model parameters indicated that IL-6 and IL-1beta were most sensitive to inhibition by DEX. TNF-alpha seemed to primarily influence edema, whereas IL-6 contributed the most to bone loss. Lower doses of corticosteroids may be sufficient to suppress the cytokines most relevant to bone erosion.