Structure-Based Virtual Screening for Ligands of G Protein-Coupled Receptors: What Can Molecular Docking Do for You?

Opportunities and challenges in drug discovery targeting the orphan receptor GPR12

G-protein-coupled receptor 12 (GPR12) is a brain-specific expression orphan G-protein-coupled receptor (oGPCR) that regulates various physiological processes. It is an emerging therapeutic target for central nervous system (CNS) disorders, including Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD), attention deficit hyperactivity disorder (ADHD), and schizophrenia, as well as other human diseases, such as cancer, obesity, and metabolic disorders. GPR12 remains a less extensively investigated oGPCR, particularly in terms of its biological functions, signaling pathways, and ligand discovery. The discovery of drug-like small-molecule modulators to probe the brain functions of GPR12 or to act as a potential drug candidates, as well as the identification of reliable biomarkers, are vital to elucidate the roles of this receptor in various human diseases and develop novel target-based therapeutics.

AI-based identification of therapeutic agents targeting GPCRs: introducing ligand type classifiers and systems biology

Identifying ligands targeting G protein coupled receptors (GPCRs) with novel chemotypes other than the physiological ligands is a challenge for in silico screening campaigns. Here we present an approach that identifies novel chemotype ligands by combining structural data with a random forest agonist/antagonist classifier and a signal-transduction kinetic model. As a test case, we apply this approach to identify novel antagonists of the human adenosine transmembrane receptor type 2A, an attractive target against Parkinson's disease and cancer. The identified antagonists were tested here in a radio ligand binding assay. Among those, we found a promising ligand whose chemotype differs significantly from all so-far reported antagonists, with a binding affinity of 310 ± 23.4 nM. Thus, our protocol emerges as a powerful approach to identify promising ligand candidates with novel chemotypes while preserving antagonistic potential and affinity in the nanomolar range.

A generalized protein-ligand scoring framework with balanced scoring, docking, ranking and screening powers

Applying machine learning algorithms to protein-ligand scoring functions has aroused widespread attention in recent years due to the high predictive accuracy and affordable computational cost. Nevertheless, most machine learning-based scoring functions are only applicable to a specific task, e.g., binding affinity prediction, binding pose prediction or virtual screening, suggesting that the development of a scoring function with balanced performance in all critical tasks remains a grand challenge. To this end, we propose a novel parameterization strategy by introducing an adjustable binding affinity term that represents the correlation between the predicted outcomes and experimental data into the training of mixture density network. The resulting residue-atom distance likelihood potential not only retains the superior docking and screening power over all the other state-of-the-art approaches, but also achieves a remarkable improvement in scoring and ranking performance. We emphatically explore the impacts of several key elements on prediction accuracy as well as the task preference, and demonstrate that the performance of scoring/ranking and docking/screening tasks of a certain model could be well balanced through an appropriate manner. Overall, our study highlights the potential utility of our innovative parameterization strategy as well as the resulting scoring framework in future structure-based drug design.

The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering

Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.

Phytoestrogen-derived multifunctional ligands for targeted therapy of breast cancer

Nano-targeted delivery systems have been widely used for breast tumor drug delivery. Estrogen receptors are considered to be significant drug delivery target receptors due to their overexpression in a variety of tumor cells. However, targeted ligands have a significant impact on the safety and effectiveness of active delivery systems, limiting the clinical transformation of nanoparticles. Phytoestrogens have shown good biosafety characteristics and some affinity with the estrogen receptor. In the present study, molecular docking was used to select tanshinone IIA (Tan IIA) among phytoestrogens as a target ligand to be used in nanodelivery systems with some modifications. Modified Tan IIA (Tan-NH2) showed a good biosafety profile and demonstrated tumor-targeting, anti-tumor and anti-tumor metastasis effects. Moreover, the ligand was utilized with the anti-tumor drug Dox-loaded mesoporous silica nanoparticles via chemical modification to generate a nanocomposite Tan-Dox-MSN. Tan-Dox-MSN had a uniform particle size, good dispersibility and high drug loading capacity. Validation experiments in vivo and in vitro showed that it also had a better targeting ability, anti-tumor effect and lower toxicity in normal organs. These results supported the idea that phytoestrogens with high affinity for the estrogen receptor could improve the therapeutic efficacy of nano-targeted delivery systems in breast tumors.

A novel signature incorporating lipid metabolism- and immune-related genes to predict the prognosis and immune landscape in hepatocellular carcinoma

Background

Liver hepatocellular carcinoma (LIHC) is a highly malignant tumor with high metastasis and recurrence rates. Due to the relation between lipid metabolism and the tumor immune microenvironment is constantly being elucidated, this work is carried out to produce a new prognostic gene signature that incorporates immune profiles and lipid metabolism of LIHC patients.

