Biochemical analysis of mutations in P450 oxidoreductase

All microsomal P450s require POR (cytochrome P450 reductase) for catalytic activity. Most of the clinically used drugs are metabolized by a small number of P450s and polymorphisms in the cytochrome P450s are known to cause changes in drug metabolism. We have recently found a number of POR missense mutations in the patients with disordered steroidogenesis. Our initial report described five missense mutations (A284P, R454H, V489E, C566Y and V605F) identified in four patients. We built bacterial expression vectors for each POR variant, purified the membranes expressing normal or variant POR and characterized their activities with cytochrome c and P450c17 assays. We have recently completed an extensive study of the range of POR mutations and characterized the mutants/polymorphisms A112V, T139A, M260V, Y456H, A500V, G536R, L562P, R613X, V628I and F643del from sequencing of patient DNA. We also studied POR variants Y179D, P225L, R313W, G410S and G501R that were available in databases or the published literature. We analysed the mutations with a three-dimensional model of human POR that was based on an essentially similar rat POR with known crystal structure. The missense mutations found in patients with disordered steroidogenesis mapped to functionally important domains of POR and the apparent polymorphisms mapped to less crucial regions. Since a variation in POR can alter the activity of all microsomal P450s, it can also affect the drug metabolism even with a normal P450. Understanding the genetic and biochemical basis of POR-mediated drug metabolism will provide valuable information about possible differences in P450-mediated reactions among the individuals carrying a variant or polymorphic form of POR.

Statin therapy induces ultrastructural damage in skeletal muscle in patients without myalgia

Muscle pain and weakness are frequent complaints in patients receiving 3-hydroxymethylglutaryl coenzymeA (HMG CoA) reductase inhibitors (statins). Many patients with myalgia have creatine kinase levels that are either normal or only marginally elevated, and no obvious structural defects have been reported in patients with myalgia only. To investigate further the mechanism that mediates statin-induced skeletal muscle damage, skeletal muscle biopsies from statin-treated and non-statin-treated patients were examined using both electron microscopy and biochemical approaches. The present paper reports clear evidence of skeletal muscle damage in statin-treated patients, despite their being asymptomatic. Though the degree of overall damage is slight, it has a characteristic pattern that includes breakdown of the T-tubular system and subsarcolemmal rupture. These characteristic structural abnormalities observed in the statin-treated patients were reproduced by extraction of cholesterol from skeletal muscle fibres in vitro. These findings support the hypothesis that statin-induced cholesterol lowering per se contributes to myocyte damage and suggest further that it is the specific lipid/protein organization of the skeletal muscle cell itself that renders it particularly vulnerable.

Neuregulin-1 beta attenuates doxorubicin-induced alterations of excitation-contraction coupling and reduces oxidative stress in adult rat cardiomyocytes

Treatment of metastatic breast cancer with doxorubicin (Doxo) in combination with trastuzumab, an antibody targeting the ErbB2 receptor, results in an increased incidence of heart failure. Doxo therapy induces reactive oxygen species (ROS) and alterations of calcium homeostasis. Therefore, we hypothesized that neuregulin-1 beta (NRG), a ligand of the cardiac ErbB receptors, reduces Doxo-induced alterations of EC coupling by triggering antioxidant mechanisms. Adult rat ventricular cardiomyocytes (ARVM) were isolated and treated for 18-48 h. SERCA protein was analyzed by Western blot, EC coupling parameters by fura-2 and video edge detection, gene expression by RT-PCR, and ROS by DCF-fluorescence microscopy. At clinically relevant doses Doxo reduced cardiomyocytes contractility, SERCA protein and SR calcium content. NRG, similarly as the antioxidant N-acetylcystein (NAC), did not affect EC coupling alone, but protected against Doxo-induced damage. NRG and Doxo showed an opposite modulation of glutathione reductase gene expression. NRG, similarly as NAC, reduced peroxide- or Doxo-induced oxidative stress. Specific inhibitors showed, that the antioxidant action of NRG depended on signaling via the ErbB2 receptor and on the Akt- and not on the MAPK-pathway. Therefore, NRG attenuates Doxo-induced alterations of EC coupling and reduces oxidative stress in ARVM. Inhibition of the ErbB2/NRG signaling pathway by trastuzumab in patients concomitantly treated with Doxo might prevent beneficial effects of NRG in the myocardium.

