Virtual Screening and Molecular Dynamics Simulation to Identify Inhibitors of the m6A-RNA Reader Protein YTHDC1

Abstract

YTHDC1 (YTH domain containing 1), a crucial reader protein of N6-methyladenosine (m6A) mRNA, plays a critical role in various cellular functions and is considered a promising target for therapeutic intervention in acute myeloid leukemia and other cancers. In this study, we identified orthosteric small-molecule ligands for YTHDC1. Using a molecular docking approach, we screened the eMolecules database and recognized 15 top-ranked ligands. Subsequently, molecular dynamics simulations and MM/PBSA analysis were used to assess the stability and binding free energy of these potential hit compounds in complex with YTHDC1. Notably, five compounds with IDs of ZINC82121447, ZINC02170552, ZINC65274016, ZINC10763862, and ZINC02412146 exhibited high binding affinities and favorable binding free energies. The results also showed that these compounds formed strong hydrogen bonds with residues SER378, ASN363, and ASN367 and interacted with the aromatic cage of the YTHDC1 reader protein through TRP377, TRP428, and hydrophobic residue LEU439. To assess their viability as lead compounds, we conducted absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies to reveal promising features for these identified small molecules, shedding light on their pharmacokinetic and safety profiles.