Ligand-Free Estrogen Receptor α (ESR1) as Master Regulator for the Expression of CYP3A4 and Other Cytochrome P450 Enzymes in the Human Liver

Quantitative analysis of the UDP-glucuronosyltransferase transcriptome in human tissues

UDP-glucuronosyltransferases (UGTs) are phase II drug metabolizing enzymes that play important roles in the detoxification of endogenous and exogenous substrates. The 22 human UGTs belong to four families (UGT1, UGT2, UGT3, and UGT8) and differ in their expression, substrate specificity, UDP-sugar preference, and physiological functions. Differential expression/activity of the UGTs contributes to interperson variability in drug responses and toxicity, hormone homeostasis, and disease/cancer risks. However, in normal tissues, the tissue-specific expression profiles and transcriptional regulation of the UGTs are still not fully understood. In this study, we comprehensively analyzed the transcriptome of 22 UGTs in 54 human tissues/regions using RNAseq data from GTEx. We then validated the findings in the liver and small intestine samples using real-time PCR. Our results showed large interindividual variability across tissues in the expression of each UGT and the overall composition of UGT pools, consisting of different UGTs and their splice isoforms. Our results also revealed coexpression of the UGTs, Cytochrome P450s, and many transcription factors in the liver, suggesting potential coregulation or functional coordination. Our results provide the groundwork for future studies to detail further the regulation of the expression and activity of the UGTs.

Transcription Factors and ncRNAs Associated with CYP3A Expression in Human Liver and Small Intestine Assessed with Weighted Gene Co-Expression Network Analysis

CYP3A4, CYP3A5, and CYP3A7, which are located in a multigene locus (CYP3A), play crucial roles in drug metabolism. To understand the highly variable hepatic expression of CYP3As, regulatory network analyses have focused on transcription factors (TFs). Since long non-coding RNAs (lncRNAs) likely contribute to such networks, we assessed the regulatory effects of both TFs and lncRNAs on CYP3A expression in the human liver and small intestine, main organs of CYP3A expression. Using weighted gene co-expression network analysis (WGCNA) of GTEx v8 RNA expression data and multiple stepwise regression analysis, we constructed TF-lncRNA-CYP3A co-expression networks. Multiple lncRNAs and TFs displayed robust associations with CYP3A expression that differed between liver and small intestines (LINC02499, HNF4A-AS1, AC027682.6, LOC102724153, and RP11-503C24.6), indicating that lncRNAs contribute to variance in CYP3A expression in both organs. Of these, HNF4A-AS1 had been experimentally demonstrated to affect CYP3A expression. Incorporating ncRNAs into CYP3A expression regulatory network revealed additional candidate TFs associated with CYP3A expression. These results serve as a guide for experimental studies on lncRNA-TF regulation of CYP3A expression in the liver and small intestines.

Regulatory variants in a novel distal enhancer regulate the expression of CYP3A4 and CYP3A5

The cytochrome P450 3As (CYP3As) are abundantly expressed in the liver and metabolize many commonly prescribed medications. Their expression is highly variable between individuals with little known genetic cause. Despite extensive investigation, cis-acting genetic elements that control the expression of the CYP3As remain uncharacterized. Using chromatin conformation capture (4C assays), we detected reciprocal interaction between a distal regulatory region (DRR) and the CYP3A4 promoter. The DRR colocalizes with a variety of enhancer marks and was found to promote transcription in reporter assays. CRISPR-mediated deletion of the DRR decreased expression of CYP3A4, CYP3A5, and CYP3A7, supporting its role as a shared enhancer regulating the expression of three CYP3A genes. Using reporter gene assays, we identified two single-nucleotide polymorphisms (rs115025140 and rs776744/rs776742) that increased DRR-driven luciferase reporter expression. In a liver cohort (n = 246), rs115025140 was associated with increased expression of CYP3A4 mRNA (1.8-fold) and protein (1.6-fold) and rs776744/rs776742 was associated with 1.39-fold increased expression of CYP3A5 mRNA. The rs115025140 is unique to the African population and in a clinical cohort of African Americans taking statins for lipid control rs115025140 carriers showed a trend toward reduced statin-mediated lipid reduction. In addition, using a published cohort of Chinese patients who underwent renal transplantation taking tacrolimus, rs776744/rs776742 carriers were associated with reduced tacrolimus concentration after adjusting for CYP3A5*3. Our results elucidate a complex regulatory network controlling expression of three CYP3A genes and identify two novel regulatory variants with potential clinical relevance for predicting CYP3A4 and CYP3A5 expression.