Differences in catalytic activity between rat testicular and ovarian carbonyl reductases are due to two amino acids

The sequences of rat testis carbonyl reductase (rCR1) and rat ovary carbonyl reductase (rCR2) are 98% identical, differing only at amino acids 140, 141, 143, 235 and 238. Despite such strong sequence identity, we find that rCR1 and rCR2 have different catalytic constants for metabolism of menadione and 4-benzoyl-pyridine. Compared to rCR1, rCR2 has a 20-fold lower K(m) and 5-fold lower k(cat) towards menadione and a 7-fold lower K(m) and 7-fold lower k(cat) towards 4-benzoyl-pyridine. We constructed hybrids of rCR1 and rCR2 that were changed at either residues 140, 141 and 143 or residues 235 and 238. rCR1 with residues 140, 141 and 143 of rCR2 has similar catalytic efficiency for menadione and 4-benzoyl-pyridine as rCR1. rCR1 with Thr-235 and Glu-238 of rCR2 has the catalytic constants of rCR2, indicating that it is this part of rCR2 that contributes to its lower K(m) for menadione and 4-benzoyl-pyridine. Comparisons of three-dimensional models of rCR1 and rCR2 show how Thr-235 and Glu-238 stabilize rCR2 binding of NADPH and menadione.

Ontogeny of mRNA abundance of nuclear receptors and nuclear receptor target genes in young cattle

After birth the development of appropriate detoxification mechanisms is important. Nuclear receptors (NR), such as constitutive androstane receptor (CAR), pregnane X receptor (PXR), peroxisome proliferator-activated receptor-alpha (PPARalpha), retinoid receptors (RAR, RXR), and NR target genes are involved in the detoxification of exogenous and endogenous substances. We quantified abundances of hepatic mRNA of NR and several NR target genes (cytochromes, CYP; cytochrome P450 reductase, CPR; UDP-glucuronosyl transferase, UDP) in calves at different ages. Gene expression was quantified by real-time RT-PCR. Abundance of mRNA of CAR and PXR increased from low levels at birth in pre-term calves (P0) and full-term calves (F0) to higher levels in 5-day-old calves (F5) and in 159-day-old veal calves (F159), whereas mRNA levels of PPARalpha did not exhibit significant ontogenetic changes. RARbeta mRNA levels were higher in F5 and F159 than in F0, whereas no age differences were observed for RARalpha levels. Levels of RXRalpha and RXRbeta mRNA were lower in F5 than in P0 and F0. Abundance of CYP2C8 and CYP3A4 increased from low levels in P0 and F0 to higher levels in F5 and to highest levels in F159. Abundance of CPR was transiently decreased in F0 and F5 calves. Levels of UGT1A1 mRNA increased from low levels in P0 and F0 to maximal level in F5 and F159. In conclusion, mRNA levels of NR and NR target genes exhibited ontogenetic changes that are likely of importance for handling of xeno- and endobiotics with increasing age.

Nuclear receptor and nuclear receptor target gene messenger ribonucleic acid levels at different sites of the gastrointestinal tract and in liver of healthy dogs

Nuclear receptors (NR) are ligand-activated transcription factors that regulate different metabolic pathways by influencing the expression of target genes. The current study examined mRNA abundance of NR and NR target genes at different sites of the gastrointestinal tract (GIT) and the liver of healthy dogs (Beagles; n = 11). Samples of GIT and liver were collected postmortem and homogenized, total RNA was extracted and reverse transcribed, and gene expression was quantified by real-time reverse-transcription PCR relative to the mean of 3 housekeeping genes (beta-actin, glyceraldehyde-3-phosphate dehydrogenase, and ubi-quitin). Differences were observed (P < or = 0.05) in the mRNA abundance among stomach (St), duodenum (Du), jejunum (Je), ileum (Il), and colon (Col) for NR [pregnane X receptor (Du, Je > Il, Col > St), peroxisome proliferator-associated receptor gamma (St, Du, Col > Je, Il), constitutive androstane receptor (Je, Du > Il, Col), and retinoid x receptor alpha (Du > Il)] and NR target genes [glutathione-S-transferase A3-3 (Du > Je > St, Il; St > Col), phenol-sulfating phenol sulfotransferase 1A1 (Du, Je > Il, St; Col > St), cytochrome P450 3A12 (Du, Je > St, Il, Col), multiple drug resistance gene 1 (Du, Je, Il, Col > St), multiple drug resistance-associated protein 2 (Je, Du > Il > St, Col), multiple drug resistance-associated protein 3 (Col > St > Il; Du > Je, Il; St > Il), NR corepressor 2 (St > Il, Col), and cytochrome P450 reductase (St, Du, Je > Il, Col)], but not for peroxisome proliferator-associated receptor alpha. Differences (P > 0.05) in mRNA abundance in the liver relative to the GIT were also observed. In conclusion, the presence of numerous differences in expression of NR and NR target genes in different parts of the GIT and in liver of healthy dogs may be associated with location-specific functions and regulation of GIT regions.

