Scaffold selection and scaffold hopping in lead generation: a medicinal chemistry perspective

AutoDesigner - Core Design, a De Novo Design Algorithm for Chemical Scaffolds: Application to the Design and Synthesis of Novel Selective Wee1 Inhibitors

The hit identification stage of a drug discovery program generally involves the design of novel chemical scaffolds with desired biological activity against the target(s) of interest. One common approach is scaffold hopping, which is the manual design of novel scaffolds based on known chemical matter. One major limitation of this approach is narrow chemical space exploration, which can lead to difficulties in maintaining or improving biological activity, selectivity, and favorable property space. Another limitation is the lack of preliminary structure-activity relationship (SAR) data around these designs, which could lead to selecting suboptimal scaffolds to advance lead optimization. To address these limitations, we propose AutoDesigner - Core Design (CoreDesign), a de novo scaffold design algorithm. Our approach is a cloud-integrated, de novo design algorithm for systematically exploring and refining chemical scaffolds against biological targets of interest. The algorithm designs, evaluates, and optimizes a vast range, from millions to billions, of molecules in silico, following defined project parameters encompassing structural novelty, physicochemical attributes, potency, and selectivity using active-learning FEP. To validate CoreDesign in a real-world drug discovery setting, we applied it to the design of novel, potent Wee1 inhibitors with improved selectivity over PLK1. Starting from a single known ligand and receptor structure, CoreDesign rapidly explored over 23 billion molecules to identify 1,342 novel chemical series with a mean of 4 compounds per scaffold. To rapidly analyze this large amount of data and prioritize chemical scaffolds for synthesis, we utilize t-Distributed Stochastic Neighbor Embedding (t-SNE) plots of in silico properties. The chemical space projections allowed us to rapidly identify a structurally novel 5-5 fused core meeting all the hit-identification requirements. Several compounds were synthesized and assayed from the scaffold, displaying good potency against Wee1 and excellent PLK1 selectivity. Our results suggest that CoreDesign can significantly speed up the hit-identification process and increase the probability of success of drug discovery campaigns by allowing teams to bring forward high-quality chemical scaffolds derisked by the availability of preliminary SAR.

Scaffold-Hopping Strategies in Aurone Optimization: A Comprehensive Review of Synthetic Procedures and Biological Activities of Nitrogen and Sulfur Analogues

Aurones, particular polyphenolic compounds belonging to the class of minor flavonoids and overlooked for a long time, have gained significative attention in medicinal chemistry in recent years. Indeed, considering their unique and outstanding biological properties, they stand out as an intriguing reservoir of new potential lead compounds in the drug discovery context. Nevertheless, several physicochemical, pharmacokinetic, and pharmacodynamic (P3) issues hinder their progression in more advanced phases of the drug discovery pipeline, making lead optimization campaigns necessary. In this context, scaffold hopping has proven to be a valuable approach in the optimization of natural products. This review provides a comprehensive and updated picture of the scaffold-hopping approaches directed at the optimization of natural and synthetic aurones. In the literature analysis, a particular focus is given to nitrogen and sulfur analogues. For each class presented, general synthetic procedures are summarized, highlighting the key advantages and potential issues. Furthermore, the biological activities of the most representative scaffold-hopped compounds are presented, emphasizing the improvements achieved and the potential for further optimization compared to the aurone class.

Synthesis and herbicidal activities of 2-((3-pyridinyl-2-oxo-2,3-dihydrobenzo[d]thiazol-6-yl)oxy)propionates

BACKGROUND

The discovery of lead compounds is fundamental to herbicide innovation, yet the limited availability of valuable lead compounds has impeded their progress in recent years. The study presents a novel molecular scaffold that exhibits remarkably potent herbicidal activity.

RESULTS

Through a scaffold-hopping strategy, a highly potent lead compound for herbicides, namely 3-(2-pyridinyl)-benzothiazol-2-one, was unexpectedly discovered during attempts to structurally modify haloxyfop, a commercial aryl-oxy-phenoxy-propionate herbicide. To investigate the structure-activity relationship (SAR) of the newly discovered herbicidal chemicals, a series of 2-(2-oxo-3-(pyridin-2-yl)-2,3-dihydrobenzo[d]thiazol-6-yloxy)propanoic acid derivatives, I-01 ~ I-27, were designed and synthesized. SAR analysis revealed that trifluoromethyl at the 5-position of pyridine is crucial for herbicidal activity, whereas additional fluorine or Cl atom at the 3-position of pyridine significantly enhances activity. Carboxylic ester derivatives exhibit superior herbicidal activity compared with amide derivatives. Moreover, the activity of carboxylic ester derivatives decreases with C chain extension, but the introduction of O atoms in the side chain benefits activity enhancement. Pot experiments conducted in a glasshouse demonstrated that I-01 and I-09 exhibited potent postemergence herbicidal activity against broadleaf weeds, and completely inhibited growth of Amaranthus retroflex, Abutilon theophrasti and Portulaca oleracea at a dosage of 75 g ha-1.

CONCLUSION