Genome-wide microRNA profiles identify miR-107 as a top miRNA associating with expression of the CYP3As and other drug metabolizing cytochrome P450 enzymes in the liver

Cytochrome P450 (CYP) drug metabolizing enzymes are responsible for the metabolism of over 70% of currently used medications with the CYP3A family being the most important CYP enzymes in the liver. Large inter-person variability in expression/activity of the CYP3As greatly affects drug exposure and treatment outcomes, yet the cause of such variability remains elusive. Micro-RNAs (miRNAs) are small noncoding RNAs that negatively regulate gene expression and are involved in diverse cellular processes including metabolism of xenobiotics and therapeutic outcomes. Target prediction and in vitro functional assays have linked several miRNAs to the control of CYP3A4 expression. Yet, their co-expression with CYP3As in the liver remain unclear. In this study, we used genome-wide miRNA profiling in liver samples to identify miRNAs associated with the expression of the CYP3As. We identified and validated both miR-107 and miR-1260 as strongly associated with the expression of CYP3A4, CYP3A5, and CYP3A43. Moreover, we found associations between miR-107 and nine transcription factors (TFs) that regulate CYP3A expression, with estrogen receptor alpha (ESR1) having the largest effect size. Including ESR1 and the other TFs in the regression model either diminished or abolished the associations between miR-107 and the CYP3As, indicating that the role of miR-107 in CYP3A expression may be indirect and occur through these key TFs. Indeed, testing the other nine CYPs previously shown to be regulated by ESR1 identified similar miR-107 associations that were dependent on the exclusion of ESR1 and other key TFs in the regression model. In addition, we found significant differences in miRNA expression profiles in liver samples between race and sex. Together, our results identify miR-107 as a potential epigenetic regulator that is strongly associated with the expression of many CYPs, likely via impacting the CYP regulatory network controlled by ESR1 and other key TFs. Therefore, both genetic and epigenetic factors that alter the expression of miR-107 may have a broad influence on drug metabolism.

Transcriptional Regulation of Carboxylesterase 1 in Human Liver: Role of the Nuclear Receptor Subfamily 1 Group H Member 3 and Its Splice Isoforms

Carboxylesterase 1 (CES1) is the predominant carboxylesterase in the human liver, involved in metabolism of both xenobiotics and endogenous substrates. Genetic or epigenetic factors that alter CES1 activity or expression are associated with changes in drug response, lipid, and glucose homeostasis. However, the transcriptional regulation of CES1 in the human liver remains uncertain. By applying both the random forest and Sobol's Sensitivity Indices (SSI) to analyze existing liver RNA expression microarray data (GSE9588), we identified nuclear receptor subfamily 1 group H member 3 (NR1H3) liver X receptor (LXR)α as a key factor regulating constitutive CES1 expression. This model prediction was validated using small interfering RNA (siRNA) knockdown and CRISPR-mediated transcriptional activation of NR1H3 in Huh7 and HepG2 cells. We found that NR1H3's activation of CES1 is splice isoform-specific, namely that increased expression of the NR1H3-211 isoform increased CES1 expression whereas NR1H3-201 did not. Also, in human liver samples, expression of NR1H3-211 and CES1 are correlated, whereas NR1H3-201 and CES1 are not. This trend also occurs during differentiation of induced pluripotent stem cells (iPSCs) to hepatocytes, where only expression of the NR1H3-211 isoform parallels expression of CES1 Moreover, we found that treatment with the NR1H3 agonist T0901317 in HepG2 cells had no effect on CES1 expression. Overall, our results demonstrate a key role of NR1H3 in maintaining the constitutive expression of CES1 in the human liver. Furthermore, our results support that the effect of NR1H3 is splice isoform-specific and appears to be ligand independent. SIGNIFICANCE STATEMENT: Despite the central role of carboxylesterase 1 (CES1) in metabolism of numerous medications, little is known about its transcriptional regulation. This study identifies nuclear receptor subfamily 1 group H member 3 as a key regulator of constitutive CES1 expression and therefore is a potential target for future studies to understand interperson variabilities in CES1 activity and drug metabolism.

