Temporally Resolved and Interpretable Machine Learning Model of GPCR Conformational Transition

Identifying target-specific drugs remains a challenge in pharmacology, especially for highly homologous proteins such as dopamine receptors D2R and D3R. Differences in target-specific cryptic druggable sites for such receptors arise from the distinct conformational ensembles underlying their dynamic behavior. While Molecular Dynamics (MD) simulations has emerged as a powerful tool for dissecting protein dynamics, the sheer volume of MD data requires scalable and unbiased data analysis strategies to pinpoint residue communities regulating conformational state ensembles. We have developed the Dynamically Resolved Universal Model for BayEsiAn network Tracking (DRUMBEAT) interpretable machine learning algorithm and validated it by identifying residue communities that enable the deactivation of the {beta}2-adrenergic receptor. Further, upon analyzing dopamine receptor dynamics we identified distinct and non-conserved residue communities around the contacts F1704.62_F172ECL2 and S1464.38_G14134.56 that are specific to D3R conformational transitions compared to D2R. This information can be tapped to design subtype-specific drugs for neuropsychiatric and substance use disorders.