Methods

We used the "DEseq2" R package and the "Venn" R package to identify differentially expressed genes related to lipid metabolism (LRDGs) in LIHC. Additionally, we performed unsupervised clustering of LIHC patients based on LRDGs to identify their subgroups and immuno-infiltration and Gene Ontology (GO) enrichment analysis on the subgroups. Next, we employed multivariate, LASSO and univariate Cox regression analyses to determine variables and to create a prognostic profile on the basis of immune- and lipid metabolism-related differential genes (IRDGs and LRDGs). We separated patients into low- and high-risk groups in accordance with the best cut-off value of risk score. We conducted Decision Curve Analysis (DCA), Receiver Operating Characteristic curve analysis as a function of time as well as Survival Analysis to evaluate this signature's prognostic value. We incorporated the clinical characteristics of patients into the risk model to obtain a nomogram prognostic model. GEO14520 and ICGC-LIRI JP datasets were employed to externally confirm the accuracy and robustness of signature. The gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied for investigating the underlying mechanisms. Immune infiltration analysis was implemented to examine the differences in immune between both risk groups. Single-cell RNA sequencing (scRNA-SEQ) was utilized to characterize the genes that were involved in the distribution of signature and expression characteristics of different LIHC cell types. The patients' sensitivity in both risk groups to commonly used chemotherapeutic agents and semi-inhibitory concentrations (IC50) of the drugs was assessed using the GDSC database. On the basis of the differentially expressed genes (DEGs) in the two groups, the CMAP database was adopted for the prediction of potential small-molecule compounds. Small-molecule compounds were molecularly docked with prognostic markers. Lastly, we investigated the prognostic gene expression levels in normal and LIHC tissues with immunohistochemistry (IHC) and quantitative reverse transcription polymerase chain reaction(qRT-PCR).

Results

We built and verified a prognostic signature with seven genes that incorporated immune profiles and lipid metabolism. Patients were classified as low- and high-risk groups depending on their prognostic profiles. The overall survival (OS) was markedly lower in the high-risk group as compared to low-risk group. Time-dependent ROC curves more precisely predicted patients' survival at 1, 3 and 5 years; the area under the ROC curve was 0.81 (1 year), 0.75 (3 years) and 0.77 (5 years). The DCA curves showed the value of the prognostic genes in this signature for clinical applications. We included the patients' clinical characteristics in the risk model for both multivariate and univariate Cox regression analyses, and the findings revealed that the risk model represents an independent factor that influences OS in LIHC patients. With immune analysis, GSVA and GSEA, we identified that there are remarkable differences between the two risk groups in immune pathways, lipid metabolism, tumor development, immune cell infiltration and immune microenvironment, response to immunotherapy, and sensitivity to chemotherapy. Moreover, those with higher risk scores presented greater sensitivity to the chemotherapeutic agents. Experiments in vitro further elucidated the roles of SPP1 and FLT3 in the LIHC immune microenvironment. Furthermore, four small-molecule drugs that could target LIHC were screened. In vitro qRT-PCR , IHC revealed that the SPP1,KIF18A expressions were raised in LIHC in tumor samples, whereas FLT3,SOCS2 showed the opposite trend.

Conclusions

We developed and verified a new signature comprising immune- and lipid metabolism-associated markers and to assess the prognosis and the immune status of LIHC patients. This signature can be applied to survival prediction, individualized chemotherapy, and immunotherapeutic guidance for patients with liver cancer. This study also provides potential targeted therapeutics and novel ideas for the immune evasion and progression of LIHC.

Metadynamics simulations of ligands binding to protein surfaces: a novel tool for rational drug design

Structure-based drug design protocols may encounter difficulties to investigate poses when the biomolecular targets do not exhibit typical binding pockets. In this study, by providing two concrete examples from our labs, we suggest that the combination of metadynamics free energy methods (validated against affinity measurements), along with experimental structural information (by X-ray crystallography and NMR), can help to identify the poses of ligands on protein surfaces. The simulation workflow proposed here was implemented in a widely used code, namely GROMACS, and it could straightforwardly be applied to various drug-design campaigns targeting ligands' binding to protein surfaces.

Structure-based discovery of conformationally selective inhibitors of the serotonin transporter

