Bioisosteric Replacement and Scaffold Hopping in Lead Generation and Optimization

Design, synthesis, and evaluation of purine and pyrimidine-based KRAS G12D inhibitors: Towards potential anticancer therapy

Oncogenic RAS mutations, commonly observed in human tumors, affect approximately 30% of cancer cases and pose a significant challenge for effective cancer treatment. Current strategies to inhibit the KRAS G12D mutation have shown limited success, emphasizing the urgent need for new therapeutic approaches. In this study, we designed and synthesized several purine and pyrimidine analogs as inhibitors for the KRAS G12D mutation. Our synthesized compounds demonstrated potent anticancer activity against cell lines with the KRAS G12D mutation, effectively impeding their growth. They also exhibited low toxicity in normal cells, indicating their selective action against cancer cells harboring the KRAS G12D mutation. Notably, the lead compound, PU1-1 induced the programmed cell death of KRAS G12D-mutated cells and reduced the levels of active KRAS and its downstream signaling proteins. Moreover, PU1-1 significantly shrunk the tumor size in a pancreatic xenograft model induced by the KRAS G12D mutation, further validating its potential as a therapeutic agent. These findings highlight the potential of purine-based KRAS G12D inhibitors as candidates for targeted cancer therapy. However, further exploration and optimization of these compounds are essential to meet the unmet clinical needs of patients with KRAS-mutant cancers.

Structure-Affinity relationships of novel σ2R/TMEM97 ligands

The sigma 2 receptor (σ2R), which was recently identified as the transmembrane protein 97 (TMEM97), is increasingly attracting interest as a possible therapeutic target for indications in neuroscience. Toward identifying novel modulators of σ2R/TMEM97, we prepared a collection of benzoxazocine, benzomorphan, and methanobenzazepine ligands related to the known bioactive norbenzomorphans DKR-1677, FEM-1689, and EES-1686 and determined their Ki values for σ2R/TMEM97 and the sigma 1 receptor (σ1R). The σ2R/TMEM97 binding affinities and selectivities relative to σ1R of these new benzoxazocine, benzomorphan, and methanobenzazepine analogs are lower, often significantly lower, than their respective norbenzomorphan counterparts, suggesting the spatial orientation of pharmacophoric substituents is critical for binding to the two proteins. The benzoxazocine, benzomorphan, and methanobenzazepine congeners of DKR-1677 and FEM-1689 tend to be weakly selective for σ2R/TMEM97 versus σ1R, whereas EES-1686 derivatives exhibit the greatest selectivity, suggesting the size and/or nature of the substituent on the nitrogen atom of the scaffold may be important for selectivity. Computational docking studies were performed for the 1S,5R-and 1R,5S-enantiomers of DKR-1677, FEM-1689, and EES-1686 and their benzoxazocine, benzomorphan, and methanobenzazepine counterparts. These computations predict that the protonated amino group of each ligand forms a highly conserved salt bridge and a H-bonding interaction with Asp29 as well as a cation-π interaction with Tyr150 of σ2R/TMEM97. These electrostatic interactions are major driving forces for binding to σ2R/TMEM97 and are similar, though not identical, for each ligand. Other interactions within the well-defined binding pocket also tend to be comparable, but there are some major differences in how the hydrophobic aryl groups of various ligands interact with the protein surface external to the binding pocket. Overall, these studies show that the orientations of aryl and N-substituents on the norbenzomorphan and related scaffolds are important determinants of binding affinity of σ2R/TMEM97 ligands, and small changes can have significant effects upon binding profiles.

Genetic Algorithm-Based Receptor Ligand: A Genetic Algorithm-Guided Generative Model to Boost the Novelty and Drug-Likeness of Molecules in a Sampling Chemical Space

Deep learning-based de novo molecular design has recently gained significant attention. While numerous DL-based generative models have been successfully developed for designing novel compounds, the majority of the generated molecules lack sufficiently novel scaffolds or high drug-like profiles. The aforementioned issues may not be fully captured by commonly used metrics for the assessment of molecular generative models, such as novelty, diversity, and quantitative estimation of the drug-likeness score. To address these limitations, we proposed a genetic algorithm-guided generative model called GARel (genetic algorithm-based receptor-ligand interaction generator), a novel framework for training a DL-based generative model to produce drug-like molecules with novel scaffolds. To efficiently train the GARel model, we utilized dense net to update the parameters based on molecules with novel scaffolds and drug-like features. To demonstrate the capability of the GARel model, we used it to design inhibitors for three targets: AA2AR, EGFR, and SARS-Cov2. The results indicate that GARel-generated molecules feature more diverse and novel scaffolds and possess more desirable physicochemical properties and favorable docking scores. Compared with other generative models, GARel makes significant progress in balancing novelty and drug-likeness, providing a promising direction for the further development of DL-based de novo design methodology with potential impacts on drug discovery.

Scaffold hopping for designing of potent and selective CYP1B1 inhibitors to overcome docetaxel resistance: synthesis and evaluation

Cytochrome P450 1B1, a tumor-specific overexpressed enzyme, significantly impairs the pharmacokinetics of several commonly used anticancer drugs including docetaxel, paclitaxel and cisplatin, leading to the problem of resistance to these drugs. Currently, there is no CYP1B1 inhibition-based adjuvant therapy available to treat this resistance problem. Hence, in the current study, exhaustive in-silico studies including scaffold hopping followed by molecular docking, three-dimensional quantitative structure-activity relationships (3D-QSAR), molecular dynamics and free energy perturbation studies were carried out to identify potent and selective CYP1B1 inhibitors. Initially, scaffold hopping analysis was performed against a well-reported potent and selective CYP1B1 inhibitor (i.e. compound 3n). A total of 200 scaffolds were identified along with their shape and field similarity scores. The top three scaffolds were further selected on the basis of these scores and their synthesis feasibility to design some potent and selective CYP1B1 inhibitors using the aforementioned in-silico techniques. Designed molecules were further synthesized to evaluate their CYP1B1 inhibitory activity and docetaxel resistance reversal potential against CYP1B1 overexpressed drug resistance MCF-7 cell line. In-vitro results indicated that compounds 2a, 2c and 2d manifested IC50 values for CYP1B1 ranging from 0.075, 0.092 to 0.088 μM with at least 10-fold selectivity. At low micromolar concentrations, compounds 1e, 1f, 2a and 2d exhibited promising cytotoxic effects in the docetaxel-resistant CYP1B1 overexpressed MCF-7 cell line. In particular, compound 2a is most effective in reversing the resistance with IC50 of 29.0 ± 3.6 μM. All of these discoveries could pave the way for the development of adjuvant therapy capable of overcoming CYP1B1-mediated resistance.Communicated by Ramaswamy H. Sarma.

