Role of DNA methylation in the tissue-specific expression of the CYP17A1 gene for steroidogenesis in rodents

The CYP17A1 gene is the qualitative regulator of steroidogenesis. Depending on the presence or absence of CYP17 activities mineralocorticoids, glucocorticoids or adrenal androgens are produced. The expression of the CYP17A1 gene is tissue as well as species-specific. In contrast to humans, adrenals of rodents do not express the CYP17A1 gene and have therefore no P450c17 enzyme for cortisol production, but produce corticosterone. DNA methylation is involved in the tissue-specific silencing of the CYP17A1 gene in human placental JEG-3 cells. We investigated the role of DNA methylation for the tissue-specific expression of the CYP17A1 gene in rodents. Rats treated with the methyltransferase inhibitor 5-aza-deoxycytidine excreted the cortisol metabolite tetrahydrocortisol in their urine suggesting that treatment induced CYP17 expression and 17alpha-hydroxylase activity through demethylation. Accordingly, bisulfite modification experiments identified a methylated CpG island in the CYP17 promoter in DNA extracted from rat adrenals but not from testes. Both methyltransferase and histone deacetylase inhibitors induced the expression of the CYP17A1 gene in mouse adrenocortical Y1 cells which normally do not express CYP17, indicating that the expression of the mouse CYP17A1 gene is epigenetically controlled. The role of DNA methylation for CYP17 expression was further underlined by the finding that a reporter construct driven by the mouse -1041 bp CYP17 promoter was active in Y1 cells, thus excluding the lack of essential transcription factors for CYP17 expression in these adrenal cells.

Differential expression and localization of lipid transporters in the bovine mammary gland during the pregnancy-lactation cycle

The transport of lipids across mammary gland epithelial cells (MEC) determines milk lipid content and composition. We investigated the expression of lipid transporters and their regulators in comparison to blood metabolites during lactation and dry period (DP) in dairy cows. Repeated mammary gland biopsies and blood samples were taken from 10 animals at 7 stages of the pregnancy-lactation cycle. Expression levels of the specific mRNAs were determined by quantitative reverse transcription-PCR, whereas ABCA1 was localized by immunohistochemistry. Blood serum metabolites were determined by common enzymatic chemistries. Elevated mRNA profiles of ABCA1 and ABCA7 were found during DP as compared with lactation and were inversely associated with blood cholesterol levels. Elevated levels of ABCG2, NPC1, SREBP1, SREBP2, LXR alpha, and PPAR gamma were found postpartum, whereas ABCG1 did not differ between the functional stages of the mammary gland. The ABCA1 protein was localized in MEC and showed differential activity between DP and lactation suggesting a role of ABCA1 in the removal of excess cellular cholesterol from MEC during the DP. The expression profiles of ABCA7 and NPC1 may reflect a role of these transporters in the clearance of apoptotic cells and the intracellular redistribution of cholesterol, respectively. Regulation of lipid transporters in the mammary gland is partially associated with transcription factors that control lipid homeostasis.

Cardiac c-kit+AT2+ cell population is increased in response to ischemic injury and supports cardiomyocyte performance

The expression pattern of angiotensin AT2 receptors with predominance during fetal life and upregulation under pathological conditions during tissue injury/repair process suggests that AT2 receptors may exert an important action in injury/repair adaptive mechanisms. Less is known about AT2 receptors in acute ischemia-induced cardiac injury. We aimed here to elucidate the role of AT2 receptors after acute myocardial infarction. Double immunofluorescence staining showed that cardiac AT2 receptors were mainly detected in clusters of small c-kit+ cells accumulating in peri-infarct zone and c-kit+AT2+ cells increased in response to acute cardiac injury. Further, we isolated cardiac c-kit+AT2+ cell population by modified magnetic activated cell sorting and fluorescence activated cell sorting. These cardiac c-kit+AT2+ cells, represented approximately 0.19% of total cardiac cells in infarcted heart, were characterized by upregulated transcription factors implicated in cardiogenic differentiation (Gata-4, Notch-2, Nkx-2.5) and genes required for self-renewal (Tbx-3, c-Myc, Akt). When adult cardiomyocytes and cardiac c-kit+AT2+ cells isolated from infarcted rat hearts were cocultured, AT2 receptor stimulation in vitro inhibited apoptosis of these cocultured cardiomyocytes. Moreover, in vivo AT2 receptor stimulation led to an increased c-kit+AT2+ cell population in the infarcted myocardium and reduced apoptosis of cardiomyocytes in rats with acute myocardial infarction. These data suggest that cardiac c-kit+AT2+ cell population exists and increases after acute ischemic injury. AT2 receptor activation supports performance of cardiomyocytes, thus contributing to cardioprotection via cardiac c-kit+AT2+ cell population.

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.

