In silico discovery of potential inhibitors targeting the MEIG1-PACRG complex for male contraceptive development

Abstract

The interaction between meiosis-expressed gene 1 (MEIG1) and Parkin co-regulated gene (PACRG) is a critical determinant of spermiogenesis, the process by which round spermatids mature into functional spermatozoa. Disruption of the MEIG1-PACRG complex can impair sperm development, highlighting its potential as a therapeutic target for addressing male infertility or for the development of non-hormonal contraceptive methods. This study used virtual screening, molecular docking, and molecular dynamics (MD) simulations to identify small molecule inhibitors targeting the MEIG1-PACRG interface. MD simulations provided representative protein conformations, which were used to virtually screen a library of over 800,000 compounds, resulting in 48 high-ranking candidates for each protein. PACRG emerged as a favorable target due to its flexible binding pockets and better docking scores compared to MEIG1. Key binding residues with compounds included W50, Y68, N70, and E74 on MEIG1, and K93, W96, E101, and H137 on PACRG. MD simulations revealed that compound stability in MEIG1 complexes is primarily maintained by hydrogen bonding with E74 and {pi}-{pi} stacking interactions with W50 and Y68. In PACRG complexes, compound stabilization is facilitated by hydrogen bonding with E101 and {pi}-{pi} interactions involving W96 and H137. These findings highlight distinct molecular determinants of ligand binding for each protein. Our work provides mechanistic insights and identifies promising compounds for further experimental validation, establishing a foundation for developing MEIG1-PACRG interaction inhibitors as male contraceptives.