BioisoIdentifier: an online free tool to investigate local structural replacements from PDB

Within the realm of contemporary medicinal chemistry, bioisosteres are empirically used to enhance potency and selectivity, improve adsorption, distribution, metabolism, excretion and toxicity profiles of drug candidates. It is believed that bioisosteric know-how may help bypass granted patents or generate novel intellectual property for commercialization. Beside the synthetic expertise, the drug discovery process also depends on efficient in silico tools. We hereby present BioisoIdentifier (BII), a web server aiming to uncover bioisosteric information for specific fragment. Using the Protein Data Bank as source, and specific substructures that the user attempt to surrogate as input, BII tries to find suitable fragments that fit well within the local protein active site. BII is a powerful computational tool that offers the ligand design ideas for bioisosteric replacing. For the validation of BII, catechol is conceived as model fragment attempted to be replaced, and many ideas are successfully offered. These outputs are hierarchically grouped according to structural similarity, and clustered based on unsupervised machine learning algorithms. In summary, we constructed a user-friendly interface to enable the viewing of top-ranking molecules for further experimental exploration. This makes BII a highly valuable tool for drug discovery. The BII web server is freely available to researchers and can be accessed at http://www.aifordrugs.cn/index/ . Scientific Contribution: By designing a more optimal computational process for mining bioisosteric replacements from the publicly accessible PDB database, then deployed on a web server for throughly free access for researchers. Additionally, machine learning methods are applied to cluster the bioisosteric replacements searched by the platform, making a scientific contribution to facilitate chemists' selection of appropriate bioisosteric replacements. The number of bioisosteric replacements obtained using BII is significantly larger than the currently available platforms, which expanding the search space for effective local structural replacements.

Molecular medicinal insights into scaffold hopping-based drug discovery success

In both academia and the pharmaceutical industry, innovative hypotheses, methodologies and technologies that can shorten the drug research and development, leading to higher success rates, are vital. In this review, we demonstrate how innovative variations of the scaffold-hopping strategy have been used to create new druggable molecular spaces, drugs, clinical candidates, preclinical candidates, and bioactive agents. We also analyze molecular modulations that enabled improvements of the pharmacodynamic (PD), physiochemical, and pharmacokinetic (PK) properties (P3 properties) of the drugs resulting from these scaffold-hopping strategies.

Application of Artificial Intelligence at All Stages of Bone Tissue Engineering

The development of artificial intelligence (AI) has revolutionized medical care in recent years and plays a vital role in a number of areas, such as diagnostics and forecasting. In this review, we discuss the most promising areas of AI application to the field of bone tissue engineering and prosthetics, which can drastically benefit from AI-assisted optimization and patient personalization of implants and scaffolds in ways ranging from visualization and real-time monitoring to the implantation cases prediction, thereby leveraging the compromise between specific architecture decisions, material choice, and synthesis procedure. With the emphasized crucial role of accuracy and robustness of developed AI algorithms, especially in bone tissue engineering, it was shown that rigorous validation and testing, demanding large datasets and extensive clinical trials, are essential, and we discuss how through developing multidisciplinary cooperation among biology, chemistry with materials science, and AI, these challenges can be addressed.

3D-QSAR assisted identification of selective CYP1B1 inhibitors: an effective bioisosteric replacement/molecular docking/electrostatic complementarity analysis

Cytochrome P450-1B1 is a majorly overexpressed drug-metabolizing enzyme in tumors and is responsible for inactivation and subsequent resistance to a variety of anti-cancer drugs, i.e., docetaxel, tamoxifen, and cisplatin. In the present study, a 3D quantitative structure-activity relationship (3D-QSAR) model has been constructed for the identification, design, and optimization of novel CYP1B1 inhibitors. The model has been built using a set of 148 selective CYP1B1 inhibitors. The developed model was evaluated based on certain statistical parameters including q2 and r2 which showed the acceptable predictive and descriptive capability of the generated model. The developed 3D-QSAR model assisted in understanding the key molecular fields which were firmly related to the selective CYP1B1 inhibition. A theoretic approach for the generation of new lead compounds with optimized CYP1B1 receptor affinity has been performed utilizing bioisosteric replacement analysis. These generated molecules were subjected to a developed 3D-QSAR model to predict the inhibitory activity potentials. Furthermore, these compounds were scrutinized through the activity atlas model, molecular docking, electrostatic complementarity, molecular dynamics, and waterswap analysis. The final hits might act as selective CYP1B1 inhibitors which could address the issue of resistance. This 3D-QSAR includes several chemically diverse selective CYP1B1 receptor ligands and well accounts for the individual ligand's inhibition affinities. These features of the developed 3D-QSAR model will ensure future prospective applications of the model to speed up the identification of new potent and selective CYP1B1 receptor ligands.

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.

3-Methoxy-2-phenylimidazo[1,2-b]pyridazines highly active against Mycobacterium tuberculosis and Mycobacterium marinum

