Comprehensive brain tissue metabolomics and biological network technology to decipher the mechanism of hydrogen-rich water on Radiation-induced cognitive impairment in rats

Associated Data

Data Availability Statement



Hydrogen-rich water (HRW) has been shown to prevent cognitive impairment caused by ionizing radiation. This study aimed to investigate the pharmacological effects and mechanisms of HRW on ionizing radiation by coupling the brain metabolomics and biological target network methods.

Methods and results

HRW significantly improves the cognitive impairment in rats exposed to ionizing radiation. Based on metabolomics and biological network results, we identified 54 differential metabolites and 93 target genes. The KEGG pathway indicates that glutathione metabolism, ascorbic acid and aldehyde acid metabolism, pentose and glucuronic acid interconversion, and glycerophospholipid metabolism play important roles in ionizing radiation therapy.


Our study has systematically elucidated the molecular mechanism of HRW against ionizing radiation, which can be mediated by modulating targets, pathways and metabolite levels. This provides a new perspective for identifying the underlying pharmacological mechanism of HRW.

Keywords: Hydrogen-rich water, Ionizing radiation, Brain tissue metabolomics, Biological network


Ionizing radiation can produce a variety of free radicals, damage nucleic acids, proteins, carbohydrates and lipid compounds, thereby destroying the structure and function of cells [, ]. Approximately 60–70% of the damage is caused by hydroxyl radicals []. The brain contains high levels of polyunsaturated fatty acids (PUFAs) and relatively low levels of antioxidant enzymes, making it highly susceptible to oxidative damage []. Radiation can affect cognitive functions such as learning and memory because it can damage neurogenesis, promote brain cell apoptosis, and cause neurodegenerative changes []. Radiation may cause long-term cognitive impairment, especially in children whose nervous system is still developing [, ].

Hydrogen has great potential in the clinical treatment of various neurological diseases []. Previous studies have shown that hydrogen treatment can effectively prevent and alleviate neonatal brain injury and cognitive dysfunction, especially cerebral ischemia-reperfusion injury []. Hydrogen can cross the blood-brain barrier and quickly cross cell membranes to interact with toxic ROS. Therefore, hydrogen has great potential in the clinical treatment of various diseases of the nervous system.

Metabolomics, as a new omics technique, study the cellular processes involving metabolites, intermediates and products of cell metabolism, and has find applications in drug metabolism, biomarkers, and environmental toxicants. Metabolomics elucidates the end points of biological processes, reflecting the physiological and pathological state through metabolic levels []. Therefore, metabolomics is a useful tool to investigate pharmacodynamics and mechanism of action of radioprotectants which reflect the end state of diseases and treatments []. However, the changes in endogenous metabolites are unknown, and using metabolomics alone as a tool may not reveal the effectiveness of radiation protection agents in the treatment of radiation-induced diseases.

Computer virtual computing based on high throughput omics data analysis and network database search may reveal complex network signal relationships among drugs, targets and diseases []. For radioprotective agents, their pharmacological effects are characterized by multi-targets and multi-levels. The biological target network has the advantage to analyze and explain the action mechanisms from a holistic perspective, promoting the in-depth investigation []. Combining the two key technologies of computer simulation and metabionomics, the molecular mechanism of agents’ action on diseases can be explained from the metabolic level and genetic levels. For example, Liu et al. linked the biological targets network with metabolomics and found that amino acid metabolism is a potential target for Xiaoyaosan anti-depressant [].

To the best of our knowledge, no studies have used metabolomics combined with biological network approaches to study the effect of hydrogen on radiation-induced cognitive impairment. This study used metabolomics combined with a biological target network to investigate the molecular biological mechanism of HRW on radiation-induced cognitive impairment in rats.


Materials and methods

Materials and reagents

Portable HRW machine (China Institute For Radiation Protection, Shanxi Zhongfu Nuclear Instruments Co.,Ltd.), Linear accelerator (Elekta Synergyt, Crawley, UK), ultra-high-performance liquid chromatography (Dionex Ultimate 3000 UHPLC, Thermo Scientific, USA) coupled to high-resolution mass spectrometry (Q-Exactive Plus, Thermo Scientific, USA). Water, methanol, acetonitrile and formic acid were purchased from CNW Technology Co., Ltd. (Dü sseldorf, Germany). L-2-chlorophenylalanine was obtained from Shanghai Hengchuang Biotechnology Co., Ltd. (Shanghai, China). All chemicals and solvents are of analytical or HPLC grade.

Animals and ethic statement

Male Sprague-Dawley rats weighing 140–170 g and, aging 6 weeks were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd., China. Animals were caged in a humidity and temperature-controlled room, with a 12/12 light-dark cycle. Water and food were provided to animals with free access. The animal study protocols were approved by the Animal Protection and Ethics Committee of the China Institute for Radiation Protection. All laboratory operations were performed in accordance with the Regulations of China Institute for Radiation Protection on the administration of laboratory animals. The study was carried out in compliance with the ARRIVE guidelines.

Establishment of the ionizing radiation model

