Identification of FDA-Approved Drugs as Potential Inhibitors of WEE2: Structure-Based Virtual Screening and Molecular Dynamics with Perspectives for Machine Learning-Assisted Prioritization

TL;DR

This study identifies FDA-approved drugs that may inhibit WEE2, a protein important for female fertility, using computer-based methods and suggests machine learning can help prioritize drug candidates.

Contribution

The study introduces an integrated in silico pipeline combining virtual screening, MD simulations, and ML for identifying and prioritizing WEE2 inhibitors.

Findings

Midostaurin and Nilotinib showed strong binding to WEE2 with affinities similar to a reference inhibitor.

MD simulations showed these drugs stabilize WEE2's catalytic domain and reduce conformational flexibility.

ML and DL methods are proposed to enhance hit prioritization and ADMET prediction for drug repurposing.

Abstract

Wee1-like protein kinase 2 (WEE2) is an oocyte-specific kinase that regulates meiotic arrest and fertilization. Its largely restricted expression in female germ cells and absence in somatic tissues make it a highly selective target for reproductive health interventions. Despite its central role in human fertility, no clinically approved WEE2 modulator is available. In this study, we employed an integrated in silico approach that combines structure-based virtual screening, molecular dynamics (MD) simulations, and MM-PBSA free-energy calculations to identify repurposed drug candidates with potential WEE2 inhibitory activity. Screening of ~3800 DrugBank compounds against the WEE2 catalytic domain yielded ten high-affinity hits, from which Midostaurin and Nilotinib emerged as the most mechanistically relevant based on kinase-targeting properties and pharmacological profiles. Docking analyses revealed strong binding affinities (−11.5 and −11.3 kcal/mol) and interaction fingerprints highly similar to the reference inhibitor MK1775, including key contacts with hinge-region residues Val220, Tyr291, and Cys292. All-atom MD simulations for 300 ns demonstrated that both compounds induce stable protein–ligand complexes with minimal conformational drift, decreased residual flexibility, preserved compactness, and stable intramolecular hydrogen-bond networks. Principal component and free-energy landscape analyses further indicate restricted conformational sampling of WEE2 upon ligand binding, supporting ligand-induced stabilization of the catalytic domain. MM-PBSA calculations confirmed favorable binding free energies for Midostaurin (−18.78 ± 2.23 kJ/mol) and Nilotinib (−17.47 ± 2.95 kJ/mol), exceeding that of MK1775. To increase the translational prioritization of candidate hits, we place our structure-based pipeline in the context of modern machine learning (ML) and deep learning (DL)-enabled virtual screening workflows. ML/DL rescoring and graph-based molecular property predictors can rapidly re-rank docking hits and estimate absorption, distribution, metabolism, excretion, and toxicity (ADMET) liabilities before in vitro evaluation.

Figures