3D Pharmacophore Research Articles

Discovery of Potential Inhibitors of CDK1 by Integrating Pharmacophore-Based Virtual Screening, Molecular Docking, Molecular Dynamics Simulation Studies, and Evaluation of Their Inhibitory Activity.

The ability of CDK1 to compensate for the absence of other cell cycle CDKs poses a great challenge to treat cancers that overexpress these proteins. Despite several studies focusing on the area, there are no FDA-approved drugs selectively targeting CDK1. Here, the study aimed to develop potential CDK1 selective inhibitors through drug repurposing and leveraging the structural insights provided by the hit molecules generated. Approximately 280,000 compounds from DrugBank, Selleckchem, Otava and an in-house library were screened initially based on fit values using 3D QSAR pharmacophores built for CDK1 and subsequently through Lipinski, ADMET, and TOPKAT filters. 10,310 hits were investigated for docking into the binding site of CDK1 determined using the crystal structure of human CDK1 in complex with NU6102. The best 55 hits with better docking scores were further analyzed, and 12 hits were selected for 100 ns MD simulations followed by binding energy calculations using the MM-PBSA method. Finally, 10 hit molecules were tested in an in vitro CDK1 Kinase inhibition assay. Out of these, 3 hits showed significant CDK1 inhibitory potential with IC50 < 5 μM. These results indicate these compounds can be used to develop subtype-selective CDK1 inhibitors with better efficacy and reduced toxicities in the future.

Identification of Novel Bromodomain-Containing Protein 4 (BRD4) Binders through 3D Pharmacophore-Based Repositioning Screening Campaign.

A 3D structure-based pharmacophore model built for bromodomain-containing protein 4 (BRD4) is reported here, specifically developed for investigating and identifying the key structural features of the (+)-JQ1 known inhibitor within the BRD4 binding site. Using this pharmacophore model, 273 synthesized and purchased compounds previously considered for other targets but yielding poor results were screened in a drug repositioning campaign. Subsequently, only six compounds showed potential as BRD4 binders and were subjected to further biophysical and biochemical assays. Compounds 2, 5, and 6 showed high affinity for BRD4, with IC50 values of 0.60 ± 0.25 µM, 3.46 ± 1.22 µM, and 4.66 ± 0.52 µM, respectively. Additionally, these compounds were tested against two other bromodomains, BRD3 and BRD9, and two of them showed high selectivity for BRD4. The reported 3D structure-based pharmacophore model proves to be a straightforward and useful tool for selecting novel BRD4 ligands.

SmartCADD: AI-QM Empowered Drug Discovery Platform with Explainability.

Artificial intelligence (AI) has emerged as a pivotal force in enhancing productivity across various sectors, with its impact being profoundly felt within the pharmaceutical and biotechnology domains. Despite AI's rapid adoption, its integration into scientific research faces resistance due to myriad challenges: the opaqueness of AI models, the intricate nature of their implementation, and the issue of data scarcity. In response to these impediments, we introduce SmartCADD, an innovative, open-source virtual screening platform that combines deep learning, computer-aided drug design (CADD), and quantum mechanics methodologies within a user-friendly Python framework. SmartCADD is engineered to streamline the construction of comprehensive virtual screening workflows that incorporate a variety of formerly independent techniques─spanning ADMET property predictions, de novo 2D and 3D pharmacophore modeling, molecular docking, to the integration of explainable AI mechanisms. This manuscript highlights the foundational principles, key functionalities, and the unique integrative approach of SmartCADD. Furthermore, we demonstrate its efficacy through a case study focused on the identification of promising lead compounds for HIV inhibition. By democratizing access to advanced AI and quantum mechanics tools, SmartCADD stands as a catalyst for progress in pharmaceutical research and development, heralding a new era of innovation and efficiency.

3D-QSAR, Pharmacophore Modeling, ADMET, and DFT Studies of Halogenated Conjugated Dienones as Potent MAO-B Inhibitors.

