Introduction to the BioChemical Library (BCL): An Application-Based Open-Source Toolkit for Integrated Cheminformatics and Machine Learning in Computer-Aided Drug Discovery

Advancements in Ligand-Based Virtual Screening through the Synergistic Integration of Graph Neural Networks and Expert-Crafted Descriptors

The fusion of traditional chemical descriptors with Graph Neural Networks (GNNs) offers a compelling strategy for enhancing ligand-based virtual screening methodologies. A comprehensive evaluation revealed that the benefits derived from this integrative strategy vary significantly among different GNNs. Specifically, while GCN and SchNet demonstrate pronounced improvements by incorporating descriptors, SphereNet exhibits only marginal enhancement. Intriguingly, despite SphereNet's modest gain, all three models-GCN, SchNet, and SphereNet-achieve comparable performance levels when leveraging this combination strategy. This observation underscores a pivotal insight: sophisticated GNN architectures may be substituted with simpler counterparts without sacrificing efficacy, provided that they are augmented with descriptors. Furthermore, our analysis reveals a set of expert-crafted descriptors' robustness in scaffold-split scenarios, frequently outperforming the combined GNN-descriptor models. Given the critical importance of scaffold splitting in accurately mimicking real-world drug discovery contexts, this finding accentuates an imperative for GNN researchers to innovate models that can adeptly navigate and predict within such frameworks. Our work not only validates the potential of integrating descriptors with GNNs in advancing ligand-based virtual screening but also illuminates pathways for future enhancements in model development and application. Our implementation can be found at https://github.com/meilerlab/gnn-descriptor.

Dirigent isoflavene-forming PsPTS2: 3D structure, stereochemical, and kinetic characterization comparison with pterocarpan-forming PsPTS1 homolog in pea

Pea phytoalexins (-)-maackiain and (+)-pisatin have opposite C6a/C11a configurations, but biosynthetically how this occurs is unknown. Pea dirigent-protein (DP) PsPTS2 generates 7,2'-dihydroxy-4',5'-methylenedioxyisoflav-3-ene (DMDIF), and stereoselectivity toward four possible 7,2'-dihydroxy-4',5'-methylenedioxyisoflavan-4-ol (DMDI) stereoisomers was investigated. Stereoisomer configurations were determined using NMR spectroscopy, electronic circular dichroism, and molecular orbital analyses. PsPTS2 efficiently converted cis-(3R,4R)-DMDI into DMDIF 20-fold faster than the trans-(3R,4S)-isomer. The 4R-configured substrate's near β-axial OH orientation significantly enhanced its leaving group abilities in generating A-ring mono-quinone methide (QM), whereas 4S-isomer's α-equatorial-OH was a poorer leaving group. Docking simulations indicated that the 4R-configured β-axial OH was closest to Asp51, whereas 4S-isomer's α-equatorial OH was further away. Neither cis-(3S,4S)- nor trans-(3S,4R)-DMDIs were substrates, even with the former having C3/C4 stereochemistry as in (+)-pisatin. PsPTS2 used cis-(3R,4R)-7,2'-dihydroxy-4'-methoxyisoflavan-4-ol [cis-(3R,4R)-DMI] and C3/C4 stereoisomers to give 2',7-dihydroxy-4'-methoxyisoflav-3-ene (DMIF). DP homologs may exist in licorice (Glycyrrhiza pallidiflora) and tree legume Bolusanthus speciosus, as DMIF occurs in both species. PsPTS1 utilized cis-(3R,4R)-DMDI to give (-)-maackiain 2200-fold more efficiently than with cis-(3R,4R)-DMI to give (-)-medicarpin. PsPTS1 also slowly converted trans-(3S,4R)-DMDI into (+)-maackiain, reflecting the better 4R configured OH leaving group. PsPTS2 and PsPTS1 provisionally provide the means to enable differing C6a and C11a configurations in (+)-pisatin and (-)-maackiain, via identical DP-engendered mono-QM bound intermediate generation, which PsPTS2 either re-aromatizes to give DMDIF or PsPTS1 intramolecularly cyclizes to afford (-)-maackiain. Substrate docking simulations using PsPTS2 and PsPTS1 indicate cis-(3R,4R)-DMDI binds in the anti-configuration in PsPTS2 to afford DMDIF, and the syn-configuration in PsPTS1 to give maackiain.

Structure-Activity Relationship Study of the High-Affinity Neuropeptide Y4 Receptor Positive Allosteric Modulator VU0506013