The serotonin transporter (SERT) removes synaptic serotonin and is the target of anti-depressant drugs. SERT adopts three conformations: outward-open, occluded, and inward-open. All known inhibitors target the outward-open state except ibogaine, which has unusual anti-depressant and substance-withdrawal effects, and stabilizes the inward-open conformation. Unfortunately, ibogaine's promiscuity and cardiotoxicity limit the understanding of inward-open state ligands. We docked over 200 million small molecules against the inward-open state of the SERT. Thirty-six top-ranking compounds were synthesized, and thirteen inhibited; further structure-based optimization led to the selection of two potent (low nanomolar) inhibitors. These stabilized an outward-closed state of the SERT with little activity against common off-targets. A cryo-EM structure of one of these bound to the SERT confirmed the predicted geometry. In mouse behavioral assays, both compounds had anxiolytic- and anti-depressant-like activity, with potencies up to 200-fold better than fluoxetine (Prozac), and one substantially reversed morphine withdrawal effects.

Allosteric modulators of solute carrier function: a theoretical framework

Large-scale drug screening is currently the basis for the identification of new chemical entities. This is a rather laborious approach, because a large number of compounds must be tested to cover the chemical space in an unbiased fashion. However, the structures of targetable proteins have become increasingly available. Thus, a new era has arguably been ushered in with the advent of methods, which allow for structure-based docking campaigns (i.e., virtual screens). Solute carriers (SLCs) are among the most promising drug targets. This claim is substantiated by the fact that a large fraction of the 400 solute carrier genes is associated with human diseases. The ability to dock large ligand libraries into selected structures of solute carriers has set the stage for rational drug design. In the present study, we show that these structure-based approaches can be refined by taking into account how solute carriers operate. We specifically address the feasibility of targeting solute carriers with allosteric modulators, because their actions differ fundamentally from those of ligands, which bind to the substrate binding site. For the pertinent analysis we used transition state theory in conjunction with the linear free energy relationship (LFER). These provide the theoretical framework to understand how allosteric modulators affect solute carrier function.

Structure-based virtual screening discovers potent and selective adenosine A1 receptor antagonists

Development of subtype-selective leads is essential in drug discovery campaigns targeting G protein-coupled receptors (GPCRs). Herein, a structure-based virtual screening approach to rationally design subtype-selective ligands was applied to the A₁ and A_2A adenosine receptors (A₁R and A_2AR). Crystal structures of these closely related subtypes revealed a non-conserved subpocket in the binding sites that could be exploited to identify A₁R selective ligands. A library of 4.6 million compounds was screened computationally against both receptors using molecular docking and 20 A₁R selective ligands were predicted. Of these, seven antagonized the A₁R with micromolar activities and several compounds displayed slight selectivity for this subtype. Twenty-seven analogs of two discovered scaffolds were designed, resulting in antagonists with nanomolar potency and up to 76-fold A₁R-selectivity. Our results show the potential of structure-based virtual screening to guide discovery and optimization of subtype-selective ligands, which could facilitate the development of safer drugs.

Integrating structure-based approaches in generative molecular design

Generative molecular design for drug discovery and development has seen a recent resurgence promising to improve the efficiency of the design-make-test-analyse cycle; by computationally exploring much larger chemical spaces than traditional virtual screening techniques. However, most generative models thus far have only utilized small-molecule information to train and condition de novo molecule generators. Here, we instead focus on recent approaches that incorporate protein structure into de novo molecule optimization in an attempt to maximize the predicted on-target binding affinity of generated molecules. We summarize these structure integration principles into either distribution learning or goal-directed optimization and for each case whether the approach is protein structure-explicit or implicit with respect to the generative model. We discuss recent approaches in the context of this categorization and provide our perspective on the future direction of the field.

Computational approaches streamlining drug discovery

Computer-aided drug discovery has been around for decades, although the past few years have seen a tectonic shift towards embracing computational technologies in both academia and pharma. This shift is largely defined by the flood of data on ligand properties and binding to therapeutic targets and their 3D structures, abundant computing capacities and the advent of on-demand virtual libraries of drug-like small molecules in their billions. Taking full advantage of these resources requires fast computational methods for effective ligand screening. This includes structure-based virtual screening of gigascale chemical spaces, further facilitated by fast iterative screening approaches. Highly synergistic are developments in deep learning predictions of ligand properties and target activities in lieu of receptor structure. Here we review recent advances in ligand discovery technologies, their potential for reshaping the whole process of drug discovery and development, as well as the challenges they encounter. We also discuss how the rapid identification of highly diverse, potent, target-selective and drug-like ligands to protein targets can democratize the drug discovery process, presenting new opportunities for the cost-effective development of safer and more effective small-molecule treatments.

Application of computational methods for class A GPCR Ligand discovery

