Case Study

Profiling dynamic severity in hypoxic ischemic encephalopathy

Severity of presentation in neonates with suspected hypoxic ischemic encephalopathy (HIE) has a dynamic metabolomic profile. Following the initation of therapeutic intervention, metabolite severity profiles evolve over time. These metabolite changes implicate different protein networks with potential implications for patient management and intervention protocols.

Executive Takeaway

Severity in HIE is dynamic. Metabolomics paired with Splash & Ripple identifies the protein networks that may be implicated.

Researchers at the University of Pretoria wanted to understand the drivers of severity in hypoxic ischemic encephalopathy (HIE). Metabolomics was used to profile the biochemical changes in these neonates at two timepoints during the application of therapeutic hypothermia. The metabolite profile revealed that severity in HIE is associated with distinct biochemical changes at different timepoints.

The Splash & Ripple application mapped these changes to protein candidates and provided mechanistic arguments for the observed perturbations. If severity is dynamic, then the appropriate intervention may have temporal implications. Metabolomics paired with Splash & Ripple helped to reveal this dynamic nature of severity in HIE, and the potential mediators involved.

Proof at a glance

What the analysis delivered

The scientific challenge

Metabolites are the readout. Proteins help explain the mechanism.

Metabolomics measures the functional end-state of biology, providing a rich readout of molecular phenotypes. That richness is exactly what makes it powerful and challenging to interpret in isolation.

A severity contrast in HIE returns dozens of significantly altered metabolites. Pathway enrichment will tell you the lipids moved, or the amino acids moved. But enrichment stops at the category. It cannot tell you which enzyme, transporter or regulator sits at the point where several of those compounds meet, and it cannot connect a metabolite in the blood to a signaling protein that never appears on a metabolic pathway diagram at all.

That is a real problem in neonatal encephalopathy. The injury is a cascade: primary energy failure, reperfusion, oxidative injury, inflammation, delayed cell death. Each phase may have different protein drivers. A metabolite list from a single timepoint, read on its own, cannot distinguish which phase you are looking at or whether what you are seeing is damage or the body’s attempt to contain it.

  • Enrichment analysis tells you the category: A categorisation of metabolite signals without any protein context is limited.
  • Correlation networks stay inside one omic: Metabolite–metabolite relationships reveal patterns within the measured profile. Connecting those patterns to proteins adds a route to mechanistic interpretation.
  • Splash & Ripple names the convergence point: Places every perturbed metabolite on a curated protein map, and looks for proteins where multiple metabolite signals converge.

The foundation

Metabolite Affiliated Proteins (M.A.P) database

Metabolon’s expertly curated Metabolite Affiliated Proteins (M.A.P) database connects named chemical compounds to documented molecular relationships. This comprehensive resource empowers the application to deliver meaningful insights.

The M.A.P database leverages the Reactome database as a foundation. Metabolon has built on this foundation to cover more chemical compounds and curate additional relationships. Of the 2,820 distinct Metabolon library chemical IDs represented in the M.A.P database, only 25.4% have a Reactome (version 92) association. Metabolon’s curated resource contains almost three times more metabolite-protein relationships than Reactome. Results from an internal validation (described in accompanying white paper) at Metabolon revealed that across 131 samples with a diagnosed mendelian disorder, Metabolon’s map recovered 72.5% of diagnostic proteins compared to 57.3% using Reactome alone.

metabolite-protein connections diagram

metabolite-protein connections diagram

Metabolite-Protein connections. In this study, 93.75% (0–6 h) and 92.79% (48–72 h) of perturbed metabolites mapped directly to at least one protein, and 100% of the top ten perturbations at 48–72 hours.

How Splash & Ripple was applied

Begin with metabolomics

Metabolomics measures the consequences of genetic and environmental insults. Splash & Ripple uses expertly curated metabolite-protein relationships to build on the metabolite measurements, leveraging a protein association network to explore neighbourhoods connected to the pattern of perturbed metabolites.

We rarely observe the insult, or the protein it acts on, directly. What we can measure are its metabolomic consequences. Splash & Ripple treats each perturbed metabolite as that measured readout and works backward across the protein network to uncover the biology behind it: the proteins whose altered behavior the metabolite reflects, and the wider neighborhood connected to them. The proteins closest to the readout form the Splash; expanding the network to more proteins with higher degrees of separation gives the Ripple.

Splash & Ripple reduces the manual work of connecting metabolites to protein biology and makes that interpretation more accessible to researchers without specialist expertise.