Studies on CYP3A activity during the menstrual cycle as measured by urinary 6β-hydroxycortisol/cortisol

The 6β-OH-cortisol/cortisol ratio (6β-OHC/C) in urine is an endogenous marker of drug-metabolizing enzyme cytochrome P450 3A (CYP3A). The primary aim of this single center, prospective, non-interventional cohort study, was to investigate the variability of 6β-OHC/C during the menstrual cycle. In addition, possible associations between the CYP3A activity and sex hormones, gut microbiota metabolite trimethylamine-N-Oxide (TMAO) and microRNA-27b, respectively, were investigated. Serum and urinary samples from healthy, regularly menstruating women followed for two menstrual cycles were analyzed. Twenty-six complete menstrual cycles including follicular, ovulatory, and luteal phase were defined based on hormone analyses in serum. 6β-OHC/C were analyzed in urine and sex hormones, TMAO and miRNA-27b were analyzed in serum at the same time points. 6β-OHC/C did not vary between the follicular, ovulatory, or luteal phases. There was a difference in the relative miRNA-27b expression between the follicular and ovulatory phase (p = .03). A significant association was found between 6β-OHC/C and progesterone during the follicular (p = .005) and ovulatory (p = .01) phases (n = 26 for each phase). In addition, a significant association was found between the ratio and TMAO during the ovulatory (p = .02) and luteal (p = .002) phases. 6β-OHC/C and gut microbiota TMAO were significantly associated (p = .003) when evaluating all values, for all phases (n = 78). Interestingly, the finding of an association between 6β-OHC/C in urine and levels of TMAO in serum suggest that gut microbiota may affect CYP3A activity.

Trans-eQTLs of the CYP3A4 and CYP3A5 associated with tacrolimus trough blood concentration in Chinese renal transplant patients

This study aimed to systematically investigate trans-eQTLs of CYP3A4 and CYP3A5 affecting tacrolimus trough blood concentrations in Chinese renal transplant patients. We used Plink v1.90 to perform data quality control and linear regression analysis on GTEx v8 data. SNPs with p-value < 0.05 were selected and the GTEx eQTL Calculator was used to further prioritize the eQTLs of CYP3A4 and CYP3A5 in the liver and small intestine. The eQTLs with a p-value < 5 × 10^-5 and MAF≥ 0.05 in the CHB population were selected as candidate eQTLs. The genotyping of candidate eQTLs was performed using high-resolution melting (HRM) assays and Sanger DNA sequencing. This study included 845 Chinese renal transplant patients who received tacrolimus as an immunosuppressive agent. Association between 103 candidate eQTLs and log-transformed tacrolimus concentration/dose ratio (log (C₀/D)) in this cohort was conducted using the SNPassoc package of R software. In the end, a total of 75,632 liver eQTLs of CYP3A4, 69,558 liver eQTLs of CYP3A5, 48,596 small intestine eQTLs of CYP3A4 and 28,616 small intestine eQTLs of CYP3A5 were obtained using the GTEx v8 eQTL Calculator. Of the 103 candidate eQTLs, rs75727207, rs181294422 and rs28522676 were significantly associated with tacrolimus log(C₀/D) in different genetic models. We discovered a substantial number of novel eQTLs of CYP3A4 and CYP3A5 in liver and small intestine, also found that rs75727207, rs181294422 and rs28522676 may affect tacrolimus trough blood concentrations in Chinese renal transplant patients.

