Both obesity and chronic inflammation are often associated with insulin resistance and type 2 diabetes. The Zucker diabetic fatty (ZDF) rat (fa/fa) is an obese animal model frequently used in type 2 diabetes research. The current study determines whether chronic administration (from 5 weeks of age through 24 weeks of age) of salsalate, a salicylate with anti-inflammatory properties, would be effective in mitigating diabetes disease progression in ZDF rats. Although a trend existed for lower blood glucose in the salsalate-treated group, significant differences were obscured by high animal-level variability. However, even in the non-drug-treated group, not all ZDF rats became diabetic as expected. Therefore, animals were parsed into two groups, regardless of drug treatment: normoglycemic ZDF rats, which maintained blood glucose profiles identical to nondiabetic Zucker lean rats (ZLRs), and hyperglycemic ZDF rats, which exhibited progressive elevation in blood glucose. To ascertain the differences between ZDF rats that became hyperglycemic and those that did not, relevant physiological indices and expression levels of adiponectin, tumor necrosis factor-α, interleukin-6, and glucocorticoid-induced leucine zipper messenger RNAs in adipose tissue were measured at sacrifice. Plasma C-reactive protein concentrations and expression levels of cytokine and glucocorticoid-induced leucine zipper messenger RNAs suggested more prevalent chronic inflammation in hyperglycemic animals. Early elevation of the insulin-sensitizing adipokine, adiponectin, was present in both ZDF groups, with the rate of its age-related decline faster in hyperglycemic animals. The most marked difference between the two groups of ZDF animals was in insulin output. Although the two ZDF populations had very similar elevated plasma insulin concentrations for the first 10 weeks, after that time, plasma insulin decreased markedly in the animals that became hyperglycemic, whereas it remained high in the normoglycemic ZDF rats. Thus, hyperglycemic ZDF animals exhibit both insulin resistance and progressive beta cell failure, whereas normoglycemic ZDF rats exhibit a lesser degree of insulin resistance that does not progress to beta cell failure. In these respects, the normoglycemic ZDF rats appear to revert back to a phenotype that strongly resembles that of nondiabetic Zucker fatty rats from which they were derived.
type 2 diabetes; ZDF rats; animal models
This study investigates the antidiabetic effects of salsalate on disease progression of diabetes in non-obese diabetic Goto-Kakizaki (GK) rats, an experimental model of type 2 diabetes. Salsalate was formulated in rat chow (1,000 ppm) and used to feed rats from 5 to 21 weeks of age. At 5 weeks of age, GK and Wistar (WIS) control rats were subdivided into four groups, each composed of six rats: GK rats with standard diet (GK-C); GK rats with salsalate-containing diet (GK-S); WIS rats with standard diet (WIS-C); and WIS rats with salsalate-containing diet (WIS-S). The GK-C rats (167.2±11.6 mg/dL) showed higher blood glucose concentrations than WIS-C rats (133.7±4.9 mg/dL, P<0.001) at the beginning of the experiment, and had substantially elevated blood glucose from an age of 15 weeks until sacrifice at 21 weeks (341.0±133.6 mg/dL). The GK-S rats showed an almost flat profile of blood glucose from 4 weeks (165.1±11.0 mg/dL) until sacrifice at 21 weeks of age (203.7±22.2 mg/dL). While this difference in blood glucose between 4 and 21 weeks in GK-S animals was significant, blood glucose at 21 weeks was significantly lower in GK-S compared to GK-C animals. At sacrifice, salsalate decreased plasma insulin (GK-S =1.0±0.3; GK-C =2.0±0.3 ng/mL, P<0.001) and increased plasma adiponectin concentrations (GK-S =15.9±0.7; GK-C =9.7±2.0 μg/mL, P<0.001). Salsalate also lowered total cholesterol in GK-S rats (96.1±8.5 mg/dL) compared with GK-C rats (128.0±11.4 mg/dL, P<0.001). Inflammation-related genes (Ifit1 and Iigp1) exhibited much higher mRNA expression in GK-C rats than WIS-C rats in liver, adipose, and muscle tissues, while salsalate decreased the Ifit1 and Iigp1 mRNA only in adipose tissue. These results suggest that salsalate acts to both increase adiponectin and decrease adipose tissue-based inflammation while preventing type 2 diabetes disease progression in GK rats.
type 2 diabetes; salicylates; inflammation; adiponectin
Development and progression of type 2 diabetes is a complex interaction between genetics and environmental influences. High dietary fat is one environmental factor that is conducive to the development of insulin-resistant diabetes. In the present report, we compare the responses of lean poly-genic, diabetic Goto-Kakizaki (GK) rats to those of control Wistar-Kyoto (WKY) rats fed a high fat diet from weaning to 20 weeks of age. This comparison included a wide array of physiological measurements along with gene expression profiling of abdominal adipose tissue using Affymetrix gene array chips. Animals of both strains fed a high fat diet or a normal diet were sacrificed at 4, 8, 12, 16, and 20 weeks for this comparison. The microarray analysis revealed that the two strains developed different adaptations to increased dietary fat. WKY rats decrease fatty acid synthesis and lipogenic processes whereas GK rats increase lipid elimination. However, on both diets the major differences between the two strains remained essentially the same. Specifically relative to the WKY strain, the GK strain showed lipoatrophy, chronic inflammation, and insulin resistance.
diabetes; high fat diet; gene expression; microarray
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.
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.
Corticosteroids; GILZ; glucocorticoids; qRT-PCR
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
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
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 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
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
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
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
Although corticosteroids (CSs) affect gene expression in multiple tissues, the array of genes that are regulated by these catabolic steroids is diverse, highly tissue specific, and depends on their functions in the tissue. Liver has many important functions in performing and regulating diverse metabolic processes. Muscle, in addition to its mechanical role, is critical in maintaining systemic energy homeostasis and accounts for about 80% of insulin-directed glucose disposal. Consequently, a better understanding of CS pharmacogenomic effects in these tissues would provide valuable information regarding the tissue-specificity of transcriptional dynamics, and would provide insights into the underlying molecular mechanisms of action for both beneficial and detrimental effects.
We performed an integrated analysis of transcriptional data from liver and muscle in response to methylprednisolone (MPL) infusion, which included clustering and functional annotation of clustered gene groups, promoter extraction and putative transcription factor (TF) identification, and finally, regulatory closeness (RC) identification.
This analysis allowed the identification of critical transcriptional responses and CS-responsive functions in liver and muscle during chronic MPL administration, the prediction of putative transcriptional regulators relevant to transcriptional responses of CS-affected genes which are also potential secondary bio-signals altering expression levels of target-genes, and the exploration of the tissue-specificity and biological significance of gene expression patterns, CS-responsive functions, and transcriptional regulation.
The analysis provided an integrated description of the genomic and functional effects of chronic MPL infusion in liver and muscle.
liver; muscle; glucocorticoids; corticosteroids; gene expression; gene regulation; promoter analysis
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.