the map is the bridge

the metabolite is the readout

The metabolite is the readout; the proteins are the candidates. At the center sits the measured metabolite. Working outward, Splash & Ripple recovers the proteins behind it and groups them by how close they are to the measured metabolites: directly mapped proteins, near neighbors at one or two interaction steps, and a separate Ripple band reached at three or four steps under stricter criteria. The groups describe network proximity, not causal certainty, and each is scored and ranked on its own terms.

Aligned themes, alongside highly ranked metabolites, can be used to filter the candidate protein list. This approach focuses analyses on discrete biological relationships that can be further explored using Sankey plots to visualize the potential protein-protein network pathways that connect metabolites of interest to a specific protein candidate.

  1. Step 1: What the metabolomic data tells usPathway enrichment on the perturbed metabolites produces themes grounded in the metabolite measurements with no prior assumption about the disease.
  2. Step 2: What the phenotype predicts
    Independently, the phenotype is used to derive the biology a domain expert might anticipate.
  3. Step 3: Aligned themes
    Where an expected theme has support from a data-driven theme it is categorised as an aligned theme. Aligned themes therefore represent signal in the metabolomic dataset that has some potential relationship to the study phenotype. It acts as a filter to narrow the candidate protein list to only those most relevant to the research question.

Study and data context

A severity contrast, held at two points in the injury cascade

The cohort was assembled at the University of Pretoria as part of the NESHIE study under the direction of Prof. Michael Pepper on behalf of the NESHIE team, with Dr. Caroline Foden as primary genomics and multi-omics researcher. Dried blood spots were collected from neonates diagnosed with HIE at two timepoints over the course of the standard-of-care intervention, therapeutic hypothermia: 0–6 hours (n = 84) and 48–72 hours (n = 102).

Critically, the statistical design was not disease versus health. At each timepoint the analysis compared the most severe presentations against the less severe. Every metabolite entering Splash & Ripple is therefore a severity signal within an already-affected, already-treated population which is far more specific and far more clinically relevant question than “what does HIE look like?”

Because the two comparisons were run independently, they act as a natural internal check on each other. Anything that appears at both timepoints is a persistent feature of severity; anything that appears at only one is time-resolved.

study design

Study design. Two independent severity contrasts, each analyzed separately, giving a time-resolved view of what distinguishes the most severely affected neonates.

Timepoint 1: 0-6 hours (n=84)

Early severity was dominated by lipid and oxidative-stress signals

Every significantly enriched pathway in the early severity contrast was a lipid pathway.

Lysoplasmalogen metabolism (adj. p = 0.0028) and omega-3/omega-6 fatty acid synthesis (adj. p = 0.0032) led the list, followed by alpha-linolenic acid metabolism, broad fatty acid metabolism and arachidonic acid oxylipin metabolism. Two of these formed aligned themes: ether lipid/plasmalogen remodeling with oxidative stress defense, and PUFA imbalance with oxylipin inflammatory signaling, because oxidative membrane injury and neuroinflammation are exactly what the HIE phenotype predicts.1,2

  • Observed: Polyunsaturated fatty acids depleted across the board: arachidonate, linoleate, linolenate, mead acid, dihomo-linolenate.
  • Observed: Plasmalogens — the brain’s endogenous membrane antioxidants — reduced at four separate species.1
  • Observed: Cysteine-glutathione disulfide reduced, alongside 14-HDoHE/17-HDoHE (rank 4), the hydroxylated DHA products.
  • Inferred: Simultaneous loss of PUFA substrate, plasmalogen antioxidant capacity and glutathione-linked species is the signature of lipid peroxidation outpacing defense.

Applied to the ranked protein table, the two aligned themes intersected with 1-(1-enyl-palmitoyl)-GPE (P-16:0) narrowed roughly 15,000 proteins to 342 and 441 candidates respectively. From those, the analysis surfaced four proteins, and three of them converge on the same axis.

the map is the bridge

Independent selections, one destination. GPX4, ALOX15 and ALOX5 were chosen from two different aligned themes and different metabolite anchors, yet all three sit on the glutathione–lipoxygenase axis that determines whether a cell undergoes ferroptosis.3–6 PPARG, the fourth candidate, sits on the opposing anti-inflammatory arm.7,8 This kind of convergence is the signal Splash & Ripple is built to expose, and it is not visible in a metabolite list.