Variation in Expression of Cytochrome P450 3A Isoforms and Toxicological Effects: Endo- and Exogenous Substances as Regulatory Factors and Substrates

The CYP3A subfamily, which includes isoforms CYP3A4, CYP3A5, and CYP3A7 in humans, plays important roles in the metabolism of various endogenous and exogenous substances. Gene and protein expression of CYP3A4, CYP3A5, and CYP3A7 show large inter-individual differences, which are caused by many endogenous and exogenous factors. Inter-individual differences can cause negative outcomes, such as adverse drug events and disease development. Therefore, it is important to understand the variations in CYP3A expression caused by endo- and exogenous factors, as well as the variation in the metabolism and kinetics of endo- and exogenous substrates. In this review, we summarize the factors regulating CYP3A expression, such as bile acids, hormones, microRNA, inflammatory cytokines, drugs, environmental chemicals, and dietary factors. In addition, variations in CYP3A expression under pathological conditions, such as coronavirus disease 2019 and liver diseases, are described as examples of the physiological effects of endogenous factors. We also summarize endogenous and exogenous substrates metabolized by CYP3A isoforms, such as cholesterol, bile acids, hormones, arachidonic acid, vitamin D, and drugs. The relationship between the changes in the kinetics of these substrates and the toxicological effects in our bodies are discussed. The usefulness of these substrates and metabolites as endogenous biomarkers for CYP3A activity is also discussed. Notably, we focused on discrimination between CYP3A4, CYP3A5, and CYP3A7 to understand inter-individual differences in CYP3A expression and function.

Beyond the X Factor: Relevance of Sex Hormones in NAFLD Pathophysiology

Non-alcoholic fatty liver disease (NAFLD) is a major health issue worldwide, being frequently associated with obesity, unbalanced dietary regimens, and reduced physical activity. Despite their greater adiposity and reduced physical activity, women show a lower risk of developing NAFLD in comparison to men, likely a consequence of a sex-specific regulation of liver metabolism. In the liver, sex differences in the uptake, synthesis, oxidation, deposition, and mobilization of lipids, as well as in the regulation of inflammation, are associated with differences in NAFLD prevalence and progression between men and women. Given the major role of sex hormones in driving hepatic sexual dimorphism, this review will focus on the role of sex hormones and their signaling in the regulation of hepatic metabolism and in the molecular mechanisms triggering NAFLD development and progression.

Cytochrome P450 3A4 (CYP3A4) protein quantification using capillary western blot technology and total protein normalization

The western blot (WB) is the predominate method for protein quantification, frequently used in pharmacological and toxicological studies. To control for technical variation, WB signals are normalized through immunodetection of an internal standard "house-keeping" gene or total protein quantification via staining of the same blot or a duplicate, sister blot. Increasing evidence suggests that house-keeping genes are subject to change after drug treatment or under disease states, causing protein quantification errors in WB. Recent advances in automated capillary-based WB technologies enable measurement of the protein of interest, internal standards, and total protein in a single capillary. Using this approach, we quantified cytochrome P450 3A4 (CYP3A4) across 179 liver samples and compared normalization by both β-actin and total protein to determine which better functions as an internal standard. CYP3A4 is responsible for metabolism of a wide array of xenobiotics and is known to exhibit large inter-person variation, making it a good candidate to evaluate protein quantification. We observed significant differences in β-actin protein levels between liver samples (~20-fold) and found better correlation between CYP3A4 protein and mRNA using total protein normalization than β-actin, indicating total protein normalization to be less error prone for estimation of CYP3A4. Furthermore, by using total protein normalization, we confirmed significant association between CYP3A4 protein expression and the functional CYP3A4 variant CYP3A4*22, which contains two linked SNPs rs35599367 and rs62471956. Our results indicate that the automatic capillary WB instrument combined with total protein normalization provides a high throughput and robust approach for protein quantification.