Gene array of primary human osteoblasts exposed to enamel matrix derivative in combination with a natural bone mineral

OBJECTIVES The application of an enamel matrix derivative (EMD) for regenerative periodontal surgery has been shown to promote formation of new cementum, periodontal ligament, and alveolar bone. In intrabony defects with a complicated anatomy, the combination of EMD with various bone grafting materials has resulted in additional clinical improvements, but the initial cellular response of osteoblasts coming in contact with these particles have not yet been fully elucidated. The objective of the present study was to evaluate the in vitro effects of EMD combined with a natural bone mineral (NBM) on a wide variety of genes, cytokines, and transcription factors and extracellular matrix proteins on primary human osteoblasts. MATERIAL AND METHODS Primary human osteoblasts were seeded on NBM particles pre-coated with versus without EMD and analyzed for gene differences using a human osteogenesis gene super-array (Applied Biosystems). Osteoblast-related genes include those transcribed during bone mineralization, ossification, bone metabolism, cell growth and differentiation, as well as gene products representing extracellular matrix molecules, transcription factors, and cell adhesion molecules. RESULTS EMD promoted gene expression of various osteoblast differentiation markers including a number of collagen types and isoforms, SMAD intracellular proteins, osteopontin, cadherin, alkaline phosphatase, and bone sialoprotein. EMD also upregulated a variety of growth factors including bone morphogenetic proteins, vascular endothelial growth factors, insulin-like growth factor, transforming growth factor, and their associated receptor proteins. CONCLUSION The results from the present study demonstrate that EMD is capable of activating a wide variety of genes, growth factors, and cytokines when pre-coated onto NBM particles. CLINICAL RELEVANCE The described in vitro effects of EMD on human primary osteoblasts provide further biologic support for the clinical application of a combination of EMD with NBM particles in periodontal and oral regenerative surgery.

Skeletal muscle mitochondria in the elderly: effects of physical fitness and exercise training.

Context: Sarcopenia is thought to be associated with mitochondrial (M) loss. It is unclear whether the decrease in M content is consequent to aging per se or to decreased physical activity. Objectives: To examine the influence of fitness on M content and function, and to assess whether exercise could improve M function in older adults. Design and subjects: Three distinct studies were conducted: 1) a cross-sectional observation comparing M content and fitness in a large heterogeneous cohort of older adults; 2) a case-control study comparing chronically endurance-trained older adults (A) and sedentary (S) subjects matched for age and gender; 3) a 4-month exercise intervention in S. Setting: University-based clinical research center Outcomes: M volume density (Mv) was assessed by electron microscopy from vastus lateralis biopsies, electron transport chain proteins (ETC) by western blotting, mRNAs for transcription factors involved in M biogenesis by qRT-PCR and in-vivo oxidative capacity (ATPmax) by (31)P-MR spectroscopy. Peak oxygen uptake (VO2peak) was measured by GXT. Results: VO2peak was strongly correlated with Mv in eighty 60-80 yo adults. Comparison of A vs. S revealed differences in Mv, ATPmax and some ETC complexes. Finally, exercise intervention confirmed that S are able to recover Mv, ATPmax and specific transcription factors. Conclusions: These data suggest that 1) aging per se is not the primary culprit leading to M dysfunction, 2) an aerobic exercise program, even at an older age, can ameliorate the loss in skeletal muscle M content and may prevent aging muscle comorbidities and 3) the improvement of M function is all about content.

Hepatic gene expression involved in glucose and lipid metabolism in transition cows: effects of fat mobilization during early lactation in relation to milk performance and metabolic changes.