A series of 3-methoxy-2-phenylimidazo[1,2-b]pyridazine derivatives which were highly active against autoluminescent Mycobacterium tuberculosis (Mtb) and Mycobacterium marinum (Mm) in an in vitro assay were identified. SAR analysis showed that the most active compounds, which included a phenyl group bearing fluoro substituent(s) at C2, a methoxy function at C3, and a benzyl-heteroatom moiety at C6, exhibited in vitro MIC90 values generally around 0.63-1.26 μM against Mtb and Mm. However, these compounds were inactive against Mtb in vivo (mice), and investigations revealed very short metabolic half-lives (<10 min) when incubated with mouse liver microsomes. Multiple observations of side products produced from oxidative cleavage of the imidazole moiety during the chemical synthesis work suggested that this is a likely metabolic pathway leading to the lack of observed activity in vivo.

Strategies for the Discovery of Oxazolidinone Antibacterial Agents: Development and Future Perspectives

Oxazolidinones represent a significant class of synthetic bacterial protein synthesis inhibitors that are primarily effective against Gram-positive bacteria. The commercial success of linezolid, the first FDA-approved oxazolidinone antibiotic, has motivated researchers to develop more potent oxazolidinones by employing various drug development strategies to fight against antimicrobial resistance, some of which have shown promising results. Thus, this Perspective aims to discuss the strategies employed in constructing oxazolidinone-based antibacterial agents and summarize recent advances in discovering oxazolidinone antibiotics to provide valuable insights for potentially developing next-generation oxazolidinone antibacterial agents or other pharmaceuticals.

Expanding the Role of Boron in New Drug Chemotypes: Properties, Chemistry, Pharmaceutical Potential of Hemiboronic Naphthoids

New chemotypes and bioisosteres can open a new chemical space in drug discovery and help meet an urgent demand for novel agents to fight infections and other diseases. With the aim of identifying new boron-containing drug chemotypes, this article details a comprehensive evaluation of the pseudoaromatic hemiboronic naphthoids, benzoxaza- and benzodiazaborines. Relevant physical properties in aqueous media (acidity, solubility, log P, and stability) of prototypic members of four subclasses were determined. Both scaffolds are amenable to common reactions used in drug discovery, such as chemoselective Suzuki-Miyaura, Chan-Lam, and amidation reactions. Small model libraries were prepared to assess the scope of these transformations, and the entire collection was screened for antifungal (Candida albicans) and antibacterial activity (MRSA, Escherichia coli), unveiling promising benzoxazaborines with low micromolar minimum inhibitory concentration values. Select DMPK assays of representative compounds suggest promising drug-like behavior for all four subclasses. Moreover, several drug isosteres were evaluated for anti-inflammatory and anticancer activity as appropriate.

Preclinical Evaluation of an Imidazole-Linked Heterocycle for Alzheimer's Disease

Humanity is facing a vast prevalence of neurodegenerative diseases, with Alzheimer's disease (AD) being the most dominant, without efficacious drugs, and with only a few therapeutic targets identified. In this scenario, we aim to find molecular entities that modulate imidazoline I2 receptors (I2-IRs) that have been pointed out as relevant targets in AD. In this work, we explored structural modifications of well-established I2-IR ligands, giving access to derivatives with an imidazole-linked heterocycle as a common key feature. We report the synthesis, the affinity in human I2-IRs, the brain penetration capabilities, the in silico ADMET studies, and the three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of this new bunch of I2-IR ligands. Selected compounds showed neuroprotective properties and beneficial effects in an in vitro model of Parkinson's disease, rescued the human dopaminergic cell line SH-SY5Y from death after treatment with 6-hydroxydopamine, and showed crucial anti-inflammatory effects in a cellular model of neuroinflammation. After a preliminary pharmacokinetic study, we explored the action of our representative 2-(benzo[b]thiophen-2-yl)-1H-imidazole LSL33 in a mouse model of AD (5xFAD). Oral administration of LSL33 at 2 mg/Kg for 4 weeks ameliorated 5XFAD cognitive impairment and synaptic plasticity, as well as reduced neuroinflammation markers. In summary, this new I2-IR ligand that promoted beneficial effects in a well-established AD mouse model should be considered a promising therapeutic strategy for neurodegeneration.

Recent Development of Novel Aminoethyl-Substituted Chalcones as Potential Drug Candidates for the Treatment of Alzheimer's Disease

No drug on the market, as a single entity, participates in different pathways involved in the pathology of Alzheimer's disease. The current study is aimed at the exploration of multifunctional chalcone derivatives which can act on multiple targets involved in Alzheimer's disease. A series of novel aminoethyl-substituted chalcones have been developed using in silico approaches (scaffold morphing, molecular docking, and ADME) and reported synthetic methods. The synthesized analogs were characterized and evaluated biologically using different in vitro assays against AChE, AGEs, and radical formation. Among all compounds, compound PS-10 was found to have potent AChE inhibitory activity (IC50 = 15.3 nM), even more than the standard drug (IC50 = 15.68 nM). Further, the in vivo evaluation of PS-10 against STZ-induced dementia in rats showed memory improvement (Morris Water Maze test) in rats. Also, PS-10 inhibited STZ-induced brain AChE activity and oxidative stress, further strengthening the observed in vitro effects. Further, the molecular dynamic simulation studies displayed the stability of the PS-10 and AChE complex. The novel aminoethyl-substituted chalcones might be considered potential multifunctional anti-Alzheimer's molecules.

Scaffold Morphing and In Silico Design of Potential BACE-1 (β-Secretase) Inhibitors: A Hope for a Newer Dawn in Anti-Alzheimer Therapeutics