It has been reported that the extension of conjugation in chalcone scaffolds considerably enhanced the potency, selectivity, reversibility, and competitive mode of MAO-B inhibition. In this study, using the experimental results of IC50 values of fifteen halogenated conjugated dienone derivatives (MK1-MK15) against MAO-B, we developed a 3DQSAR model. Further, we created a 3D pharmacophore model in active compounds in the series. The built model selected three variables (G2U, RDF115m, RDF155m) among the 653 AlvaDesc molecular descriptors, with a r2 value of 0.87 and a Q2 cv for cross-validation equal to 0.82. The three variables were mostly associated with the direction of symmetry and the likelihood of discovering massive atoms at great distances. The evaluated molecules exhibited a good correlation between experimental and predicted data, indicating that the IC50 value of the structure MK2 was related to the interatomic distances of 15.5 Å between bromine and chloro substituents. Furthermore, the molecules in the series with the highest activity were those with enhanced second component symmetry directional index from the 3D representation, which included the structures MK5 and MK6. Additionally, a pharmacophore hypothesis was developed and validated using the decoy Schrodinger dataset, with an ROC score of 0.87 and an HHRR 1 fitness score that ranged from 2.783 to 3.00. The MK series exhibited a significant blood-brain barrier (BBB) permeability, according to exploratory analyses and in silico projections, and almost all analogues were expected to have strong BBB permeability. Further DFT research revealed that electrostatics were important in the interactions with MAO-B.

Lead Identification Through In Silico Studies: Targeting Acetylcholinesterase Enzyme Against Alzheimer's Disease.

In this work, we aimed to acquire the best potential small molecule for Alzheimer's disease (AD) treatment using different models in Biovia Discovery Studio to identify new potential inhibitors of acetylcholinesterase (AChE) via in silico studies. The prevalence of cognitive impairment-related neurodegenerative disorders, such as AD, has been observed to escalate rapidly. However, we still know little about the underlying functions, outcome predictors, or intervention targets causing AD. The objective of the study was to optimize and identify the lead compound to target AChE against Alzheimer's disease. Different in silico studies were employed, including the pharmacophore model, virtual screening, molecular docking, de novo evolution model, and molecular dynamics. The pharmacophoric features of AChE inhibitors were determined by ligand-based pharmacophore models and 3D QSAR pharmacophore generation. Further validation of the best pharmacophore model was done using the cost analysis method, Fischer's randomization method, and test set. The molecules that harmonized the best pharmacophore model with the estimated activity < 1 nM and ADMET parameters were filtered, and 12 molecules were subjected to molecular docking studies to obtain binding energy. 3vsp_EK8_1 secured the highest binding energy of 65.60 kcal/mol. Further optimization led to a 3v_Evo_4 molecule with a better binding energy of 70.17 kcal/mol. The molecule 3v_evo_4 was subjected to 100 ns molecular simulation compared to donepezil, which showed better stability at the binding site. A lead compound, 3v_Evo_4 molecule, was identified to inhibit AChE, and it could be further studied to develop as a drug with better efficacy than the existing available drugs for treating AD.

Discovery of Novel, Selective, and Nonbasic Agonists for the Kappa-Opioid Receptor Determined by Salvinorin A-Based Virtual Screening.

Modulating the kappa-opioid receptor (KOR) is a promising strategy for treating various human diseases. KOR agonists show potential for treating pain, pruritus, and epilepsy, while KOR antagonists show potential for treating depression, anxiety, and addiction. The diterpenoid Salvinorin A (SalA), a secondary metabolite of Salvia divinorum, is a potent and selective KOR agonist. Unlike typical opioids, SalA lacks a basic nitrogen, which encouraged us to search for nonbasic KOR ligands. Through structure-based virtual screening using 3D pharmacophore models based on the binding mode of SalA, we identified novel, nonbasic, potent, and selective KOR agonists. In vitro studies confirmed two virtual hits, SalA-VS-07 and SalA-VS-08, as highly selective for the KOR and showing G protein-biased KOR agonist activity. Both KOR ligands share a novel spiro-moiety and a nonbasic scaffold. Our findings provide novel starting points for developing therapeutics aimed at treating pain and other conditions in which KOR is a central player.

