1. Academic Validation
  2. Discovery of PAR2 antagonists through fragment molecular orbital-based hotspot analysis and allosteric sodium site dynamics simulations

Discovery of PAR2 antagonists through fragment molecular orbital-based hotspot analysis and allosteric sodium site dynamics simulations

  • Bioorg Chem. 2025 Dec:167:109285. doi: 10.1016/j.bioorg.2025.109285.
Junyoung Park 1 Hocheol Lim 2 Yechan Lee 3 Taegun Kim 4 Heejeong Song 5 Suah Lee 5 Yunkyung Heo 3 Jinhyuk Kim 6 Wan Namkung 3 Gyoonhee Han 7
Affiliations

Affiliations

  • 1 The Interdisciplinary Program in Integrative Biotechnology, Yonsei University, Incheon 21983, Republic of Korea.
  • 2 Bioinformatics and Molecular Design Research Center (BMDRC), Incheon 21983, Republic of Korea.
  • 3 College of Pharmacy and Yonsei Institute of Pharmaceutical Sciences, Yonsei University, Incheon 21983, Republic of Korea.
  • 4 Department of Biotechnology, Yonsei University, Seoul 03722, Republic of Korea.
  • 5 Graduate Program of Industrial Pharmaceutical Science, College of Pharmacy, Yonsei University, Incheon 21983, Republic of Korea.
  • 6 Department of Integrative Biotechnology, Yonsei University, Incheon 21983, Republic of Korea.
  • 7 The Interdisciplinary Program in Integrative Biotechnology, Yonsei University, Incheon 21983, Republic of Korea; Bioinformatics and Molecular Design Research Center (BMDRC), Incheon 21983, Republic of Korea; Department of Biotechnology, Yonsei University, Seoul 03722, Republic of Korea; Postech Biotech Center, Pohang, Gyeongbuk 37673, Republic of Korea. Electronic address: gyoonhee@yonsei.ac.kr.
Abstract

Protease-activated Receptor 2 (PAR2) is a G protein-coupled receptor (GPCR) implicated in inflammatory diseases and Cancer, representing a validated therapeutic target. In our previous study, we combined molecular dynamics (MD) and fragment molecular orbital (FMO) analyses to rationalize binding mode differences between PAR2 agonists and antagonists. Here, we advanced this framework into a practical ligand discovery platform that integrates FMO-based hotspot mapping, pharmacophore modeling, structure-based virtual screening, and functional validation. FMO analysis of PAR2 complexes with a peptide agonist, antibody, and antagonist identified residue-level binding hotspots that were translated into a pharmacophore model for screening ∼50,000 in-house compounds. Hit compounds were virtually screened by docking and MD simulations, with Na+ ion displacement at the conserved allosteric site serving as a functional filter. Application of a 1.4 Å cutoff excluded agonist-like perturbations while retaining antagonist-like Na+ displacement profiles. The initial hit, an indazole derivative S757-0124, engaged key hotspots and exhibited antagonist-like Na+ dynamics but suffered from poor solubility and moderate potency. Scaffold replacement with an sp3-rich tetrahydropyrazolopyridine core yielded P2L-003, which showed 35-fold improved potency, 11.6-fold higher solubility, and maintained reduced Na+ displacement characteristic of antagonist binding. Functional assays further confirmed that P2L-003 selectively antagonized PAR2 without affecting PAR1, PAR4, or ATP-mediated signaling, while Western blotting demonstrated dose-dependent inhibition of ERK1/2 and p38 phosphorylation. Together, these results highlight how integrating electronic precision from FMO, functional filtering from MD, and medicinal chemistry optimization provides a robust and generalizable workflow for GPCR antagonist discovery.

Keywords

Antagonist; Fragment molecular orbital; GPCR; Molecular dynamics; Na(+) ion displacement; PAR2; Pharmacophore.

Figures
Products