Despite the initial goal of scaffold-hopping not being achieved, we have successfully identified a novel molecular scaffold exhibiting exceptional herbicidal activity, thereby presenting innovative prospects for herbicide development. © 2024 Society of Chemical Industry.

A dual diffusion model enables 3D molecule generation and lead optimization based on target pockets

Structure-based generative chemistry is essential in computer-aided drug discovery by exploring a vast chemical space to design ligands with high binding affinity for targets. However, traditional in silico methods are limited by computational inefficiency, while machine learning approaches face bottlenecks due to auto-regressive sampling. To address these concerns, we have developed a conditional deep generative model, PMDM, for 3D molecule generation fitting specified targets. PMDM consists of a conditional equivariant diffusion model with both local and global molecular dynamics, enabling PMDM to consider the conditioned protein information to generate molecules efficiently. The comprehensive experiments indicate that PMDM outperforms baseline models across multiple evaluation metrics. To evaluate the applications of PMDM under real drug design scenarios, we conduct lead compound optimization for SARS-CoV-2 main protease (Mpro) and Cyclin-dependent Kinase 2 (CDK2), respectively. The selected lead optimization molecules are synthesized and evaluated for their in-vitro activities against CDK2, displaying improved CDK2 activity.

Optimization of binding affinities in chemical space with generative pre-trained transformer and deep reinforcement learning

Background

The key challenge in drug discovery is to discover novel compounds with desirable properties. Among the properties, binding affinity to a target is one of the prerequisites and usually evaluated by molecular docking or quantitative structure activity relationship (QSAR) models.

Methods

In this study, we developed SGPT-RL, which uses a generative pre-trained transformer (GPT) as the policy network of the reinforcement learning (RL) agent to optimize the binding affinity to a target. SGPT-RL was evaluated on the Moses distribution learning benchmark and two goal-directed generation tasks, with Dopamine Receptor D2 (DRD2) and Angiotensin-Converting Enzyme 2 (ACE2) as the targets. Both QSAR model and molecular docking were implemented as the optimization goals in the tasks. The popular Reinvent method was used as the baseline for comparison.

Results

The results on the Moses benchmark showed that SGPT-RL learned good property distributions and generated molecules with high validity and novelty. On the two goal-directed generation tasks, both SGPT-RL and Reinvent were able to generate valid molecules with improved target scores. The SGPT-RL method achieved better results than Reinvent on the ACE2 task, where molecular docking was used as the optimization goal. Further analysis shows that SGPT-RL learned conserved scaffold patterns during exploration.

Conclusions

The superior performance of SGPT-RL in the ACE2 task indicates that it can be applied to the virtual screening process where molecular docking is widely used as the criteria. Besides, the scaffold patterns learned by SGPT-RL during the exploration process can assist chemists to better design and discover novel lead candidates.

Curation and cheminformatics analysis of a Ugi-reaction derived library (URDL) of synthetically tractable small molecules for virtual screening application

Virtual screening (VS) is an important approach in drug discovery and relies on the availability of a virtual library of synthetically tractable molecules. Ugi reaction (UR) represents an important multi-component reaction (MCR) that reliably produces a peptidomimetic scaffold. Recent literature shows that a tactically assembled Ugi adduct can be subjected to further chemical modifications to yield a variety of rings and scaffolds, thus, renewing the interest in this old reaction. Given the reliability and efficiency of UR, we collated an UR derived library (URDL) of small molecules (total = 5773) for VS. The synthesis of the majority of URDL molecules may be carried out in 1-2 pots in a time and cost-effective manner. The detailed analysis of the average property and chemical space of URDL was also carried out using the open-source Datawarrior program. The comparison with FDA-approved oral drugs and inhibitors of protein-protein interactions (iPPIs) suggests URDL molecules are 'clean', drug-like, and conform to a structurally distinct space from the other two categories. The average physicochemical properties of compounds in the URDL library lie closer to iPPI molecules than oral drugs thus suggesting that the URDL resource can be applied to discover novel iPPI molecules. The URDL molecules consist of diverse ring systems, many of which have not been exploited yet for drug design. Thus, URDL represents a small virtual library of drug-like molecules with unexplored chemical space designed for VS. The structures of all molecules of URDL, oral drugs, and iPPI compounds are being made freely accessible as supplementary information for broader application.

Systematical Evaluation of the Structure-Cardiotoxicity Relationship of 7-Azaindazole-based PI3K Inhibitors Designed by Bioisosteric Approach