Insufficient feed intake during early lactation results in elevated body fat mobilization to meet energy demands for milk production. Hepatic energy metabolism is involved by increasing endogenous glucose production and hepatic glucose output for milk synthesis and by adaptation of postcalving fuel oxidation. Given that cows differ in their degree of fat mobilization around parturition, indicated by variable total liver fat concentration (LFC), the study investigated the influence of peripartum fat mobilization on hepatic gene expression involved in gluconeogenesis, fatty acid oxidation, ketogenesis, and cholesterol synthesis, as well as transcriptional factors referring to energy metabolism. German Holstein cows were grouped according to mean total LFC on d 1, 14, and 28 after parturition as low [<200mg of total fat/g of dry matter (DM); n=10], medium (200-300 mg of total fat/g of DM; n=10), and high (>300 mg of total fat/g of DM; n=7), indicating fat mobilization during early lactation. Cows were fed total mixed rations ad libitum and held under equal conditions. Liver biopsies were taken at d 56 and 15 before and d 1, 14, 28, and 49 after parturition to measure mRNA abundances of pyruvate carboxylase (PC); phosphoenolpyruvate carboxykinase; glucose-6-phosphatase; propionyl-coenzyme A (CoA) carboxylase α; carnitine palmitoyl-transferase 1A (CPT1A); acyl-CoA synthetase, long chain 1 (ASCL1); acyl-CoA dehydrogenase, very long chain; 3-hydroxy-3-methylglutaryl-CoA synthase 1 and 2; sterol regulatory element-binding factor 1; and peroxisome proliferator-activated factor α. Total LFC postpartum differed greatly among cows, and the mRNA abundance of most enzymes and transcription factors changed with time during the experimental period. Abundance of PC mRNA increased at parturition to a greater extent in high- and medium-LFC groups than in the low-LFC group. Significant LFC × time interactions for ACSL1 and CPT1A during the experimental period indicated variable gene expression depending on LFC after parturition. Correlations between hepatic gene expression and performance data and plasma concentrations of metabolites and hormones showed time-specific relations during the transition period. Elevated body fat mobilization during early lactation affected gene expression involved in gluconeogenesis to a greater extent than gene expression involved in lipid metabolism, indicating the dependence of hepatic glucose metabolism on hepatic lipid status and fat mobilization during early lactation.

Differential regulation of human 3β-hydroxysteroid dehydrogenase type 2 for steroid hormone biosynthesis by starvation and cyclic AMP stimulation: studies in the human adrenal NCI-H295R cell model

Human steroid biosynthesis depends on a specifically regulated cascade of enzymes including 3β-hydroxysteroid dehydrogenases (HSD3Bs). Type 2 HSD3B catalyzes the conversion of pregnenolone, 17α-hydroxypregnenolone and dehydroepiandrosterone to progesterone, 17α-hydroxyprogesterone and androstenedione in the human adrenal cortex and the gonads but the exact regulation of this enzyme is unknown. Therefore, specific downregulation of HSD3B2 at adrenarche around age 6-8 years and characteristic upregulation of HSD3B2 in the ovaries of women suffering from the polycystic ovary syndrome remain unexplained prompting us to study the regulation of HSD3B2 in adrenal NCI-H295R cells. Our studies confirm that the HSD3B2 promoter is regulated by transcription factors GATA, Nur77 and SF1/LRH1 in concert and that the NBRE/Nur77 site is crucial for hormonal stimulation with cAMP. In fact, these three transcription factors together were able to transactivate the HSD3B2 promoter in placental JEG3 cells which normally do not express HSD3B2. By contrast, epigenetic mechanisms such as methylation and acetylation seem not involved in controlling HSD3B2 expression. Cyclic AMP was found to exert differential effects on HSD3B2 when comparing short (acute) versus long-term (chronic) stimulation. Short cAMP stimulation inhibited HSD3B2 activity directly possibly due to regulation at co-factor or substrate level or posttranslational modification of the protein. Long cAMP stimulation attenuated HSD3B2 inhibition and increased HSD3B2 expression through transcriptional regulation. Although PKA and MAPK pathways are obvious candidates for possibly transmitting the cAMP signal to HSD3B2, our studies using PKA and MEK1/2 inhibitors revealed no such downstream signaling of cAMP. However, both signaling pathways were clearly regulating HSD3B2 expression.

"Hot spot" in the PROP1 gene responsible for combined pituitary hormone deficiency.

As pituitary function depends on the integrity of the hypothalamic-pituitary axis, any defect in the development and organogenesis of this gland may account for a form of combined pituitary hormone deficiency (CPHD). Although pit-1 was 1 of the first factors identified as a cause of CPHD in mice, many other homeodomain and transcription factors have been characterized as being involved in different developmental stages of pituitary gland development, such as prophet of pit-1 (prop-1), P-Lim, ETS-1, and Brn 4. The aims of the present study were first to screen families and patients suffering from different forms of CPHD for PROP1 gene alterations, and second to define possible hot spots and the frequency of the different gene alterations found. Of 73 subjects (36 families) analyzed, we found 35 patients, belonging to 18 unrelated families, with CPHD caused by a PROP1 gene defect. The PROP1 gene alterations included 3 missense mutations, 2 frameshift mutations, and 1 splice site mutation. The 2 reported frameshift mutations could be caused by any 2-bp GA or AG deletion at either the 148-GGA-GGG-153 or 295-CGA-GAG-AGT-303 position. As any combination of a GA or AG deletion yields the same sequencing data, the frameshift mutations were called 149delGA and 296delGA, respectively. All but 1 mutation were located in the PROP1 gene encoding the homeodomain. Importantly, 3 tandem repeats of the dinucleotides GA at location 296-302 in the PROP1 gene represent a hot spot for CPHD. In conclusion, the PROP1 gene seems to be a major candidate gene for CPHD; however, further studies are needed to evaluate other genetic defects involved in pituitary development.