The interpretation is coherent with the metabolomics: lipoxygenases consume polyunsaturated fatty acids to generate hydroperoxides; GPX4 uses glutathione to reduce them.4,5 If GPX4 capacity is exceeded, peroxidation propagates and the cell dies by ferroptosis.3,6 The measured profile shows substrate falling, plasmalogen antioxidant capacity falling, and the glutathione-linked pool falling all at once, and all more strongly in the most severe neonates.

Timepoint 2: 48-72 hours (n=102)

The severity signal shifted to nicotinamide metabolism.

Not one lipid pathway reached significance in the later contrast. The enrichment profile had moved to amino acids, nucleotides, nicotinate and nicotinamide metabolism, tryptophan, and ascorbate.

Three aligned themes emerged, each mapping onto a recognized phase of the encephalopathy cascade: excitotoxicity and amino acid neurotransmitter imbalance,9 mitochondrial energy failure with NAD⁺/ATP disruption,10,21 and oxidative stress with antioxidant depletion.11 And the analysis anchored all three on a single compound.

A path chosen and interrogated

A path chosen and interrogated

A path chosen and interrogated. Having ranked BDNF in Splash Indirect and 1-MNA as the top key metabolite, the analyst selected this specific protein–metabolite connection in the Sankey view and asked Splash & Ripple to reconstruct the intervening path. Each step is returned with an explicit evidence grade — published, emerging, or speculative — so the hypothesis can be evaluated rather than simply accepted.

1-methylnicotinamide (1-MNA)

Absent from the key metabolite profile at 0–6 hours. At 48–72 hours it is the highest-ranked compound to satisfy every key-metabolite criterion: significantly pathway-associated, not a xenobiotic, tractably connected, and it is significantly elevated in the most severe cases.

Nicotinamide

Nicotinamide itself is down. Its methylated product 1-MNA is up. 1-MNA’s own downstream catabolite, N1-methyl-2-pyridone-5-carboxamide, is up. NAD⁺ and NMN are up. Cortisol is elevated at rank 5.

Together, these changes are consistent with increased engagement of the nicotinamide methylation and clearance pathway.

Why interrogate BDNF in the first place?

Brain-derived neurotrophic factor is the brain’s principal pro-survival, pro-repair signal, and it is directly relevant to the biology flagged here: it counters glutamate excitotoxicity, the exact mechanism behind this timepoint’s excitotoxicity aligned theme, and enhancing BDNF/TrkB signaling is neuroprotective in neonatal hypoxia–ischemia models.12,13

BDNF is a promising biomarker in this clinical space, but it is not yet clear how to read it. Studies have linked both elevated and reduced circulating BDNF to worse encephalopathy severity or outcome, and the direction appears to depend on sampling time, matrix and treatment context.14–17 If the relationship to prognosis remains unsettled when BDNF is considered alone, placing it inside a network with other severity-linked metabolites and proteins may help clarify how the signal should be interpreted, or point toward alternative markers that travel with the same biology.

That is what makes the BDNF–1-MNA connection interesting. The link is not one-directional: 1-MNA has been reported to restore hippocampal BDNF signaling and produce TrkB-dependent effects in an independent model.18 So the elevated 1-MNA seen in the most severe neonates and the brain’s neurotrophic repair program may be coupled rather than separate stories. This remains a conjecture, but importantly it is a specific, testable one.

Compensation, or compensation that is failing?

NNMT transfers a methyl group from S-adenosylmethionine to nicotinamide, producing 1-MNA.19 Elevated 1-MNA is therefore a readout of NNMT flux: a recognized stress-responsive adaptation with anti-inflammatory and vasoprotective properties.19,20 Read on its own, the elevation looks protective.

But the samples were taken 48 to 72 hours after therapeutic hypothermia, and it is the most severe neonates in whom 1-MNA is highest. A compensatory response that is still active may imply that the insult is ongoing, or that the compensation is insufficient to prevent further damage.

(Scroll left/right to see full table)

