# Identification of Novel Extracellular-Signal-Regulated Kinase 2 Inhibitors Through Machine Learning-Driven De Novo Design, Molecular Docking, and Free-Energy Perturbation

**Authors:** Ibrahim A. Alsarra, Mahima Sudhir Kolpe, Md Ataul Islam

PMC · DOI: 10.3390/ph19020337 · 2026-02-20

## TL;DR

This study uses machine learning and computational methods to design new ERK2 inhibitors, which could lead to better drugs for targeting cell signaling pathways.

## Contribution

A novel ML-driven de novo design approach identifies ERK2 inhibitors with high binding affinity and drug-like properties.

## Key findings

- Four molecules (Ek1-Ek4) showed high ERK2 binding affinity (-9.50 to -10.50 kcal/mol) via molecular docking.
- MD simulations confirmed strong ERK2 binding and low backbone deviation for the proposed molecules.
- Ek1 exhibited higher FEP energy (-26.85 kJ/mol) than a standard molecule, indicating stronger binding.

## Abstract

Background: The extracellular-signal-regulated kinase (ERK) cascade regulates cell proliferation, differentiation, and survival, and ERK2 mediates substrate phosphorylation, influencing gene expression and cellular functions. Methods: In the current study, a pool of new molecules was generated using the DeLA-Drug, a machine learning (ML)-assisted de novo design tool. The chemical space was reduced through a similarity search against active ERK2 inhibitors and molecular docking with AutoDock vina, followed by pharmacokinetic assessment in DeepPK. Poses of the final selected molecules were refined in DiffDock, and dynamicity was assessed through molecular dynamics (MD) simulation. Finally, the free-energy perturbation (FEP)-based binding affinity was explored in Gromacs2023.4. Results: From the above approaches, four molecules (Ek1, Ek2, Ek3, and Ek4) were identified as promising candidates with favorable binding interactions. Molecular docking revealed that the selected molecules exhibited higher binding affinity for ERK2, ranging from −9.50 to −10.50 kcal/mol. The dynamics assessment via MD simulation clearly revealed their strong association with ERK2, corroborated by the lower deviation of the ERK2 backbone in dynamic states. All four screened molecules have satisfactory pharmacokinetic properties, medicinal chemistry properties, and good synthetic accessibility scores, indicating their potential as drug-like compounds under Lipinski’s rule of five to inhibit or modulate ERK2 activity. The FEP energy of Ek1 was found to be −26.85 kJ/mol, which is higher than the standard molecule (−22.77 kJ/mol) and indicates its strong affinity toward ERK2. Conclusions: These results suggest that all proposed ERK2 modulators are potential avenues for future drug discovery targeting ERK2, subject to experimental validation.

## Linked entities

- **Proteins:** MAPK1 (mitogen-activated protein kinase 1)

## Full-text entities

- **Genes:** MAPK3 (mitogen-activated protein kinase 3) [NCBI Gene 5595] {aka ERK-1, ERK1, ERT2, HS44KDAP, HUMKER1A, P44ERK1}, MAPK1 (mitogen-activated protein kinase 1) [NCBI Gene 5594] {aka ERK, ERK-2, ERK2, ERT1, MAPK2, NS13}, EPHA3 (EPH receptor A3) [NCBI Gene 2042] {aka EK4, ETK, ETK1, HEK, HEK4, TYRO4}, Mapk1 (mitogen activated protein kinase 1) [NCBI Gene 116590] {aka ERK-2, ERT1, Erk2, p42-MAPK}, EPHB3 (EPH receptor B3) [NCBI Gene 2049] {aka EK2, ETK2, HEK2, TYRO6}, MAP2K7 (mitogen-activated protein kinase kinase 7) [NCBI Gene 5609] {aka JNKK2, MAPKK7, MEK, MEK 7, MKK7, PRKMK7}, EPHA8 (EPH receptor A8) [NCBI Gene 2046] {aka EEK, EK3, HEK3}, BRAF (B-Raf proto-oncogene, serine/threonine kinase) [NCBI Gene 673] {aka B-RAF1, B-raf, BRAF-1, BRAF1, NS7, RAFB1}, Ephb1 (Eph receptor B1) [NCBI Gene 24338] {aka Ephb2, Erk, elk}
- **Diseases:** solid (MESH:D018250), colorectal, lung, and pancreatic carcinomas (MESH:D015179), Toxicity (MESH:D064420), 19A (OMIM:615528), autoimmune disorders (MESH:D001327), ADV (MESH:D049932), Rat sarcoma (MESH:D012509), skin sensitization (MESH:D012871), injury to (MESH:D014947), inflammation (MESH:D007249), melanoma (MESH:D008545), cancer (MESH:D009369)
- **Chemicals:** 19A. (MESH:C039200), ATP (MESH:D000255), MK-8353 (MESH:C000632607), Hydrogen (MESH:D006859), Na+ (MESH:D012964), DeLA (-), LY3214996 (MESH:C000719760), Cl- (MESH:D002713), abemaciclib (MESH:C000590451), benzene (MESH:D001554), water (MESH:D014867), BVD-523 (MESH:C000618314), FR180204 (MESH:C505241), pyrazole (MESH:C031280), SCH772984 (MESH:C587178)
- **Species:** Homo sapiens (human, species) [taxon 9606]
- **Mutations:** serine/threonine, V600E

## Figures

12 figures with captions in the complete paper: https://tomesphere.com/paper/PMC12943845/full.md

---
Source: https://tomesphere.com/paper/PMC12943845