Scaffold hopping from synthetic RXR modulators by virtual screening and de novo design

Artifical intelligence: a virtual chemist for natural product drug discovery

Nature is full of a bundle of medicinal substances and its product perceived as a prerogative structure to collaborate with protein drug targets. The natural product's (NPs) structure heterogeneity and eccentric characteristics inspired scientists to work on natural product-inspired medicine. To gear NP drug-finding artificial intelligence (AI) to confront and excavate unexplored opportunities. Natural product-inspired drug discoveries based on AI to act as an innovative tool for molecular design and lead discovery. Various models of machine learning produce quickly synthesizable mimetics of the natural products templates. The invention of novel natural products mimetics by computer-assisted technology provides a feasible strategy to get the natural product with defined bio-activities. AI's hit rate makes its high importance by improving trail patterns such as dose selection, trail life span, efficacy parameters, and biomarkers. Along these lines, AI methods can be a successful tool in a targeted way to formulate advanced medicinal applications for natural products. 'Prediction of future of natural product based drug discovery is not magic, actually its artificial intelligence'Communicated by Ramaswamy H. Sarma.

Rational design and virtual screening identify mimetics of the RXR agonist valerenic acid

The ligand-sensing transcription factor retinoid X receptor (RXR) is the universal heterodimer partner of nuclear receptors and involved in multiple physiological processes. Its pharmacological modulation holds therapeutic potential in cancer and neurodegeneration but many available RXR ligands lack specificity. The sesquiterpenoid valerenic acid has been identified as RXR agonist with unprecedented subtype and homodimer preference. Here, we identified simplified mimetics of the complex natural product by rational design and virtual screening that exhibited similar activity profiles on RXR and informed about structural elements contributing to the favorable activity.

Pharmacoinformatics studies of coenzyme Q10 and potassium polyacrylate on angiotensin-converting enzyme associated with hypertension

Coenzyme Q10's (CoQ10) favorable impact on cardiovascular diseases risk factors like hypertension and atherosclerosis is linked to the antioxidant action of CoQ10 in these conditions. This study showed the possible effects of CoQ10, potassium polyacrylate (PCK), and valsartan, a reference drug, on the angiotensin-converting enzyme (ACE), a crucial component of the renin-angiotensin system. The Glide tool on Maestro 11.1 was used to calculate the respective binding affinity and binding energy of these compounds towards ACE. The Schrödinger suite was used to run molecular dynamic simulations for 100 ns. The pkCSM tool was used to forecast the pharmacokinetic characteristics and toxicological effects. The SwissADME server was used to estimate the drug-like properties of these compounds. Based on their corresponding scoring values and the negative values of the binding free energies, molecular docking analysis of CoQ10 and PCK revealed that both exhibited favorable binding affinities towards the ACE, with CoQ10 having the highest binding scores. The results showed that both CoQ10 and PCK and the reference drug, valsartan, have some amino acids in common (at the pocket site of ACE) as the key residues for binding to ACE. Both CoQ10 and PCK demonstrated drug-like qualities and were not harmful, according to the predicted pharmacokinetics and toxicology studies. The results of this study suggest that because of its inhibitory interactions with ACE, CoQ10 in particular could be useful in regulating and reducing hypertension.Communicated by Ramaswamy H. Sarma.

Discovery of novel RARα agonists using pharmacophore-based virtual screening, molecular docking, and molecular dynamics simulation studies

Nuclear retinoic acid receptors (RARs) are ligand-dependent transcription factors involved in various biological processes, such as embryogenesis, cell proliferation, differentiation, reproduction, and apoptosis. These receptors are regulated by retinoids, i.e., retinoic acid (RA) and its analogs, as receptor agonists. RAR agonists are promising therapeutic agents for the treatment of serious dermatological disorders, including some malignant conditions. By inducing apoptosis, they are able to inhibit the proliferation of diverse cancer cell lines. Also, RAR agonists have recently been identified as therapeutic options for some neurodegenerative diseases. These features make retinoids very attractive molecules for medical purposes. Synthetic selective RAR agonists have several advantages over endogenous ones, but they suffer poor pharmacokinetic properties. These compounds are normally lipophilic acids with unfavorable drug-like features such as poor oral bioavailability. Recently, highly selective, potent, and less toxic RAR agonists with proper lipophilicity, thus, good oral bioavailability have been developed for some therapeutic applications. In the present study, ligand and structure-based virtual screening technique was exploited to introduce some novel RARα agonists. Pharmacokinetic assessment was also performed in silico to suggest those compounds which have optimized drug-like features. Finally, two compounds with the best in silico pharmacological features are proposed as lead molecules for future development of RARα agonists.