G protein-coupled receptors (GPCR) are integral membrane proteins of considerable interest as targets for drug development due to their role in transmitting cellular signals in a multitude of biological processes. Of the six classes categorizing GPCR (A, B, C, D, E, and F), class A contains the largest number of therapeutically relevant GPCR. Despite their importance as drug targets, many challenges exist for the discovery of novel class A GPCR ligands serving as drug precursors. Though knowledge of the structural and functional characteristics of GPCR has grown significantly over the past 20 years, a large portion of GPCR lack reported, experimentally determined structures. Furthermore, many GPCR have no known endogenous and/or synthetic ligands, limiting further exploration of their biochemical, cellular, and physiological roles. While many successes in GPCR ligand discovery have resulted from experimental high-throughput screening, computational methods have played an increasingly important role in GPCR ligand identification in the past decade. Here we discuss computational techniques applied to GPCR ligand discovery. This review summarizes class A GPCR structure/function and provides an overview of many obstacles currently faced in GPCR ligand discovery. Furthermore, we discuss applications and recent successes of computational techniques used to predict GPCR structure as well as present a summary of ligand- and structure-based methods used to identify potential GPCR ligands. Finally, we discuss computational hit list generation and refinement and provide comprehensive workflows for GPCR ligand identification.

Modeling the Orthosteric Binding Site of the G Protein-Coupled Odorant Receptor OR5K1

With approximately 400 encoding genes in humans, odorant receptors (ORs) are the largest subfamily of class A G protein-coupled receptors (GPCRs). Despite its high relevance and representation, the odorant-GPCRome is structurally poorly characterized: no experimental structures are available, and the low sequence identity of ORs to experimentally solved GPCRs is a significant challenge for their modeling. Moreover, the receptive range of most ORs is unknown. The odorant receptor OR5K1 was recently and comprehensively characterized in terms of cognate agonists. Here, we report two additional agonists and functional data of the most potent compound on two mutants, L104^3.32 and L255^6.51. Experimental data was used to guide the investigation of the binding modes of OR5K1 ligands into the orthosteric binding site using structural information from AI-driven modeling, as recently released in the AlphaFold Protein Structure Database, and from homology modeling. Induced-fit docking simulations were used to sample the binding site conformational space for ensemble docking. Mutagenesis data guided side chain residue sampling and model selection. We obtained models that could better rationalize the different activity of active (agonist) versus inactive molecules with respect to starting models and also capture differences in activity related to minor structural differences. Therefore, we provide a model refinement protocol that can be applied to model the orthosteric binding site of ORs as well as that of GPCRs with low sequence identity to available templates.

A novel strategy for designing the magic shotguns for distantly related target pairs

Due to its promising capacity in improving drug efficacy, polypharmacology has emerged to be a new theme in the drug discovery of complex disease. In the process of novel multi-target drugs (MTDs) discovery, in silico strategies come to be quite essential for the advantage of high throughput and low cost. However, current researchers mostly aim at typical closely related target pairs. Because of the intricate pathogenesis networks of complex diseases, many distantly related targets are found to play crucial role in synergistic treatment. Therefore, an innovational method to develop drugs which could simultaneously target distantly related target pairs is of utmost importance. At the same time, reducing the false discovery rate in the design of MTDs remains to be the daunting technological difficulty. In this research, effective small molecule clustering in the positive dataset, together with a putative negative dataset generation strategy, was adopted in the process of model constructions. Through comprehensive assessment on 10 target pairs with hierarchical similarity-levels, the proposed strategy turned out to reduce the false discovery rate successfully. Constructed model types with much smaller numbers of inhibitor molecules gained considerable yields and showed better false-hit controllability than before. To further evaluate the generalization ability, an in-depth assessment of high-throughput virtual screening on ChEMBL database was conducted. As a result, this novel strategy could hierarchically improve the enrichment factors for each target pair (especially for those distantly related/unrelated target pairs), corresponding to target pair similarity-levels.

Characterization of dietary and herbal sourced natural compounds that modulate SEL1L-HRD1 ERAD activity and alleviate protein misfolding in the ER

Dysregulated production of peptide hormones is the key pathogenic factor of various endocrine diseases. Endoplasmic reticulum (ER) associated degradation (ERAD) is a critical machinery in maintaining ER proteostasis in mammalian cells by degrading misfolded proteins. Dysfunction of ERAD leads to maturation defect of many peptide hormones, such as provasopressin (proAVP), which results in the occurrence of Central Diabetes Insipidus. However, drugs targeting ERAD to regulate the production of peptide hormones are very limited. Herbal products provide not only nutritional sources, but also alternative therapeutics for chronic diseases. Virtual screening provides an effective and high-throughput strategy for identifying protein structure-based interacting compounds extracted from a variety of dietary or herbal sources, which could be served as (pro)drugs for preventing or treating endocrine diseases. Here, we performed a virtual screening by directly targeting SEL1L of the most conserved SEL1L-HRD1 ERAD machinery. Further, we analyzed 58 top-ranked compounds and demonstrated that Cryptochlorogenic acid (CCA) showed strong affinity with the binding pocket of SEL1L with HRD1. Through structure-based docking, protein expression assays, and FACS analysis, we revealed that CCA enhanced ERAD activity and promoted the degradation of misfolded proAVP, thus facilitated the secretion of well-folded proAVP. These results provide us with insights into drug discovery strategies targeting ER protein homeostasis, as well as candidate compounds for treating hormone-related diseases.

