Mechanistic Insights into Ligand-Mediated Modulation of 14-3-3η/Parkin Interaction through Molecular Dynamics and...

Parkinson's disease (PD) is a complex neurological disorder driven by the dysfunction of key protein components. In this regard, 14-3-3η, a ubiquitous protein-protein interaction (PPI) facilitator, complexes with Parkin, an E3 ligase, inactivating the latter. This interaction promotes oxidative stress in dysfunctional mitochondria, ultimately causing neuronal cell death. The current study employs various in silico methods to investigate this complexation, identifying the orthosteric pocket and predicting small-molecule inhibitors to modulate their interaction. To achieve this, multiple independent 1 μs classical molecular dynamics (MD) simulations, combined with enhanced sampling protocols, were performed alongside contact map and alanine-scanning studies to identify the probable orthosteric pocket. A virtual screening of 446,736 small molecules was conducted, from which five candidates were shortlisted based on their docking scores and their Adsorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties. Structural stability analyses and Molecular Mechanics Generalized Born Surface Area (MM/GBSA) calculations with entropy correction identified four potential PPI inhibitors with their absolute binding free energies (ABFE) estimated using the Streamlined Alchemical Free Energy Perturbation (SAFEP) protocol. Notably, Lig-2 and Lig-4 triggered the opening of the 14-3-3η binding pocket, as supported by shortest path map (SPM) analysis, where reduced edge weights between terminal residues (E16 and L222) indicated weakened information transfer across the protein network. Umbrella sampling demonstrated a significant reduction in protein-protein unbinding potential of mean force (PMF) for Lig-1, Lig-2, and Lig-4, compared to that in the apo state. Uniform manifold approximation and projection (UMAP) and Density-Based Spatial Clustering of Applications with Noise(DBSCAN) clustering identified intermediate conformations linking the pocket opening to an increased 14-3-3η/Parkin center of mass (COM) distance. Finally, dynamic cross correlation maps (DCCMs) indicated diminished intra- and interprotein communication in the presence of orthosteric ligands, suggesting rigidification of both proteins. Overall, this study offers comprehensive insight into the 14-3-3η/Parkin complexation mechanism and identifies Lig-2 and Lig-4 as promising orthosteric inhibitors capable of reducing their binding affinity, presenting a potential therapeutic strategy for PD.