In Silico Exploration of Molecular Mechanisms for Inhibiting Inflammatory Responses by 3Н-Thiazolo[4,5-b]pyridin-2-one Derivatives.

Combined in silico strategy for molecular mechanisms exploration of a series 3H-thiazolo[4,5-b]pyridin-2-ones exhibiting strong anti-exudative action through QSAR analysis, molecular docking and pharmacophore modelling is reported. GA-ML technique was used for QSAR models generation with 2D autocorrelation descriptors. One- and two-parameter regressions revealed that certain structural patterns or heteroatoms contribute mutually to the anti-exudative activity potentiation. Possible action mechanisms were discovered through flexible docking simulations with cyclooxygenase pathway enzymes (COX-1, COX-2, mPGES-1). Docking results indicated the possibility of stable complexes formation with the effective docking scores and proper orientation of ligands within the enzymes active sites. Pharmacophore modelling was carried out using protein-ligand interaction fingerprints methodology. Two- and three-centre 3D pharmacophore queries were constructed. Their analysis indicated the functionality of bicyclic thiazolopyridine scaffold proved by the steric placement of heteroatoms in the corresponding pharmacophore centres.

DMXP: A De Novo Small-Molecule 3D Structure Predictor with Graph Attention Networks.

Generating the three-dimensional (3D) structure of small molecules is crucial in both structure- and ligand-based drug design. Structure-based drug design needs bioactive conformations of compounds for lead identification and optimization. Ligand-based drug design techniques, such as 3D shape similarity search, 3D pharmacophore model, 3D-QSAR, etc., all require high-quality small-molecule ligand conformations to obtain reliable results. Although predicting a small molecular bioactive conformer requires information from the receptor, a crystal structure of the molecule is a proper approximation to its bioactive conformer in a specific receptor because the binding pose of a small molecule in its receptor's binding pockets should be energetically close to the crystal structures. This study presents a de novo small molecular structure predictor (dMXP) with graph attention networks based on crystal data derived from the Cambridge Structural Database (CSD) combined with molecular electrostatic information calculated by density-functional theory (DFT). Two featuring strategies (topological and atomic partial change features) were employed to explore the relation between these features and the 3D crystal structure of a small molecule. These features were then assembled to construct the holistic 3D crystal structure of a molecule. Molecular graphs were encoded using a graph attention mechanism to deal with the issues of the inconsistencies of local substructures contributing to the entire molecular structure. The root-mean-square deviation (RMSDs) of approximately 80% dMXP predicted structures and the native binding poses within receptors are less than 2.0 Å.

Ensemble-Based Virtual Screening Led to the Discovery of Novel Lead Molecules as Potential NMBAs.

Neuromuscular blocking agents (NMBAs) are routinely used during anesthesia to relax skeletal muscle. Nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels; NMBAs can induce muscle paralysis by preventing the neurotransmitter acetylcholine (ACh) from binding to nAChRs situated on the postsynaptic membranes. Despite widespread efforts, it is still a great challenge to find new NMBAs since the introduction of cisatracurium in 1995. In this work, an effective ensemble-based virtual screening method, including molecular property filters, 3D pharmacophore model, and molecular docking, was applied to discover potential NMBAs from the ZINC15 database. The results showed that screened hit compounds had better docking scores than the reference compound d-tubocurarine. In order to further investigate the binding modes between the hit compounds and nAChRs at simulated physiological conditions, the molecular dynamics simulation was performed. Deep analysis of the simulation results revealed that ZINC257459695 can stably bind to nAChRs' active sites and interact with the key residue Asp165. The binding free energies were also calculated for the obtained hits using the MM/GBSA method. In silico ADMET calculations were performed to assess the pharmacokinetic properties of hit compounds in the human body. Overall, the identified ZINC257459695 may be a promising lead compound for developing new NMBAs as an adjunct to general anesthesia, necessitating further investigations.