An Sp1/Sp3 binding polymorphism confers methylation protection.

Hundreds of genes show aberrant DNA hypermethylation in cancer, yet little is known about the causes of this hypermethylation. We identified RIL as a frequent methylation target in cancer. In search for factors that influence RIL hypermethylation, we found a 12-bp polymorphic sequence around its transcription start site that creates a long allele. Pyrosequencing of homozygous tumors revealed a 2.1-fold higher methylation for the short alleles (P<0.001). Bisulfite sequencing of cancers heterozygous for RIL showed that the short alleles are 3.1-fold more methylated than the long (P<0.001). The comparison of expression levels between unmethylated long and short EBV-transformed cell lines showed no difference in expression in vivo. Electrophorectic mobility shift assay showed that the inserted region of the long allele binds Sp1 and Sp3 transcription factors, a binding that is absent in the short allele. Transient transfection of RIL allele-specific transgenes showed no effects of the additional Sp1 site on transcription early on. However, stable transfection of methylation-seeded constructs showed gradually decreasing transcription levels from the short allele with eventual spreading of de novo methylation. In contrast, the long allele showed stable levels of expression over time as measured by luciferase and approximately 2-3-fold lower levels of methylation by bisulfite sequencing (P<0.001), suggesting that the polymorphic Sp1 site protects against time-dependent silencing. Our finding demonstrates that, in some genes, hypermethylation in cancer is dictated by protein-DNA interactions at the promoters and provides a novel mechanism by which genetic polymorphisms can influence an epigenetic state.

Bifurcation and singularity analysis of a molecular network for the induction of long-term memory.

Withdrawal reflexes of the mollusk Aplysia exhibit sensitization, a simple form of long-term memory (LTM). Sensitization is due, in part, to long-term facilitation (LTF) of sensorimotor neuron synapses. LTF is induced by the modulatory actions of serotonin (5-HT). Pettigrew et al. developed a computational model of the nonlinear intracellular signaling and gene network that underlies the induction of 5-HT-induced LTF. The model simulated empirical observations that repeated applications of 5-HT induce persistent activation of protein kinase A (PKA) and that this persistent activation requires a suprathreshold exposure of 5-HT. This study extends the analysis of the Pettigrew model by applying bifurcation analysis, singularity theory, and numerical simulation. Using singularity theory, classification diagrams of parameter space were constructed, identifying regions with qualitatively different steady-state behaviors. The graphical representation of these regions illustrates the robustness of these regions to changes in model parameters. Because persistent protein kinase A (PKA) activity correlates with Aplysia LTM, the analysis focuses on a positive feedback loop in the model that tends to maintain PKA activity. In this loop, PKA phosphorylates a transcription factor (TF-1), thereby increasing the expression of an ubiquitin hydrolase (Ap-Uch). Ap-Uch then acts to increase PKA activity, closing the loop. This positive feedback loop manifests multiple, coexisting steady states, or multiplicity, which provides a mechanism for a bistable switch in PKA activity. After the removal of 5-HT, the PKA activity either returns to its basal level (reversible switch) or remains at a high level (irreversible switch). Such an irreversible switch might be a mechanism that contributes to the persistence of LTM. The classification diagrams also identify parameters and processes that might be manipulated, perhaps pharmacologically, to enhance the induction of memory. Rational drug design, to affect complex processes such as memory formation, can benefit from this type of analysis.

Simulation of Drosophila circadian oscillations, mutations, and light responses by a model with VRI, PDP-1, and CLK.

A model of Drosophila circadian rhythm generation was developed to represent feedback loops based on transcriptional regulation of per, Clk (dclock), Pdp-1, and vri (vrille). The model postulates that histone acetylation kinetics make transcriptional activation a nonlinear function of [CLK]. Such a nonlinearity is essential to simulate robust circadian oscillations of transcription in our model and in previous models. Simulations suggest that two positive feedback loops involving Clk are not essential for oscillations, because oscillations of [PER] were preserved when Clk, vri, or Pdp-1 expression was fixed. However, eliminating positive feedback by fixing vri expression altered the oscillation period. Eliminating the negative feedback loop in which PER represses per expression abolished oscillations. Simulations of per or Clk null mutations, of per overexpression, and of vri, Clk, or Pdp-1 heterozygous null mutations altered model behavior in ways similar to experimental data. The model simulated a photic phase-response curve resembling experimental curves, and oscillations entrained to simulated light-dark cycles. Temperature compensation of oscillation period could be simulated if temperature elevation slowed PER nuclear entry or PER phosphorylation. The model makes experimental predictions, some of which could be tested in transgenic Drosophila.