Impact of Estrogen Withdrawal and Replacement in Female Mice along the Intestinal Tract. Comparison of E2 Replacement with the Effect of a Mixture of Low Dose Pollutants

Postmenopausal women represent a vulnerable population towards endocrine disruptors due to hormonal deficit. We previously demonstrated that chronic exposure of ovariectomized C57Bl6/J mice fed a high-fat, high-sucrose diet to a low-dose mixture of chemicals with one dioxin, one polychlorobiphenyl, one phthalate, and bisphenol A triggered metabolic alterations in the liver but the intestine was not explored. Yet, the gastrointestinal tract is the main route by which pollutants enter the body. In the present study, we investigated the metabolic consequences of ovarian withdrawal and E2 replacement on the various gut segments along with investigating the impact of the mixture of pollutants. We showed that genes encoding estrogen receptors (Esr1, Gper1 not Esr2), xenobiotic processing genes (e.g., Cyp3a11, Cyp2b10), and genes related to gut homeostasis in the jejunum (e.g., Cd36, Got2, Mmp7) and to bile acid biosynthesis in the gut (e.g., Fgf15, Slc10a2) and liver (e.g., Abcb11, Slc10a1) were under estrogen regulation. Exposure to pollutants mimicked some of the effects of E2 replacement, particularly in the ileum (e.g., Esr1, Nr1c1) suggesting that the mixture had estrogen-mimetic activities. The present findings have important implications for the understanding of estrogen-dependent metabolic alterations with regards to situations of loss of estrogens as observed after menopause.

Characterization of a Novel CYP1A2 Knockout Rat Model Constructed by CRISPR/Cas9

CYP1A2, as one of the most important cytochrome P450 isoforms, is involved in the biotransformation of many important endogenous and exogenous substances. CYP1A2 also plays an important role in the development of many diseases because it is involved in the biotransformation of precancerous substances and poisons. Although the generation of Cyp1a2 knockout (KO) mouse model has been reported, there are still no relevant rat models for the study of CYP1A2-mediated pharmacokinetics and diseases. In this report, CYP1A2 KO rat model was established successfully by CRISPR/Cas9 without any detectable off-target effect. Compared with wild-type rats, this model showed a loss of CYP1A2 protein expression in the liver. The results of pharmacokinetics in vivo and incubation in vitro of specific substrates of CYP1A2 confirmed the lack of function of CYP1A2 in KO rats. In further studies of potential compensatory effects, we found that CYP1A1 was significantly upregulated, and CYP2E1, CYP3A2, and liver X receptor β were downregulated in KO rats. In addition, CYP1A2 KO rats exhibited a significant increase in serum cholesterol and free testosterone accompanied by mild liver damage and lipid deposition, suggesting that CYP1A2 deficiency affects lipid metabolism and liver function to a certain extent. In summary, we successfully constructed the CYP1A2 KO rat model, which provides a useful tool for studying the metabolic function and physiologic function of CYP1A2. SIGNIFICANCE STATEMENT: Human CYP1A2 not only metabolizes clinical drugs and pollutants but also mediates the biotransformation of endogenous substances and plays an important role in the development of many diseases. However, there are no relevant CYP1A2 rat models for the research of pharmacokinetics and diseases. This study successfully established CYP1A2 knockout rat model by using CRISPR/Cas9. This rat model provides a powerful tool to study the function of CYP1A2 in drug metabolism and diseases.

Co-expression of drug metabolizing cytochrome P450 enzymes and estrogen receptor alpha (ESR1) in human liver: racial differences and the regulatory role of ESR1

OBJECTIVES

