Exploring potential targets and molecular mechanisms of traumatic brain injury exacerbated by Benzo(a)pyrene via network toxicology and molecular dynamics simulation.

Benzo(a)pyrene (BaP) is a common environmental pollutant from combustion sources that promotes oxidative stress, neuroinflammation and disruption of blood-brain barrier (BBB). However, its contribution to worsening traumatic brain injury (TBI) remains unclear. In this study, we aimed to assess the contribution of BaP to secondary injury in TBI. By integrating data from e.g., the Comparative Toxicogenomics Database, GeneCards, and Online Mendelian Inheritance in Man, 121 overlapping core targets were identified between BaP and TBI. Enrichment analyses via Gene Ontology and Kyoto Encyclopedia of Genes and Genomes, combined with protein-protein interaction networks and topological algorithms (degree, closeness centrality, betweenness centrality, average shortest path length, topological coefficient and partner of multi-edged node pairs), highlighted five hub genes (TP53, EGFR, AKT1, ACTB, and TNF) implicated in mitogen-activated protein kinase signaling, oxidative stress, and neuroinflammation. Molecular docking showed strong binding affinities of BaP to these hub proteins, with energies from -9.3 to -12.1 kcal/mol, tighter than co-crystal ligands and existing protein-binding drugs. Molecular dynamics simulations confirmed interaction stability through low root-mean-square deviation (< 0.5 nm), fluctuation, and radius of gyration values. Calculation of binding free energies using MM-PBSA validated the strong binding affinity between BaP and binding pockets of each hub genes. Toxicity prediction analysis revealed an oral LD