Annotation of Allosteric Compounds to Enhance Bioactivity Modeling for Class A GPCRs

Comparative Study of Allosteric GPCR Binding Sites and Their Ligandability Potential

The steadily growing number of experimental G-protein-coupled receptor (GPCR) structures has revealed diverse locations of allosteric modulation, and yet few drugs target them. This gap highlights the need for a deeper understanding of allosteric modulation in GPCR drug discovery. The current work introduces a systematic annotation scheme to structurally classify GPCR binding sites based on receptor class, transmembrane helix contacts, and, for membrane-facing sites, membrane sublocation. This GPCR specific annotation scheme was applied to 107 GPCR structures bound by small molecules contributing to 24 distinct allosteric binding sites for comparative evaluation of three binding site detection methods (BioGPS, SiteMap, and FTMap). BioGPS identified the most in 22 of 24 sites. In addition, our property analysis showed that extrahelical allosteric ligands and binding sites represent a distinct chemical space characterized by shallow pockets with low volume, and the corresponding allosteric ligands showed an enrichment of halogens. Furthermore, we demonstrated that combining receptor and ligand similarity can be a viable method for ligandability assessment. One challenge regarding site prediction is the ligand shaping effect on the observed binding site, especially for extrahelical sites where the ligand-induced effect was most pronounced. To our knowledge, this is the first study presenting a binding site annotation scheme standardized for GPCRs, and it allows a comparison of allosteric binding sites across different receptors in an objective way. The insight from this study provides a framework for future GPCR binding site studies and highlights the potential of targeting allosteric sites for drug development.

Integrating Molecular Dynamics and Machine Learning Algorithms to Predict the Functional Profile of Kinase Ligands

The modulation of protein function via designed small molecules is providing new opportunities in chemical biology and medicinal chemistry. While drugs have traditionally been developed to block enzymatic activities through active site occupation, a growing number of strategies now aim to control protein functions in an allosteric fashion, allowing for the tuning of a target's activation or deactivation via the modulation of the populations of conformational ensembles that underlie its function. In the context of the discovery of new active leads, it would be very useful to generate hypotheses for the functional impact of new ligands. Since the discovery and design of allosteric modulators (inhibitors/activators) is still a challenging and often serendipitous target, the development of a rapid and robust approach to predict the functional profile of a new ligand would significantly speed up candidate selection. Herein, we present different machine learning (ML) classifiers to distinguish between potential orthosteric and allosteric binders. Our approach integrates information on the chemical fingerprints of the ligands with descriptors that recapitulate ligand effects on protein functional motions. The latter are derived from molecular dynamics (MD) simulations of the target protein in complex with orthosteric or allosteric ligands. In this framework, we train and test different ML architectures, which are initially probed on the classification of orthosteric versus allosteric ligands for cyclin-dependent kinases (CDKs). The results demonstrate that different ML methods can successfully partition allosteric versus orthosteric effectors (although to different degrees). Next, we further test the models with FDA-approved CDK drugs, not included in the original dataset, as well as ligands that target other kinases, to test the range of applicability of these models outside of the domain on which they were developed. Overall, the results show that enriching the training dataset with chemical physics-based information on the protein-ligand dynamic cross-talk can significantly expand the reach and applicability of approaches for the prediction and classification of the mode of action of small molecules.

Allostery Illuminated: Harnessing AI and Machine Learning for Drug Discovery

In the past several years there has been rapid adoption of artificial intelligence (AI) and machine learning (ML) tools for drug discovery. In this Microperspective, we comment on recent AI/ML applications to the discovery of allosteric modulators, focusing on breakthroughs with AlphaFold, structure-based drug discovery (SBDD), and medicinal chemistry applications. We discuss how these technologies are facilitating drug discovery and the remaining challenges to identify allosteric binding sites and ligands.

Prospective de novo drug design with deep interactome learning

De novo drug design aims to generate molecules from scratch that possess specific chemical and pharmacological properties. We present a computational approach utilizing interactome-based deep learning for ligand- and structure-based generation of drug-like molecules. This method capitalizes on the unique strengths of both graph neural networks and chemical language models, offering an alternative to the need for application-specific reinforcement, transfer, or few-shot learning. It enables the "zero-shot" construction of compound libraries tailored to possess specific bioactivity, synthesizability, and structural novelty. In order to proactively evaluate the deep interactome learning framework for protein structure-based drug design, potential new ligands targeting the binding site of the human peroxisome proliferator-activated receptor (PPAR) subtype gamma are generated. The top-ranking designs are chemically synthesized and computationally, biophysically, and biochemically characterized. Potent PPAR partial agonists are identified, demonstrating favorable activity and the desired selectivity profiles for both nuclear receptors and off-target interactions. Crystal structure determination of the ligand-receptor complex confirms the anticipated binding mode. This successful outcome positively advocates interactome-based de novo design for application in bioorganic and medicinal chemistry, enabling the creation of innovative bioactive molecules.

Ligand-Receptor Interactions and Machine Learning in GCGR and GLP-1R Drug Discovery

The large amount of data that has been collected so far for G protein-coupled receptors requires machine learning (ML) approaches to fully exploit its potential. Our previous ML model based on gradient boosting used for prediction of drug affinity and selectivity for a receptor subtype was compared with explicit information on ligand-receptor interactions from induced-fit docking. Both methods have proved their usefulness in drug response predictions. Yet, their successful combination still requires allosteric/orthosteric assignment of ligands from datasets. Our ligand datasets included activities of two members of the secretin receptor family: GCGR and GLP-1R. Simultaneous activation of two or three receptors of this family by dual or triple agonists is not a typical kind of information included in compound databases. A precise allosteric/orthosteric ligand assignment requires a continuous update based on new structural and biological data. This data incompleteness remains the main obstacle for current ML methods applied to class B GPCR drug discovery. Even so, for these two class B receptors, our ligand-based ML model demonstrated high accuracy (5-fold cross-validation Q² > 0.63 and Q² > 0.67 for GLP-1R and GCGR, respectively). In addition, we performed a ligand annotation using recent cryogenic-electron microscopy (cryo-EM) and X-ray crystallographic data on small-molecule complexes of GCGR and GLP-1R. As a result, we assigned GLP-1R and GCGR actives deposited in ChEMBL to four small-molecule binding sites occupied by positive and negative allosteric modulators and a full agonist. Annotated compounds were added to our recently released repository of GPCR data.