A growing concern of cardiotoxicity induced by PI3K inhibitors has raised the requirements to evaluate the structure-cardiotoxicity relationship (SCR) in the development process of novel inhibitors. Based on three bioisosteric 7-azaindazole-based candidate inhibitors namely FD269, FD268 and FD274 that give same order of inhibitory concentration 50% (IC50) magnitude against PI3Ks, in this work, we proposed to systematically evaluate the SCR of 7-azaindazole-based PI3K inhibitors designed by bioisosteric approach. The 24-h lethal concentrations 50% (LC50) of FD269, FD268 and FD274 against zebrafish embryos were 0.35, 4.82 and above 50 μM (not detected), respectively. Determination of the heart rate, pericardial and yolk-sac areas and vascular malformation confirmed the remarkable reduction in the cardiotoxicity of from FD269 to FD268 and to FD274. The IC50s of all three compounds against the hERG channel were tested on the CHO cell line that constitutively expressing hERG channel, which were all higher than 20 μM. The transcriptomic analysis revealed that FD269 and FD268 induced the up-regulation of noxo1b, which encodes a subunit of an NADPH oxidase evoking the oxidative stress. Furthermore, immunohistochemistry tests confirmed the structure-dependent attenuation of the overproduction of ROS and cardiac apoptosis. Our results verified the feasibility of bioisosteric replacement to attenuate the cardiotoxicity of 7-azaindazole-based PI3K inhibitors, suggesting that the screening for PI3K inhibitors with both high potency and low cardiotoxicity from bioisosteres would be a beneficial trial.

Synthesis and Biological Activity of Silicon-Containing Ethylsulfonylpyridine Insecticides

Ethylsulfonylpyridines are a novel chemical class of insecticides with excellent broad-spectrum activity and an unprecedented mode of action. With the objective of discovering novel ethylsulfonylpyridines with a broader spectrum, wider range of uses, and/or improved properties, we have started a research program aimed at introducing silicon motifs and studying their biological effects. We designed a series of Oxazosulfyl analogues where the hydrogen atom at the 5-position of the pyridyl moiety is replaced by a trialkylsilyl group and prepared these compounds applying denovo synthetic methodology. Our novel ethylsulfonylpyridines exhibit excellent insecticidal activities. The best compound, A18, resulting from our research exhibited an LC50 value of 0.30 mg/L against Plutella xylostella and reached the activity level of the commercial standard Oxazosulfyl. Our findings confirmed our working hypothesis that at the 5-position of the pyridyl moiety larger groups with different hydrophobic, electronic, and steric properties are tolerated.

Strategies in the Design and Development of Non-Nucleoside Reverse Transcriptase Inhibitors (NNRTIs)

AIDS (acquired immunodeficiency syndrome) is a potentially life-threatening infectious disease caused by human immunodeficiency virus (HIV). To date, thousands of people have lost their lives annually due to HIV infection, and it continues to be a big public health issue globally. Since the discovery of the first drug, Zidovudine (AZT), a nucleoside reverse transcriptase inhibitor (NRTI), to date, 30 drugs have been approved by the FDA, primarily targeting reverse transcriptase, integrase, and/or protease enzymes. The majority of these drugs target the catalytic and allosteric sites of the HIV enzyme reverse transcriptase. Compared to the NRTI family of drugs, the diverse chemical class of non-nucleoside reverse transcriptase inhibitors (NNRTIs) has special anti-HIV activity with high specificity and low toxicity. However, current clinical usage of NRTI and NNRTI drugs has limited therapeutic value due to their adverse drug reactions and the emergence of multidrug-resistant (MDR) strains. To overcome drug resistance and efficacy issues, combination therapy is widely prescribed for HIV patients. Combination antiretroviral therapy (cART) includes more than one antiretroviral agent targeting two or more enzymes in the life cycle of the virus. Medicinal chemistry researchers apply different optimization strategies including structure- and fragment-based drug design, prodrug approach, scaffold hopping, molecular/fragment hybridization, bioisosterism, high-throughput screening, covalent-binding, targeting highly hydrophobic channel, targeting dual site, and multi-target-directed ligand to identify and develop novel NNRTIs with high antiviral activity against wild-type (WT) and mutant strains. The formulation experts design various delivery systems with single or combination therapies and long-acting regimens of NNRTIs to improve pharmacokinetic profiles and provide sustained therapeutic effects.

Oxetane-, Azetidine-, and Bicyclopentane-Bearing N-Heterocycles from Ynones: Scaffold Diversification via Ruthenium-Catalyzed Oxidative Alkynylation

A process for 3-fold scaffold diversification is achieved via ruthenium-catalyzed oxidative alkynylation of commercially available oxetanols, azetidinols and bicyclopentanols to form α,β-acetylenic ketones (ynones), which are subsequently converted to oxetane-, azetidine- and bicyclopentane-bearing pyrazoles, isoxazoles and pyrimidines. A one-pot oxidative alkynylation-condensation protocol that directly converts azetidinols to azetidine-substituted pyrazoles or pyrimidines is demonstrated.

Data-Driven Multi-Objective Optimization Tactics for Catalytic Asymmetric Reactions Using Bisphosphine Ligands

Optimization of the catalyst structure to simultaneously improve multiple reaction objectives (e.g., yield, enantioselectivity, and regioselectivity) remains a formidable challenge. Herein, we describe a machine learning workflow for the multi-objective optimization of catalytic reactions that employ chiral bisphosphine ligands. This was demonstrated through the optimization of two sequential reactions required in the asymmetric synthesis of an active pharmaceutical ingredient. To accomplish this, a density functional theory-derived database of >550 bisphosphine ligands was constructed, and a designer chemical space mapping technique was established. The protocol used classification methods to identify active catalysts, followed by linear regression to model reaction selectivity. This led to the prediction and validation of significantly improved ligands for all reaction outputs, suggesting a general strategy that can be readily implemented for reaction optimizations where performance is controlled by bisphosphine ligands.