Alzheimer's disease (AD) is the prime cause of 65-80% of dementia cases and is caused by plaque and tangle deposition in the brain neurons leading to brain cell degeneration. β-secretase (BACE-1) is a key enzyme responsible for depositing extracellular plaques made of β-amyloid protein. Therefore, efforts are being applied to develop novel BACE-1 enzyme inhibitors to halt plaque build-up. In our study, we analyzed some Elenbecestat analogues (a BACE-1 inhibitor currently in clinical trials) using a structure-based drug design and scaffold morphing approach to achieve a superior therapeutic profile, followed by in silico studies, including molecular docking and pharmacokinetics methodologies. Among all the designed compounds, SB306 and SB12 showed good interactions with the catalytic dyad motifs (Asp228 and Asp32) of the BACE-1 enzyme with drug-likeliness properties and a high degree of thermodynamic stability confirmed by the molecular dynamic and stability of the simulated system indicating the inhibitory nature of the SB306 and SB12 on BACE 1.

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.

Cryo-EM structures of actin binding proteins as tool for drug discovery

Cellular actin dynamic is controlled by a plethora of actin binding proteins (ABPs), including actin nucleating, bundling, cross-linking, capping, and severing proteins. In this review, regulation of actin dynamics by ABPs will be introduced, and the role of the F-actin severing protein cofilin-1 and the F-actin bundling protein L-plastin in actin dynamics discussed in more detail. Since up-regulation of these proteins in different kinds of cancers is associated with malignant progression of cancer cells, we suggest the cryogenic electron microscopy (Cryo-EM) structure of F- actin with the respective ABP as template for in silico drug design to specifically disrupt the interaction of these ABPs with F-actin.

Recent developments of imidazo[1,2-a]pyridine analogues as antituberculosis agents

Over the past 2000 years, tuberculosis (TB) has killed more people than any other infectious disease. In 2021, TB claimed 1.6 million lives worldwide, making it the second leading cause of death from an infectious disease after COVID-19. Unfortunately, TB drug discovery research was neglected in the last few decades of the twentieth century. Recently, the World Health Organization has taken the initiative to develop new TB drugs. Imidazopyridine, an important fused bicyclic 5,6 heterocycle has been recognized as a "drug prejudice" scaffold for its wide range of applications in medicinal chemistry. A few examples of imidazo[1,2-a]pyridine exhibit significant activity against multidrug-resistant TB (MDR-TB) and extensively drug-resistant TB (XDR-TB). Here, we critically review anti-TB compounds of the imidazo[1,2-a]pyridine class by discussing their development based on the structure-activity relationship, mode-of-action, and various scaffold hopping strategies over the last decade, which is identified as a renaissance era of TB drug discovery research.

Molecular Docking and Dynamics Simulation Studies Predict Potential Anti-ADAR2 Inhibitors: Implications for the Treatment of Cancer, Neurological, Immunological and Infectious Diseases

Altered RNA editing has been linked to several neurodevelopmental disorders, including autism spectrum disorder (ASD) and intellectual disability, in addition to depression, schizophrenia, some cancers, viral infections and autoimmune disorders. The human ADAR2 is a potential therapeutic target for managing these various disorders due to its crucial role in adenosine to inosine editing. This study applied consensus scoring to rank potential ADAR2 inhibitors after performing molecular docking with AutoDock Vina and Glide (Maestro), using a library of 35,161 compounds obtained from traditional Chinese medicine. A total of 47 compounds were predicted to be good binders of the human ADAR2 and had insignificant toxicity concerns. Molecular dynamics (MD) simulations, including the molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) procedure, also emphasized the binding of the shortlisted compounds. The potential compounds had plausible binding free energies ranging from −81.304 to −1068.26 kJ/mol from the MM/PBSA calculations. ZINC000085511995, a naphthoquinone had more negative binding free energy (−1068.26 kJ/mol) than inositol hexakisphosphate (IHP) [−873.873 kJ/mol], an agonist and a strong binder of ADAR2. The potential displacement of IHP by ZINC000085511995 in the IHP binding site of ADAR2 could be explored for possible deactivation of ADAR2. Bayesian-based biological activity prediction corroborates the neuropharmacological, antineoplastic and antiviral activity of the potential lead compounds. All the potential lead compounds, except ZINC000014612330 and ZINC000013462928, were predicted to be inhibitors of various deaminases. The potential lead compounds also had probability of activity (Pa) > 0.442 and probability of inactivity (Pi) < 0.116 values for treating acute neurologic disorders, except for ZINC000085996580 and ZINC000013462928. Pursuing these compounds for their anti-ADAR2 activities holds a promising future, especially against neurological disorders, some cancers and viral infections caused by RNA viruses. Molecular interaction, hydrogen bond and per-residue decomposition analyses predicted Arg400, Arg401, Lys519, Trp687, Glu689, and Lys690 as hot-spot residues in the ADAR2 IHP binding site. Most of the top compounds were observed to have naphthoquinone, indole, furanocoumarin or benzofuran moieties. Serotonin and tryptophan, which are beneficial in digestive regulation, improving sleep cycle and mood, are indole derivatives. These chemical series may have the potential to treat neurological disorders, prion diseases, some cancers, specific viral infections, metabolic disorders and eating disorders through the disruption of ADAR2 pathways. A total of nine potential lead compounds were shortlisted as plausible modulators of ADAR2.

A Perspective on Explanations of Molecular Prediction Models

Chemists can be skeptical in using deep learning (DL) in decision making, due to the lack of interpretability in "black-box" models. Explainable artificial intelligence (XAI) is a branch of artificial intelligence (AI) which addresses this drawback by providing tools to interpret DL models and their predictions. We review the principles of XAI in the domain of chemistry and emerging methods for creating and evaluating explanations. Then, we focus on methods developed by our group and their applications in predicting solubility, blood-brain barrier permeability, and the scent of molecules. We show that XAI methods like chemical counterfactuals and descriptor explanations can explain DL predictions while giving insight into structure-property relationships. Finally, we discuss how a two-step process of developing a black-box model and explaining predictions can uncover structure-property relationships.

Synthesis of the 1,5-Benzothiazepane Scaffold - Established Methods and New Developments