Rational Design of a New RXR Agonist Scaffold Enabling Single-Subtype Preference for RXRα, RXRβ, and RXRγ

The three retinoid X receptor subtypes (RXRα, RXRβ, RXRγ) exhibit critical regulatory roles in cell proliferation and differentiation, metabolism, and inflammation. Due to their importance in nuclear receptor signaling, RXRs are widely distributed and pan-RXR agonists cause adverse effects, but the three highly conserved RXR ligand binding sites render the development of subtype-selective ligands a major challenge. We have fused elements of known RXR ligands to obtain a new RXR agonist chemotype on which minor structural modifications enabled the development of tools with single-subtype preference for RXRα, RXRβ, and RXRγ. Molecular modeling indicated different binding conformations and interaction patterns with the RXR LBDs as factors of preferential binding. In a phenotypic adipocyte differentiation experiment, only the RXRα preferential tool enhanced the adipogenic effects of pioglitazone, suggesting this subtype as particularly relevant in adipogenesis and highlighting the set of subtype-preferential RXR agonist tools as suitable for functional cellular studies.

Validation strategies for target prediction methods

Computational methods for target prediction, based on molecular similarity and network-based approaches, machine learning, docking and others, have evolved as valuable and powerful tools to aid the challenging task of mode of action identification for bioactive small molecules such as drugs and drug-like compounds. Critical to discerning the scope and limitations of a target prediction method is understanding how its performance was evaluated and reported. Ideally, large-scale prospective experiments are conducted to validate the performance of a model; however, this expensive and time-consuming endeavor is often not feasible. Therefore, to estimate the predictive power of a method, statistical validation based on retrospective knowledge is commonly used. There are multiple statistical validation techniques that vary in rigor. In this review we discuss the validation strategies employed, highlighting the usefulness and constraints of the validation schemes and metrics that are employed to measure and describe performance. We address the limitations of measuring only generalized performance, given that the underlying bioactivity and structural data are biased towards certain small-molecule scaffolds and target families, and suggest additional aspects of performance to consider in order to produce more detailed and realistic estimates of predictive power. Finally, we describe the validation strategies that were employed by some of the most thoroughly validated and accessible target prediction methods.

Targeting Nuclear Receptors in Neurodegeneration and Neuroinflammation

Nuclear receptors, also known as ligand-activated transcription factors, regulate gene expression upon ligand signals and present as attractive therapeutic targets especially in chronic diseases. Despite the therapeutic relevance of some nuclear receptors in various pathologies, their potential in neurodegeneration and neuroinflammation is insufficiently established. This perspective gathers preclinical and clinical data for a potential role of individual nuclear receptors as future targets in Alzheimer's disease, Parkinson's disease, and multiple sclerosis, and concomitantly evaluates the level of medicinal chemistry targeting these proteins. Considerable evidence suggests the high promise of ligand-activated transcription factors to counteract neurodegenerative diseases with a particularly high potential of several orphan nuclear receptors. However, potent tools are lacking for orphan receptors, and limited central nervous system exposure or insufficient selectivity also compromises the suitability of well-studied nuclear receptor ligands for functional studies. Medicinal chemistry efforts are needed to develop dedicated high-quality tool compounds for the therapeutic validation of nuclear receptors in neurodegenerative pathologies.

Molecular Scaffold Hopping via Holistic Molecular Representation

Molecular descriptors encode a variety of molecular representations for computer-assisted drug discovery. Here, we focus on the Weighted Holistic Atom Localization and Entity Shape (WHALES) descriptors, which were originally designed for scaffold hopping from natural products to synthetic molecules. WHALES descriptors capture molecular shape and partial charges simultaneously. We introduce the key aspects of the WHALES concept and provide a step-by-step guide on how to use these descriptors for virtual compound screening and scaffold hopping. The results presented can be reproduced by using the code freely available from URL: github.com/ETHmodlab/scaffold_hopping_whales .

Comprehensive Set of Tertiary Complex Structures and Palmitic Acid Binding Provide Molecular Insights into Ligand Design for RXR Isoforms