Construction of a novel choline metabolism-related signature to predict prognosis, immune landscape, and chemotherapy response in colon adenocarcinoma

BACKGROUND

Colon adenocarcinoma (COAD) is a common digestive system malignancy with high mortality and poor prognosis. Accumulating evidence indicates that choline metabolism is closely related to tumorigenesis and development. However, the efficacy of choline metabolism-related signature in predicting patient prognosis, immune microenvironment and chemotherapy response has not been fully clarified.

METHODS

Choline metabolism-related differentially expressed genes (DEGs) between normal and COAD tissues were screened using datasets from The Cancer Genome Atlas (TCGA), Kyoto Encyclopedia of Genes and Genomes (KEGG), AmiGO2 and Reactome Pathway databases. Two choline metabolism-related genes (CHKB and PEMT) were identified by univariate and multivariate Cox regression analyses. TCGA-COAD was the training cohort, and GSE17536 was the validation cohort. Patients in the high- and low-risk groups were distinguished according to the optimal cutoff value of the risk score. A nomogram was used to assess the prognostic accuracy of the choline metabolism-related signature. Calibration curves, decision curve analysis (DCA), and clinical impact curve (CIC) were used to improve the clinical applicability of the prognostic signature. Gene Ontology (GO) and KEGG pathway enrichment analyses of DEGs in the high- and low-risk groups were performed. KEGG cluster analysis was conducted by the KOBAS-i database. The distribution and expression of CHKB and PEMT in various types of immune cells were analyzed based on single-cell RNA sequencing (scRNA-seq). The CIBERSORT and ESTIMATE algorithms evaluated tumor immune cell infiltration in the high- and low-risk groups. Evaluation of the half maximal inhibitory concentration (IC50) of common chemotherapeutic drugs based on the choline metabolism-related signature was performed. Small molecule compounds were predicted using the Connectivity Map (CMap) database. Molecular docking is used to simulate the binding conformation of small molecule compounds and key targets. By immunohistochemistry (IHC), Western blot, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) experiments, the expression levels of CHKB and PEMT in human, mouse, and cell lines were detected.

RESULTS

We constructed and validated a choline metabolism-related signature containing two genes (CHKB and PEMT). The overall survival (OS) of patients in the high-risk group was significantly worse than that of patients in the low-risk group. The nomogram could effectively and accurately predict the OS of COAD patients at 1, 3, and 5 years. The DCA curve and CIC demonstrate the clinical utility of the nomogram. scRNA-seq showed that CHKB was mainly distributed in endothelial cells, while PEMT was mainly distributed in CD4+ T cells and CD8+ T cells. In addition, multiple types of immune cells expressing CHKB and PEMT differed significantly. There were significant differences in the immune microenvironment, immune checkpoint expression and chemotherapy response between the two risk groups. In addition, we screened five potential small molecule drugs that targeted treatment for COAD. Finally, the results of IHC, Western blot, and qRT-PCR consistently showed that the expression of CHKB in human, mouse, and cell lines was elevated in normal samples, while PMET showed the opposite trend.

CONCLUSION

In conclusion, we constructed a choline metabolism-related signature in COAD and revealed its potential application value in predicting the prognosis, immune microenvironment, and chemotherapy response of patients, which may lay an important theoretical basis for future personalized precision therapy.

Identification of Active Compounds against Melanoma Growth by Virtual Screening for Non-Classical Human DHFR Inhibitors

Antifolates such as methotrexate (MTX) have been largely known as anticancer agents because of their role in blocking nucleic acid synthesis and cell proliferation. Their mechanism of action lies in their ability to inhibit enzymes involved in the folic acid cycle, especially human dihydrofolate reductase (hDHFR). However, most of them have a classical structure that has proven ineffective against melanoma, and, therefore, inhibitors with a non-classical lipophilic structure are increasingly becoming an attractive alternative to circumvent this clinical resistance. In this study, we conducted a protocol combining virtual screening (VS) and cell-based assays to identify new potential non-classical hDHFR inhibitors. Among 173 hit compounds identified (average logP = 3.68; average MW = 378.34 Da), two-herein, called C1 and C2-exhibited activity against melanoma cell lines B16 and A375 by MTT and Trypan-Blue assays. C1 showed cell growth arrest (39% and 56%) and C2 showed potent cytotoxic activity (77% and 51%) in a dose-dependent manner. The effects of C2 on A375 cell viability were greater than MTX (98% vs 60%) at equivalent concentrations and times. Our results indicate that the integrated in silico/in vitro approach provided a benchmark to identify novel promising non-classical DHFR inhibitors showing activity against melanoma cells.