The 1,5-benzothiazepane structure is an important heterocyclic moiety present in a variety of commercial drugs and pharmaceuticals. This privileged scaffold exhibits a diversity of biological activities, including antimicrobial, antibacterial, anti-epileptic, anti-HIV, antidepressant, antithrombotic and anticancer properties. Its important pharmacological potential renders research into the development of new and efficient synthetic methods of high relevance. In the first part of this review, an overview of different synthetic approaches toward 1,5-benzothiazepane and its derivatives is provided, ranging from established protocols to recent (enantioselective) methods that promote sustainability. In the second part, several structural characteristics influencing biological activity are briefly explored, providing a few insights into the structure-activity relationships of these compounds.

In silico studies on recreational drugs: 3D quantitative structure activity relationship prediction of classified and de novo designer benzodiazepines

Currently, increasing availability and popularity of designer benzodiazepines (DBZDs) constitutes a primary threat to public health. To assess this threat, the biological activity/potency of DBZDs was investigated using in silico studies. Specific Quantitative Structure Activity Relationship (QSAR) models were developed in Forge™ for the prediction of biological activity (IC50 ) on the γ-aminobutyric acid A receptor (GABA-AR) of previously identified classified and unclassified DBDZs. A set of new potential ligands resulting from scaffold hopping studies conducted with MOE® was also evaluated. Two generated QSAR models (i.e. 3D-field QSAR and RVM) returned very good performance statistics (r2  = 0.98 [both] and q2  = 0.75 and 0.72, respectively). The DBZDs predicted to be the most active were flubrotizolam, clonazolam, pynazolam and flucotizolam, consistently with what reported in literature and/or drug discussion fora. The scaffold hopping studies strongly suggest that replacement of the pendant phenyl moiety with a five-membered ring could increase biological activity and highlight the existence of a still unexplored chemical space for DBZDs. QSAR could be of use as a preliminary risk assessment model for (newly) identified DBZDs, as well as scaffold hopping for the creation of computational libraries that could be used by regulatory bodies as support tools for scheduling procedures.

ReMODE: a deep learning-based web server for target-specific drug design

Deep learning (DL) and machine learning contribute significantly to basic biology research and drug discovery in the past few decades. Recent advances in DL-based generative models have led to superior developments in de novo drug design. However, data availability, deep data processing, and the lack of user-friendly DL tools and interfaces make it difficult to apply these DL techniques to drug design. We hereby present ReMODE (Receptor-based MOlecular DEsign), a new web server based on DL algorithm for target-specific ligand design, which integrates different functional modules to enable users to develop customizable drug design tasks. As designed, the ReMODE sever can construct the target-specific tasks toward the protein targets selected by users. Meanwhile, the server also provides some extensions: users can optimize the drug-likeness or synthetic accessibility of the generated molecules, and control other physicochemical properties; users can also choose a sub-structure/scaffold as a starting point for fragment-based drug design. The ReMODE server also enables users to optimize the pharmacophore matching and docking conformations of the generated molecules. We believe that the ReMODE server will benefit researchers for drug discovery. ReMODE is publicly available at http://cadd.zju.edu.cn/relation/remode/ .

Indoline-5-Sulfonamides: A Role of the Core in Inhibition of Cancer-Related Carbonic Anhydrases, Antiproliferative Activity and Circumventing of Multidrug Resistance

The overexpression and activity of carbonic anhydrase (CA, EC 4.2.1.1) isoforms CA IX and CA XII promote the accumulation of exceeding protons and acidosis in the extracellular tumor environment. Sulfonamides are effective inhibitors of most families of CAs. In this study, using scaffold-hopping, indoline-5-sulfonamide analogs 4a-u of the CA IX-selective inhibitor 3 were designed and synthesized to evaluate their biological properties. 1-Acylated indoline-5-sulfonamides demonstrated inhibitory activity against tumor-associated CA IX and XII with K_I values up to 132.8 nM and 41.3 nM. Compound 4f, as one of the most potent inhibitors of CA IX and XII, exhibits hypoxic selectivity, suppressing the growth of MCF7 cells at 12.9 µM, and causes partial inhibition of hypoxia-induced CA IX expression in A431 skin cancer cells. 4e and 4f reverse chemoresistance to doxorubicin of K562/4 with overexpression of P-gp.

Papaya Leaf Extracts as Potential Dengue Treatment: An In-Silico Study

Dengue fever (DF), dengue hemorrhagic fever (DHF), and dengue shock syndrome (DSS) cause serious public health problems, with nearly 390 million people affected and 20,000 deaths per year in tropical and subtropical countries. Despite numerous attempts, no antiviral drug or vaccine is currently available to combat the manifestation. The challenge of discovering an efficient vaccine is enhanced by the surplus presence of efficient vectors and drug resistance from the virus. For centuries, papaya (Carica papaya) extracts have been traditionally used to treat DF, DHF, and DSS. In the present study, we systematically investigated seven compounds isolated from papaya leaf extract with regard to their potential as inhibitors for non-structural (NS) proteins, NS3 and NS5, which play a crucial role in viral RNA replication. The computational tools applied stretched across classical molecular docking, molecular dynamics (MD) simulations and SwissADME used to calculate binding affinities; binding free energies; Absorption, Distribution, Metabolism, and Excretion (ADME); and drug-likeness properties, thus, identifying Kaempferol, Chlorogenic acid, and Quercetin as potential candidates, with Kaempferol and Quercetin scoring best. Therefore, for the Kaempferol and Quercetin complexes, hybrid quantum mechanical/molecular mechanical (QM/MM) geometry and frequency calculations were performed, followed by the local mode analysis developed in our group to quantify Kaempferol-NS and Quercetin-NS hydrogen bonding. Given the non-toxic nature and the wide availability of the Kaempferol and Quercetin papaya extract in almost all of the susceptible regions, and our results showing high NS3 and NS5 binding affinities and energies, strong hydrogen bonding with both NS3 and NS5, and excellent ADME properties, we suggest Kaempferol and Quercetin as a strong NS3 and NS5 inhibitor to be further investigated in vitro.