Integration of the Butina algorithm and ensemble learning strategies for the advancement of a pharmacophore ligand-based model: an in silico investigation of apelin agonists.

Introduction: 3D pharmacophore models describe the ligand's chemical interactions in their bioactive conformation. They offer a simple but sophisticated approach to decipher the chemically encoded ligand information, making them a valuable tool in drug design. Methods: Our research summarized the key studies for applying 3D pharmacophore models in virtual screening for 6,944 compounds of APJ receptor agonists. Recent advances in clustering algorithms and ensemble methods have enabled classical pharmacophore modeling to evolve into more flexible and knowledge-driven techniques. Butina clustering categorizes molecules based on their structural similarity (indicated by the Tanimoto coefficient) to create a structurally diverse training dataset. The learning method combines various individual pharmacophore models into a set of pharmacophore models for pharmacophore space optimization in virtual screening. Results: This approach was evaluated on Apelin datasets and afforded good screening performance, as proven by Receiver Operating Characteristic (AUC score of 0.994 ± 0.007), enrichment factor of (EF1% of 50.07 ± 0.211), Güner-Henry score of 0.956 ± 0.015, and F-measure of 0.911 ± 0.031. Discussion: Although one of the high-scoring models achieved statistically superior results in each dataset (AUC of 0.82; an EF1% of 19.466; GH of 0.131 and F1-score of 0.071), the ensemble learning method including voting and stacking method balanced the shortcomings of each model and passed with close performance measures.

The discovery of a novel IκB kinase β inhibitor based on pharmacophore modeling, virtual screening and biological evaluation.

Background: IκB kinase β (IKKβ) plays a pivotal role in the NF-κB signaling pathway and is considered a promising therapeutic target for various diseases. Materials & methods: The authors developed and validated a 3D pharmacophore model of IKKβ inhibitors via the HypoGen algorithm in Discovery Studio 2019, then performed virtual screening, molecular docking and kinase assays to identify hit compounds from the ChemDiv database. The compound with the highest inhibitory activity was further evaluated in adjuvant-induced arthritis rat models. Results: Among the four hit compounds, Hit 4 had the highest IKKβ inhibitory activity (IC50=30.4±3.8), and it could significantly ameliorate joint inflammation and damage in vivo. Conclusion: The identified compound, Hit 4, can be optimized as a therapeutic agent for inflammatory diseases.

Discovery of dual kinase inhibitors targeting VEGFR2 and FAK: structure-based pharmacophore modeling, virtual screening, and molecular docking studies

VEGFR2 and FAK signaling pathways are interconnected and have synergistic effects on tumor angiogenesis, growth, and metastasis. Thus, instead of the conventional targeting of each of these proteins individually with a specific inhibitor, the present work aimed to discover novel dual inhibitors targeting both VEGFR2 and FAK exploiting their association. To this end, receptor-based pharmacophore modeling technique was opted to generate 3D pharmacophore models for VEGFR2 and FAK type II kinase inhibitors. The generated pharmacophore models were validated by assessing their ability to discriminate between active and decoy compounds in a pre-compiled test set of VEGFR2 and FAK active compounds and decoys. ZINCPharmer web tool was then used to screen the ZINC database purchasable subset using the validated pharmacophore models retrieving 42,616 hits for VEGFR2 and 28,475 hits for FAK. Subsequently, they were filtered using various filters leaving 13,023 and 6,832 survived compounds for VEGFR2 and FAK, respectively, with 124 common compounds. Based on molecular docking simulations, thirteen compounds were found to satisfy all necessary interactions with VEGFR2 and FAK kinase domains. Thus, they are predicted to have a possible dual VEGFR2/FAK inhibitory activity. Finally, SwissADME web tool showed that compound ZINC09875266 is not only promising in terms of binding pattern to our target kinases, but also in terms of pharmacokinetic properties.