Effects of dexamethasone on mRNA abundance of nuclear receptors and hepatic nuclear receptor target genes in neonatal calves

Hepatic nuclear receptors (NR), particularly constitutive androstane receptor (CAR) and pregnane X receptor (PXR), are involved in the coordinated transcriptional control of genes that encode proteins involved in the metabolism and detoxification of xeno- and endobiotics. A broad spectrum of metabolic processes are mediated by NR acting in concert with ligands such as glucocorticoids. This study examined the role of dexamethasone on hepatic mRNA expression of CAR, PXR and several NR target genes. Twenty-eight male calves were allotted to one of four treatment groups in a 2 x 2 arrangement of treatments: feed source (colostrum or milk-based formula) and glucocorticoid administration (twice daily intramuscular dexamethasone). Liver biopsies were obtained at 5 days of age. Real-time reverse transcription polymerase chain reaction was used to quantify mRNA abundances. No effects of feed source on mRNA abundances were observed. For the NR examined, mRNA abundance of both CAR and PXR in dexamethasone-treated calves was lower (p < 0.05) by 39% and 40%, respectively, than in control calves. Abundance of NR target genes exhibited a mixed response. SULT1A1 mRNA abundance was 39% higher (p < 0.05) in dexamethasone-treated calves compared with control calves. mRNA abundance of CYP2C8 tended also to be higher (+44%; p = 0.053) after dexamethasone treatment. No significant treatment effects (p > 0.10) were observed for mRNA abundances of CYP3A4, CYP2E1, SULT2A1, UGT1A1 or cytochrome P450 reductase (CPR). In conclusion, an enhanced glucocorticoid status, induced by pharmacological amounts of dexamethasone, had differential and in part unexpected effects on NR and NR target systems in 5-day-old calves. Part of the unexpected responses may be due the immaturity of NR and NR receptor target systems.

Hexose-6-phosphate dehydrogenase modulates 11beta-hydroxysteroid dehydrogenase type 1-dependent metabolism of 7-keto- and 7beta-hydroxy-neurosteroids

BACKGROUND: The role of 11beta-hydroxysteroid dehydrogenase type 1 (11beta-HSD1) in the regulation of energy metabolism and immune system by locally reactivating glucocorticoids has been extensively studied. Experiments determining initial rates of enzyme activity revealed that 11beta-HSD1 can catalyze both the reductase and the dehydrogenase reaction in cell lysates, whereas it predominantly catalyzes the reduction of cortisone to cortisol in intact cells that also express hexose-6-phosphate dehydrogenase (H6PDH), which provides cofactor NADPH. Besides its role in glucocorticoid metabolism, there is evidence that 11beta-HSD1 is involved in the metabolism of 7-keto- and 7-hydroxy-steroids; however the impact of H6PDH on this alternative function of 11beta-HSD1 has not been assessed. METHODOLOGY: We investigated the 11beta-HSD1-dependent metabolism of the neurosteroids 7-keto-, 7alpha-hydroxy- and 7beta-hydroxy-dehydroepiandrosterone (DHEA) and 7-keto- and 7beta-hydroxy-pregnenolone, respectively, in the absence or presence of H6PDH in intact cells. 3D-structural modeling was applied to study the binding of ligands in 11beta-HSD1. PRINCIPAL FINDINGS: We demonstrated that 11beta-HSD1 functions in a reversible way and efficiently catalyzed the interconversion of these 7-keto- and 7-hydroxy-neurosteroids in intact cells. In the presence of H6PDH, 11beta-HSD1 predominantly converted 7-keto-DHEA and 7-ketopregnenolone into their corresponding 7beta-hydroxy metabolites, indicating a role for H6PDH and 11beta-HSD1 in the local generation of 7beta-hydroxy-neurosteroids. 3D-structural modeling offered an explanation for the preferred formation of 7beta-hydroxy-neurosteroids. CONCLUSIONS: Our results from experiments determining the steady state concentrations of glucocorticoids or 7-oxygenated neurosteroids suggested that the equilibrium between cortisone and cortisol and between 7-keto- and 7-hydroxy-neurosteroids is regulated by 11beta-HSD1 and greatly depends on the coexpression with H6PDH. Thus, the impact of H6PDH on 11beta-HSD1 activity has to be considered for understanding both glucocorticoid and neurosteroid action in different tissues.