Candidate Tier Rank / percentile Metabolite coverage Why it was selected
CXCL12 Splash Direct 6 · top 0.04% 36% CXCL12/CXCR4 signaling limits ROS-mediated damage; transiently upregulated after neonatal hypoxic–ischemic injury
APP Splash Direct 14 · top 0.09% 50% Regulates glutamatergic signaling; rapidly induced after neonatal HI injury
CCR2 Splash Direct 18 · top 0.12% 31% Induced by excitotoxic injury; mediates monocyte-driven inflammatory damage
CXCL8 Splash Direct 24 · top 0.16% 42% IL-8 is elevated in newborns with HIE and correlates with injury severity
ENPP1 Splash Indirect 2,518 · top 16.5% 42% Hydrolyzes extracellular ATP and NAD⁺-derived ADPR, modulating adenosine signaling
BDNF Splash Indirect 2,526 · top 16.57% 55% Counters glutamate excitotoxicity; promising neurotrophic biomarker in HIE whose direction remains unsettled — analyst-selected for path reconstruction to 1-MNA

CXCL8 acts like an internal control: interleukin-8 is an established plasma biomarker of HIE severity in newborns16. The modest rank of BDNF demonstrates that multiple features, apart from raw rank, are considered during lead selection.

The Results

Severity in HIE is fluid

a time-resolved mechanism

a time-resolved mechanism

Where is the value?

If severity in HIE is fluid, i.e., it looks different at different timepoints, implicating different systems, proteins, pathways, then this has repercussions for intervention management. The first step towards addressing this challenge is a full understanding of this temporal evolution of severity. Metabolomics paired with the Splash & Ripple application is one method to profile and contextualise how severity in HIE evolves over time as a step towards better patient management protocols.

The value of the approach

Search-space reduction reveals hypotheses to consider

Splash & Ripple helps researchers move from metabolite measurements to network-informed hypotheses.

Metabolomics revealed two sets of statistically significant compounds implicating divergent drivers of severity at different timepoints. Splash & Ripple provides greater context to each severity signal, and proposes mechanistic hypotheses for a selection of the observed changes.

By placing metabolomic data within the larger context of protein biology, researchers have an increased opportunity for idenitfying targetable biology. For example, GPX4 is a key regulator of ferroptosis, and understanding its activity in the early severity profile and any pressures acting on its capacity could ultimately inform novel or revised intervention strategies.

search-space reduction

Search-space reduction. Nothing is thrown away. Every protein retains a rank and can be revisited, but the researcher’s attention is directed to a handful of explainable starting points instead of a whole proteome.

  • Observed | Measured evidence: The metabolomic profile.
  • Inferred | Network position: The selected protein candidates and their relative rank based on their relationship to the metabolomic profile.
  • Hypothesis | Mechanistic reading: The relationship between individual protein-metabolite pairs is rationalised within the context of the study phenotype. Reasonable arguments are proposed for consideration.

Limitations and next steps

Informed decision making

Splash & Ripple generates hypotheses. It does not generate proof, and the outputs are explicitly built to make that distinction legible.

Limitations: What still requires validation
  • Protein-level measurements were not made. Every protein here is a network inference, not an observation.
  • Dried blood spots are peripheral. A circulating metabolite cannot be assigned to brain production. 1-MNA in particular has substantial renal contribution.
  • Metabolite abundance is an indirect measurement of pathway flux or protein activity. The ferroptosis and NNMT interpretations remain testable mechanistic hypotheses.
  • Severity associations cannot distinguish drivers of injury from downstream damage, treatment effects or compensatory responses.
Next steps: Where the hypotheses go
  • Targeted quantification of 1-MNA, nicotinamide and NAD⁺ intermediates against graded outcome measures.
  • Ferroptosis markers and GPX4/lipoxygenase activity in an appropriately matched model at the acute timepoint.
  • Replication of the severity contrast in an independent HIE cohort.
  • Longitudinal sampling to determine whether the 1-MNA elevation resolves in neonates with better outcomes.
Splash & Ripple case study featured image

Download the Case Study

Download the full case study PDF.

Acknowledgements

University of Pretoria

We acknowledge Prof. Michael Pepper, Dr. Caroline Foden, and the University of Pretoria research team for their work on the NESHIE study supporting this case study.

For inquiries about the NESHIE study, please contact the University of Pretoria team: Prof. Michael Pepper, email: michael.pepper@up.ac.za.