One class classification for the detection of β2 adrenergic receptor agonists using single-ligand dynamic interaction data

G protein-coupled receptors are involved in many biological processes, relaying the extracellular signal inside the cell. Signaling is regulated by the interactions between receptors and their ligands, it can be stimulated by agonists, or inhibited by antagonists or inverse agonists. The development of a new drug targeting a member of this family requires to take into account the pharmacological profile of the designed ligands in order to elicit the desired response. The structure-based virtual screening of chemical libraries may prioritize a specific class of ligands by combining docking results and ligand binding information provided by crystallographic structures. The performance of the method depends on the relevance of the structural data, in particular the conformation of the targeted site, the binding mode of the reference ligand, and the approach used to compare the interactions formed by the docked ligand with those formed by the reference ligand in the crystallographic structure. Here, we propose a new method based on the conformational dynamics of a single protein–ligand reference complex to improve the biased selection of ligands with specific pharmacological properties in a structure-based virtual screening exercise. Interactions patterns between a reference agonist and the receptor, here exemplified on the β2 adrenergic receptor, were extracted from molecular dynamics simulations of the agonist/receptor complex and encoded in graphs used to train a one-class machine learning classifier. Different conditions were tested: low to high affinity agonists, varying simulation duration, considering or ignoring hydrophobic contacts, and tuning of the classifier parametrization. The best models applied to post-process raw data from retrospective virtual screening obtained by docking of test libraries effectively filtered out irrelevant poses, discarding inactive and non-agonist ligands while identifying agonists. Taken together, our results suggest that consistency of the binding mode during the simulation is a key to the success of the method.

Mind the Gap-Deciphering GPCR Pharmacology Using 3D Pharmacophores and Artificial Intelligence

G protein-coupled receptors (GPCRs) are amongst the most pharmaceutically relevant and well-studied protein targets, yet unanswered questions in the field leave significant gaps in our understanding of their nuanced structure and function. Three-dimensional pharmacophore models are powerful computational tools in in silico drug discovery, presenting myriad opportunities for the integration of GPCR structural biology and cheminformatics. This review highlights success stories in the application of 3D pharmacophore modeling to de novo drug design, the discovery of biased and allosteric ligands, scaffold hopping, QSAR analysis, hit-to-lead optimization, GPCR de-orphanization, mechanistic understanding of GPCR pharmacology and the elucidation of ligand-receptor interactions. Furthermore, advances in the incorporation of dynamics and machine learning are highlighted. The review will analyze challenges in the field of GPCR drug discovery, detailing how 3D pharmacophore modeling can be used to address them. Finally, we will present opportunities afforded by 3D pharmacophore modeling in the advancement of our understanding and targeting of GPCRs.

Structure-Based Discovery of Negative Allosteric Modulators of the Metabotropic Glutamate Receptor 5

Recently determined structures of class C G protein-coupled receptors (GPCRs) revealed the location of allosteric binding sites and opened new opportunities for the discovery of novel modulators. In this work, molecular docking screens for allosteric modulators targeting the metabotropic glutamate receptor 5 (mGlu₅) were performed. The mGlu₅ receptor is activated by the main excitatory neurotransmitter of the nervous central system, L-glutamate, and mGlu₅ receptor activity can be allosterically modulated by negative or positive allosteric modulators. The mGlu₅ receptor is a promising target for the treatment of psychiatric and neurodegenerative diseases, and several allosteric modulators of this GPCR have been evaluated in clinical trials. Chemical libraries containing fragment- (1.6 million molecules) and lead-like (4.6 million molecules) compounds were docked to an allosteric binding site of mGlu₅ identified in X-ray crystal structures. Among the top-ranked compounds, 59 fragments and 59 lead-like compounds were selected for experimental evaluation. Of these, four fragment- and seven lead-like compounds were confirmed to bind to the allosteric site with affinities ranging from 0.43 to 8.6 μM, corresponding to a hit rate of 9%. The four compounds with the highest affinities were demonstrated to be negative allosteric modulators of mGlu₅ signaling in functional assays. The results demonstrate that virtual screens of fragment- and lead-like chemical libraries have complementary advantages and illustrate how access to high-resolution structures of GPCRs in complex with allosteric modulators can accelerate lead discovery.