Modulation of human CYP19A1 activity by mutant NADPH P450 oxidoreductase

Mutations in NADPH P450 oxidoreductase (POR) cause a broad spectrum of human disease with abnormalities in steroidogenesis. We have studied the impact of P450 reductase mutations on the activity of CYP19A1. POR supported CYP19A1 activity with a calculated Km of 126 nm for androstenedione and a Vmax of 1.7 pmol/min. Mutations R457H and V492E located in the FAD domain of POR that disrupt electron transfer caused a complete loss of CYP19A1 activity. The A287P mutation of POR decreased the activities of CYP17A1 by 60-80% but had normal CYP19A1 activity. Molecular modeling and protein docking studies suggested that A287P is involved in the interaction of POR:CYP17A1 but not in the POR:CYP19A1 interaction. Mutations C569Y and V608F in the NADPH binding domain of POR had 49 and 28% of activity of CYP19A1 compared with normal reductase and were more sensitive to the amount of NADPH available for supporting CYP19A1 activity. Substitution of NADH for NADPH had a higher impact on C569Y and V608F mutants of POR. Similar effects were obtained at low/high (5.5/8.5) pH, but using octanol to limit the flux of electrons from POR to CYP19A1 inhibited activity supported by all variants. High molar ratios of KCl also reduced the CYP19A1 supporting activities of C569Y and V608F mutants of POR to a greater extent compared to normal POR and A287P mutant. Because POR supports many P450s involved in steroidogenesis, bone formation, and drug metabolism, variations in the effects of POR mutations on specific enzyme activities may explain the broad clinical spectrum of POR deficiency.

Fructose-1,6-bisphosphate preserves intracellular glutathione and protects cortical neurons against oxidative stress

Fructose-1,6-bisphosphate (FBP), an endogenous intermediate of glycolysis, protects the brain against ischemia-reperfusion injury. The mechanisms of FBP protection after cerebral ischemia are not well understood. The current study was undertaken to determine whether FBP protects primary neurons against hypoxia and oxidative stress by preserving reduced glutathione (GSH). Cultures of pure cortical neurons were subjected to oxygen deprivation, a donor of nitric oxide and superoxide radicals (3-morpholinosydnonimine), an inhibitor of glutathione synthesis (L-buthionine-sulfoximine) or glutathione reductase (1,3-bis(2-chloroethyl)-1-nitrosourea) in the presence or absence of FBP (3.5 mM). Neuronal viability was determined using an 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide assay. FBP protected neurons against hypoxia-reoxygenation and oxidative stress under conditions of compromised GSH metabolism. The efficacy of FBP depended on duration of hypoxia and was associated with higher intracellular GSH concentration, an effect partly mediated via increased glutathione reductase activity.

Oxidative stress in lungs of mice infected with influenza A virus

As oxidative stress has been implicated in the pathogenesis of certain viral diseases we determined antioxidant and prooxidant parameters in lungs and bronchoalveolar lavage fluid (BALF) of mice infected with a lethal dose of influenza A/PR8/34 virus. Viral infection was characterized by massive infiltration of leukocytes, mainly polymorphonuclear leukocytes, into the alveolar space. The total number of BALF cells increased up to 8-fold (day 3 post-infection) and these cells appeared activated as judged by their increased rates of superoxide anion radical (O2-.) generation upon stimulation. Maximal rates of radical generation by BALF cells during the early stages of infection were 15- or 70-fold higher than those of cells from control animals when expressed per cell or total BALF cells, respectively. At the terminal stages of infection the total capacity of BALF cells to release O2-. declined to approximately 35-fold the control values. Infection also resulted in increased in vivo formation of hydrogen peroxide (H2O2) within the lungs at a time that coincided with the maximal capacity of BALF cells to release O2-.. Whereas pulmonary activities of glutathione peroxidase and reductase remained unaltered, levels of ascorbate in the cell-free BALF decreased significantly during the early stages of the infection and then returned to normal levels and above, late in infection. The oxidation state of the dehydroascorbic acid/ascorbate couple increased concomitantly with the decrease in ascorbate concentrations early in infection and remained elevated throughout the infection. As assessed by the prevention of peroxyl radical-induced loss of phycoerythrin fluorescence, the total antioxidant capacity present in lung tissue homogenate from terminally ill animals was not diminished when compared to that prepared from lungs of control mice. We conclude that although early stages of influenza infection are associated with the presence of oxidative stress in the lung tissue and alveolar fluid lining the epithelial cells, this stress does not appear to overwhelm local antioxidant defenses. The results therefore do not support a direct causative role of oxidative tissue damage in the pathogenesis of influenza virus infection.