The function and expression of cytochrome P450 (CYP) drug metabolizing enzymes is highly variable, greatly affecting drug exposure, and therapeutic outcomes. The expression of these enzymes is known to be controlled by many transcription factors (TFs), including ligand-free estrogen receptor alpha (ESR1, in the absence of estrogen). However, the relationship between the expression of ESR1, other TFs, and CYP enzymes in human liver is still unclear.

METHODS

Using real-time PCR, we quantified the mRNA levels of 12 CYP enzymes and nine TFs in 246 human liver samples from European American (EA, n = 133) and African American (AA, n = 113) donors.

RESULTS

Our results showed higher expression levels of ESR1 and six CYP enzymes in EA than in AA. Partial least square regression analysis showed that ESR1 is the top-ranking TF associating with the expression of eight CYP enzymes, six of which showed racial difference in expression. Conversely, four CYP enzymes without racial difference in expression did not have ESR1 as a top-ranking TF. These results indicate that ESR1 may contribute to variation in CYP enzyme expression between these two ancestral backgrounds.

CONCLUSIONS

These results are consistent with our previous study showing ESR1 as a master regulator for the expression of several CYP enzymes. Therefore, factors affecting ESR1 expression may have broad influence on drug metabolism through altered expression of CYP enzymes.

Hepcidin-induced reduction in iron content and PGC-1β expression negatively regulates osteoclast differentiation to play a protective role in postmenopausal osteoporosis

As a necessary trace element, iron is involved in many physiological processes. Clinical and basic studies have found that disturbances in iron metabolism, especially iron overload, might lead to bone loss and even be involved in postmenopausal osteoporosis. Hepcidin is a key regulator of iron homeostasis. However, the exact role of hepcidin in bone metabolism and the underlying mechanism remain unknown. In this study, we found that in postmenopausal osteoporosis cohort, the concentration of hepcidin in the serum was significantly reduced and positively correlated with bone mineral density. Ovariectomized (OVX) mice were then used to construct an osteoporosis model. Hepcidin overexpression in these mice significantly improved bone mass and rescued the phenotype of bone loss. Additionally, overexpression of hepcidin in OVX mice greatly reduced the number and differentiation of osteoclasts in vivo and in vitro. This study found that overexpression of hepcidin significantly inhibited ROS production, mitochondrial biogenesis, and PGC-1β expression. These data showed that hepcidin protected osteoporosis by reducing iron levels in bone tissue, and in conjunction with PGC-1β, reduced ROS production and the number of mitochondria, thus inhibiting osteoclast differentiation and bone absorption. Hepcidin could provide new targets for the clinical treatment of postmenopausal osteoporosis.

DNA methylation and expression of estrogen receptor alpha in fathead minnows exposed to 17α-ethynylestradiol

The gene expression response thought to underlie the negative apical effects resulting from estrogen exposure have been thoroughly described in fish. Although epigenetics are believed to play a critical role translating environmental exposures into the development of adverse apical effects, they remain poorly characterized in fish species. This study investigated alterations of DNA methylation of estrogen receptor alpha (esr1) in brain and liver tissues from 8 to 10 month old male fathead minnows (Pimephales promelas) after a 2d exposure to either 2.5 ng/L or 10 ng/L 17α-ethynylestradiol (EE2). Changes in the patterns of methylation were evaluated using targeted deep sequencing of bisulfite treated DNA in the 5' region of esr1. Methylation and gene expression were assessed at 2d of exposure and after a 7 and 14d depuration period. After 2d EE2 exposure, males exhibited significant demethylation in the 5' upstream region of esr1 in liver tissue, which was inversely correlated to gene expression. This methylation pattern reflected what was seen in females. No gene body methylation (GBM) was observed for liver of exposed males. Differential methylation was observed for a single upstream CpG site in the liver after the 14d depuration. A less pronounced methylation response was observed in the upstream region in brain tissue, however, several CpGs were necessarily excluded from the analysis. In contrast to the liver, a significant GBM response was observed across the entire gene body, which was sustained until at least 7d post-exposure. No differential expression was observed in the brain, limiting functional interpretation of methylation changes. The identification of EE2-dependent changes in methylation levels strongly suggests the importance of epigenetic mechanisms as a mediator of the organismal response to environmental exposures and the need for further characterization of the epigenome. Further, differential methylation following depuration indicates estrogenic effects persist well after the active exposure, which has implications for the risk posed by repeated exposures..