Silicon-Containing Complex II Acaricides─Design, Synthesis, and Pharmacological Optimization

Bioisosteric replacement has been proven to be a powerful strategy in life science research. In this review, general aspects of carbon-silicon bioisosteric substitution and its applications in pharmaceutical and crop protection research are described. Carbon and their silicon analogues possess similar intrinsic properties. Replacing carbon with silicon in pharmaceuticals and pesticides has shown to result in positive effects on efficacy and selectivity, physicochemical properties, and bioavailability and also to eliminate or improve human or environmental safety properties as well as to provide novelty and new intellectual property in many cases. Furthermore, the application of carbon-silicon substitution in the search for new complex II acaricides is highlighted. This research led to the discovery of sila-cyflumetofen 23a and other silicon-containing analogues of cyflumetofen that match or exceed the acaricidal activity of cyflumetofen. The molecular design strategy, synthetic aspects, biological activity, computational modeling work, and structure-activity relationships will be discussed.

Application of Molecular-Modeling, Scaffold-Hopping, and Bioisosteric Approaches to the Discovery of New Heterocyclic Picolinamides

Macrocyclic natural products and their derivatives are a valuable source for biologically active crop protection products and have had significant impact on the development of conventional agrochemicals. However, they can be challenging starting points for lead-generation efforts because of their size, structural complexity, and developability. Using molecular modeling and electrostatic analysis, alternative bicyclic isosteres were identified as replacements for the antifungal nine-membered macrocycle UK-2A. By application of a structure-based conformational approach, a series of heterocyclic replacements were derivatized to deliver promising fungicidal activity and scaffold bioisosteres were further diversified to investigate structure-activity relationships.

Bioisosterism and Scaffold Hopping in Modern Nematicide Research

The application of agrochemicals is critical to global food safety. Nowadays, environmentally friendly green agrochemicals are the trend in field crop protection. The research and development of nematicides absorbed more attention as a typical representation of agrochemicals. This review describes the origin of recently commercialized nematicides, the application of bioisosterism and scaffold hopping in the discovery and optimization of agrochemicals, especially nematicides, and novel bioisosteric design strategies for the identification of fluensulfone analogues. Pesticide repurposing, high-throughput screening, computer-aided drug design, and incorporation of known pharmacophoric fragments have been the most successful approach for the discovery of new nematicides. As outlined, the strategies of bioisosteric replacements and scaffold hopping have been very successful approaches in the search for new nematicides for sustainable crop protection. In the exploration of novel fluensulfone analogues with nematicidal activity, bioisosteric replacement of sulfone by amide, chain extension by insertion of a methylene group, and reversal of the amide group have proven to be successful approaches and yielded new and highly active fluensulfone analogues. These attempts might result in compounds with an optimal balance of steric, hydrophobic, electronic, and hydrogen-bonding properties and contribute to deal with the complex problem during the research and development of new nematicides. Further ideas are also put forward to provide new approaches for the molecular design of nematicides.

Scaffold-Hopping Approach To Identify New Chemotypes of Dimpropyridaz

Dimpropyridaz is a pyrazole carboxamide insecticide with a novel mode of action, currently under worldwide development by BASF, providing excellent activity against sucking pests. A series of dimpropyridaz analogues were designed to investigate the impact of bioisosteric heterocyclic replacements on the biological activity and molecular properties. Focus was given to prepare analogues where the 4-pyridazinyl moiety was replaced by 5-pyrimidinyl, 2-pyrimidinyl, 3-pyridazinyl, and 2-pyrazinyl groups. Five different synthetic routes were developed for the preparation of these analogues, delivering the target compounds in moderate to good yields. We explained some aspects of the observed structure-activity relationship by a density functional theory (DFT) calculation and DFT-derived Multiwfn and VMD program models. These findings provide first insights into the important role of the 4-pyridazinyl heterocyclic moiety in the pyrazole carboxamide insecticide chemical class and the mechanism of action of dimpropyridaz.

Bioisosteric-Replacement-Driven Lead Optimization of Tyclopyrazoflor

A heterocyclic replacement strategy has been applied attempting to further optimize the biological properties of the new insecticide tyclopyrazoflor, a pyridinyl-pyrazole compound with excellent activity against sap-feeding insects recently discovered by Corteva. Our research program on this novel chemical class of insecticide was targeted to evaluate the importance of both the pyridyl and pyrazole moieties on the biological activity, properties, and structure-activity. Synthetic methodologies have been developed delivering our designed target compounds in moderate to good yields. In our attempt to rationalize the biological activity differences of tyclopyrazoflor analogues, a potential surface analysis based on a density functional theory (DFT) calculation and DFT-derived Multiwfn and VMD program model was performed. Furthermore, the SwissBioisostere online database was applied as a great source for inspiration for bioisosteric replacements.

Bridging Known and Unknown Unknowns: From Natural Products and Their Mimics to Unmet Needs in Neuroscience