Twenty-four rats were randomized into four groups (6 in each group), as follows: HRW group received HRW(hydrogen concentration: 0.8–0.9 ppm) without radiation treatment; Radiotherapy group received 30 Gy whole brain irradiation of purified water; Radiation -HRW group received 30 Gy whole brain irradiation with HRW; Control group received purified water without radiation treatment. Rats were dosed with purified water or HRW gavage at 20ml/kg per day 10 min before irradiation and 30 days after irradiation. Immobilize the animal directly with clamps before preparing the animal for radiotherapy with deionized water using a portable HRW machine. Then, at room temperature, a 6-MeV electron beam was emitted at a rate of 3 Gy/min by a linear accelerator to irradiate the whole brain. On the 36th day after irradiation, rats were sacrificed by intraperitoneal injection of pentobarbital and blood was taken from the abdominal aorta, the brains of the animals were collected.

Morris Water Maze Test

On the 30th day after irradiation, the rats in each group were subjected to a cognitive test in the form of a position navigation test and a spatial search test. In the site navigation task, rats acclimated themselves with the test environment by swimming for 2 min in the circular pool without a platform the day before the formal test. Rats were then trained for five consecutive days. To start training, rats were placed in the pool facing the pool wall, the swimming route of the rats were tracked, the escape latency (the time from entering the water to finding platform–all limbs on the platform) was measured. Rats were allowed to rest on the platform for a few seconds before returning them to their cages. The escape latency was recorded as 120 s if the rat did not find the platform within 2 min and were guided to the platform. The arithmetic mean of the escape latencies was calculated four times per day. The position navigation test was followed by a spatial probe test in which the platform was withdrawn and the rats were placed at any entry point in the pool on day six. The swimming distance of rats were recorded and analyzed within 1 min. Record the passage time they cross the original platform, the time the rat retained in the original platform quadrant, their swimming distance in the platform quadrant, and the number of times the rat crossed the platform.

Brain tissue samples collection and preparation

After the brain tissue was collected, 30 mg brain tissue sample were accurately weighed and transferred to a 1.5 mL Eppendorf tube. Two small steel spheres and 20 µL of the internal standard (2-chloro-l-phenylalanine in methanol, 0.3 mg/mL) and 400 µL of methanol/water (4/1, v/v) extraction solvent were added to the sample tube. The samples were stored at -80 ℃ for 2 min, ground at 60 HZ for 2 min, sonicated at ambient temperature (25 ℃ to 28 ℃) for 10 min and stored at -20 ℃ for 30 min Centrifuge the extract at 13,000 rpm, 4 ℃ for 15 min. Dry 300 mL of the supernatant in a brown glass vial in a freeze concentration centrifugal dryer. Then 200 µL mixture of methanol and water (1/4, vol/vol) were added to each sample, the sample was vortexed for 30 s and, placed at 4 ℃ for 2 min. Samples were centrifuged at 13,000 rpm, 4 ℃ for 5 min. The supernatant (150 µL) was collected from each tube using a crystal syringe and was filtered through a 0.22 μm microfilter and transferred to an LC vial. Before LC -MS analysis, the vials were stored at -80 ℃ until. QC samples were prepared through mixing aliquots of all samples into a pooled sample.

Metabolomics analysis


UHPLC-MS parameters

Samples were analyzed by ultra-high-performance liquid chromatography coupled to high-resolution mass spectrometry in both ESI positive and ESI negative ion modes. In both positive and negative ionization modes, an ACQUITY UPLC HSS T3 column (100 mm×2.1 mm, 1.8 μm) were employed. Water and acetonitrile, both containing 0.1% formic acid, were used as mobile phases A and B, respectively. The linear gradient was: 0 min, 5% B; 2 min, 5% B; 4 min, 25%B; 8 min, 50% B; 10 min, 80% B; 14 min, 100% B; 15 min, 100% B; 15.1 min, 5%B and 16 min, 5%B. The column temperature was 45 ℃ and the flow rate was 0.35 mL/min. The injection volume was 2 µL. All samples were kept at 4℃ during analysis.

Data acquisition was performed in full scan mode (m/z ranges from 100 to 1000) combined with MSE mode, consisting of alternating acquisition of 2 independent scans with different collision energies (CE) during the run. The mass spectrometry parameters were as follows: a low-energy scan (CE 4 eV) and a high-energy scan (CE ramp 20-45 eV) to fragment the ions. Argon (99.999%) was used as collision-induced dissociation gas; scan rate was 0.2 s/scan; capillary voltages were − 3 kV (negative mode) and 3.8 kV (positive mode); capillary and heater temperatures were 320℃ and 300℃. The auxiliary gas volume flow was 8 arb, sheath gas volume flow was 35 arb. The stepped normalized collision energy (NCE) was set to 10, 20, and 40 eV respectively. The QCs were injected periodically (every 10 samples) throughout the analysis to provide a data set from which repeatability can be assessed.


Data processing