References

  1. Paul S, Lancaster GI, Meikle PJ. Plasmalogens: A potential therapeutic target for neurodegenerative and cardiometabolic disease. Prog Lipid Res. 2019. PMID: 30974122
  2. Zhang W, et al. Omega-3 polyunsaturated fatty acid supplementation confers long-term neuroprotection against neonatal hypoxic-ischemic brain injury through anti-inflammatory actions. Stroke. 2010. PMID: 20705927
  3. Peeples ES, Genaro-Mattos TC. Ferroptosis: A Promising Therapeutic Target for Neonatal Hypoxic-Ischemic Brain Injury. Int J Mol Sci. 2022. PMID: 35806425
  4. Cozza G, et al. Glutathione peroxidase 4-catalyzed reduction of lipid hydroperoxides in membranes: The polar head of membrane phospholipids binds the enzyme and addresses the fatty acid hydroperoxide group toward the redox center. Free Radic Biol Med. 2017. PMID: 28709976
  5. Haynes RL, van Leyen K. 12/15-lipoxygenase expression is increased in oligodendrocytes and microglia of periventricular leukomalacia. Dev Neurosci. 2013. PMID: 23838566
  6. Li C, et al. Ferroptosis contributes to hypoxic-ischemic brain injury in neonatal rats: Role of the SIRT1/Nrf2/GPx4 signaling pathway. CNS Neurosci Ther. 2022. PMID: 36184790
  7. Yi JH, et al. PPARgamma agonist rosiglitazone is neuroprotective after traumatic brain injury via anti-inflammatory and anti-oxidative mechanisms. Brain Res. 2008. PMID: 18948087
  8. Fang M, et al. The Neuroprotective Mechanisms of PPAR-γ: Inhibition of Microglia-Mediated Neuroinflammation and Oxidative Stress in a Neonatal Mouse Model of Hypoxic-Ischemic White Matter Injury. CNS Neurosci Ther. 2024. PMID: 39496476
  9. Zhang X, Peng K, Zhang X. The Function of the NMDA Receptor in Hypoxic-Ischemic Encephalopathy. Front Neurosci. 2020. PMID: 33117117
  10. Kawamura T, et al. Therapeutic Effect of Nicotinamide Mononucleotide for Hypoxic-Ischemic Brain Injury in Neonatal Mice. ASN Neuro. 2023. PMID: 37787108
  11. Zhao M, et al. Oxidative Stress in Hypoxic-Ischemic Encephalopathy: Molecular Mechanisms and Therapeutic Strategies. Int J Mol Sci. 2016. PMID: 27973415
  12. Almeida RD, et al. Neuroprotection by BDNF against glutamate-induced apoptotic cell death is mediated by ERK and PI3-kinase pathways. Cell Death Differ. 2005. PMID: 15905876
  13. Almli CR, et al. BDNF protects against spatial memory deficits following neonatal hypoxia-ischemia. Exp Neurol. 2000. PMID: 11031087
  14. Imam SS, et al. Cord blood brain derived neurotrophic factor: diagnostic and prognostic marker in fullterm newborns with perinatal asphyxia. Pak J Biol Sci. 2009. PMID: 20180326
  15. Parry SM, Peeples ES. The impact of hypoxic-ischemic brain injury on stem cell mobilization, migration, adhesion, and proliferation. Neural Regen Res. 2018. PMID: 30028311
  16. Massaro AN, et al. Plasma biomarkers of brain injury in neonatal hypoxic-ischemic encephalopathy. J Pediatr. 2018. PMID: 29478510
  17. Dietrick B, et al. Plasma and cerebrospinal fluid candidate biomarkers of neonatal encephalopathy severity and neurodevelopmental outcomes. J Pediatr. 2020. PMID: 32610169
  18. Zhao J, et al. Antidepressant-like effects of 1-methylnicotinamide in a chronic unpredictable mild stress model of depression. Neurosci Lett. 2021. PMID: 33248165
  19. Pissios P. Nicotinamide N-Methyltransferase: More Than a Vitamin B3 Clearance Enzyme. Trends Endocrinol Metab. 2017. PMID: 28291578
  20. Domagala TB, et al. Nitric oxide production and endothelium-dependent vasorelaxation ameliorated by N1-methylnicotinamide in human blood vessels. Hypertension. 2012. PMID: 22353616
  21. Kong W, Lu C. Role of mitochondria in neonatal hypoxic-ischemic encephalopathy. Histol Histopathol. 2024. PMID: 38314617
Splash & Ripple computer with softwarae screenshot

Splash & Ripple

Identify the Proteins Driving Your Biochemical Profile

Metabolon’s Splash & Ripple workflow works backward from metabolomic readouts, connecting perturbed metabolites to the proteins most likely to explain observed biological responses.

Get Started Today

Turn your metabolomic results into prioritized protein candidates, explicit mechanistic hypotheses, and a clearer path toward experimental validation.

Fill out the form to speak to an expert about leveraging Splash & Ripple for your next study.