The retinoid X receptor (RXR) is a ligand-sensing transcription factor acting mainly as a universal heterodimer partner for other nuclear receptors. Despite presenting as a potential therapeutic target for cancer and neurodegeneration, adverse effects typically observed for RXR agonists, likely due to the lack of isoform selectivity, limit chemotherapeutic application of currently available RXR ligands. The three human RXR isoforms exhibit different expression patterns; however, they share high sequence similarity, presenting a major obstacle toward the development of subtype-selective ligands. Here, we report the discovery of the saturated fatty acid, palmitic acid, as an RXR ligand and disclose a uniform set of crystal structures of all three RXR isoforms in an active conformation induced by palmitic acid. A structural comparison revealed subtle differences among the RXR subtypes. We also observed an ability of palmitic acid as well as myristic acid and stearic acid to induce recruitment of steroid receptor co-activator 1 to the RXR ligand-binding domain with low micromolar potencies. With the high, millimolar endogenous concentrations of these highly abundant lipids, our results suggest their potential involvement in RXR signaling.

A Novel Biphenyl-based Chemotype of Retinoid X Receptor Ligands Enables Subtype and Heterodimer Preferences

The nuclear retinoid X receptors (RXRs) are key ligand sensing transcription factors that serve as universal nuclear receptor heterodimer partners and are thus involved in numerous physiological processes. Therapeutic targeting of RXRs holds high potential but available RXR activators suffer from limited safety. Selectivity for RXR subtypes or for certain RXR heterodimers are promising strategies for safer RXR modulation. Here, we report systematic structure-activity relationship studies on biphenyl carboxylates as new RXR ligand chemotype. We discovered specific structural modifications that enhance potency on RXRs, govern subtype preference, and vary modulation of different RXR heterodimers. Fusion of these structural motifs enabled specific tuning of subtype preferential profiles with markedly improved potency. Our results provide further evidence that RXR subtype selective ligands can be designed and present a novel chemotype of RXR modulators that can be tuned for subtype and heterodimer preferences.

Alternative retinoid X receptor (RXR) ligands

Retinoid X receptors (RXRs) control a wide variety of functions by virtue of their dimerization with other nuclear hormone receptors (NRs), contributing thereby to activities of different signaling pathways. We review known RXR ligands as transcriptional modulators of specific RXR-dimers and the associated biological processes. We also discuss the physiological relevance of such ligands, which remains frequently a matter of debate and which at present is best met by member(s) of a novel family of retinoids, postulated as Vitamin A5. Through comparison with other natural, but also with synthetic ligands, we discuss high diversity in the modes of ligand binding to RXRs resulting in agonistic or antagonistic profiles and selectivity towards specific subtypes of permissive heterodimers. Despite such diversity, direct ligand binding to the ligand binding pocket resulting in agonistic activity was preferentially preserved in the course of animal evolution pointing to its functional relevance, and potential for existence of other, species-specific endogenous RXR ligands sharing the same mode of function.

Design of Natural-Product-Inspired Multitarget Ligands by Machine Learning

A virtual screening protocol based on machine learning models was used to identify mimetics of the natural product (-)-galantamine. This fully automated approach identified eight compounds with bioactivities on at least one of the macromolecular targets of (-)-galantamine, with different polypharmacological profiles. Two of the computer-generated hits possess an expanded spectrum of bioactivity on targets relevant to the treatment of Alzheimer's disease and are suitable for hit-to-lead expansion. These results advocate multitarget drug design by advanced virtual screening protocols based on chemically informed machine learning models.

Machine Learning Consensus To Predict the Binding to the Androgen Receptor within the CoMPARA Project

The nuclear androgen receptor (AR) is one of the most relevant biological targets of Endocrine Disrupting Chemicals (EDCs), which produce adverse effects by interfering with hormonal regulation and endocrine system functioning. This paper describes novel in silico models to identify organic AR modulators in the context of the Collaborative Modeling Project of Androgen Receptor Activity (CoMPARA), coordinated by the National Center of Computational Toxicology (U.S. Environmental Protection Agency). The collaborative project involved 35 international research groups to prioritize the experimental tests of approximatively 40k compounds, based on the predictions provided by each participant. In this paper, we describe our machine learning approach to predict the binding to AR, which is based on a consensus of a multivariate Bernoulli Naive Bayes, a Random Forest, and N-Nearest Neighbor classification models. The approach was developed in compliance with the Organization of Economic Cooperation and Development (OECD) principles, trained on 1687 ToxCast molecules classified according to 11 in vitro assays, and further validated on a set of 3,882 external compounds. The models provided robust and reliable predictions and were used to gather novel data-driven insights on the structural features related to AR binding, agonism, and antagonism.