The Gup1 homologue of Trypanosoma brucei is a GPI glycosylphosphatidylinositol remodelase

Glycosylphosphatidylinositol (GPI) lipids of Trypanosoma brucei undergo lipid remodelling, whereby longer fatty acids on the glycerol are replaced by myristate (C14:0). A similar process occurs on GPI proteins of Saccharomyces cerevisiae where Per1p first deacylates, Gup1p subsequently reacylates the anchor lipid, thus replacing a shorter fatty acid by C26:0. Heterologous expression of the GUP1 homologue of T. brucei in gup1Delta yeast cells partially normalizes the gup1Delta phenotype and restores the transfer of labelled fatty acids from Coenzyme A to lyso-GPI proteins in a newly developed microsomal assay. In this assay, the Gup1p from T. brucei (tbGup1p) strongly prefers C14:0 and C12:0 over C16:0 and C18:0, whereas yeast Gup1p strongly prefers C16:0 and C18:0. This acyl specificity of tbGup1p closely matches the reported specificity of the reacylation of free lyso-GPI lipids in microsomes of T. brucei. Depletion of tbGup1p in trypanosomes by RNAi drastically reduces the rate of myristate incorporation into the sn-2 position of lyso-GPI lipids. Thus, tbGup1p is involved in the addition of myristate to sn-2 during GPI remodelling in T. brucei and can account for the fatty acid specificity of this process. tbGup1p can act on GPI proteins as well as on GPI lipids.

DENITRIFICATION IN SOILS: FROM GENES TO ENVIRONMENTAL OUTCOMES

Denitrification is an important process of global nitrogen cycle as it removes reactive nitrogen from the biosphere, and acts as the primary source of nitrous oxide (N2O). This thesis seeks to gain better understanding of the biogeochemistry of denitrification by investigating the process from four different aspects: genetic basis, enzymatic kinetics, environmental interactions, and environmental consequences. Laboratory and field experiments were combined with modeling efforts to unravel the complexity of denitrification process under microbiological and environmental controls. Dynamics of denitrification products observed in laboratory experiments revealed an important role of constitutive denitrification enzymes, whose presence were further confirmed with quantitative analysis of functional genes encoding nitrite reductase and nitrous oxide reductase. A metabolic model of denitrification developed with explicit denitrification enzyme kinetics and representation of constitutive enzymes successfully reproduced the dynamics of N2O and N2 accumulation observed in the incubation experiments, revealing important regulatory effect of denitrification enzyme kinetics on the accumulation of denitrification products. Field studies demonstrated complex interaction of belowground N2O production, consumption and transport, resulting in two pulse pattern in the surface flux. Coupled soil gas diffusion/denitrification model showed great potential in simulating the dynamics of N2O below ground, with explicit representation of the activity of constitutive denitrification enzymes. A complete survey of environmental variables showed distinct regulation regimes on the denitrification activity from constitutive enzymes and new synthesized enzymes. Uncertainties in N2O estimation with current biogeochemical models may be reduced as accurate simulation of the dynamics of N2O in soil and surface fluxes is possible with a coupled diffusion/denitrification model that includes explicit representation of denitrification enzyme kinetics. In conclusion, denitrification is a complex ecological function regulated at cellular level. To assess the environmental consequences of denitrification and develop useful tools to mitigate N2O emissions require a comprehensive understanding of the regulatory network of denitrification with respect to microbial physiology and environmental interactions.

Role of CYP27A1 in progesterone metabolism in vitro and in vivo