Scientific excursions into the unknown are often characterized by unanticipated twists and turns that may lead in directions that never could have been predicted. Decisions made during the course of these explorations determine what we discover. This Account chronicles one such journey that began with a challenge encountered during the synthesis of a natural product and then unfolded over more than 30 years to focus on unmet needs in neuroscience. Specifically, while developing a concise approach to tetrahydroalstonine, a heteroyohimboid alkaloid having α-adrenergic activity, we faced the predicament of assembling a key intermediate. Solving this problem resulted in the serendipitous discovery of the vinylogous Mannich reaction and a productive program wherein we used this powerful construction as a key step in the syntheses of numerous alkaloids. However, we also realized that lessons learned from the synthesis of tetrahydroalstonine could be generalized to invent a new strategy for preparing diverse collections of substituted nitrogen heterocycles that could be screened against biological targets. The approach featured the combination of several reactants in a multicomponent assembly process to give a functionalized intermediate that could be elaborated by various ring-forming reactions to give heterocyclic scaffolds that could be further diversified. Screening these compound sets against a broad range of biological targets revealed some intriguing hits, but none of them led to a productive collaboration in translational research. Notwithstanding this setback, we screened curated members of our collections against proteins in the central nervous system and discovered some substituted B-norbenzomorphans that were selective for the enigmatic sigma-2 receptor (σ₂R), an understudied protein that had been primarily associated with cancer. With scant knowledge of its role in neuroscience, we posited that small molecules that bind to σ₂R might be neuroprotective, thus launching a new venture. In parallel investigations we prepared analogues of the initial hits, explored their effects in animal models of neurodegenerative and neurological conditions, and identified σ₂R as transmembrane protein 97 (TMEM97). After first establishing the neuroprotective effects of several σ₂R/TMEM97 ligands in a transgenic Caenorhabditis elegans model of neurodegeneration, we showed that one of these has procognitive effects and reduces levels of proinflammatory cytokines in a transgenic mouse model of Alzheimer's disease. We then identified a closely related σ₂R/TMEM97 ligand that mitigates hippocampal-dependent memory deficits, prevents axon degeneration, and protects neurons and oligodendrocytes after traumatic brain injury. In a recent study, this compound was shown to protect retinal ganglion cells from retinal ischemia/reperfusion injury. In other collaborative investigations, we have shown that related, but structurally distinct, σ₂R/TMEM97 ligands alleviate neuropathic pain, while a σ₂R/TMEM97 ligand representing yet another chemotype reduces impairments associated with alcohol withdrawal. More recently, we have shown that σ₂R/TMEM97 ligands enhance survival of cortical neurons in a neuronal model of Huntington's disease. Translational and mechanistic studies in these and other areas are in progress. Solving a problem we faced in natural product synthesis thus served as an unexpected gateway to discoveries that could lead to entirely new approaches to treat neurodegenerative and neurological conditions by targeting σ₂R/TMEM97, a protein that has never been associated with these afflictions.

Scaffold hopping from indoles to indazoles yields dual MCL-1/BCL-2 inhibitors from MCL-1 selective leads

Overexpression of the anti-apoptotic BCL-2 proteins is associated with the development and progression of a range of cancers. Venetoclax, an FDA-approved BCL-2 inhibitor, is fast becoming the standard-of-care for acute myeloid leukemia and chronic lymphocytic leukemia. However, the median survival offered by venetoclax is only 18 months (as part of a combination therapy regimen), and one of the primary culprits for this is the concomitant upregulation of sister anti-apoptotic proteins, in particular MCL-1 (and BCL-xL), which provides an escape route that manifests as venetoclax resistance. Since inhibition of BCL-xL leads to thrombocytopenia, we believe that a dual MCL-1/BCL-2 inhibitor may provide an enhanced therapeutic effect relative to a selective BCL-2 inhibitor. Beginning with a carboxylic acid-containing literature compound that is a potent inhibitor of MCL-1 and a moderate inhibitor of BCL-2, we herein describe our efforts to develop dual inhibitors of MCL-1 and BCL-2 by scaffold hopping from an indole core to an indazole framework. Subsequently, further elaboration of our novel N2-substituted, indazole-3-carboxylic acid lead into a family of indazole-3-acylsulfonamides resulted in improved inhibition of both MCL-1 and BCL-2, possibly through occupation of the p4 pocket, with minimal or no inhibition of BCL-xL.

RELATION: A Deep Generative Model for Structure-Based De Novo Drug Design

Deep learning (DL)-based de novo molecular design has recently gained considerable traction. Many DL-based generative models have been successfully developed to design novel molecules, but most of them are ligand-centric and the role of the 3D geometries of target binding pockets in molecular generation has not been well-exploited. Here, we proposed a new 3D-based generative model called RELATION. In the RELATION model, the BiTL algorithm was specifically designed to extract and transfer the desired geometric features of the protein-ligand complexes to a latent space for generation. The pharmacophore conditioning and docking-based Bayesian sampling were applied to efficiently navigate the vast chemical space for the design of molecules with desired geometric properties and pharmacophore features. As a proof of concept, the RELATION model was used to design inhibitors for two targets, AKT1 and CDK2. The calculation results demonstrated that the RELATION model could efficiently generate novel molecules with favorable binding affinity and pharmacophore features.

Free energy perturbation (FEP)-guided scaffold hopping

Scaffold hopping refers to computer-aided screening for active compounds with different structures against the same receptor to enrich privileged scaffolds, which is a topic of high interest in organic and medicinal chemistry. However, most approaches cannot efficiently predict the potency level of candidates after scaffold hopping. Herein, we identified potent PDE5 inhibitors with a novel scaffold via a free energy perturbation (FEP)-guided scaffold-hopping strategy, and FEP shows great advantages to precisely predict the theoretical binding potencies ΔG FEP between ligands and their target, which were more consistent with the experimental binding potencies ΔG EXP (the mean absolute deviations| Δ G FEP - Δ G EXP |  < 2 kcal/mol) than those ΔG MM-PBSA or ΔG MM-GBSA predicted by the MM-PBSA or MM-GBSA method. Lead L12 had an IC50 of 8.7 nmol/L and exhibited a different binding pattern in its crystal structure with PDE5 from the famous starting drug tadalafil. Our work provides the first report via the FEP-guided scaffold hopping strategy for potent inhibitor discovery with a novel scaffold, implying that it will have a variety of future applications in rational molecular design and drug discovery.

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.

Discovery of Novel TRPM8 Blockers Suitable for the Treatment of Somatic and Ocular Painful Conditions: A Journey through pKa and LogD Modulation