Contemporary Medicinal Chemistry Strategies for the Discovery and Development of Novel HIV-1 Non-nucleoside Reverse Transcriptase Inhibitors

Currently, HIV-1 non-nucleoside reverse transcriptase inhibitors (NNRTIs) are a major component of the highly active anti-retroviral therapy (HAART) regimen. However, the occurrence of drug-resistant strains and adverse reactions after long-term usage have inevitably compromised the clinical application of NNRTIs. Therefore, the development of novel inhibitors with distinct anti-resistance profiles and better pharmacological properties is still an enormous challenge. Herein, we summarize state-of-the-art medicinal chemistry strategies for the discovery of potent NNRTIs, such as structure-based design strategies, contemporary computer-aided drug design, covalent-binding strategies, and the application of multi-target-directed ligands. The strategies described here will facilitate the identification of promising HIV-1 NNRTIs.

Synthesis and In Silico Docking Studies of Ethyl 2-(2-Arylidene-1-Alkylhydrazinyl)Thiazole-4-Carboxylates as Antiglycating agents

Ethyl 2-(2-arylidene-1-alkylhydrazinyl)thiazole-4-carboxylates ( 1a-k ) were synthesized by alkylation on HN- of ethyl 2-(2-arylidenehydrazinyl)thiazole-4-carboxylates. The proposed structures ( 1a-k ) are corroborated by spectro-analytical techniques like UV, FT-IR, 1 H-, 13 C-NMR and HRMS. The compounds ( 1a-k ) were screened for their antiglycation and antioxidant assays. The in vitro antiglycation results revealed promising activity of compounds 1a , 1b , 1d , 1e , 1f , 1g , 1j and 1k with IC 50 values 0.0004 ± 1.097-17.22 ± 0.538 µM/mL when compared to standard, aminoguanidine (IC 50 = 25.50 ± 0.337 µM/mL). Among all tested compounds 1j and 1k are the best antiglycating agents with IC 50 values 1.848 ± 0.646 and 0.0004 ± 1.097 µM/mL, respectively. The in-silico studies also agree with these results where binding energy for 1j and 1k was found to be -9.25 and -8.42 kcal/mol with calculated dissociation constants of 0.16 and 0.67 µM. The antiglycation results demonstrate the application of these compounds in reducing diabetic complications.

Single-Pot Preparation of 4-Amino-2-(het)aryl-5-Substituted Thiazoles Employing Functionalized Dithioesters as Thiocarbonyl Precursors

An effective, diversity oriented, one-pot reaction of 4-amino-2-(het)aryl/alkyl-5-functionalized thiazoles has been disclosed, utilizing aryl/heteroaryl/alkyl dithioesters as thiocarbonyl coupling partners in a modified Thorpe-Ziegler type cyclization. The reaction proceeds at room temperature, under mild conditions, in excellent yields, displaying broad functional group compatibility at 2 and 5 positions of thiazoles. This synthetic strategy has been further expanded for the one-pot construction of two highly potent tubulin polymerization inhibitors, i.e., 2-(het)aryl-4-amino-5-(3,4,5-trimethoxyaroyl) thiazoles, in high yields.

Heterocyclic Compounds: Pharmacology of Pyrazole Analogs From Rational Structural Considerations

Low quality of life and life-threatening conditions often demand pharmacological screening of lead compounds. A spectrum of pharmacological activities has been attributed to pyrazole analogs. The substitution, replacement, or removal of functional groups on a pyrazole ring appears consistent with diverse molecular interactions, efficacy, and potency of these analogs. This mini-review explores cytotoxic, cytoprotective, antinociceptive, anti-inflammatory, and antidepressant activities of some pyrazole analogs to advance structure-related pharmacological profiles and rational design of new analogs. Numerous interactions of these derivatives at their targets could impact future research considerations and prospects while offering opportunities for optimizing therapeutic activity with fewer adverse effects.

Global Assessment of Substituents on the Basis of Analogue Series

While bioisosteric replacements have been extensively investigated, comprehensive analyses of R-/functional groups have thus far been rare in medicinal chemistry. We introduce a new analysis concept for the exploration of chemical substituent space that is based upon bioactive analogue series as a source. From ∼24,000 analogue series, more than 19,000 substituents were isolated that were differently distributed. A subset of ∼400 substituent fragments occurred most frequently in different structural contexts. These substituents contained well-known R-groups as well as novel structures. Substitution site-specific replacement and network analysis revealed that chemically similar substituents preferentially occurred at given sites and identified intuitive substitution pathways that can be explored for compound design. Taken together, the results of our analysis provide new insights into substituent space and identify preferred substituents on the basis of analogue series. As a part of our study, all the data reported are made freely available.

Structural Similarities between Some Common Fluorophores Used in Biology, Marketed Drugs, Endogenous Metabolites, and Natural Products

