A proteome-level time-series study
of drug effects (i.e., pharmacodynamics)
is critical for understanding mechanisms of action and systems pharmacology,
but is challenging, because of the requirement of a proteomics method
for reliable quantification of many biological samples. Here, we describe a highly reproducible strategy, enabling a global,
large-scale investigation of the expression dynamics of corticosteroid-regulated
proteins in livers from adrenalectomized rats over 11 time points
after drug dosing (0.5–66 h, N = 5/point).
The analytical advances include (i) exhaustive tissue extraction with
a Polytron/sonication procedure in a detergent cocktail buffer, and
a cleanup/digestion procedure providing very consistent protein yields
(relative standard deviation (RSD%) of 2.7%–6.4%) and peptide
recoveries (4.1–9.0%) across the 60 animals; (ii) an ultrahigh-pressure
nano-LC setup with substantially improved temperature stabilization,
pump-noise suppression, and programmed interface cleaning, enabling
excellent reproducibility for continuous analyses of numerous samples;
(iii) separation on a 100-cm-long column (2-μm particles) with
high reproducibility for days to enable both in-depth profiling and
accurate peptide ion-current match; and (iv) well-controlled ion-current-based
quantification. To obtain high-quality quantitative data necessary
to describe the 11 time-points protein expression temporal profiles,
strict criteria were used to define “quantifiable proteins”.
A total of 323 drug-responsive proteins were revealed with confidence,
and the time profiles of these proteins provided new insights into
the diverse temporal changes of biological cascades associated with
hepatic metabolism, response to hormone stimuli, gluconeogenesis,
inflammatory responses, and protein translation processes. Most profile
changes persisted well after the drug was eliminated. The developed
strategy can also be broadly applied in preclinical and clinical research,
where the analysis of numerous biological replicates is crucial.
The PK / PD of abatacept, a selective T-cell co-stimulation modulator, was examined in rats with collagen-induced arthritis (CIA) using a nonlinear mixed effect modeling approach. Male Lewis rats underwent collagen induction to produce rheumatoid arthritis. Two single-dose groups received either 10 mg/kg intravenous (IV) or 20 mg/kg subcutaneous (SC) abatacept, and one multiple-dose group received one 20 mg/kg SC abatacept dose and four additional 10 mg/kg SC doses. Effects on disease progression (DIS) were measured by paw swelling. Plasma concentrations of abatacept were assayed by enzyme-linked immunosorbent assay (ELISA). The PK / PD data were sequentially fitted using NONMEM VI. Goodness-of-fit was assessed by objective functions and visual inspection of diagnostic plots. The PK of abatacept followed a two-compartment model with linear elimination. For SC doses, short-term zero-order absorption was assumed with F = 59.2 %. The disease progression component was an indirect response model with a time-dependent change in paw edema production rate constant (kin) that was inhibited by abatacept. Variation in the PK data could be explained by inter-individual variability in clearance (CL) and central compartment volume (V1), while the large variability of the PD data may be the result of paw edema production (kin0) and loss rate constant (kout). Abatacept has modest effects on paw swelling in CIA rats. The PK / PD profiles were well described by the proposed model and allowed evaluation of inter-individual variability on drug- and DIS-related parameters.
Abatacept; arthritis; model; pharmacokinetics; pharmacodynamics; disease progression
It was hypothesized that expression profiling using gene arrays can be used to distinguish temporal patterns of changes in gene expression in response to a drug in vivo, and that these patterns can be used to identify groups of genes regulated by common mechanisms. A corticosteroid, methylprednisolone (MPL), was administered intravenously to a group of 47 rats (Rattus rattus) that were sacrificed at 17 timepoints over 72 h after MPL administration. Plasma drug concentrations and hepatic glucocorticoid receptors were measured from each animal. In addition, RNAs prepared from individual livers were used to query Affymetrix genechips for mRNA expression patterns. Statistical analyses using Affymetrix and GeneSpring software were applied to the results. Cluster analysis revealed six major temporal patterns containing 196 corticosteroid-responsive probe sets representing 153 different genes. Four clusters showed increased expression with differences in lag-time, onset rate, and/or duration of transcriptional effect. A fifth cluster showed rapid reduction persisting for 18 h. The final cluster identified showed decreased expression followed by an extended period of increased expression. These results lend new insights into the diverse hepatic genes involved in the physiologic, therapeutic, and adverse effects of corticosteroids and suggest that a limited array of control processes account for the dynamics of their pharmacogenomic effects.
Corticosteroids; Glucocorticoids; Expression profiling; Cluster analysis
A fifth-generation model for receptor/gene-mediated corticosteroid effects was proposed based on results from a 50 mg/kg IV bolus dose of methylprednisolone (MPL) in male adrenalectomized rats, and confirmed using data from other acute dosage regimens. Steady-state equations for receptor down-regulation and tyrosine aminotransferase (TAT) enzyme induction patterns were derived. Five groups of male Wistar rats (n=5/group) were subcutaneously implanted with Alzet mini-pumps primed to release saline or 0.05, 0.1, 0.2, and 0.3 mg/kg/hr of MPL for 7 days. Rats were sacrificed at the end of the infusion. Plasma MPL concentrations, blood lymphocyte counts, and hepatic cytosolic free receptor density, receptor mRNA, TAT mRNA, and TAT enzyme levels were quantitated. The pronounced steroid effects were evidenced by marked losses in body weights and changes in organ weights. All four treatments caused a dose-dependent reduction in hepatic receptor levels, which correlated with the induction of TAT mRNA and TAT enzyme levels. The 7 day receptor mRNA and free receptor density correlated well with the model predicted steady-state levels. However, the extent of enzyme induction was markedly higher than that predicted by the model suggesting that the usual receptor/gene-mediated effects observed upon single/intermittent dosing of MPL may be countered by alterations in other aspects of the system. A mean IC50 of 6.1 ng/mL was estimated for the immunosuppressive effects of methylprednisolone on blood lymphocytes. The extent and duration of steroid exposure play a critical role in mediating steroid effects and advanced PK/PD models provide unique insights into controlling factors.
pharmacodynamics; pharmacogenomics; methylprednisolone; tyrosine amino-transferase
The Goto-Kakizaki (GK) rat, a polygenic non-obese model of type 2 diabetes, is a useful surrogate for study of diabetes-related changes independent of obesity. GK rats and appropriate controls were killed at 4, 8, 12, 16 and 20 weeks post-weaning and differential muscle gene expression along with body and muscle weights, plasma hormones and lipids, and blood cell measurements were carried out. Gene expression analysis identified 204 genes showing 2-fold or greater differences between GK and controls in at least 3 ages. Array results suggested increased oxidative capacity in GK muscles, as well as differential gene expression related to insulin resistance, which was also indicated by HOMA-IR measurements. In addition, potential new biomarkers in muscle gene expression were identified that could be either a cause or consequence of T2DM. Furthermore, we demonstrate here the presence of chronic inflammation evident both systemically and in the musculature, despite the absence of obesity.
type 2 diabetes; skeletal muscle; inflammation; microarrays; gene expression
The pharmacogenomic effects of a corticosteroid (CS) were assessed in rat skeletal muscle using microarrays. Adrenalectomized (ADX) rats were treated with methylprednisolone (MPL) by either 50 mg/kg intravenous injection or 7-day 0.3 mg/kg/h infusion through subcutaneously implanted pumps. RNAs extracted from individual rat muscles were hybridized to Affymetrix Rat Genome Genechips. Data mining yielded 653 and 2316 CS-responsive probe sets following MPL bolus and infusion treatments. Of these, 196 genes were controlled by MPL under both dosing conditions. Cluster analysis revealed that 124 probe sets exhibited three typical expression dynamic profiles following acute dosing. Cluster A consisted of up-regulated probe sets which were grouped into five subclusters each exhibiting unique temporal patterns during the infusion. Cluster B comprised down-regulated probe sets which were divided into two subclusters with distinct dynamics during the infusion. Cluster C probe sets exhibited delayed down-regulation under both bolus and infusion conditions. Among those, 104 probe sets were further grouped into subclusters based on their profiles following chronic MPL dosing. Several mathematical models were proposed and adequately captured the temporal patterns for each subcluster. Multiple types of dosing regimens are needed to resolve common determinants of gene regulation as chronic exposure results in unexpected differences in gene expression compared to acute dosing. Pharmacokinetic/pharmacodynamic (PK/PD) modeling provides a quantitative tool for elucidating the complexities of CS pharmacogenomics in skeletal muscle.
Microarray studies; pharmacokinetics; pharmacodynamics; mathematical models; computational biology
This study examines methylprednisolone (MPL) effects on the dynamics of hepatic low-density lipoprotein receptor (LDLR) mRNA and plasma lipids associated with increased risks for atherosclerosis.
Materials and methods
Normal male Wistar rats were given 50 mg/kg MPL intramuscularly (IM) and sacrificed at various times. Measurements included plasma MPL and CST, hepatic glucocorticoid receptor (GR) mRNA, cytosolic GR density and hepatic LDLR mRNA, and plasma total cholesterol (TC), low-density lipoprotein cholesterol (LDLC), high density lipoprotein cholesterol (HDLC), and triglycerides (TG).
MPL showed bi-exponential disposition with two first-order absorption components. Hepatic GR and LDLR mRNA exhibited circadian patterns which were disrupted by MPL. Down-regulation in GR mRNA (40–50%) was followed by a delayed rebound phase. LDLR mRNA exhibited transient down-regulation (60–70%). Cytosolic GR density was significantly suppressed but returned to baseline by 72 h. Plasma TC and LDLC showed increases (55 and 142%) at 12 h. A mechanistic receptor/gene pharmacokinetic/pharmacodynamic model was developed to describe CS effects on hepatic LDLR mRNA and plasma cholesterols.
Our PK/PD model was able to satisfactorily capture the MPL effects on hepatic LDLR, its relationship to various plasma cholesterols, and builds the foundation to explore this area in the future.
cholesterol; corticosteroids; glucocorticoid receptors; LDL receptors; lipids; pharmacodynamics
The transcriptional response of skeletal muscle to chronic corticosteroid exposure was examined over 168 h and compared with the response profiles observed following a single dose of corticosteroid. Male adrenalectomized Wistar rats were given a constant-rate infusion of 0.3 mg•kg−1•h−1 methylprednisolone for up to 7 days via subcutaneously implanted minipumps. Four control and forty drug-treated animals were killed at ten different time points during infusion. Liver total RNAs were hybridized to 44 individual Affymetrix REA230A gene chips. Previously, we described a filtration approach for identifying genes of interest in microarray data sets developed from tissues of rats treated with methylprednisolone (MPL) following acute dosing. Here, a similar approach involving a series of three filters was applied sequentially to identify genes of interest. These filters were designed to eliminate probe sets that were not expressed in the tissue, not regulated by the drug, or did not meet defined quality control standards. Filtering eliminated 86% of probe sets, leaving a remainder of 2,316 for further consideration. In a previous study, 653 probe sets were identified as MPL regulated following administration of a single (acute) dose of the drug. Comparison of the two data sets yielded 196 genes identified as regulated by MPL in both dosing regimens. Because of receptor downregulation, it was predicted that genes regulated by receptor-glucocorticoid response element interactions would exhibit tolerance in chronic profiles. However, many genes did not exhibit steroid tolerance, indicating that present perspectives on the mechanism of glucocorticoid action cannot entirely explain all temporal profiles.
glucocorticoids; corticosteroids; Affymetrix gene chips; gene expression; time series
A retrospective analysis was performed to modify our fourth-generation pharmacodynamic model for glucocorticoid receptor (GR) dynamics with incorporation of more physiological features. This modified model was developed by integrating previously reported free cytosolic GR and GR mRNA data following single (10, 50 mg/kg) and dual (50 mg/kg at 0 and 24 hr) intravenous doses of methylprednisolone (MPL) in adrenalectomized (ADX) male Wistar rats with several in vitro studies describing real-time kinetics of the transfer of rat steroid-receptor complex from the cell cytosol to the nucleus. Additionally, free hepatic cytosolic GR and its mRNA data from a chronic infusion dosing study of MPL (0.1 and 0.3 mg/kg/hr) in male ADX Wistar rats were used to verify the predictability of the model. Incorporation of information regarding in vitro receptor kinetics allowed us to describe the receptor-mediated pharmacogenomic effects of MPL for a larger variety of genes in rat liver from microarray studies. These included early responsive gene like CCAAT/enhancer binding protein-β (CEBP-β), a transcription factor, as well as the later responsive gene for tyrosine aminotransferase (TAT), a classical biomarker of glucocorticoid (GC) genomic effects. This more mechanistic model of GR dynamics can be applied to characterize profiles for a greater number of genes in liver.
glucocorticoids; glucocorticoid receptor; nuclear localization; pharmacodynamics; methylprednisolone; pharmacogenomics
One of the goals of systems biology is the identification of regulatory mechanisms that govern an organism’s response to external stimuli. Transcription factors have been hypothesized as a major contributor to an organism’s response to various outside stimuli, and a great deal of work has been done to predict the set of transcription factors which regulate a given gene. Most of the current methods seek to identify possible binding sites from genomic sequence. Initial attempts at predicting transcription factors from genomic sequences suffered from the problem of false positives. Making the problem more difficult, it has also been shown that while predicted binding sites might be false positives, they can be shown to bind to their corresponding sequences in vitro. One method for rectifying this is through the use of phylogenetic analysis in which only regions which show high evolutionary conservation are analyzed. However such an approach may be too stringent because of the level of degeneracy shown in transcription factor binding site position weight matrices. Due to the degeneracy, there may be only a few bases that need to be conserved across species. Therefore, while a sequence may not show a high level of evolutionary conservation, these sequences may still show high affinity for the same transcription factor. In predicting transcription factor binding we explore the notion that “Co-expression implies co-regulation” [Allocco et al. BMC Bioinformatics 5:18, 2004]. With multiple genes requiring similar transcription factors binding sites, there exists a basis for eliminating false positives. This method allows for the selection of transcription factors binding sites that are active under a given experimental paradigm, thereby allowing us to indirectly incorporate the effects of chromosome and recognition site presentation upon transcription factor binding prediction. Rather than having to rationalize that a few transcription factors binding sites are over-represented in a cluster of genes, one can show that a few transcription factors are active in the cluster of genes that have been grouped together. Although the method focuses on predicting experiment-specific transcription factor binding sites, it is possible that if such a methodology were used in an iterative process where different experiments were analyzed, one could obtain a comprehensive set of transcription factors binding sites which regulate the various dynamic responses shown by biological systems under a variety of conditions hence building a more comprehensive model of transcriptional regulation.
Corticosteroids; Gene expression; Transcription factor binding site; Phylogenetics
Microarray analyses were performed on livers from adrenalectomized male Wistar rats chronically infused with methylprednisolone (MPL) (0.3 mg/kg·h) using Alzet mini-osmotic pumps for periods ranging from 6 h to 7 d. Four control and 40 drug-treated animals were killed at 10 different times during drug infusion. Total RNA preparations from the livers of these animals were hybridized to 44 individual Affymetrix REA230A gene chips, generating data for 15,967 different probe sets for each chip. A series of three filters were applied sequentially. These filters were designed to eliminate probe sets that were not expressed in the tissue, were not regulated by the drug, or did not meet defined quality control standards. These filters eliminated 13,978 probe sets (87.5%) leaving a remainder of 1989 probe sets for further consideration. We previously described a similar dataset obtained from animals after administration of a single dose of MPL (50 mg/kg given iv). That study involved 16 time points over a 72-h period. A similar filtering schema applied to the single-bolus-dose data-set identified 1519 probe sets as being regulated by MPL. A comparison of datasets from the two different dosing regimens identified 358 genes that were regulated by MPL in response to both dosing regimens. Regulated genes were grouped into 13 categories, mainly on gene product function. The temporal profiles of these common genes were subjected to detailed scrutiny. Examination of temporal profiles demonstrates that current perspectives on the mechanism of glucocorticoid action cannot entirely explain the temporal profiles of these regulated genes.
Methylprednisolone (MPL) pharmacokinetics was examined in adrenalectomized (ADX) and normal rats to assess the feasibility of intramuscular (i.m.) dosing for use in pharmacodynamic studies. Several study phases were pursued. Parallel group studies were performed in normal and ADX rats given 50 mg/kg MPL (i.v. or i.m.) and blood samples were collected up to 6 h. Data from studies where normal rats were dosed with 50 mg/kg MPL i.m. and killed over either 6 or 96 h were combined to determine muscle site and plasma MPL concentrations. Lastly, ADX rats were dosed with 50 mg/kg MPL i.m. and killed over 18 h to assess hepatic tyrosine aminotransferase (TAT) dynamics. MPL exhibited bi-exponential kinetics after i.v. dosing with a terminal slope of 2.1 h−1. The i.m. drug was absorbed slowly with two first-order absorption rate constants, 1.26 and 0.219 h−1 indicating flip-flop kinetics with overall 50% bioavailability. The kinetics of MPL at the injection site exhibited slow, dual absorption rates. Although i.m. MPL showed lower bioavailability compared with other corticosteroids in rats, TAT dynamics revealed similar i.m. and i.v. response profiles. The more convenient intramuscular dosing can replace the i.v. route without causing marked differences in pharmacodynamics.
methylprednisolone; corticosteroids; pharmacokinetics; intramuscular injection; tyrosine aminotransferase
Microarrays have been utilized in many biological, physiological and pharmacological studies as a high-throughput genomic technique. Several generations of Affymetrix GeneChip® microarrays are widely used in gene expression studies. However, differences in intensities of signals for different probe sets that represent the same gene on various types of Affymetrix chips make comparison of datasets complicated.
Materials and Methods
A power coefficient scaling factor was applied in the pharmacokinetic/ pharmacodynamic (PK/PD) modeling to account for differences in probe set sensitivities (i.e., signal intensities). Microarray data from muscle and liver following methylprednisolone 50 mg/kg i.v. bolus and 0.3 mg/kg/h infusion regimens were taken as an exemplar.
The scaling factor applied to the pharmacodynamic output function was used to solve the problem of intensity differences between probe sets. This approach yielded consistent pharmacodynamic parameters for the applied models.
Modeling of pharmacodynamic/pharmacogenomic (PD/PG) data from diverse chips should be performed with caution due to differential probe set intensities. In such circumstances, a power scaling factor can be applied in the modeling.
bioinformatics; computational biology; pharmacodynamics; pharmacogenomics; pharmacokinetics
Receptor/gene-mediated effects of corticosteroids on hepatic tyrosine aminotransferase (TAT) were evaluated in normal rats. A group of normal male Wistar rats were injected with 50 mg/kg methylprednisolone (MPL) intramuscularly at the nadir of their plasma corticosterone (CST) rhythm (early light cycle) and sacrificed at various time points up to 96 h post-treatment. Blood and livers were collected to measure plasma MPL, CST, hepatic glucocorticoid receptor (GR) mRNA, cytosolic GR density, TAT mRNA, and TAT activity. The pharmacokinetics of MPL showed bi-exponential disposition with two first-order absorption components from the injection site and bioavailability was 21%. Plasma CST was reduced after MPL dosing, but resumed its daily circadian pattern within 36 h. Cytosolic receptor density was significantly suppressed (90%) and returned to baseline by 72 h resuming its biphasic pattern. Hepatic GR mRNA follows a circadian pattern which was disrupted by MPL and did not return during the study. MPL caused significant down-regulation (50%) in GR mRNA which was followed by a delayed rebound phase (60–70 h). Hepatic TAT mRNA and activity showed up-regulation as a consequence of MPL, and returned to their circadian baseline within 72 and 24 h of treatment. A mechanistic receptor/gene-mediated pharmacokinetic/pharmacodynamic model was able to satisfactorily describe the complex interplay of exogenous and endogenous corticosteroid effects on hepatic GR mRNA, cytosolic free GR, TAT mRNA, and TAT activity in normal rats.
Methylprednisolone; Corticosteroids; Pharmacokinetics; Pharmacodynamics; Tyrosine aminotransferase; Glucocorticoid receptors
The circadian rhythm of endogenous corticosterone (CS) may produce fluctuations of downstream gene expression in normal rats. This study examined changes in glucocorticoid receptor (GR) and glutamine synthetase (GS) expression in rat skeletal muscle in relation to plasma CS over a 24-h period.
Fifty-four normal male Wistar rats were sacrificed at 18 time points (n = 3) over 24 h. Plasma CS concentrations and gastrocnemius muscle GR and GS mRNA and GS activity were measured.
The circadian rhythm of plasma CS was captured by a two-harmonic function. The expression of GR and GS mRNA and GS activity follow a circadian rhythm in normal rat skeletal muscle. GR mRNA reaches a trough at 4 h after the peak of plasma CS and it fluctuates between 0.55 and 0.9 fmol g tissue−1. GS mRNA and activity reach peaks at 6 and 12 h after the endogenous CS peak. GS mRNA oscillates between 3 and 6 fmol g tissue−1, whereas GS activity fluctuates between 17 and 23 µmol min−1 g protein−1. Mechanistic receptor/gene-mediated pharmacodynamic models were applied to describe the temporal patterns of GR mRNA, GS mRNA, and GS activity within the circadian cycle.
The integrated models were able to capture the circadian expression patterns of plasma CS, and GR and GS in normal rat skeletal muscle showing a dependence of tissue gene expression on plasma CS.
biological rhythm; mathematical model; pharmacodynamics; pharmacokinetics
Pyruvate dehydrogenase kinase 4 (PDK4) is a lipid status responsive gene involved in muscle fuel selection. Evidence is mounting in support of the therapeutic potential of PDK4 inhibitors to treat diabetes. Factors that regulate PDK4 mRNA expression include plasma corticosterone, insulin and free fatty acids. Our objective was to determine the impact of those plasma factors on PDK4 mRNA and to develop and validate a population mathematical model to differentiate aging, diet and disease effects on muscle PDK4 expression. The Goto-Kakizaki (GK) rat, a polygenic non-obese model of type 2 diabetes, was used as the diabetic animal model. We examined muscle PDK4 mRNA expression by real-time QRTPCR. Groups of GK rats along with controls fed with either a normal or high fat diet were sacrificed at 4, 8, 12, 16, and 20 weeks of age. Plasma corticosterone, insulin and free fatty acid were measured. The proposed mechanism-based model successfully described the age, disease and diet effects and the relative contribution of these plasma regulators on PDK4 mRNA expression. Muscle growth reduced the PDK4 mRNA production rate by 14% per gram increase. High fat diet increased the initial production rate constant in GK rats by 2.19-fold. The model indicated that corticosterone had a moderate effect and PDK4 was more sensitive to free fatty acid than insulin fluxes, which was in good agreement with the literature data.
population model; type 2 diabetes; disease progression; PDK4; Goto-Kakizaki rats
Methotrexate (MTX) is an anchor drug used to treat rheumatoid arthritis (RA), but responsiveness is variable in effectiveness and toxicity. Methotrexate and its polyglutamate conjugates (MTXPGn) in red blood cells (RBC) have been associated with patient response. In the current study, 13 collagen-induced arthritic (CIA) rats and 12 healthy rats were given subcutaneous doses of either saline or 0.3 or 1.5 mg/kg per 2 days of MTX from day 21 to 43 post-induction. Blood samples were obtained at various times to measure MTX in plasma, and MTX and MTXPGn in RBC. Effects on disease progression were indicated by body weight and paw size. After multiple-doses, RBC MTX reached steady-state (82.4 nM) within 4 days. The MTXPG2 and MTXPG3 in RBC kept increasing until the end of the study attaining 12.5 and 17.7 nM. Significant weight loss was observed after dosing of 1.5 mg/kg/2 days, whereas moderate effectiveness was observed after dosing of 0.3 mg/kg/2 days. A pharmacokinetic/ pharmacodynamic/disease (PK/PD/DIS) model with indirect mechanisms and transduction components incorporating plasma MTX, RBC MTX, and RBC MTXPGn concentrations, and paw size was developed using naïve data pooling and ADAPT 5. The PK/PD in CIA rats dosed at 0.3 mg/kg/2 days were captured well by our proposed model. MTX showed modest (Imaxd = 0.16) but sensitive (IC50d = 0.712 nM) effectiveness on paw edema. The higher dose produced toxicity. The proposed model offers improved understanding of MTX effects on rheumatoid arthritis.
Methotrexate; rheumatoid arthritis; pharmacokinetics; pharmacodynamics; disease progression
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.
glucocorticoids; glucocorticoid receptor; GRβ; qRTPCR
Chronopharmacology is an important but under-explored aspect of therapeutics. Rhythmic variations in biological processes can influence drug action, including pharmacodynamic responses, due to circadian variations in the availability or functioning of drug targets. We hypothesized that global gene expression analysis can be useful in the identification of circadian regulated genes involved in drug action. Circadian variations in gene expression in rat liver were explored using Affymetrix gene arrays. A rich time series involving animals analyzed at 18 time points within the 24 hour cycle was generated. Of the more than 15,000 probe sets on these arrays, 265 exhibited oscillations with a 24 hour frequency. Cluster analysis yielded 5 distinct circadian clusters, with approximately two-thirds of the transcripts reaching maximum expression during the animal’s dark/active period. Of the 265 probe sets, 107 of potential therapeutic importance were identified. The expression levels of clock genes were also investigated in this study. Five clock genes exhibited circadian variation in liver, and data suggest that these genes may also be regulated by corticosteroids.
A retrospective meta-modeling analysis was performed to integrate previously reported data of glucocorticoid (GC) effects on glucose regulation following a single intramuscular dose (50 mg/kg), single intravenous doses (10, 50 mg/kg), and intravenous infusions (0.1, 0.2, 0.3 and 0.4 mg/kg/h) of methylprednisolone (MPL) in normal and adrenalectomized (ADX) male Wistar rats. A mechanistic pharmacodynamic (PD) model was developed based on the receptor/gene/protein-mediated GC effects on glucose regulation. Three major target organs (liver, white adipose tissue and skeletal muscle) together with some selected intermediate controlling factors were designated as important regulators involved in the pathogenesis of GC-induced glucose dysregulation. Assessed were dynamic changes of food intake and systemic factors (plasma glucose, insulin, free fatty acids (FFA) and leptin) and tissue-specific biomarkers (cAMP, phosphoenolpyruvate carboxykinase (PEPCK) mRNA and enzyme activity, leptin mRNA, interleukin 6 receptor type 1 (IL6R1) mRNA and Insulin receptor substrate-1 (IRS-1) mRNA) after acute and chronic dosing with MPL along with the GC receptor (GR) dynamics in each target organ. Upon binding to GR in liver, MPL dosing caused increased glucose production by stimulating hepatic cAMP and PEPCK activity. In adipose tissue, the rise in leptin mRNA and plasma leptin caused reduction of food intake, the exogenous source of glucose input. Down-regulation of IRS-1 mRNA expression in skeletal muscle inhibited the stimulatory effect of insulin on glucose utilization further contributing to hyperglycemia. The nuclear drug-receptor complex served as the driving force for stimulation or inhibition of downstream target gene expression within different tissues. Incorporating information such as receptor dynamics, as well as the gene and protein induction, allowed us to describe the receptor-mediated effects of MPL on glucose regulation in each important tissue. This advanced mechanistic model provides unique insights into the contributions of major tissues and quantitative hypotheses for the multi-factor control of a complex metabolic system.
The pharmacokinetics (PK) of salsalate (SS) and salicylic acid (SA) was assessed in normal Wistar and diabetic Goto-Kakizaki rats. Three PK studies were conducted: 1) PK of SA in normal rats after intravenous dosing of SA at 20, 40, 80 mg/kg. 2) PK of SS and SA in normal rats after oral dosing of SS at 28, 56, 112 mg/kg. 3) PK during 4 months feeding of SS-containing diet in both normal and diabetic rats. The disposition of SS and SA were simultaneously evaluated using a pharmacokinetic model comprised of several transit absorption steps and linear and nonlinear dual elimination pathways for SA. The results indicated that the nonlinear elimination pathway of SA only accounted for a small fraction of the total clearance (< 12%) at therapeutic concentrations. A flat profile of SA was observed after oral dosing SS, particularly at a high dose. The possible reasons for this flat profile were posed. During the SS-diet feeding, diabetic rats achieved lower blood concentrations of SA than normal rats with a higher apparent clearance (CL/F) possibly due to incomplete (47%) bioavailability. Such CL/F decreased with age in both diabetic and normal rats. The effect of diabetes on SA pharmacokinetics may necessitate increased dosing in future usage of SS in diabetes.
salsalate; salicylic acid; pharmacokinetics; diabetes
Kidney is a major target for adverse effects associated with corticosteroids. A microarray dataset was generated to examine changes in gene expression in rat kidney in response to methylprednisolone. Four control and 48 drug-treated animals were killed at 16 times after drug administration. Kidney RNA was used to query 52 individual Affymetrix chips, generating data for 15,967 different probe sets for each chip. Mining techniques applicable to time series data that identify drug-regulated changes in gene expression were applied. Four sequential filters eliminated probe sets that were not expressed in the tissue, not regulated by drug, or did not meet defined quality control standards. These filters eliminated 14,890 probe sets (94%) from further consideration. Application of judiciously chosen filters is an effective tool for data mining of time series datasets. The remaining data can then be further analyzed by clustering and mathematical modeling. Initial analysis of this filtered dataset identified a group of genes whose pattern of regulation was highly correlated with prototype corticosteroid enhanced genes. Twenty genes in this group, as well as selected genes exhibiting either downregulation or no regulation, were analyzed for 5′ GRE half-sites conserved across species. In general, the results support the hypothesis that the existence of conserved DNA binding sites can serve as an important adjunct to purely analytic approaches to clustering genes into groups with common mechanisms of regulation. This dataset, as well as similar datasets on liver and muscle, are available online in a format amenable to further analysis by others.
data mining; gene arrays; glucocorticoids; pharmacogenomics; evolutionary conservation
Increased expression of inducible nitric oxide synthase (iNOS) resulting in nitric oxide elevation represents an important component of inflammatory responses. We assess the effects of methylprednisolone (MPL) on these processes during endotoxin-induced acute inflammation and provide a mechanism-based model to quantitatively describe them.
Male Lewis rats were dosed with lipopolysaccharide (50 μg/kg LPS) alone or with methylprednisolone (10 and 50 mg/kg) and sacrificed at different time points. Plasma MPL, lung iNOS mRNA expression, plasma nitric oxide (NO) and other physiological factors were measured. Sodium nitrate (750 μmole/kg) was given to a separate cohort of rats to assess NO disposition kinetics. PK-PD modeling was performed with ADAPT 5.
Disposition kinetics of plasma MPL and NO showed bi-exponential decline and were described by two-compartment models. LPS increased expression of iNOS mRNA in lung and increased plasma NO, while MPL dosing palliated this increase in a dose-dependent manner. These effects were well captured using tandem indirect response and precursor-pool models.
The model provides a quantitative assessment of the suppression of NO production by MPL and shows that the major effects are at the transcriptional level by reducing expression of iNOS mRNA.
corticosteroids; inflammation; iNOS; nitric oxide; PK-PD modeling
The dynamics of aging and type 2 diabetes (T2D) disease progression were investigated in normal [Wistar-Kyoto (WKY)] and diabetic [Goto-Kakizaki (GK)] rats and a mechanistic disease progression model was developed for glucose, insulin, and glycosylated hemoglobin (HbA1c) changes over time. The study included 30 WKY and 30 GK rats. Plasma glucose and insulin, blood glucose and HbA1c concentrations and hematological measurements were taken at ages 4, 8, 12, 16 and 20 weeks. A mathematical model described the development of insulin resistance (IR) and β-cell function with age/growth and diabetes progression. The model utilized transit compartments and an indirect response model to quantitate biomarker changes over time. Glucose, insulin and HbA1c concentrations in WKY rats increased to a steady-state at 8 weeks due to developmental changes. Glucose concentrations at 4 weeks in GK rats were almost twice those of controls, and increased to a steady-state after 8 weeks. Insulin concentrations at 4 weeks in GK rats were similar to controls, and then hyperinsulinemia occurred until 12–16 weeks of age indicating IR. Subsequently, insulin concentrations in GK rats declined to slightly below WKY controls due to β-cell failure. HbA1c showed a delayed increase relative to glucose. Modeling of HbA1c was complicated by age-related changes in hematology in rats. The diabetes model quantitatively described the glucose/insulin inter-regulation and HbA1c production and reflected the underlying pathogenic factors of T2D—IR and β-cell dysfunction. The model could be extended to incorporate other biomarkers and effects of various anti-diabetic drugs.
Type 2 diabetes; Disease progression modeling; Insulin resistance; β-cell function
This comparison employs mathematical disease progression models to identify a rat model of arthritis with the least inter-animal variability and features lending to better study designs.
Arthritis was induced with either collagen (CIA) or mycobacterium (AIA) in either Lewis or Dark Agouti (DA) rats. Disease progression was monitored by paw edema and body weight. Models with production, loss, and feedback components were constructed and population analysis using NONMEM software was employed to identify inter-animal variability in the various disease progression parameters.
Onset time was the only parameter different within all four groups (DA–AIA 11.5 days, DA–CIA 16.5 days, Lewis–AIA 11.9 days, Lewis–CIA 13.9 days). The loss-of-edema rate constant was 20% slower in DA (0.362 h−1) than Lewis (0.466 h−1) rats. Most models exhibited peak paw edema 20 days post-induction. Edema in CIA returned to 150% of the initial value after the disease peaked. DA rats displayed more severe overall responses.
No statistical differences between groups were observed for inter-animal variation in disease onset, progression and severity parameters. Onset time varies and should be noted in the design of future studies. DA rats may offer a more dynamic range of edema response than Lewis rats.
arthritis; disease; model; progression; rat