Ethyl Acetate Extract of Dracocephalum heterophyllum Benth Ameliorates Nonalcoholic Steatohepatitis and Fibrosis via Regulating Bile Acid Metabolism, Oxidative Stress and Inhibiting Inflammation

Dracocephalum heterophyllum Benth is well-known for its ability to alleviate liver heat. In this study, Dracocephalum heterophyllum Benth ethyl acetate extracts were evaluated on mouse models of nonalcoholic steatohepatitis and liver fibrosis. After 6 and 8 weeks of treatment, serum parameters and gene expressions in tissue samples, as well as stained tissue sections, demonstrated that the ethyl acetate extracts were effective in treating these liver diseases. The principal bioactive constituent (rosmarinic acid) was identified and screened by high pressure liquid chromatography-1,1-diphenyl-2-picrylhydrazyl and affinity ultrafiltration-HPLC. The rosmarinic acid was separated from extracts with high purity by medium- and high-pressure liquid chromatography. Finally, the interactions between rosmarinic acid and the key targets of lipid metabolism, oxidative stress and inflammation were verified by molecular docking. Thereby, an indirect regulation of lipid and cholesterol metabolism and inhibition of liver inflammation and liver fibrosis by the studied extract has been observed. This study demonstrated that Dracocephalum heterophyllum Benth ethyl acetate extracts have the potential to treat nonalcoholic steatohepatitis and liver fibrosis, revealing their multi-target and multi-pathway therapeutic characteristics.

Therapeutic and prognostic potential of GPCRs in prostate cancer from multi-omics landscape

Prostate cancer (PRAD) is a common and fatal malignancy. It is difficult to manage clinically due to drug resistance and poor prognosis, thus creating an urgent need for novel therapeutic targets and prognostic biomarkers. Although G protein-coupled receptors (GPCRs) have been most attractive for drug development, there have been lack of an exhaustive assessment on GPCRs in PRAD like their molecular features, prognostic and therapeutic values. To close this gap, we herein systematically investigate multi-omics profiling for GPCRs in the primary PRAD by analyzing somatic mutations, somatic copy-number alterations (SCNAs), DNA methylation and mRNA expression. GPCRs exhibit low expression levels and mutation frequencies while SCNAs are more prevalent. 46 and 255 disease-related GPCRs are identified by the mRNA expression and DNA methylation analysis, respectively, complementing information lack in the genome analysis. In addition, the genomic alterations do not exhibit an observable correlation with the GPCR expression, reflecting the complex regulatory processes from DNA to RNA. Conversely, a tight association is observed between the DNA methylation and mRNA expression. The virtual screening and molecular dynamics simulation further identify four potential drugs in repositioning to PRAD. The combination of 3 clinical characteristics and 26 GPCR molecular features revealed by the transcriptome and genome exhibit good performance in predicting progression-free survival in patients with the primary PRAD, providing candidates as new biomarkers. These observations from the multi-omics analysis on GPCRs provide new insights into the underlying mechanism of primary PRAD and potential of GPCRs in developing therapeutic strategies on PRAD.

Boosting Protein-Ligand Binding Pose Prediction and Virtual Screening Based on Residue-Atom Distance Likelihood Potential and Graph Transformer

The past few years have witnessed enormous progress toward applying machine learning approaches to the development of protein-ligand scoring functions. However, the robust performance and wide applicability of scoring functions remain a big challenge for increasing the success rate of docking-based virtual screening. Herein, a novel scoring function named RTMScore was developed by introducing a tailored residue-based graph representation strategy and several graph transformer layers for the learning of protein and ligand representations, followed by a mixture density network to obtain residue-atom distance likelihood potential. Our approach was resolutely validated on the CASF-2016 benchmark, and the results indicate that RTMScore can outperform almost all of the other state-of-the-art methods in terms of both the docking and screening powers. Further evaluation confirms the robustness of our approach that can not only retain its docking power on cross-docked poses but also achieve improved performance as a rescoring tool in larger-scale virtual screening.

Molecular Pharmacology of P2X Receptors: Exploring Druggable Domains Revealed by Structural Biology