It is known that at least some fluorophores can act as 'surrogate' substrates for solute carriers (SLCs) involved in pharmaceutical drug uptake, and this promiscuity is taken to reflect at least a certain structural similarity. As part of a comprehensive study seeking the 'natural' substrates of 'orphan' transporters that also serve to take up pharmaceutical drugs into cells, we have noted that many drugs bear structural similarities to natural products. A cursory inspection of common fluorophores indicates that they too are surprisingly 'drug-like', and they also enter at least some cells. Some are also known to be substrates of efflux transporters. Consequently, we sought to assess the structural similarity of common fluorophores to marketed drugs, endogenous mammalian metabolites, and natural products. We used a set of some 150 fluorophores along with standard fingerprinting methods and the Tanimoto similarity metric. Results: The great majority of fluorophores tested exhibited significant similarity (Tanimoto similarity > 0.75) to at least one drug, as judged via descriptor properties (especially their aromaticity, for identifiable reasons that we explain), by molecular fingerprints, by visual inspection, and via the "quantitative estimate of drug likeness" technique. It is concluded that this set of fluorophores does overlap with a significant part of both the drug space and natural products space. Consequently, fluorophores do indeed offer a much wider opportunity than had possibly been realised to be used as surrogate uptake molecules in the competitive or trans-stimulation assay of membrane transporter activities.

Streptomyces-Derived Metabolites with Potential Photoprotective Properties-A Systematic Literature Review and Meta-Analysis on the Reported Chemodiversity

Sun overexposure is associated with the development of diseases that primarily affect the skin, which can lead to skin cancer. Among the main measures of photoprotection is the use of sunscreens. However, there is currently concern about the reported harmful effects to both humans and the environment due to several of the sunscreen ingredients available on the market. For this reason, the search for and development of new agents with photoprotective properties is required. In searching for these metabolites, researchers have turned their attention to microbial sources, especially the microbiota in unusual hostile environments. Among the diverse microorganisms available in nature, Actinobacteria and specifically Streptomyces, have been shown to be a source of metabolites with various biological activities of interest, such as antimicrobial, antitumor and immunomodulator activities. Herein, we present the results of a systematic review of the literature in which Streptomyces isolates were studied as a source of compounds with photoprotective properties. A meta-analysis of the structure-property and structure-activity relationships of those metabolites identified in the qualitative analysis phase was also carried out. These findings indicate that Streptomyces are a source of metabolites with potential applications in the development of new, safe and more eco-friendly sunscreens.

Auto In Silico Ligand Directing Evolution to Facilitate the Rapid and Efficient Discovery of Drug Lead