ESR1 ChIP-Seq Identifies Distinct Ligand-Free ESR1 Genomic Binding Sites in Human Hepatocytes and Liver Tissue

The estrogen receptor alpha (ESR1) is an important gene transcriptional regulator, known to mediate the effects of estrogen. Canonically, ESR1 is activated by its ligand estrogen. However, the role of unliganded ESR1 in transcriptional regulation has been gaining attention. We have recently shown that ligand-free ESR1 is a key regulator of several cytochrome P450 (CYP) genes in the liver, however ligand-free ESR1 has not been characterized genome-wide in the human liver. To address this, ESR1 ChIP-Seq was conducted in human liver samples and in hepatocytes with or without 17beta-estradiol (E2) treatment. We identified both ligand-dependent and ligand-independent binding sites throughout the genome. These two ESR1 binding categories showed different genomic localization, pathway enrichment, and cofactor colocalization, indicating different ESR1 regulatory function depending on ligand availability. By analyzing existing ESR1 data from additional human cell lines, we uncovered a potential ligand-independent ESR1 activity, namely its co-enrichment with the zinc finger protein 143 (ZNF143). Furthermore, we identified ESR1 binding sites near many gene loci related to drug therapy, including the CYPs. Overall, this study shows distinct ligand-free and ligand-bound ESR1 chromatin binding profiles in the liver and suggests the potential broad influence of ESR1 in drug metabolism and drug therapy.

Licorice Ameliorates Cisplatin-Induced Hepatotoxicity Through Antiapoptosis, Antioxidative Stress, Anti-Inflammation, and Acceleration of Metabolism

Cisplatin (CP) is one of the most effective antitumor drugs in the clinic, but has serious adverse reactions, and its hepatotoxicity has not been fully investigated. Licorice (GC), a traditional herbal medicine, has been commonly used as a detoxifier for poisons and drugs, and may be an effective drug for CP-induced hepatotoxicity. However, its mechanism and the effector molecules remain ambiguous. Therefore, in this study, a network pharmacology and proteomics-based approach was established, and a panoramic view of the detoxification of GC on CP-induced hepatotoxicity was provided. The experimental results indicated that GC can recover functional indices and pathological liver injury, inhibit hepatocyte apoptosis, upregulate B-cell lymphoma/leukemia 2 (Bcl-2) and superoxide dismutase (SOD) levels, and downregulate cellular tumor antigen p53 (p53), caspase-3, malondialdehyde high mobility group protein B1 (HMGB1), tumor necrosis factor alpha (TNF-α), and interleukin 1β (IL-1β) levels. Proteomics indicated that GC regulates phosphatidylcholine translocator ABCB1 (ABCB1B), canalicular multispecific organic anion transporter 1 (ABCC2), cytochrome P450 4A2 (CYP4A2), cytochrome P450 1A1 (CYP1A1), cytochrome P450 1A2 (CYP1A2), estrogen receptor (ESR1), and DNA topoisomerase 2-alpha (TOP2A), inhibits oxidative stress, apoptosis, and inflammatory responses, and accelerates drug metabolism. In this study, we provide the investigation of the efficacy of GC against CP-induced hepatotoxicity, and offer a promising alternative for the clinic.

Uncovering the Mechanism of the Effects of Pien-Tze-Huang on Liver Cancer Using Network Pharmacology and Molecular Docking