Single-stranded DNA-binding proteins regulate the abundance of LIM domain and LIM domain-binding proteins.

The LIM domain-binding protein Ldb1 is an essential cofactor of LIM-homeodomain (LIM-HD) and LIM-only (LMO) proteins in development. The stoichiometry of Ldb1, LIM-HD, and LMO proteins is tightly controlled in the cell and is likely a critical determinant of their biological actions. Single-stranded DNA-binding proteins (SSBPs) were recently shown to interact with Ldb1 and are also important in developmental programs. We establish here that two mammalian SSBPs, SSBP2 and SSBP3, contribute to an erythroid DNA-binding complex that contains the transcription factors Tal1 and GATA-1, the LIM domain protein Lmo2, and Ldb1 and binds a bipartite E-box-GATA DNA sequence motif. In addition, SSBP2 was found to augment transcription of the Protein 4.2 (P4.2) gene, a direct target of the E-box-GATA-binding complex, in an Ldb1-dependent manner and to increase endogenous Ldb1 and Lmo2 protein levels, E-box-GATA DNA-binding activity, and P4.2 and beta-globin expression in erythroid progenitors. Finally, SSBP2 was demonstrated to inhibit Ldb1 and Lmo2 interaction with the E3 ubiquitin ligase RLIM, prevent RLIM-mediated Ldb1 ubiquitination, and protect Ldb1 and Lmo2 from proteasomal degradation. These results define a novel biochemical function for SSBPs in regulating the abundance of LIM domain and LIM domain-binding proteins.

Conversion of myoblasts to physiologically active neuronal phenotype.

Repressor element 1 (RE1)-silencing transcription factor (REST)/neuron-restrictive silencer factor (NRSF) can repress several terminal neuronal differentiation genes by binding to a specific DNA sequence (RE1/neuron-restrictive silencer element [NRSE]) present in their regulatory regions. REST-VP16 binds to the same RE1/NRSE, but activates these REST/NRSF target genes. However, it is unclear whether REST-VP16 expression is sufficient to cause formation of functional neurons either from neural stem cells or from heterologous stem cells. Here we show that the expression of REST-VP16 in myoblasts grown under muscle differentiation conditions blocked entry into the muscle differentiation pathway, countered endogenous REST/NRSF-dependent repression, activated the REST/NRSF target genes, and, surprisingly, activated other neuronal differentiation genes and converted the myoblasts to a physiologically active neuronal phenotype. Furthermore, in vitro differentiated neurons produced by REST-VP16-expressing myoblasts, when injected into mouse brain, survived, incorporated into the normal brain, and did not form tumors. This is the first instance in which myoblasts were converted to a neuronal phenotype. Our results suggest that direct activation of REST/NRSF target genes with a single transgene, REST-VP16, is sufficient to activate other terminal neuronal differentiation genes and to override the muscle differentiation pathways, and they suggest that this approach provides an efficient way of triggering neuronal differentiation in myoblasts and possibly other stem cells.

Extended kalman filter for estimation of parameters in nonlinear state-space models of biochemical networks.

It is system dynamics that determines the function of cells, tissues and organisms. To develop mathematical models and estimate their parameters are an essential issue for studying dynamic behaviors of biological systems which include metabolic networks, genetic regulatory networks and signal transduction pathways, under perturbation of external stimuli. In general, biological dynamic systems are partially observed. Therefore, a natural way to model dynamic biological systems is to employ nonlinear state-space equations. Although statistical methods for parameter estimation of linear models in biological dynamic systems have been developed intensively in the recent years, the estimation of both states and parameters of nonlinear dynamic systems remains a challenging task. In this report, we apply extended Kalman Filter (EKF) to the estimation of both states and parameters of nonlinear state-space models. To evaluate the performance of the EKF for parameter estimation, we apply the EKF to a simulation dataset and two real datasets: JAK-STAT signal transduction pathway and Ras/Raf/MEK/ERK signaling transduction pathways datasets. The preliminary results show that EKF can accurately estimate the parameters and predict states in nonlinear state-space equations for modeling dynamic biochemical networks.