Discovery of Novel Molecular Frameworks of Farnesoid X Receptor Modulators by Ensemble Machine Learning

The bile acid activated transcription factor farnesoid X receptor (FXR) has revealed therapeutic potential as a molecular drug target for the treatment of hepatic and metabolic disorders. Despite strong efforts in FXR ligand development, the structural diversity among the known FXR modulators is limited. Only four molecular frameworks account for more than 50 % of the FXR modulators annotated in ChEMBL. Here, we leverage machine learning methods to expand the chemical space of FXR-targeting small molecules by employing an ensemble of three complementary machine learning approaches. A counter-propagation artificial neural network, a k-nearest neighbor learner, and a three-dimensional pharmacophore descriptor were combined to retrieve novel FXR ligands from a collection of more than 3 million compounds. The ensemble machine learning model identified six new FXR modulators among ten top-ranked candidates. These active hits comprise both FXR activators and antagonists with micromolar potencies. With four novel FXR ligand scaffolds, these computationally identified bioactive compounds appreciably expand the chemical space of known FXR modulators and may serve as starting points for hit-to-lead expansion.

Biaryl scaffold-focused virtual screening for anti-aggregatory and neuroprotective effects in Alzheimer's disease

Background

Alzheimer’s disease (AD) is a primary cause of dementia in ageing population affecting more than 35 million people around the globe. It is a chronic neurodegenerative disease caused by defected folding and aggregation of amyloid beta (Aβ) protein. Aβ is formed by the cleavage of membrane embedded amyloid precursor protein (APP) by using enzyme ‘transmembrane aspartyl protease, β-secretase’. Inhibition of β-secretase is a viable strategy to prevent neurotoxicity in AD. Another strategy in the treatment of AD is inhibition of acetylcholinesterase. This inhibition reduces the degradation of acetylcholine and temporarily restores the cholinergic function of neurons and improves cognitive function. Monoamine oxidase and higher glutamate levels are also found to be linked with Aβ peptide related oxidative stress. Oxidative stress leads to reduced activity of glutamate synthase resulting in significantly higher level of glutamate in brain. The aim of this study is to perform in silico screening of a virtual library of biaryl scaffold containing compounds potentially used for the treatment of AD. Screening was done against the primary targets of AD therapeutics, acetylcholinesterase, β-secretase (BACE1), Monoamine oxidases (MAO) and N-Methyl-D-aspartate (NMDA) receptor. Compounds were screened for their inhibitory potential by employing molecular docking approach using AutoDock vina. Binding energy scores were embodied in the heatmap to display varies strengths of interactions of the ligands targeting AD.

Results

Several ligands showed notable interaction with at least two targets, but the strong interaction with all the targets is shown by very few ligands. The pharmacokinetics of the interacting ligands was also predicted. The interacting ligands have good drug-likeness and brain availability essential for drugs with intracranial targets.

Conclusion

These results suggest that biaryl scaffold may be pliable to drug development for neuroprotection in AD and that the synthesis of further analogues to optimize these properties should be considered.

Electronic supplementary material

The online version of this article (10.1186/s12868-018-0472-6) contains supplementary material, which is available to authorized users.

Scaffold-Hopping from Synthetic Drugs by Holistic Molecular Representation

The discovery of novel ligand chemotypes allows to explore uncharted regions in chemical space, thereby potentially improving synthetic accessibility, potency, and the drug-likeness of molecules. Here, we demonstrate the scaffold-hopping ability of the new Weighted Holistic Atom Localization and Entity Shape (WHALES) molecular descriptors compared to seven state-of-the-art molecular representations on 30,000 compounds and 182 biological targets. In a prospective application, we apply WHALES to the discovery of novel retinoid X receptor (RXR) modulators. WHALES descriptors identified four agonists with innovative molecular scaffolds, populating uncharted regions of the chemical space. One of the agonists, possessing a rare non-acidic chemotype, revealed high selectivity on 12 nuclear receptors and comparable efficacy as bexarotene on induction of ATP-binding cassette transporter A1, angiopoietin like protein 4 and apolipoprotein E. The outcome of this research supports WHALES as an innovative tool to explore novel regions of the chemical space and to detect novel bioactive chemotypes by straightforward similarity searching.