Pien-Tze-Huang (PTH) has a long history in the treatment of liver cancer. However, its molecular mechanism of action remains unclear. TCMSP and TCM were used to collect the active ingredients. Bioactive compounds targets were predicted by reverse pharmacophore models. The antiliver cancer targets of PTH were selected by gene comparison of liver cancer in the GEO database. Molecular docking was used to verify the binding activity of the targets and the active ingredients. The DAVID was used to analyze the gene function and signal pathway. A model was built with Cytoscape. The core genes were obtained by PPI network. We screened the 4 main medicinal ingredients of PTH to obtain 16 active ingredient, 190 potential targets, and 6 core genes. We found that active small molecules exert anticancer effects by multiple pathways. The core genes were involved in multiple biological processes. We also found that eight chemical components play a greater role in inhibiting liver cancer. PTH achieves the effect of inhibiting liver cancer through the synergistic effect of multiple components, multiple targets, and multiple pathways. This study provides a potential scientific basis for further elucidating the molecular mechanism of action of PTH against liver cancer.

Prognostic Score-based Clinical Factors and Metabolism-related Biomarkers for Predicting the Progression of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a common malignant tumor representing more than 90% of primary liver cancer. This study aimed to identify metabolism-related biomarkers with prognostic value by developing the novel prognostic score (PS) model. Transcriptomic profiles derived from TCGA and EBIArray databases were analyzed to identify differentially expressed genes (DEGs) in HCC tumor samples compared with normal samples. The overlapped genes between DEGs and metabolism-related genes (crucial genes) were screened and functionally analyzed. A novel PS model was constructed to identify optimal signature genes. Cox regression analysis was performed to identify independent clinical factors related to prognosis. Nomogram model was constructed to estimate the predictability of clinical factors. Finally, protein expression of crucial genes was explored in different cancer tissues and cell types from the Human Protein Atlas (HPA). We screened a total of 305 overlapped genes (differentially expressed metabolism-related genes). These genes were mainly involved in "oxidation reduction," "steroid hormone biosynthesis," "fatty acid metabolic process," and "linoleic acid metabolism." Furthermore, we screened ten optimal DEGs (CYP2C9, CYP3A4, and TKT, among others) by using the PS model. Two clinical factors of pathologic stage (P < .001, HR: 1.512 [1.219-1.875]) and PS status (P <.001, HR: 2.259 [1.522-3.354]) were independent prognostic predictors by cox regression analysis. Nomogram model showed a high predicted probability of overall survival time, and the AUC value was 0.837. The expression status of 7 proteins was frequently altered in normal or differential tumor tissues, such as liver cancer and stomach cancer samples.We have identified several metabolism-related biomarkers for prognosis prediction of HCC based on the PS model. Two clinical factors were independent prognostic predictors of pathologic stage and PS status (high/low risk). The prognosis prediction model described in this study is a useful and stable method for novel biomarker identification.

Motifs enable communication efficiency and fault-tolerance in transcriptional networks

Analysis of the topology of transcriptional regulatory networks (TRNs) is an effective way to study the regulatory interactions between the transcription factors (TFs) and the target genes. TRNs are characterized by the abundance of motifs such as feed forward loops (FFLs), which contribute to their structural and functional properties. In this paper, we focus on the role of motifs (specifically, FFLs) in signal propagation in TRNs and the organization of the TRN topology with FFLs as building blocks. To this end, we classify nodes participating in FFLs (termed motif central nodes) into three distinct roles (namely, roles A, B and C), and contrast them with TRN nodes having high connectivity on the basis of their potential for information dissemination, using metrics such as network efficiency, path enumeration, epidemic models and standard graph centrality measures. We also present the notion of a three tier architecture and how it can help study the structural properties of TRN based on connectivity and clustering tendency of motif central nodes. Finally, we motivate the potential implication of the structural properties of motif centrality in design of efficient protocols of information routing in communication networks as well as their functional properties in global regulation and stress response to study specific disease conditions and identification of drug targets.