Motivated by the growing demand for reducing the chemical optimization burden of H2L, we developed auto in silico ligand directing evolution (AILDE, http://chemyang.ccnu.edu.cn/ccb/server/AILDE), an efficient and general approach for the rapid identification of drug leads in accessible chemical space. This computational strategy relies on minor chemical modifications on the scaffold of a hit compound, and it is primarily intended for identifying new lead compounds with minimal losses or, in some cases, even increases in ligand efficiency. We also described how AILDE greatly reduces the chemical optimization burden in the design of mesenchymal-epithelial transition factor (c-Met) kinase inhibitors. We only synthesized eight compounds and found highly efficient compound 5g, which showed an ∼1,000-fold improvement in in vitro activity compared with the hit compound. 5g also displayed excellent in vivo antitumor efficacy as a drug lead. We believe that AILDE may be applied to a large number of studies for rapid design and identification of drug leads.

•

AILDE was developed for the rapid identification of drug leads

•

A potent drug lead targeted to c-Met was found by synthesizing only eight compounds

Development of the triazole-fused pyrimidine derivatives as highly potent and reversible inhibitors of histone lysine specific demethylase 1 (LSD1/KDM1A)

Histone lysine specific demethylase 1 (LSD1) has been recognized as an important modulator in post-translational process in epigenetics. Dysregulation of LSD1 has been implicated in the development of various cancers. Herein, we report the discovery of the hit compound 8a (IC₅₀ = 3.93 μmol/L) and further medicinal chemistry efforts, leading to the generation of compound 15u (IC₅₀ = 49 nmol/L, and K_i = 16 nmol/L), which inhibited LSD1 reversibly and competitively with H3K4me2, and was selective to LSD1 over MAO-A/B. Docking studies were performed to rationalize the potency of compound 15u. Compound 15u also showed strong antiproliferative activity against four leukemia cell lines (OCL-AML3, K562, THP-1 and U937) as well as the lymphoma cell line Raji with the IC₅₀ values of 1.79, 1.30, 0.45, 1.22 and 1.40 μmol/L, respectively. In THP-1 cell line, 15u significantly inhibited colony formation and caused remarkable morphological changes. Compound 15u induced expression of CD86 and CD11b in THP-1 cells, confirming its cellular activity and ability of inducing differentiation. The findings further indicate that targeting LSD1 is a promising strategy for AML treatment, the triazole-fused pyrimidine derivatives are new scaffolds for the development of LSD1/KDM1A inhibitors.

QligFEP: an automated workflow for small molecule free energy calculations in Q

The process of ligand binding to a biological target can be represented as the equilibrium between the relevant solvated and bound states of the ligand. This which is the basis of structure-based, rigorous methods such as the estimation of relative binding affinities by free energy perturbation (FEP). Despite the growing capacity of computing power and the development of more accurate force fields, a high throughput application of FEP is currently hampered due to the need, in the current schemes, of an expert user definition of the "alchemical" transformations between molecules in the series explored. Here, we present QligFEP, a solution to this problem using an automated workflow for FEP calculations based on a dual topology approach. In this scheme, the starting poses of each of the two ligands, for which the relative affinity is to be calculated, are explicitly present in the MD simulations associated with the (dual topology) FEP transformation, making the perturbation pathway between the two ligands univocal. We show that this generalized method can be applied to accurately estimate solvation free energies for amino acid sidechain mimics, as well as the binding affinity shifts due to the chemical changes typical of lead optimization processes. This is illustrated in a number of protein systems extracted from other FEP studies in the literature: inhibitors of CDK2 kinase and a series of A2A adenosine G protein-coupled receptor antagonists, where the results obtained with QligFEP are in excellent agreement with experimental data. In addition, our protocol allows for scaffold hopping perturbations to identify the binding affinities between different core scaffolds, which we illustrate with a series of Chk1 kinase inhibitors. QligFEP is implemented in the open-source MD package Q, and works with the most common family of force fields: OPLS, CHARMM and AMBER.

Role of Hydrophobicity and Charge of Amyloid-Beta Oligomer Eliminating d-Peptides in the Interaction with Amyloid-Beta Monomers

Inhibition of the self-assembly process of amyloid-beta and even more the removal of already existing toxic amyloid-beta assemblies represent promising therapeutic strategies against Alzheimer's disease. To approach this aim, we selected a d-enantiomeric peptide by phage-display based on the interaction with amyloid-beta monomers. This lead compound was successfully optimized by peptide microarrays with respect to its affinity and specificity to the target resulting in d-peptides with both increased hydrophobicity and net charge. Here, we present a detailed biophysical characterization of the interactions between these optimized d-peptides and amyloid-beta monomers in comparison to the original lead compound in order to obtain a more thorough understanding of the physicochemical determinants of the interactions. Kinetics and apparent stoichiometry of complex formation were studied using surface plasmon resonance. Potential modes of binding to amyloid-beta were identified, and the influences of ionic strength on complex stability, as well as on the inhibitory effect on amyloid-beta aggregation were investigated. The results reveal a very different mode of interaction of the optimized d-peptides based on a combination of electrostatic and hydrophobic interactions as compared to the mostly electrostatically driven interaction of the lead compound. These conclusions were supported by the thermodynamic profiles of the interaction between optimized d-peptides and Aβ monomers, which indicate an increase in binding entropy with respect to the lead compound.

Synthesis of 1,2,4-Trisubstituted-(1H)-imidazoles through Cu(OTf)2-/I2-Catalyzed C-C Bond Cleavage of Chalcones and Benzylamines

1,2,4-Trisubstituted-(1H)-imidazoles have been synthesized by the Cu(OTf)₂- and I₂-catalyzed unusual C-C bond cleavage of chalcones and benzylamines. After the α,β-unsaturated C-C bond cleavage, the β-portion is eliminated from the reaction. Various aryl- and heteroaryl-substituted chalcones and benzylamines were well tolerated in this unusual transformation to yield the trisubstituted-(1H)-imidazoles.

Current insights into anti-HIV drug discovery and development: a review of recent patent literature (2014-2017)

INTRODUCTION

To deal with the rapid emergence of drug resistance challenges, together with the difficulty to eradicate the virus, off-target effects and significant cumulative drug toxicities, it is still imperative to develop next-generation anti-HIV agents with novel chemical classes or new mechanisms of action.

AREAS COVERED

We primarily focused on current strategies to discover novel anti-HIV agents. Moreover, examples of anti-HIV lead compounds were mainly selected from recently patented publications (reported between 2014 and 2017). In particular, 'privileged structure'-focused substituents decorating approach, scaffold hopping, natural-product diversification and prodrug are focused on. Furthermore, exploitation of new compounds with unexplored mechanisms of action and medicinal chemistry strategies to deplete the HIV reservoir were also described. Perspectives that could inspire future anti-HIV drug discovery are delineated.

EXPERT OPINION

Even if a large number of patents have been disclosed recently, additional HIV inhibitors are still required, especially novel chemical skeletons displaying a unexploited mechanism of action. Current medicinal chemistry strategies are inadequate, and appropriate and new methodologies and technologies should be exploited to identify novel anti-HIV drug candidates in a time- and cost- effective manner.

Agrochemical lead optimization by scaffold hopping

Scaffold hopping, the exchange of a specific portion of a potential active ingredient with another substructure with the aim of finding isofunctional molecular structures with significantly different molecular backbones, often offers the chance in lead discovery or optimization to mitigate problems related to toxicity, intellectual property, and insufficient potency or stability. Scaffold hopping tools such as isosteric ring replacement including 1,3 nitrogen shift and cyclic imine-amide isosterism, but also ring opening and ring closure approaches, functional group isosterism, reversion of functional groups, chain shortening, chain lengthening, and scaffolds delivered by natural products, have become a permanent fixture of the innovation and optimization process in crop protection research. Their appropriate use will be explained through examples of success stories in the field of agrochemistry. Analogies to, but also differences from, the main categories of scaffold hopping in medicinal drug discovery are discussed. © 2017 Society of Chemical Industry.

Pharmacological relationships and ligand discovery of G protein-coupled receptors revealed by simultaneous ligand and receptor clustering

Conventional ligand and receptor similarity methods have been extensively used for exposing pharmacological relationships and drug lead discovery. They may in some cases neglect minor relationships useful for target hopping particularly against the remote family members. To complement the conventional methods for capturing these minor relationships, we developed a new method that uses a SLARC (Simultaneous Ligand And Receptor Clustering) 2D map to simultaneously characterize the ligand structural and receptor binding-site sequence relationships of a receptor family. The SLARC maps of the rhodopsin-like GPCR family comprehensively revealed scaffold hopping, target hopping, and multi-target relationships for the ligands of both homologous and remote family members. Their usefulness in new ligand discovery was validated by guiding the prospective discovery of novel indole piperazinylpyrimidine dual-targeting adenosine A2A receptor antagonist and dopamine D2 agonist compounds. The SLARC approach is useful for revealing pharmacological relationships and discovering new ligands at target family levels.

Indole-based allosteric inhibitors of HIV-1 integrase

Employing a scaffold hopping approach, a series of allosteric HIV-1 integrase (IN) inhibitors (ALLINIs) have been synthesized based on an indole scaffold. These compounds incorporate the key elements utilized in quinoline-based ALLINIs for binding to the IN dimer interface at the principal LEDGF/p75 binding pocket. The most potent of these compounds displayed good activity in the LEDGF/p75 dependent integration assay (IC50=4.5μM) and, as predicted based on the geometry of the five- versus six-membered ring, retained activity against the A128T IN mutant that confers resistance to many quinoline-based ALLINIs.

High-Throughput Synthesis of Diverse Compound Collections for Lead Discovery and Optimization

Small-molecule intervention of protein function is one central dogma of drug discovery. The generation of small-molecule libraries fuels the discovery pipeline at many stages and thereby resembles a key aspect of this endeavor. High-throughput synthesis is a major source for compound libraries utilized in academia and industry, seeking new chemical modulators of pharmacological targets. Here, we discuss the crucial factors of library design strategies from the perspective of synthetic chemistry, giving a brief historic background and a summary of current approaches. Simple measures of success of a high-throughput synthesis such as quantity or diversity have long been discarded and replaced by more integrated measures. Case studies are presented and put into context to highlight the cross-connectivity of the various stages of the drug discovery process.

Creation of Novel Cores for β-Secretase (BACE-1) Inhibitors: A Multiparameter Lead Generation Strategy

In order to find optimal core structures as starting points for lead optimization, a multiparameter lead generation workflow was designed with the goal of finding BACE-1 inhibitors as a treatment for Alzheimer's disease. De novo design of core fragments was connected with three predictive in silico models addressing target affinity, permeability, and hERG activity, in order to guide synthesis. Taking advantage of an additive SAR, the prioritized cores were decorated with a few, well-characterized substituents from known BACE-1 inhibitors in order to allow for core-to-core comparisons. Prediction methods and analyses of how physicochemical properties of the core structures correlate to in vitro data are described. The syntheses and in vitro data of the test compounds are reported in a separate paper by Ginman et al. [J. Med. Chem. 2013, 56, 4181-4205]. The affinity predictions are described in detail by Roos et al. [J. Chem. Inf. 2014, DOI: 10.1021/ci400374z].

Potency prediction of β-secretase (BACE-1) inhibitors using density functional methods

Scoring potency is a main challenge for structure based drug design. Inductive effects of subtle variations in the ligand are not possible to accurately predict by classical computational chemistry methods. In this study, the problem of predicting potency of ligands with electronic variations participating in key interactions with the protein was addressed. The potency was predicted for a large set of cyclic amidine and guanidine cores extracted from β-secretase (BACE-1) inhibitors. All cores were of similar size and had equal interaction motifs but were diverse with respect to electronic substitutions. A density functional theory approach, in combination with a representation of the active site of a protein using only key residues, was shown to be predictive. This computational approach was used to guide and support drug design, within the time frame of a normal drug discovery design cycle.

Cofactor-independent phosphoglycerate mutase from nematodes has limited druggability, as revealed by two high-throughput screens

Cofactor-independent phosphoglycerate mutase (iPGAM) is essential for the growth of C. elegans but is absent from humans, suggesting its potential as a drug target in parasitic nematodes such as Brugia malayi, a cause of lymphatic filariasis (LF). iPGAM's active site is small and hydrophilic, implying that it may not be druggable, but another binding site might permit allosteric inhibition. As a comprehensive assessment of iPGAM's druggability, high-throughput screening (HTS) was conducted at two different locations: ∼220,000 compounds were tested against the C. elegans iPGAM by Genzyme Corporation, and ∼160,000 compounds were screened against the B. malayi iPGAM at the National Center for Drug Screening in Shanghai. iPGAM's catalytic activity was coupled to downstream glycolytic enzymes, resulting in NADH consumption, as monitored by a decline in visible-light absorbance at 340 nm. This assay performed well in both screens (Z′-factor >0.50) and identified two novel inhibitors that may be useful as chemical probes. However, these compounds have very modest potency against the B. malayi iPGAM (IC₅₀ >10 µM) and represent isolated singleton hits rather than members of a common scaffold. Thus, despite the other appealing properties of the nematode iPGAMs, their low druggability makes them challenging to pursue as drug targets. This study illustrates a “druggability paradox” of target-based drug discovery: proteins are generally unsuitable for resource-intensive HTS unless they are considered druggable, yet druggability is often difficult to predict in the absence of HTS data.

Parasitic worms like Brugia malayi cause widespread lymphatic filariasis (LF) in southeast Asia and sub-Saharan Africa. The adult worms causing most of the symptoms of LF are difficult to treat with existing drugs. As a possible step toward new LF drugs, we searched for inhibitors of the B. malayi cofactor-independent phosphoglycerate mutase (iPGAM), an enzyme thought to be critical to survival and development of this parasite. Despite testing over 100,000 compounds at each of two screening centers, we found only two compounds that consistently inhibited the B. malayi enzyme more strongly than the cofactor-dependent enzyme found in humans. These compounds have limited potency and are not especially great starting points for drug development. The 3-dimensional structure of iPGAM suggests that the active site is difficult to access from the surrounding solvent, which may partly explain our very low yield of inhibitors. We conclude that iPGAM may not be an ideal drug target in B. malayi or related organisms because it is difficult to inhibit with druglike compounds.

Rational hopping of a peptidic scaffold into non-peptidic scaffolds: structurally novel potent proteasome inhibitors derived from a natural product, belactosin A

Rational scaffold hopping of a natural product belactosin A derivative based on the pharmacophore model constructed resulted in the identification of the significantly simplified highly potent non-peptide derivatives.

Core refinement toward permeable β-secretase (BACE-1) inhibitors with low hERG activity

By use of iterative design aided by predictive models for target affinity, brain permeability, and hERG activity, novel and diverse compounds based on cyclic amidine and guanidine cores were synthesized with the goal of finding BACE-1 inhibitors as a treatment for Alzheimer's disease. Since synthesis feasibility had low priority in the design of the cores, an extensive synthesis effort was needed to make the relevant compounds. Syntheses of these compounds are reported, together with physicochemical properties and structure-activity relationships based on in vitro data. Four crystal structures of diverse amidines binding in the active site are deposited and discussed. Inhibitors of BACE-1 with 3 μM to 32 nM potencies in cells are shown, together with data on in vivo brain exposure levels for four compounds. The results presented show the importance of the core structure for the profile of the final compounds.

Exploring inhibition of Pdx1, a component of the PLP synthase complex of the human malaria parasite Plasmodium falciparum

Malaria tropica is a devastating infectious disease caused by Plasmodium falciparum. This parasite synthesizes vitamin B6 de novo via the PLP (pyridoxal 5'-phosphate) synthase enzymatic complex consisting of PfPdx1 and PfPdx2 proteins. Biosynthesis of PLP is largely performed by PfPdx1, ammonia provided by PfPdx2 subunits is condensed together with R5P (D-ribose 5-phosphate) and G3P (DL-glyceraldehyde 3-phosphate). PfPdx1 accommodates both the R5P and G3P substrates and intricately co-ordinates the reaction mechanism, which is composed of a series of imine bond formations, leading to the production of PLP. We demonstrate that E4P (D-erythrose 4-phosphate) inhibits PfPdx1 in a dose-dependent manner. We propose that the acyclic phospho-sugar E4P, with a C1 aldehyde group similar to acyclic R5P, could interfere with R5P imine bond formations in the PfPdx1 reaction mechanism. Molecular docking and subsequent screening identified the E4P hydrazide analogue 4PEHz (4-phospho-D-erythronhydrazide), which selectively inhibited PfPdx1 with an IC50 of 43 μM. PfPdx1 contained in the heteromeric PLP synthase complex was shown to be more sensitive to 4PEHz and was inhibited with an IC50 of 16 μM. Moreover, the compound had an IC50 value of 10 μM against cultured P. falciparum intraerythrocytic parasites. To analyse further the selectivity of 4PEHz, transgenic cell lines overexpressing PfPdx1 and PfPdx2 showed that additional copies of the protein complex conferred protection against 4PEHz, indicating that the PLP synthase is directly affected by 4PEHz in vivo. These PfPdx1 inhibitors represent novel lead scaffolds which are capable of targeting PLP biosynthesis, and we propose this as a viable strategy for the development of new therapeutics against malaria.