The raw data were exported using Xcalibur workstation and imported into Composite Discoverer 2.0 software to obtain matched and aligned peak data. The parameters are as follows: the scanning range is 100–1000 m/z. Ass deviation is 5 ppm. The retention time is 0.05 min. The signal-to-noise ratio threshold is 1.5. The peak data containing retention time, molecular formula, accurate molecular weight, and peak area information were imported into Excel for peak area normalization. Finally, the peak area normalization data were imported into SIMCA-P 14.0 for partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA). VIP > 1 in the S-plot and P < 0.05 in the independent sample t-test were used to screen for the most contributing differential metabolites. Differential metabolites were matched and screened using Metlin, HMDB, Pubchem, KEGG, m/zcloud, and others online databases. Finally, the identified differential metabolites were directed to MetaboAnalyst 5.0 for pathway enrichment analysis.

Construction of biological network of HRW for ionizing radiation therapy

The GeneCards database ( and the DisGeNET database ( were used to identify potential genes associated with ionizing radiation. Metabolites were introduced into the Cytoscape 3.9.2 plug-in MetScape for metabolic enzyme analysis. Finally, Cytoscape 3.9.2 was used to construct a “gene-metabolite” regulatory network for ionizing radiation therapy.

Biological function and pathway analysis

To elucidate the gene functions and their role in signal transduction, DAVID ( was used to evaluate the KEGG and GO enrichment profiles of HRW against ionizing radiation (P < 0.05).

ROC analysis

ROC curve analysis and calculation of AUC were performed using the Omicstudio tool ( to screen for potential therapeutic biomarkers and assesse the efficacy of central metabolites in the HRW treatment of ionizing radiation. These metabolites were combined by logistic regression analysis. P < 0.05 with a 95% confidence interval was considered statistically significant.

Statistical analysis

All data were presented as mean ± standard deviation. SPSS 22.0 and GraphPad Prism 7.0 were used for statistical analysis. The independent samples t test and one-way analysis of variance were used for comparison between two groups and multiple groups. P < 0.05 was considered statistically significant.


Effect of HRW on learning and memory in rats with ionizing radiation

Morris water maze test was used to assess the functional consequences of radiation exposure on the brain and whether HRW produced any neurocognitive benefits. During the training period of the first five days, the radiation group had significantly longer escape latency than the control group, the HRW group and the radiation -HRW treatment group (P < 0.05, Fig. 1A.). The spatial search test was performed from day six. The radiation-HRW group had significantly longer retention time in the original platform quadrant, greater platform crossings, and longer swimming distance than the radiation group (P < 0.05 Fig. 1B, C, D). These results indicated that brain radiation damaged the spatial memory and spatial search ability of rats, and HRW had a protective effect on radiation-induced cognitive impairment.

Effects of hydrogen-rich water (HRW) on cognitive function by Morris water maze test. (A) The escape latencies of different groups during the five-day training period. *P < 0.05 vs. radiation-HRW, one-way ANOVA. (B) The retention time in the original platform quadrant. (C) The swimming distance in the original platform quadrant. (D) The number of platform crossings. *P < 0.05


HRW regulated metabolic disorders in ionizing radiation rats

PLS-DA was used for metabolic profiling to reduce the dimensionality of complex data obtained from brain tissue samples from the control and radiation groups. The model was clearly separated from the control group, indicating that this radiation had a significant effect on the endogenous metabolic spectrum of rat brain tissue (Fig. 2A). The PLS-DA model was then validated with an arrangement test (Fig. 2B). All the permutation values of R2 and Q2 were smaller than the original values, indicating that the established model was stable and reliable. We further analyzed the differences between groups using PLS-DA (Fig. 2E). The control group, the radiation group, and radiation-HRW group could be clearly distinguished, indicating that HRW was effective for radiation group rats. The QC samples were collected together in PLS-DA to demonstrate the robustness and applicability of the system and method.

Multivariate statistical analysis of LC-MS data. (A) PLS-DA score chart of control and radiation group. (B) PLS-DA model verification diagram. (C) OPLS-DA score chart of control rats and model rats. (D) S-plot of control rats and model rats. (E) PLS-DA scores of brain tissue from different groups


Brain tissue sample data from control and model groups were analyzed with OPLS- DA to identify radiation-related differential metabolites. The combination of VIP > 1 (in S-plot) and P < 0.05 (in T-test) revealed a significant contribution of metabolites to homeostasis (Fig. 2 C and and2D).2D). Compared with the control group, a total of 80 metabolites including LysoPE(18:1(11Z)/0:0), pantothenic acid, gamma-glutamylleucine, 1-phosphate ribose, stearoylethanolamide, etc., changed significantly in the model group. After the intervention of HRW, such changes were significantly adjusted (Table 1), indicating that HRW dramatically reversed the metabolic disorders of ionizing radiation rats.

Table 1

Statistics of differential metabolites recovered by HRW

No Metabolites m/z Retention time (min) Ion mode Formula Radiation vs. Control Radiation-HRW vs. Radiation
1 Bisdiphosphoinositol tetrakisphosphate 704.8074441 0.74 Neg C6H20O24P8
2 PC(20:1(11Z)/24:0) 880.7169796 0.74 Neg C52H102NO8P
3 Choline 104.1073217 0.78 Pos C5H13NO
4 alpha-methyl butyric Acid 120.1020555 0.78 Pos C5H10O2
5 Ribothymidine 276.1187823 0.8 Pos C10H14N2O6
6 2-N,6-N-Bis(2,3-dihydroxybenzoyl)-L-lysine 399.1192105 0.81 Neg C20H22N2O8
7 Ribose 1-phosphate 229.0114577 0.82 Neg C5H11O8P
8 Norophthalmic acid 274.1046339 0.82 Neg C10H17N3O6
9 Caproic acid 134.1174357 0.84 Pos C6H12O2
10 Drazoxolon 238.0377343 0.86 Pos C10H8ClN3O2
11  S-Cysteinosuccinic acid 236.0230962 0.87 Neg C7H11NO6S
12 Succinylproline 233.1129797 0.88 Pos C9H13NO5
13 Glutathionate(1-) 308.090742 0.88 Pos C10H17N3O6S
14 gamma-Glutamylglutamic acid 275.0886761 0.89 Neg C10H16N2O7
15 Maleic hydrazide 113.0348169 0.89 Pos C4H4N2O2
16 (4 S,5 S)-4,5-dihydroxy-2,6-dioxohexanoic acid 351.0571659 0.89 Neg C6H8O6
17 Itaconic acid 175.0237955 0.89 Neg C5H6O4
18 Aristolodione 306.0765161 0.89 Neg C18H13NO4
19 Isobutyric acid 175.098179 0.91 Neg C4H8O2
20 D-Lysine 191.1036811 0.94 Neg C6H14N2O2
21 Propaphos 343.0535487 0.95 Pos C13H21O4PS
22 ADP-ribose 604.0705895 1.05 Neg C15H23N5O14P2
23 a-L-threo-4-Hex-4-enopyranuronosyl-D-galacturonic acid 351.0571152 1.17 Neg C12H16O12
24 D-Glucuronic acid 175.0238612 1.17 Neg C6H10O7
25 Glutathione 306.0763654 1.23 Neg C10H17N3O6S
26 Letrozole 308.0905898 1.24 Pos C17H11N5
27 3,4,5-trihydroxy-6-(2-oxo-1,3-diphenylpropoxy)oxane-2-carboxylic acid 385.1282845 1.24 Pos C21H22O8
28 Pseudouridine 245.0766217 1.27 Pos C9H12N2O6
29 (R)-3-hydroxybutyrylcarnitine 248.1490203 1.28 Pos C11H21NO5
30 [2-hydroxy-5-(2-hydroxy-3-methoxy-3-oxopropyl)phenyl]oxidanesulfonic acid 310.0586886 1.38 Pos C10H12O8S
31 2-(alpha-D-Galactosyl)-sn-glycerol 3-phosphate 143.0702116 1.39 Pos C7H12O4
32 Cysteine-Homocysteine disulfide 299.037998 1.39 Neg C7H14N2O4S2
33 Isoguanosine 284.0986875 1.42 Pos C10H13N5O5
34 Guanosine 282.0844741 1.44 Neg C10H13N5O5
35  L-Ascorbic acid-2-glucoside 383.0846362 1.44 Neg C12H18O11
36 1,4-beta-D-Glucan 537.1685822 1.45 Pos C18H32O18
37 Allopurinol-1-ribonucleoside 269.0878499 1.45 Pos C10H12N4O5
38 Arabinosylhypoxanthine 267.0736398 1.45 Neg C10H12N4O5
39 Allopurinol 137.0457685 1.45 Pos C5H4N4O
40 5-fluoroindole-2-carboxylic Acid 357.0689298 1.45 Neg C9H6FNO2
41 Elephantorrhizol 303.0502085 1.45 Neg C15H14O8
42 Saphenamycin 803.2361443 1.45 Neg C23H18N2O5
43 Pyrocatechol glucuronide 571.1311958 1.46 Neg C12H14O8
44 Isosorbide Dinitrate 281.0272222 1.49 Neg C6H8N2O8
45 Carotamine 310.0588562 1.8 Pos C14H13N3O3
46 Glutamyltyrosine 311.1235281 2.82 Pos C14H18N2O6
47 Pantothenic Acid 220.1178617 3.09 Pos C9H17NO5
48  S-(4,5-Dihydro-2-methyl-3-furanyl) ethanethioate 176.0739269 4.07 Pos C7H10O2S
49 5β-cholanic Acid-3α, 12α-diol N-(2-sulphoethyl)-amide 538.2619276 4.20 Pos C26H45NO6S
50 gamma-Glutamylleucine 261.1444168 4.55 Pos C11H20N2O5
51 LysoPE(0:0/18:1(11Z)) 480.3084895 11.46 Pos C23H46NO7P
52 LysoPE(18:1(11Z)/0:0) 480.3083533 11.61 Pos C23H46NO7P
53 Stearoylethanolamide 328.3207028 14.15 Pos C20H41NO2
54 4-(beta-Acetylaminoethyl)imidazole 154.0973617 15.53 Pos C7H11N3O


Regulation of key metabolite-related metabolic pathways by HRW

To explore the metabolic pathways of ionizing radiation rats, differential metabolites were introduced into MetaboAnalyst for pathway enrichment analysis. Pathway effect > 0.1 was considered to be a key metabolic pathway (the results are shown in Fig. 3), and five metabolic pathways related to ionizing radiation were identified, including: glutathione metabolism, ascorbic acid and aldonic acid metabolism, pentose and glucuronic acid interconversion, glycerophospholipid metabolism, alanine, aspartic acid and glutamic acid metabolism. Four metabolic pathways (glutathione metabolism, ascorbic acid and aldonic acid metabolism, pentose and glucuronic acid interconversion, and glycerophospholipid metabolism) are specifically regulated by HRW. The above results suggest that HRW may improve ionizing radiation by regulating multiple metabolic pathways. It reflects the multi-objective and multi-pathway mechanism of HRW.

Metabolic pathway analysis of MetPA. (A) Analysis of ionizing radiation-related metabolic pathways. (B) Analysis of HRW-regulated metabolic pathways


Integrative analysis of metabolomics and network pharmacology

An interaction network based on metabolomics and network pharmacology was constructed to fully understand the mechanism of HRW against ionizing radiation (Fig. 4). Differential metabolites were introduced into the Cytoscape’s MetScape plug-in to collect a compound-response-enzyme-gene network. By matching potential targets identified in network pharmacology to genes from the MetScape analysis, we identified one key target. The key relevant metabolites are phosphatidylcholine, choline and glucuronide salts. The pathways affected are glutathione metabolism, ascorbic acid and aldehyde acid metabolism, pentose and glucuronic acid interconversion, and glycerophospholipid metabolism. They may play an important role in the therapeutic effects of FWR on ionizing radiation.

An external file that holds a picture, illustration, etc.
Object name is 12860_2023_491_Fig4_HTML.jpg

The “active ingredient-target/gene-metabolite” biological network constructed by Cytoscape. (Blue nodes: metabolite-related target genes, red nodes: differential metabolites.)


Pathway analysis

The KEGG enrichment analysis was performed on the targets and genes associated with recovery of differential metabolites, and the functions of these 93 key target genes were predicted and analyzed. The KEGG enrichment results showed that key signaling pathways sensitive to ionizing radiation, such as glutathione metabolism (C2), arachidonic acid metabolism(C1), glycerophospholipid metabolism(C1), and ether lipid metabolism(C1) were significantly enriched []. There were interactions between these signal pathways that can synergistically affect pathways that in turn, synergistically affected the HRW treatment with ionizing radiation (Fig. 5).

The KEGG enrichment analysis of HRW regulation-related targets


Therapeutic performance evaluation of biomarkers

Brain tissue biomarkers are the cornerstone of ionizing radiation diagnosis, and their function was assessed by ROC analysis. A total of three metabolites were selected as potential biomarkers including phosphatidylcholine, choline, and D-glucuronate. Figure 6 A shows that the AUC for phosphatidylcholine, choline, and D- glucuronide were 0.8147, 0.6843, and 0.6503, respectively. In addition, using a biomarker combination that provided greater predictive power than a single biomarker, the combined result for phosphatidylcholine, choline, and D-glucuronide was 0.8771 by logistic regression analysis, as shown in Fig. 6B. These metabolites may be potential biomarkers of HRW in ionizing radiation therapy.

ROC curves of the combined markers used to evaluate the effect of the HRW treatment. (A) Phosphatidylcholine; (B) Choline; (C) D-Glucuronate. The combination of phosphatidylcholine, choline, and D-Glucuronate, and had the highest AUC = 0.8771



Radiation remains a health threat and safety hazards, and. radiation protection agents are attracting increasing attention. WR-2721 is the only FDA-approved radioprotective agent that can prevent radiation-induced damage, especially hydroxyl radical damage. WR-2721 can cause dizziness, vomiting and other adverse reactions. Therefore, there is an urgent need to develop an effective and safe radioprotectants. Hydrogen is a safe, convenient, non-toxic substance with little or no side effects and selective antioxidant function. Hydrogen scavenging prevents oxidative stress-induced damage [, ]. Hydrogen has a protective effect on various tissues from radiation damage. For example, intake of hydrogen-rich water after radiation can alleviate damages to the immune system, skin, hematopoietic system, small intestine, and heart by reducing reactive oxygen species (ROS) []. Their selective neutralization of harmful oxygen radicals has attracted increasing attention in the field of radiation protection. Hydrogen can cross the blood-brain barrier and quickly cross cell membranes to interact with toxic ROS. However, the material basis and potential mechanisms of hydrogen utilization in many neurological diseases remain unclear.

In our previous study, a single partial whole-brain irradiation with a 6-MeV electron beam of 30 Gy could successfully induce cognitive dysfunction in rats []. Compared with control rats, the number of new neurons decreased by 67% on the 7th day after 30 Gy whole brain irradiation []; and on the 30th day after irradiation, there were no surviving new neurons, indicating that irradiation with 30 Gy successfully induced acute cognitive dysfunction. In this study, a single 30 Gy electron beam whole brain irradiation was used to establish a model of cognitive dysfunction in rats. Then, we investigated the effects of HRW on cognitive function and its possible mechanism using Metabolomics combined with biological networks. The results showed that after HRW treatment, the cognitive dysfunction of rats was alleviated, and the protective effect was related to the regulation of phosphatidylcholine, choline, and D-Glucuronate. The results of brain metabolomics showed that HRW can significantly regulate the levels of 54 different metabolites, including glutathione metabolism, ascorbate and aldarate metabolism, pentose and glucuronate interconversions, and glycerophospholipid metabolism. These pathways were discussed below.

Pyroglutamic acid is a cyclized derivative of L- glutamic acid produced by non-enzymatic reaction of glutamic acid, glutamine and γ-glutamine peptide. It can also be generated by the reaction of γ-glutamyl cyclotransferase with amino acids. Acquired pyroglutamic acid deficiency (penicillins) and glutathione deficiency (e.g., malnutrition or sepsis) often result in moderate to severe encephalopathy clinically []. In our study, decreased levels of relevant metabolites were observed to decrease in ionizing radiation rats, and these metabolites could be restored by HRW, indicating that glutathione metallurgy may be the mechanism by which HRW treats ionizing radiation-induced brain injury.

Cognitive impairment seriously affects people’s physical and mental health and also causes a substantial economic burden []. This study showed that HRW can regulate the glycerophospholipid metabolism, which is the most abundant phospholipid in the body, and can activate the phospholipase A2 pathway to release arachidonic acid, thereby regulating the metabolism of glycerophospholipids and improve depression []. In addition, ascorbate and aldarate metabolism, pentose and gluconate interconversions are closely related to cognitive impairment [, ]. Our study found that HRW can improve ionizing radiation-induced brain cognitive impairment by modulating these pathways.

Although this study systematically explained the mechanism by which HRW improves cognitive impairment, there are limitations that require further study: This study only used non-targeted brain Metabolomics and biological target network at the level of metabolites and metabolic pathways. This model could be easily exploited in future works for the elucidation of pharmacological effects and mechanisms using a multi-omics approach based on integrated transcriptomics, proteomics and metabolomics studies.Further studies warrant adding corresponding pathway inhibitors at the cellular and animal levels for verification. In depth analysis though might be worth being conducted in future studies.


In summary, the mechanisms of action of HRW on ionizing radiation were predicted through integrating biological networks and metabolomics. HRW may exert its potential effects by regulating 54 metabolites. In addition, glutathione metabolism, ascorbate and aldonic acid metabolism, pentose and glucuronate interconversions, and glycerophospholipid metabolism as well as the levels of phosphatidylcholine, choline, and D-Glucuronate are key metabolic pathways/metabolites for HRW function. These findings will help reveal the mechanism of HRW on ionizing radiation and provide a rational way to clarify this mechanism.


We thank the support of the Department of Radiology and Environmental Medicine, China Institute for Radiation Protection,Shanxi,China.

Authors’ contributions

LXM conceived and performed the experiment. LMY wrote the manuscript. LH and WY analyzed the data, YH and CXM assisted in research. QXJ and LJG designed the study and protocol. All authors have understood, agreed and approved the final edition of this manuscript.


The authors of this study thank their supporters and financial support. This work was supported by the Natural Science Foundation of Shanxi Province (grant no. 202203021211005 and 201801D1213).

Data Availability

All data generated or analyzed during this study are included in this published article.


Ethics approval and consent to participate

All procedures performed in studies were approved by the Animal Protection and Ethics Committee of the China Institute for Radiation Protection, and was reported in accordance with ARRIVE guidelines. All methods are performed in accordance with the Regulations of China Institute for Radiation Protection on the administration of laboratory animals.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.



Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.



Xiaoming Liu and Mengya Liu contributed equally to this work.



1. Jafari M, Salehi M, Asgari A, Ahmadi S, Abbasnezhad M, Hajihoosani R, et al. Effects of paraoxon on sbrain tissue biochemical parameters and oxidative stress induction in various tissues of Wistar and Norway rats. Environ Toxicol Pharmacol. 2012;34:876–87. doi: 10.1016/j.etap.2012.08.011. [PubMed] [CrossRef] []
2. Khazaie S, Jafari M, Heydari J, Asgari A, Tahmasebi K, Salehi M et al. Modulatory effects of vitamin C on biochemical and oxidative changes induced by acute exposure to diazinon in rat various tissues: prophylactic and therapeutic roles. J Anim Physiol Anim Nutr. 2019, 1–10. [PubMed]
3. Rola R, Raber J, Rizk A, Otsuka S, VandenBerg SR, Morhardt DR, et al. Radiation-induced impairment of hippocampal neurogenesis is associated with cognitive deficits in young mice. Exp Neurol. 2004;188:316–30. doi: 10.1016/j.expneurol.2004.05.005. [PubMed] [CrossRef] []
4. Kalm M, Karlsson N, Nilsson MK, Blomgren K. Loss of hippocampal neurogenesis, increased novelty-induced activity, decreased home cage activity, and impaired reversal learning one year after irradiation of the young mouse brain. Exp Neurol. 2013;247:402–9. doi: 10.1016/j.expneurol.2013.01.006. [PubMed] [CrossRef] []
5. Ohta S. Molecular hydrogen as a novel antioxidant: overview of the advantages of hydrogen for medical applications. Methods Enzymol. 2015;555:289–317. doi: 10.1016/bs.mie.2014.11.038. [PubMed] [CrossRef] []
6. Cai J, Kang Z, Liu WW, Luo X, Qiang S, Zhang JH, et al. Hydrogen therapy reduces apoptosis in neonatal hypoxia-ischemia rat model. Neurosci Lett. 2008;441:167–72. doi: 10.1016/j.neulet.2008.05.077. [PubMed] [CrossRef] []
7. Zhou J, Chen Y, Huang GQ, Li J, Wu GM, Liu L, et al. Hydrogenrich saline reverses oxidative stress, cognitive impairment, and mortality in rats submitted to sepsis by cecal ligation and puncture. J Surg Res. 2012;178:390–400. doi: 10.1016/j.jss.2012.01.041. [PubMed] [CrossRef] []
8. Ichihara M, Sobue S, Ito M, Ito M, Hirayama M, Ohno K. Beneficial biological effects and the underlying mechanisms of molecular hydrogen - comprehensive review of 321 original articles. Med Gas Res. 2015;5:12. doi: 10.1186/s13618-015-0035-1. [PMC free article] [PubMed] [CrossRef] []
9. Collino S, Martin FP, Rezzi S. Clinical Metabolomics paves the way towards Future Healthcare Strategies. Br J Clin Pharmacol. 2013;75:619–29. doi: 10.1111/j.1365-2125.2012.04216.x. [PMC free article] [PubMed] [CrossRef] []
10. Liu XJ, Wang YZ, Lv M, Zhao SJ, Chen SJ, Li SY, Qin XM. Brain tissue Metabolomics reveals compatibility rules of the antidepressant Effects of Xiaoyaosan and its efficacy groups. Psychiatry Res. 2021;299:113827. doi: 10.1016/j.psychres.2021.113827. [PubMed] [CrossRef] []
11. Zhang WN, Yang L, He SS, Qin XM, Li AP. Metabolomics coupled with integrative pharmacology reveal the Protective Effect of FangjiHuangqi Decoction against Adriamycin-induced rat nephropathy model. J Pharm Biomed Anal. 2019;174:525–33. doi: 10.1016/j.jpba.2019.05.023. [PubMed] [CrossRef] []
12. Fu C, Wu Q, Zhang Z, Xia Z, Ji H, Lu H, Wang Y, UPLC-ESI-IT -TOF-MS metabolomic study of the therapeutic effect of Xuefu Zhuyu decoction on rats with traumatic brain Injury. J Ethnopharmacol. 2019;245:112149. doi: 10.1016/j.jep.2019.112149. [PubMed] [CrossRef] []
13. Simón-Manso Y, Lowenthal MS, Kilpatrick LE, Sampson ML, Telu KH, Rudnick PA, Mallard WG, Bearden DW, Schock TB, Tchekhovskoi DV, et al. Metabolite profiling of a NIST Standard Reference Material for Human plasma (SRM 1950): GC-MS, LC-MS, NMR, and clinical laboratory analyses, libraries, and web-based Resources. Anal Chem. 2013;85:11725–31. doi: 10.1021/ac402503m. [PubMed] [CrossRef] []
14. Singh K, Trivedi R, Verma A, D’souza MM, Koundal S, Rana P, Baishya B, Khushu S. Altered Metabolites of the Rat Hippocampus After Mild and Moderate Traumatic Brain Injury - a Combined in Vivo and in Vitro 1H-MRS Study. NMR. Biomed. 2017, 30. [PubMed]
15. Zhou WX, Cheng XR, Zhang YX. Network Pharmacology - a New Philosophy for understanding of Drug Action and Discovery of New Drugs. Chin J Pharmacolo Toxicolo. 2012;26:4–9. []
16. Niu M, Zhang SQ, Zhang B, Yang K, Li S. Interpretation of network pharmacology evaluation Method Guidance. Chin Tradit Herbal Drugs. 2021;52:4119–29. []
17. Liu XJ, Wei FX, Liu HL, Zhao SJ, Du GH, Qin XM. Integrating hippocampal metabolomics and Network Pharmacology deciphers the antidepressant mechanisms of Xiaoyaosan. J Ethnopharmacol. 2021;268:113549. doi: 10.1016/j.jep.2020.113549. [PubMed] [CrossRef] []
18. Stirling ER, Cook KL, Roberts DD, Soto-Pantoja DR. Metabolomic analysis reveals unique biochemical Signatures Associated with Protection from Radiation Induced Lung Injury by lack of cd47 receptor gene expression. Metabolites. 2019;9:218. doi: 10.3390/metabo9100218. [PMC free article] [PubMed] [CrossRef] []
19. Lin CL, Huang WN, Li HH, Huang CN, Hsieh S, Lai C, et al. Hydrogen-rich water attenuates amyloid beta-induced cytotoxicity through upregulation of Sirt1-FoxO3a by stimulation of AMPactivated protein kinase in SK-N-MC cells. Chem Biol Interact. 2015;240:12–21. doi: 10.1016/j.cbi.2015.07.013. [PubMed] [CrossRef] []
20. Wu Q, Su N, Cai J, Shen Z, Cui J. Hydrogen-rich water enhances cadmium tolerance in chinese cabbage by reducing cadmium uptake and increasing antioxidant capacities. J Plant Physiol. 2015;175:174–82. doi: 10.1016/j.jplph.2014.09.017. [PubMed] [CrossRef] []
21. DeCarolis NA, Rivera PD, Ahn F, Amaral WZ, LeBlanc JA, Malhotra S, et al. (56)fe particle exposure results in a long-lasting increase in a cellular index of genomic instability and transiently suppresses adult hippocampal neurogenesis in vivo. Life Sci Space Res (Amst) 2014;2:70–9. doi: 10.1016/j.lssr.2014.06.004. [PMC free article] [PubMed] [CrossRef] []
22. Qian L, Cao F, Cui J, Huang Y, Zhou X, Liu S, et al. Radioprotective effect of hydrogen in cultured cells and mice. Free Radic Res. 2010;44:275–82. doi: 10.3109/10715760903468758. [PubMed] [CrossRef] []
23. Verheyde J, Benotmane MA. Unraveling the fundamental molecular mechanisms of morphological and cognitive defects in the irradiated brain. Brain Res Rev. 2007;53:312–20. doi: 10.1016/j.brainresrev.2006.09.004. [PubMed] [CrossRef] []
24. Zhao S, Yang Y, Liu W, Xuan Z, Wu S, Yu S, et al. Protective effect of hydrogen-rich saline against radiation-induced immune dysfunction. J Cell Mol Med. 2014;18:938–46. doi: 10.1111/jcmm.12245. [PMC free article] [PubMed] [CrossRef] []
25. Zhou P, Lin B, Wang P, Pan T, Wang S, Chen W, et al. The healing effect of hydrogen-rich water on acute radiation-induced skin injury in rats. J Radiat Res. 2019;60:17–22. doi: 10.1093/jrr/rry074. [PMC free article] [PubMed] [CrossRef] []
26. Liu M, Yuan H, Yin J, Wang R, Song J, Hu B, Li J, Qin X. Effect of Hydrogen-Rich Water on Radiation-Induced Cognitive Dysfunction in rats. Radiat Res. 2020;193:16–23. doi: 10.1667/RR15464.1. [PubMed] [CrossRef] []
27. Ji S, Tian Y, Lu Y, Sun R, Ji J, Zhang L, Duan S. Irradiation-induced hippocampal neurogenesis impairment is associated with epigenetic regulation of bdnf gene transcription. Brain Res. 2014;1577:77–88. doi: 10.1016/j.brainres.2014.06.035. [PubMed] [CrossRef] []
28. Luyasu S, Wamelink MM, Galanti L, Dive A. Pyroglutamic acid-induced metabolic acidosis: a case report. Acta Clin Belg. 2014;69:221–3. doi: 10.1179/2295333714Y.0000000022. [PubMed] [CrossRef] []
29. Dharshini SAP, Jemimah S, Taguchi YH, Gromiha MM. Exploring common therapeutic targets for neurodegenerative Disorders using Transcriptome Study. Front Genet. 2021;12:639160. doi: 10.3389/fgene.2021.639160. [PMC free article][PubMed] [CrossRef] []
30. Wang J, Li J, Liu K, Wang S, Su Q, Cheng Y, Wang Y, Wang Y. Integrated lipidomics and network pharmacology analysis of the protective effects and mechanism of Yuanzhi San on rats with cognitive impairment. Bioorg Med Chem. 2022;58:116651. doi: 10.1016/j.bmc.2022.116651. [PubMed] [CrossRef] []
31. Chen KD, Chang PT, Ping YH, Lee HC, Yeh CW, Wang PN. Gene expression profiling of peripheral blood leukocytes identifies and validates ABCB1 as a novel biomarker for Alzheimer’s disease. Neurobiol Dis. 2011;43:698–705. doi: 10.1016/j.nbd.2011.05.023. [PubMed] [CrossRef] []
32. Wise EL, Yew WS, Akana J, Gerlt JA, Rayment I. Evolution of enzymatic activities in the orotidine 5’-monophosphate decarboxylase suprafamily: structural basis for catalytic promiscuity in wild-type and designed mutants of 3-keto-L-gulonate 6-phosphate decarboxylase. Biochemistry. 2005;44:1816–23. doi: 10.1021/bi0478143. [PubMed] [CrossRef] []
Back to blog

Leave a comment