In the kidney, progesterone is inactivated to 20alpha-dihydro-progesterone (20alpha-DH-progesterone) to protect the mineralocorticoid receptor from progesterone excess. In an attempt to clone the enzyme with 20alpha-hydroxysteroid activity using expression cloning in CHOP cells and a human kidney expression library, serendipitously cDNA encoding CYP27A1 was isolated. Overexpression of CYP27A1 in CHOP cells decreased progesterone conversion to 20alpha-DH-progesterone in a dose-dependent manner, an effect enhanced by cotransfection with adrenodoxin and adrenodoxin reductase. Incubation of CHOP cells with 27-hydroxycholesterol, a product of CYP27A1, increased the ratio of progesterone/20alpha-DH-progesterone in a concentration-dependent manner, indicating that the effect of CYP27A1 overexpression was mediated by 27-hydroxycholesterol. In order to analyze whether these observations are relevant in vivo, progesterone and 20alpha-DH-progesterone were measured by GC-MS in 24-h urine of CYP27A1 gene knock out (ko) mice and their control wild type (wt) and heterozygote (hz) littermates. In CYP27A1 ko mice, urinary progesterone concentrations were decreased, 20alpha-DH-progesterone increased and the progesterone/20alpha-DH-progesterone ratio decreased threefold (p<0.001). Thus, CYP27A1 modulates progesterone concentrations. The underlying mechanism is inhibition of 20alpha-hydroxysteroid dehydrogenase by 27-hydroxycholesterol. Key words: Progesterone, sterol 27-hydroxylase, 27-hydroxycholesterol, 20a-steroid dehydrogenase, 20a-DH-progesterone.

Variation in hepatic regulation of metabolism during the dry period and in early lactation in dairy cows

The purpose of this study was to investigate variations in hepatic regulation of metabolism during the dry period, after parturition, and in early lactation in dairy cows. For this evaluation, cows were divided into 2 groups based on the plasma concentration of beta-hydroxybutyric acid (BHBA) in wk 4 postpartum (PP; group HB, BHBA >0.75 mmol/L; group LB, BHBA <0.75 mmol/L, respectively). Liver biopsies were obtained from 28 cows at drying off (mean 59 +/- 8 d antepartum), on d 1, and in wk 4 and 14 PP. Blood samples were collected every 2 wk during this entire period. Liver samples were analyzed for mRNA abundance of genes related to carbohydrate metabolism (pyruvate carboxylase, PC; phosphoenolpyruvate carboxykinase, PEPCK; citrate synthase, CS), fatty acid biosynthesis (ATP citrate lyase, ACLY) and oxidation (acyl-CoA synthetase long-chain, ACSL; carnitine palmitoyltransferase 1A, CPT 1A; carnitine palmitoyltransferase 2, CPT 2; acyl-coenzyme A dehydrogenase very long chain, ACADVL), cholesterol biosynthesis (3-hydroxy-3-methylglutaryl-coenzyme A synthase 1, HMGCS1), ketogenesis (3-hydroxy-3-methylglutaryl-coenzyme A synthase 2, HMGCS2), and of genes encoding the transcription factors peroxisome proliferator-activated receptor alpha (PPARalpha), peroxisome proliferator-activated receptor gamma (PPARgamma), and sterol regulatory element binding factor 1 (SREBF1). Blood plasma was assayed for concentrations of glucose, BHBA, nonesterified fatty acids, cholesterol, triglycerides, insulin, insulin-like growth factor-I, and thyroid hormones. In both groups, plasma parameters followed a pattern usually observed in dairy cows. However, changes were moderate and the energy balance in cows turned positive in wk 7 PP for both groups. Additionally, the energy balance and milk yield were similar for both groups after parturition onwards. Significant group effects were found at drying off, when plasma concentrations of triglycerides were higher in LB than in HB, and in wk 4 PP, when plasma concentrations of glucose and IGF-I were lower in HB than in LB. Similarly, moderate changes in mRNA expression of hepatic genes between the different time points were observed, although HB cows showed more adaptive performance than LB cows based on changes in mRNA expression of PEPCKc, PEPCKm, CS, CPT 1A, CPT 2, and PPARalpha. Part of the variation measured in this study was explained by parity. Significant Spearman rank correlation coefficients between the variables were not similar at each time point and were not similar between the groups at each time point, suggesting that metabolic regulation differs between cows. In conclusion, metabolic regulation in dairy cows is a dynamic system, and differs obviously between cows at different metabolic stages related to parturition.