Integrating pharmacophore model and deep learning for activity prediction of molecules with BRCA1 gene.

In this paper, we propose a novel approach for predicting the activity/inactivity of molecules with the BRCA1 gene by combining pharmacophore modeling and deep learning techniques. Initially, we generated 3D pharmacophore fingerprints using a pharmacophore model, which captures the essential features and spatial arrangements critical for biological activity. These fingerprints served as informative representations of the molecular structures. Next, we employed deep learning algorithms to train a predictive model using the generated pharmacophore fingerprints. The deep learning model was designed to learn complex patterns and relationships between the pharmacophore features and the corresponding activity/inactivity labels of the molecules. By utilizing this integrated approach, we aimed to enhance the accuracy and efficiency of activity prediction. To validate the effectiveness of our approach, we conducted experiments using a dataset of known molecules with BRCA1 gene activity/inactivity from diverse sources. Our results demonstrated promising predictive performance, indicating the successful integration of pharmacophore modeling and deep learning. Furthermore, we utilized the trained model to predict the activity/inactivity of unknown molecules extracted from the ChEMBL database. The predictions obtained from the ChEMBL database were assessed and compared against experimentally determined values to evaluate the reliability and generalizability of our model. Overall, our proposed approach showcased significant potential in accurately predicting the activity/inactivity of molecules with the BRCA1 gene, thus enabling the identification of potential candidates for further investigation in drug discovery and development processes.

Clustering of atoms relative to vector space in the Z-matrix coordinate system and 'graphical fingerprint' analysis of 3D pharmacophore structure.

The behavior of a molecule within its environment is governed by chemical fields present in 3D space. However, beyond local descriptors in 3D, the conformations a molecule assumes, and the resulting clusters also play a role in influencing structure-activity models. This study focuses on the clustering of atoms according to the vector space of four atoms aligned in the Z-Matrix Reference system for molecular similarity. Using 3D-QSAR analysis, it was aimed to determine the pharmacophore groups as interaction points in the binding region of the β2-adrenoceptor target of fenoterol stereoisomers. Different types of local reactive descriptors of ligands have been used to elucidate points of interaction with the target. Activity values for ligand-receptor interaction energy were determined using the Levenberg-Marquardt algorithm. Using the Molecular Comparative Electron Topology method, the 3D pharmacophore model (3D-PhaM) was obtained after aligning and superimposing the molecules and was further validated by the molecular docking method. Best guesses were calculated with a non-output validation (LOO-CV) method. Finally, the data were calculated using the 'graphic fingerprint' technique. Based on the eLKlopman (Electrostatic LUMO Klopman) descriptor, the Q2 value of this derivative set was calculated as 0.981 and the R2ext value is calculated as 0.998.

Lead generation of UPPS inhibitors targeting MRSA: Using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking, and molecular dynamic simulations

Undecaprenyl Pyrophosphate Synthase (UPPS) is a vital target enzyme in the early stages of bacterial cell wall biosynthesis. UPPS inhibitors have antibacterial activity against resistant strains such as MRSA and VRE. In this study, we used several consecutive computer-based protocols to identify novel UPPS inhibitors. The 3D QSAR pharmacophore model generation (HypoGen algorithm) protocol was used to generate a valid predictive pharmacophore model using a set of UPPS inhibitors with known reported activity. The developed model consists of four pharmacophoric features: one hydrogen bond acceptor, two hydrophobic, and one aromatic ring. It had a correlation coefficient of 0.86 and a null cost difference of 191.39, reflecting its high predictive power. Hypo1 was proven to be statistically significant using Fischer’s randomization at a 95% confidence level. The validated pharmacophore model was used for the virtual screening of several databases. The resulting hits were filtered using SMART and Lipinski filters. The hits were docked into the binding site of the UPPS protein, affording 70 hits with higher docking affinities than the reference compound (6TC, − 21.17 kcal/mol). The top five hits were selected through extensive docking analysis and visual inspection based on docking affinities, fit values, and key residue interactions with the UPPS receptor. Moreover, molecular dynamic simulations of the top hits were performed to confirm the stability of the protein–ligand complexes, yielding five promising novel UPPS inhibitors.Graphical