Extracellular ATP is a critical signaling molecule that is found in a wide range of concentrations across cellular environments. The family of nonselective cation channels that sense extracellular ATP, termed P2X receptors (P2XRs), is composed of seven subtypes (P2X1-P2X7) that assemble as functional homotrimeric and heterotrimeric ion channels. Each P2XR is activated by a distinct concentration of extracellular ATP, spanning from high nanomolar to low millimolar. P2XRs are implicated in a variety of physiological and pathophysiological processes in the cardiovascular, immune, and central nervous systems, corresponding to the spatiotemporal expression, regulation, and activation of each subtype. The therapeutic potential of P2XRs is an emerging area of research in which structural biology has seemingly exceeded medicinal chemistry, as there are several published P2XR structures but currently no FDA-approved drugs targeting these ion channels. Cryogenic electron microscopy is ideally suited to facilitate structure-based drug design for P2XRs by revealing and characterizing novel ligand-binding sites. This review covers structural elements in P2XRs including the extracellular orthosteric ATP-binding site, extracellular allosteric modulator sites, channel pore, and cytoplasmic substructures, with an emphasis on potential therapeutic ligand development.

Chenodeoxycholic acid inhibits lung adenocarcinoma progression via the integrin α5β1/FAK/p53 signaling pathway

BACKGROUND

Lung cancer is the leading cause of cancer-associated death worldwide and is classified into non-small cell lung cancer (NSCLC) and small-cell lung cancer (SCLC). NSCLC accounts for approximately 80%-85% of all lung cancer cases. Chenodeoxycholic acid (CDCA), a primary bile acid, has been reported to inhibit carcinoma cell proliferation. Here, we aimed to determine the effects and mechanism of action of CDCA against lung adenocarcinoma (LUAD).

METHODS

Western blotting and quantitative real-time polymerase chain reaction were used to evaluate the protein and mRNA expression levels in LUAD cell lines, respectively. Cell Counting Kit-8 and clone formation assays were performed to evaluate the proliferation ability of different cell types in vitro. Tumor cell motility was evaluated using Transwell assays. The transcriptional profile of A549 cells treated with CDCA was determined through RNA sequencing analysis. A xenograft model was established to evaluate the effects of CDCA on LUAD progression in vivo.

RESULTS

CDCA inhibited LUAD cell proliferation, migration, and invasion. Furthermore, it promoted apoptosis in LUAD cells. Mechanistically, CDCA inhibited the integrin α5β1 signaling pathway in LUAD cells by inhibiting the expression of the α5 and β1 subunits of integrin and phosphorylated FAK. Moreover, CDCA induced an increase in the levels of p53, a downstream gene of the integrin α5β1/FAK pathway. In addition, CDCA significantly decreased tumor volume in mice without inducing significant toxicity.

CONCLUSIONS

Our findings indicate that CDCA attenuates LUAD pathogenesis in vitro and in vivo via the integrin α5β1/FAK/p53 axis.

Sodium or Not Sodium: Should Its Presence Affect the Accuracy of Pose Prediction in Docking GPCR Antagonists?

The function of the allosteric sodium ion in stabilizing the inactive form of GPCRs has been extensively described in the past decades. Its presence has been reported to be essential for the binding of antagonist molecules in the orthosteric site of these very important therapeutical targets. Among the GPCR–antagonist crystal structures available, in most cases, the sodium ion could not be experimentally resolved, obliging computational scientists using GPCRs as targets for virtual screening to ask: “Should the sodium ion affect the accuracy of pose prediction in docking GPCR antagonists?” In the present study, we examined the performance of three orthogonal docking programs in the self-docking of GPCR antagonists to try to answer this question. The results of the present work highlight that if the sodium ion is resolved in the crystal structure used as the target, it should also be taken into account during the docking calculations. If the crystallographic studies were not able to resolve the sodium ion then no advantage would be obtained if this is manually inserted in the virtual target. The outcomes of the present analysis are useful for researchers exploiting molecular docking-based virtual screening to efficiently identify novel GPCR antagonists.

The Impact of the Secondary Binding Pocket on the Pharmacology of Class A GPCRs

G-protein coupled receptors (GPCRs) are considered important therapeutic targets due to their pathophysiological significance and pharmacological relevance. Class A receptors represent the largest group of GPCRs that gives the highest number of validated drug targets. Endogenous ligands bind to the orthosteric binding pocket (OBP) embedded in the intrahelical space of the receptor. During the last 10 years, however, it has been turned out that in many receptors there is secondary binding pocket (SBP) located in the extracellular vestibule that is much less conserved. In some cases, it serves as a stable allosteric site harbouring allosteric ligands that modulate the pharmacology of orthosteric binders. In other cases it is used by bitopic compounds occupying both the OBP and SBP. In these terms, SBP binding moieties might influence the pharmacology of the bitopic ligands. Together with others, our research group showed that SBP binders contribute significantly to the affinity, selectivity, functional activity, functional selectivity and binding kinetics of bitopic ligands. Based on these observations we developed a structure-based protocol for designing bitopic compounds with desired pharmacological profile.