Transient receptor potential melastatin 8 (TRPM8) is crucially involved in pain modulation and perception, and TRPM8 antagonists have been proposed as potential therapeutic approaches for pain treatment. Previously, we developed two TRPM8 antagonists and proposed them as drug candidates for topical and systemic pain treatment. Here, we describe the design and synthesis of these two TRPM8 antagonists (27 and 45) and the rational approach of modulation/replacement of bioisosteric chemical groups, which allowed us to identify a combination of narrow ranges of pK_a and LogD values that were crucial to ultimately optimize their potency and metabolic stability. Following the same approach, we then pursued the development of new TRPM8 antagonists suitable for the topical treatment of ocular painful conditions and identified two new compounds (51 and 59), N-alkoxy amide derivatives, that can permeate across ocular tissue and reduce the behavioral responses induced by the topical ocular menthol challenge in vivo.

SwissBioisostere 2021: updated structural, bioactivity and physicochemical data delivered by a reshaped web interface

At several stages of drug discovery, bioisosteric replacement is a common and efficient practice to find new bioactive chemotypes or to optimize series of molecules toward drug candidates. The critical steps consisting in selecting which molecular moiety should be replaced by which other chemical fragment is often relying on the expertise of specialists. Nowadays, valuable support can be obtained through the wealth of dedicated structural and knowledge data. The present article details the update of SwissBioisostere, a database of >25 millions of unique molecular replacements with data on bioactivity, physicochemistry, chemical and biological contexts extracted from the literature and related resources. The content of the database together with analysis and visualization capacities is freely available at www.swissbioisostere.ch.

The role of DMPK science in improving pharmaceutical research and development efficiency

The successful regulatory authority approval rate of drug candidates in the drug development pipeline is crucial for determining pharmaceutical research and development (R&D) efficiency. Regulatory authorities include the US Food and Drug Administration (FDA), European Medicines Agency (EMA), and Pharmaceutical and Food Safety Bureau Japan (PFSB), among others. Optimal drug metabolism and pharmacokinetics (DMPK) properties influence the progression of a drug candidate from the preclinical to the clinical phase. In this review, we provide a comprehensive assessment of essential concepts, methods, improvements, and challenges in DMPK science and its significance in drug development. This information provides insights into the association of DMPK science with pharmaceutical R&D efficiency.

Deep generative design with 3D pharmacophoric constraints

Generative models have increasingly been proposed as a solution to the molecular design problem. However, it has proved challenging to control the design process or incorporate prior knowledge, limiting their practical use in drug discovery. In particular, generative methods have made limited use of three-dimensional (3D) structural information even though this is critical to binding. This work describes a method to incorporate such information and demonstrates the benefit of doing so. We combine an existing graph-based deep generative model, DeLinker, with a convolutional neural network to utilise physically-meaningful 3D representations of molecules and target pharmacophores. We apply our model, DEVELOP, to both linker and R-group design, demonstrating its suitability for both hit-to-lead and lead optimisation. The 3D pharmacophoric information results in improved generation and allows greater control of the design process. In multiple large-scale evaluations, we show that including 3D pharmacophoric constraints results in substantial improvements in the quality of generated molecules. On a challenging test set derived from PDBbind, our model improves the proportion of generated molecules with high 3D similarity to the original molecule by over 300%. In addition, DEVELOP recovers 10× more of the original molecules compared to the baseline DeLinker method. Our approach is general-purpose, readily modifiable to alternate 3D representations, and can be incorporated into other generative frameworks. Code is available at https://github.com/oxpig/DEVELOP.

Guided rational design with scaffold hopping leading to novel histamine H3 receptor ligands

During the past decades, histamine H₃ receptors have received widespread attention in pharmaceutical research due to their involvement in pathophysiology of several diseases such as neurodegenerative disorders. In this context, blocking of these receptors is of paramount importance in progression of such diseases. In the current investigation, novel histamine H₃ receptor ligands were designed by exploiting scaffold-hopping drug-design strategy. We inspected the designed molecules in terms of ADME properties, drug-likeness, as well as toxicity profiles. Additionally molecular docking and dynamics simulation studies were performed to predict binding mode and binding free energy calculations, respectively. Among the designed structures, we selected compound d2 and its demethylated derivative as examples for synthesis and affinity measurement. In vitro binding assays of the synthesized molecules demonstrated that d2 has lower binding affinity (K_i = 2.61 μM) in radioligand displacement assay to hH₃R than that of demethylated form (K_i = 12.53 μM). The newly designed compounds avoid of any toxicity predictors resulted from extended in silico and experimental studies, can offer another scaffold for histamine H₃R antagonists for further structure-activity relationship studies.

Novel Pyridine Bioisostere of Cabozantinib as a Potent c-Met Kinase Inhibitor: Synthesis and Anti-Tumor Activity against Hepatocellular Carcinoma

Two novel bioisosteres of cabozantinib, 3 and 4, were designed and synthesized. The benzene ring in the center of the cabozantinib structure was replaced by trimethylpyridine (3) and pyridine (4), respectively. Surprisingly, the two compounds showed extremely contrasting mesenchymal-epithelial transition factor (c-Met) inhibitory activities at 1 μM concentration (4% inhibition of 3 vs. 94% inhibition of 4). The IC₅₀ value of compound 4 was 4.9 nM, similar to that of cabozantinib (5.4 nM). A ligand-based docking study suggested that 4 includes the preferred conformation for the binding to c-Met in the conformational ensemble, but 3 does not. The anti-proliferative activity of compound 4 against hepatocellular carcinoma (Hep3B and Huh7) and non-small-cell lung cancer (A549 and H1299) cell lines was better than that of cabozantinib, whereas 3 did not show a significant anti-proliferative activity. Moreover, the tumor selectivity of compound 4 toward hepatocellular carcinoma cell lines was higher than that of cabozantinib. In the xenograft chick tumor model, compound 4 inhibited Hep3B tumor growth to a much greater extent than cabozantinib. The present study suggests that compound 4 may be a good therapeutic candidate against hepatocellular carcinoma.