Positive allosteric modulators targeting the Y₄ receptor (Y₄R), a G protein-coupled receptor (GPCR) involved in the regulation of satiety, offer great potential in anti-obesity research. In this study, we selected 603 compounds by using quantitative structure-activity relationship (QSAR) models and tested them in high-throughput screening (HTS). Here, the novel positive allosteric modulator (PAM) VU0506013 was identified, which exhibits nanomolar affinity and pronounced selectivity toward the Y₄R in engineered cell lines and mouse descending colon mucosa natively expressing the Y₄R. Based on this lead structure, we conducted a systematic SAR study in two regions of the scaffold and presented a series of 27 analogues with modifications in the N- and C-terminal heterocycles of the molecule to obtain insight into functionally relevant positions. By mutagenesis and computational docking, we present a potential binding mode of VU0506013 in the transmembrane core of the Y₄R. VU0506013 presents a promising scaffold for developing in vivo tools to move toward anti-obesity drug research focused on the Y₄R.

Opening of capsaicin receptor TRPV1 is stabilized equally by its four subunits

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules, such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here, we report an approach to control the number of bound RTX molecules (0-4) in functional rat TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.70 to 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

Improving the Modeling of Extracellular Ligand Binding Pockets in RosettaGPCR for Conformational Selection

G protein-coupled receptors (GPCRs) are the largest class of drug targets and undergo substantial conformational changes in response to ligand binding. Despite recent progress in GPCR structure determination, static snapshots fail to reflect the conformational space of putative binding pocket geometries to which small molecule ligands can bind. In comparative modeling of GPCRs in the absence of a ligand, often a shrinking of the orthosteric binding pocket is observed. However, the exact prediction of the flexible orthosteric binding site is crucial for adequate structure-based drug discovery. In order to improve ligand docking and guide virtual screening experiments in computer-aided drug discovery, we developed RosettaGPCRPocketSize. The algorithm creates a conformational ensemble of biophysically realistic conformations of the GPCR binding pocket between the TM bundle, which is consistent with a knowledge base of expected pocket geometries. Specifically, tetrahedral volume restraints are defined based on information about critical residues in the orthosteric binding site and their experimentally observed range of Cα-Cα-distances. The output of RosettaGPCRPocketSize is an ensemble of binding pocket geometries that are filtered by energy to ensure biophysically probable arrangements, which can be used for docking simulations. In a benchmark set, pocket shrinkage observed in the default RosettaGPCR was reduced by up to 80% and the binding pocket volume range and geometric diversity were increased. Compared to models from four different GPCR homology model databases (RosettaGPCR, GPCR-Tasser, GPCR-SSFE, and GPCRdb), the here-created models showed more accurate volumes of the orthosteric pocket when evaluated with respect to the crystallographic reference structure. Furthermore, RosettaGPCRPocketSize was able to generate an improved realistic pocket distribution. However, while being superior to other homology models, the accuracy of generated model pockets was comparable to AlphaFold2 models. Furthermore, in a docking benchmark using small-molecule ligands with a higher molecular weight between 400 and 700 Da, a higher success rate in creating native-like binding poses was observed. In summary, RosettaGPCRPocketSize can generate GPCR models with realistic orthosteric pocket volumes, which are useful for structure-based drug discovery applications.

TRPV1 Opening is Stabilized Equally by Its Four Subunits

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here we report an approach to control the number of bound RTX molecules (0-to-4) in functional mouse TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

In-Silico Drug Toxicity and Interaction Prediction for Plant Complexes Based on Virtual Screening and Text Mining

Potential drug toxicities and drug interactions of redundant compounds of plant complexes may cause unexpected clinical responses or even severe adverse events. On the other hand, super-additivity of drug interactions between natural products and synthetic drugs may be utilized to gain better performance in disease management. Although without enough datasets for prediction model training, based on the SwissSimilarity and PubChem platforms, for the first time, a feasible workflow of prediction of both toxicity and drug interaction of plant complexes was built in this study. The optimal similarity score threshold for toxicity prediction of this system is 0.6171, based on an analysis of 20 different herbal medicines. From the PubChem database, 31 different sections of toxicity information such as "Acute Effects", "NIOSH Toxicity Data", "Interactions", "Hepatotoxicity", "Carcinogenicity", "Symptoms", and "Human Toxicity Values" sections have been retrieved, with dozens of active compounds predicted to exert potential toxicities. In Spatholobus suberectus Dunn (SSD), there are 9 out of 24 active compounds predicted to play synergistic effects on cancer management with various drugs or factors. The synergism between SSD, luteolin and docetaxel in the management of triple-negative breast cancer was proved by the combination index assay, synergy score detection assay, and xenograft model.