Identification of novel inhibitors against VP40 protein of Marburg virus by integrating molecular modeling and dynamics approaches

Marburg virus (MV) is a highly etiological agent of haemorrhagic fever in humans and has spread across the world. Its outbreaks caused a 23–90% human death rate. However, there are currently no authorized preventive or curative measures yet. VP40 is the MV matrix protein, which builds protein shell underneath the viral envelope and confers hallmark filamentous. VP40 alone is able to induce assembly and budding of filamentous virus-like particles (VLPs), which resemble authentic virions. As a result, this research is credited with clarifying the function of VP40 and leading to the discovery of new therapeutic targets effective in combating MV disease (MVD). Virtual screening, molecular docking and molecular dynamics (MD) simulation were used to find the putative active chemicals based on a 3D pharmacophore model of the protein’s active site cavity. Initially, andrographidine-C, a potent inhibitor was selected for the development of the pharmacophore model. Later, a library of 30,000 compounds along with the andrographidine-C was docked against VP40 protein. Three best hits including avanafil, diuvaretin and macrourone were subjected to further MD simulation analysis, as these compounds had better binding affinities as compared to andrographidine-C. Furthermore, throughout the 100 ns simulations, the back bone of VP40 protein in presence of avanafil, diuvaretin and macrourone remained stable which was further validated by MM-PBSA analysis. Additionally, all of these compounds depict maximum drug-like properties. The predicted drugs based on the ligand, avanafil, diuvaretin and macrourone could be exploited and developed as an alternative or complementary therapy for the treatment of MVD. Communicated by Ramaswamy H. Sarma

Effect of the Electronic Structure of Para-Substituted Benzaldehyde Benzohydrazone on Its Antimicrobial Activity: A DFT Analysis

Bacillus subtilis is a bacterium that has demonstrated its efficacy across various domains, including industry, agriculture, and commerce, owing to its protective, inhibitory, and biological mechanisms against specific microbes. However, at high concentrations, it can lead to food poisoning and severe infections, resulting in symptoms such as diarrhea and vomiting. Bacterial spores produced by Bacillus subtilis can induce conditions like gas gangrene and tetanus. In this context, benzohydrazones are recognized for their antimicrobial activity, particularly against Bacillus subtilis. This study aims to elucidate the relationship between the electronic structure of para-substituted benzaldehyde benzohydrazone derivatives and their antimicrobial activity. This leads to the proposal of a 2D pharmacophore for predicting the antibacterial activity of these derivatives. The quantitative structure-activity relationship (QSAR) approach employed is the KPG method. The electronic structures were optimized using the density functional theory (DFT) method with the B3LYP functional and the 6-31G (d,p) basis set. Charge and local molecular orbitals were considered in the optimization process. The resulting prediction equation (R=98.95%, R²=97.91%, Adjusted R²=96.76%, F(5,9)=84.52) derived from multiple linear regression provides the basis for the proposed 2D pharmacophore. This equation shows that antimicrobial activity of benzohydrazone derivative is on charge and orbital controlled. This pharmacophore holds potential utility in designing new molecular structures with enhanced activity against Bacillus subtilis.

Virtual screening approach for the discovery of selective 5α-reductase type II inhibitors for benign prostatic hyperplasia treatment.

Background: 5α-Reductase type II (5αR2) inhibition is a promising strategy for benign prostatic hyperplasia treatment. A computational approach including virtual screening, ligand-based 3Dpharmacophore modeling, 2D quantitative structure-activity relationship and molecular docking simulations were adopted to develop novel inhibitors. Results: Hits were first filtered via the validated pharmacophore and 2D quantitative structure-activity relationship models. Docking on the recently determined cocrystallized structure of 5αR2 showed three promising hits. Visual inspection results were compared with finasteride ligand and dihydrotestosterone as reference, to explain the role of binding to Glu57 and Tyr91 for 5αR2 selective inhibition. Conclusion: Alignment between Hit 2 and finasteride in the binding pocket showed similar binding modes. The biological activity prediction showed antitumor and androgen targeting activity of the new hits.

Distinct inhibitory strength of bisphenol A analogues on human and rat 11β-hydroxysteroid dehydrogenase 1: 3D quantitative structure-activity relationship and in silico molecular docking analysis

Bisphenol A (BPA) analogues are developed to replace BPA usage. However, their effects on 11β-hydroxysteroid dehydrogenase 1 (11β-HSD1) are largely unknown. The inhibitory effects of BPA and 10 BPA analogues with the substituents on the bridge moiety on human and rat 11β-HSD1 were explored in human and rat liver microsomes. The strength of inhibiting human 11β-HSD1 was bisphenol FL (IC50, 3.87 μM) > bisphenol Z (6.86 μM) > bisphenol AF (9.42 μM) > bisphenol C (16.14 μM) > bisphenol AP (32.14 μM) = bisphenol B (32.34 μM) > 4,4′-thiodiphenol (67.35 μM) > BPA (297.35 μM) > other BPA analogues (ineffective at 100 μM). The strength of inhibiting rat 11β-HSD1 was bisphenol Z (IC50, 14.44 μM) > 4,4′-thiodiphenol (19.01 μM) > bisphenol B (20.13 μM) > bisphenol F (22.10 μM) > bisphenol E (33.04 μM) > bisphenol AF (49.67 μM) > bisphenol C > (56.97 μM) > bisphenol AP (62.71 μM) >bisphenol FL (96.31 μM) > other BPA analogues (ineffective at 100 μM). Bisphenol A, AF, AP, B, C, F, FL, Z, and 4,4′-thiodiphenol bind to the active sites of human and rat 11β-HSD1. Regression of LogP and molecular weight with IC50 values revealed distinct inhibitory pattern (negative correlation for human 11β-HSD1 vs. positive correlation for rat enzyme). Regression of the lowest binding energy with IC50 values revealed a significant positive regression. 3D QSAR pharmacophore analysis showed one hydrogen bond acceptor and two hydrogen bond donors for human 11β-HSD1. In conclusion, most BPA analogues are more potent inhibitors of human and rat 11β-HSD1 enzymes and there is structure-dependent and species-dependent inhibition.

Discovery of ICOS-Targeted Small Molecules Using Pharmacophore-Based Screening.

There are currently no small molecules clinically approved as immune checkpoint modulators. Besides possessing oral bioavailability, cell-penetrating capabilities and enhanced tumor penetration compared to monoclonal antibodies (mAbs), small molecules are amenable to pharmacokinetic optimization, which allows adopting flexible dosage regimens that may avoid immune-related adverse events associated with mAbs. The interaction of inducible co-stimulator (ICOS) with its ligand (ICOS-L) plays key roles in T-cell differentiation and activation of T-cell to B-cell functions. This study represents the development and validation of a virtual screening strategy to identify small molecules that bind a novel druggable binding pocket in human ICOS. We used a lipophilic canyon in the apo-structure of ICOS and the ICOS/ICOS-L interface individually as templates for molecular dynamics simulation to generate 3D pharmacophores subsequently used for virtual screening campaigns. Our strategy was successful finding a first-in-class small molecule ICOS binder (5P, KD value=108.08±26.76 μM) and validating biophysical screening platforms for ICOS-targeted small molecules. We anticipate that future structural optimization of 5P will result in the discovery of high affinity chemical ligands for ICOS.