Discovery of Novel Inhibitors Targeting Human O-GlcNAcase: Docking-Based Virtual Screening, Biological Evaluation, Structural Modification, and Molecular Dynamics Simulation

Insights into the computer-aided drug design and discovery based on anthraquinone scaffold for cancer treatment: A systematic review

BACKGROUND

In the search for better anticancer drugs, computer-aided drug design (CADD) techniques play an indispensable role in facilitating the lengthy and costly drug discovery process especially when natural products are involved. Anthraquinone is one of the most widely-recognized natural products with anticancer properties. This review aimed to systematically assess and synthesize evidence on the utilization of CADD techniques centered on the anthraquinone scaffold for cancer treatment.

METHODS

The conduct and reporting of this review were done in accordance to the Preferred Reporting Items for Systematic Reviews and Meta-analysis (PRISMA) 2020 guideline. The protocol was registered in the "International prospective register of systematic reviews" database (PROSPERO: CRD42023432904) and also published recently. The search strategy was designed based on the combination of concept 1 "CADD or virtual screening", concept 2 "anthraquinone" and concept 3 "cancer". The search was executed in PubMed, Scopus, Web of Science and MedRxiv on 30 June 2023.

RESULTS

Databases searching retrieved a total of 317 records. After deduplication and applying the eligibility criteria, the final review ended up with 32 articles in which 3 articles were found by citation searching. The CADD methods used in the studies were either structure-based alone (69%) or combined with ligand-based methods via parallel (9%) or sequential (22%) approaches. Molecular docking was performed in all studies, with Glide and AutoDock being the most popular commercial and public software used respectively. Protein data bank was used in most studies to retrieve the crystal structure of the targets of interest while the main ligand databases were PubChem and Zinc. The utilization of in-silico techniques has enabled a deeper dive into the structural, biological and pharmacological properties of anthraquinone derivatives, revealing their remarkable anticancer properties in an all-rounded fashion.

CONCLUSION

By harnessing the power of computational tools and leveraging the natural diversity of anthraquinone compounds, researchers can expedite the development of better drugs to address the unmet medical needs in cancer treatment by improving the treatment outcome for cancer patients.

Biochemical and structural characterization of chlorhexidine as an ATP-assisted inhibitor against type 1 methionyl-tRNA synthetase from Gram-positive bacteria

Methionyl-tRNA synthetase (MetRS) catalyzes the attachment of l-methionine (l-Met) to tRNAMet to generate methionyl-tRNAMet, an essential substrate for protein translation within ribosome. Owing to its indispensable biological function and the structural discrepancies with human counterpart, bacterial MetRS is considered an ideal target for developing antibacterials. Herein, chlorhexidine (CHX) was identified as a potent binder of Staphylococcus aureus MetRS (SaMetRS) through an ATP-aided affinity screening. The co-crystal structure showed that CHX simultaneously occupies the enlarged l-Met pocket (EMP) and the auxiliary pocket (AP) of SaMetRS with its two chlorophenyl groups, while its central hexyl linker swings upwards to interact with some conserved hydrophobic residues. ATP adopts alternative conformations in the active site cavity, and forms ionic bonds and water-mediated hydrogen bonds with CHX. Consistent with this synergistic binding mode, ATP concentration-dependently enhanced the binding affinity of CHX to SaMetRS from 10.2 μM (no ATP) to 0.45 μM (1 mM ATP). While it selectively inhibited two representative type 1 MetRSs from S. aureus and Enterococcus faecalis, CHX did not show significant interactions with three tested type 2 MetRSs, including human cytoplasmic MetRS, in the enzyme inhibition and biophysical binding assays, probably due to the conformational differences between two types of MetRSs at their EMP and AP. Our findings on CHX may inspire the design of MetRS-directed antimicrobials in future.

Identification of Triazolo-quinazolinone Derivatives as Novel and Potent Chitinase OfChi-h Inhibitors Based on Structure-Based Virtual Screening

Insect chitinase, OfChi-h, from Ostrinia furnacalis, is considered as a promising target for the development of green pesticides. On the basis of the crystal structure of OfChi-h, we successfully obtained a triazolo-quinazolinone scaffold as the novel class of OfChi-h inhibitor via a structure-based virtual screening approach. Rational compound screening enabled us to acquire a potent OfChi-h inhibitor TQ19 with a Ki value of 0.33 μM. Furthermore, the in vivo biological activity of target compounds was assayed. The results showed that compounds TQ8 and TQ19 could dramatically inhibit the growth and development of Ostrinia nubilalis larvae, and most of the compounds showed higher insecticidal activity than hexaflumuron. This present work reveals that triazolo-quinazolinone derivatives can serve as novel candidates for insect growth regulators.

Identification of novel nematode succinate dehydrogenase inhibitors: Virtual screening based on ligand-pocket interactions

To discover new nematicidal succinate dehydrogenase (SDH) inhibitors with novel structures, we conducted a virtual screening of the ChemBridge library with 1.7 million compounds based on ligand-pocket interactions. The homology model of Caenorhabditis elegans SDH was established, along with a pharmacophore model based on ligand-pocket interactions. After the pharmacophore-based and docking-based screening, 19 compounds were selected for the subsequent enzymatic assays. The results showed that compound 1 (ID: 7607321) exhibited inhibitory activity against SDH with a determined IC₅₀ value of 19.6 μM. Structural modifications and nematicidal activity studies were then carried out, which provided further evidence that compound 1 exhibited excellent nematicidal activity. Molecular dynamics simulations were then conducted to investigate the underlying molecular basis for the potency of these inhibitors against SDH. This work provides a reliable strategy and useful information for the future design of nematode SDH inhibitors.

Discovery of novel inhibitors targeting nematode chitinase CeCht1: Virtual screening, biological evaluation, and molecular dynamics simulation

Plant-parasitic nematodes are a main limiting factor for worldwide agriculture. To reduce the global burden of nematode infections, chemical nematicides are still the most effective methods to manage nematodes. With the increasing resistance of nematodes, the development of new anti-nematicides drug is urgent. Nematode chitinases are found to play important roles in various physiological functions, such as larva moulting, hatching from eggshell, and host infection. Inhibition of nematode chitinase is considered a promising strategy for the development of eco-friendly nematicides. In this study, to develop novel nematode chitinase CeCht1 inhibitors, virtual screening of the ZINC database was performed using the pesticide-likeness rules, pharmacophore-based and docking-based approach in turn. Compounds HAU-4 and HAU-7 were identified as potent CeCht1 inhibitors with the IC50 values of 4.2 μM and 10.0 μM, respectively. Moreover, molecular dynamics simulations combined with binding free energy and free energy decomposition calculations were conducted to investigate the basis for the potency of the two inhibitors toward CeCht1. This work gives an insight into the future rational development of novel and potent nematode chitinase inhibitors.

Discovery of Azo-Aminopyrimidines as Novel and Potent Chitinase OfChi-h Inhibitors via Structure-Based Virtual Screening and Rational Lead Optimization

Chitinase OfChi-h, from the destructive agricultural pest Ostrinia furnacalis, is considered as a promising target for green pest control and management. In this study, structure-based virtual screening and rational molecular optimization led to the synthesis of a series of azo-aminopyrimidine derivatives as a novel class of OfChi-h inhibitors. Among them, the most potent compound 8f, with a benzyl on the amino group at the 4-position of pyrimidine, exhibited a K_i value of 64.7 nM against OfChi-h. In addition, molecular docking studies were carried out to investigate the basis for the potency of the aminopyrimidines against OfChi-h. Furthermore, the insecticidal activity of the target compounds against Plutella xylostella and Ostrinia nubilalis was assessed, and the potent OfChi-h inhibitors 8f and 8i showed higher insecticidal activity than the control pesticide hexaflumuron. The present work revealed that the azo-aminopyrimidine skeletons characterized by concise chemical structure and high efficiency could be further developed as potential pesticides for the control of lepidopteran pests.

Rational Design, Synthesis, and Biological Investigations of N-Methylcarbamoylguanidinyl Azamacrolides as a Novel Chitinase Inhibitor

Chitinase is one of the most important glycoside hydrolyases, widely existing in bacteria, fungi, insects, and plants. It is involved in fungal cell wall remodeling and insect molting. Chitinase inhibitors are an effective means of controlling pathogens and pests. Natural product argifin is a 17-membered pentapeptide that exhibits efficient chitinase inhibitory activity. However, the complexity of the synthetic process results in a lot of restrictions for wide range of applications. In this work, we designed a series of azamacrolide chitinase inhibitors based on the structural features of argifin that have high inhibitory activities against bacterial and insectile chitinase. The most potent chitinase inhibitor compound 19c exhibited IC₅₀ values of 56 nM and 110 nM against OfChi-h and SmChiB, respectively. The molecular docking and molecular dynamics simulations revealed that all inhibitors were bound to the -1 subsite of chitinases via N-methylcarbamoylguanidinyl as well as argifin. Finally, a bioactivity assay against pests was carried out. Compound 18a showed 80% mortality for Mythimna separata at a concentration of 50 mg/L. Besides, insecticides 19b and 19c exhibited high mortality against Plutella xylostella (76 and 73% mortalities at 50 mg/L, respectively).

Exploration and augmentation of pharmacological space via adversarial auto-encoder model for facilitating kinase-centric drug development

Predicting compound-protein interactions (CPIs) is of great importance for drug discovery and repositioning, yet still challenging mainly due to the sparse nature of CPI matrixes, resulting in poor generalization performance. Hence, unlike typical CPI prediction models focused on representation learning or model selection, we propose a deep neural network-based strategy, PCM-AAE, that re-explores and augments the pharmacological space of kinase inhibitors by introducing the adversarial auto-encoder model (AAE) to improve the generalization of the prediction model. To complete the data space, we constructed Ensemble of PCM-AAE (EPA), an ensemble model that quickly and accurately yields quantitative predictions of binding affinity between any human kinase and inhibitor. In rigorous internal validation, EPA showed excellent performance, consistently outperforming the model trained with the imbalanced set, especially for targets with relatively fewer training data points. Improved prediction accuracy of EPA for external datasets enhances its generalization ability, making it possible to gracefully handle previously unseen kinases and inhibitors. EPA showed promising potential when directly applied to virtual screening and off-target prediction, exhibiting its practicality in hit prediction. Our strategy is expected to facilitate kinase-centric drug development, as well as to solve more challenging prediction problems with insufficient data points.

Synthesis, Anticancer Activity, Structure-Activity Relationship, and Molecular Modeling Studies of α-Mangostin Derivatives as hERα Inhibitor

α-Mangostin is one of the secondary metabolites in mangosteen pericarp, which has been reported to have anti-breast cancer activity. In our previous study, three α-mangostin derivatives were computationally designed as hERα antagonists. In this present study, the designed compounds were synthesized undergoing a benzoylation reaction between α-mangostin with three benzoyl chloride derivatives to produce three derivatives, namely, AMB-1, AMB-2, and AMB-10. The synthesized compounds were then evaluated for their antiproliferative activity against the MCF-7 breast cancer cell model with hERα as the protein target. The in vitro assay shows moderate activity (57-126 μM) for all derivatives. The dynamic behaviors of all ligands, including α-mangostin and 4-hydroxytamoxifen (4-OHT), were studied with 100 ns of MD simulation. The structure-activity relationship shows that although it does not entirely concord with the expected design, it can explain the trend of α-mangostin and its derivatives antiproliferative activities against MCF-7, which associates with hERα antagonism.

Molecular Mechanism Exploration of Potent Fluorinated PI3K Inhibitors with a Triazine Scaffold: Unveiling the Unusual Synergistic Effect of Pyridine-to-Pyrimidine Ring Interconversion and CF3 Defluorination

The phosphatidylinostitol-3-kinase (PI3K)/AKT/mammalian target of rapamycin signaling pathway is a vital regulator of cell proliferation, growth, and survival, which is frequently overactivated in many human cancers. To this effect, PI3K, which is an important mediator of this pathway, has been pinpointed as a crucial target in cancer therapy and hence the importance of PI3K inhibitors. It was recently reported that defluorination and pyridine-to-pyrimidine ring interconversion increase the potency of specific small-molecule inhibitors of PI3K. Compound 4, an inhibitor with the difluorinated pyrimidine motif, was found to be eight times more potent against PI3K than compound 1, an inhibitor with the trifluorinated pyridine motif. This observation presents the need to rationally resolve the differential inhibitory mechanisms exhibited by both compounds. In this present work, we employed multiple computational approaches to investigate and distinguish the binding modes of 1 and 4 in addition to the effects they mediate on the secondary structure of PI3K. Likewise, we evaluated two other derivatives, compounds 2 with the difluorinated pyridine motif and 3 with the trifluorinated pyrimidine motif, to investigate the cooperativity effect between the defluorination of CF₃ and pyridine-to-pyrimidine ring interconversion. Findings revealed that PI3K, upon interaction with 4, exhibited a series of structural changes that favored the binding of the inhibitor at the active-site region. Furthermore, a positive (synergistic) cooperativity effect was observed between CF₃ defluorination and pyridine-to-pyrimidine ring interconversion. Moreover, there was a good correlation between the binding free energy estimated and the biological activity reported experimentally. Energy decomposition analysis revealed that the major contributing force to binding affinity variations between 1 and 4 is the electrostatic energy. Per-residue energy-based hierarchical clustering analysis further identified four hot-spot residues ASP841, TYR867, ASP964, and LYS833 and four warm-spot residues ASP836, SER806, ASP837, and LYS808, which essentially mediate the optimal and higher-affinity binding of compound 4 to PI3K relative to 1. This study therefore provides rational insights into the mechanisms by which 4 exhibited superior PI3K-inhibitory activities over 1, which is vital for future structure-based drug discovery efforts in PI3K targeting.

Preliminary Insight of Pyrrolylated-Chalcones as New Anti-Methicillin-Resistant Staphylococcus aureus (Anti-MRSA) Agents

Bacterial infections are regarded as one of the leading causes of fatal morbidity and death in patients infected with diseases. The ability of microorganisms, particularly methicillin-resistant Staphylococcus aureus (MRSA), to develop resistance to current drugs has evoked the need for a continuous search for new drugs with better efficacies. Hence, a series of non-PAINS associated pyrrolylated-chalcones (1–15) were synthesized and evaluated for their potency against MRSA. The hydroxyl-containing compounds (8, 9, and 10) showed the most significant anti-MRSA efficiency, with the MIC and MBC values ranging from 0.08 to 0.70 mg/mL and 0.16 to 1.88 mg/mL, respectively. The time-kill curve and SEM analyses exhibited bacterial cell death within four hours after exposure to 9, suggesting its bactericidal properties. Furthermore, the docking simulation between 9 and penicillin-binding protein 2a (PBP2a, PDB ID: 6Q9N) suggests a relatively similar bonding interaction to the standard drug with a binding affinity score of −7.0 kcal/mol. Moreover, the zebrafish model showed no toxic effects in the normal embryonic development, blood vessel formation, and apoptosis when exposed to up to 40 µM of compound 9. The overall results suggest that the pyrrolylated-chalcones may be considered as a potential inhibitor in the design of new anti-MRSA agents.

Analyses of structural dynamics revealed flexible binding mechanism for the Agrilus mali odorant binding protein 8 towards plant volatiles

BACKGROUND

Volatiles from host plants are an important source of insect pest attractants and repellents. Insect odorant binding proteins (OBPs) have been widely characterized, but the molecular binding dynamics and underlying mechanisms are still not well understood. Thus, we characterized binding characteristics of AmalOBP8 from the apple buprestid beetle (Agrilus mali Matsumura), an unprecedented serious threat to rare apple germplasm resources and local ecosystems.

RESULTS

Fluorescence studies demonstrated that the quenching mechanism was clearly static. AmalOBP8 was found to bind with both volatiles at single independent sites. Negative thermodynamic parameters suggested that binding interactions between AmalOBP8 and both volatiles could occur spontaneously. Hydrogen bonding was the key force in AmalOBP8's binding to geranyl formate, for which the amino acid residue Trp106 played a critical role in the binding pocket. Multiple Leu residues in AmalOBP8 created a strong hydrophobic environment, and formed the binding pocket for (Z)-3-hexenyl hexanoate. Compared to classic OBPs, in addition to lack of one disulfide bridge, AmalOBP8 had a small α-helix (α7) at the C-terminus, resulting in greater flexibility and adaptability for this protein to bind with different compound molecules.

CONCLUSION

Key residues of AmalOBP8 in binding interactions with plant volatiles were clarified. AmalOPB8 had a large ligand binding spectrum and great flexibility in binding with plant volatiles, providing good molecular targets for screening insect attractants and repellents. Our results can promote understanding of insects' perception of various odorants, and establish a foundation for discovery of new pest control agents. © 2020 Society of Chemical Industry.

Region- and voxel-based quantification in human brain of [18F]LSN3316612, a radioligand for O-GlcNAcase

Background

Previous studies found that the positron emission tomography (PET) radioligand [¹⁸F]LSN3316612 accurately quantified O-GlcNAcase in human brain using a two-tissue compartment model (2TCM). This study sought to assess kinetic model(s) as an alternative to 2TCM for quantifying [¹⁸F]LSN3316612 binding, particularly in order to generate good-quality parametric images.

Methods

The current study reanalyzed data from a previous study of 10 healthy volunteers who underwent both test and retest PET scans with [¹⁸F]LSN3316612. Kinetic analysis was performed at the region level with 2TCM using 120-min PET data and arterial input function, which was considered as the gold standard. Quantification was then obtained at both the region and voxel levels using Logan plot, Ichise's multilinear analysis-1 (MA1), standard spectral analysis (SA), and impulse response function at 120 min (IRF₁₂₀). To avoid arterial sampling, a noninvasive relative quantification (standardized uptake value ratio (SUVR)) was also tested using the corpus callosum as a pseudo-reference region. Venous samples were also assessed to see whether they could substitute for arterial ones.

Results

Logan and MA1 generated parametric images of good visual quality and their total distribution volume (V_T) values at both the region and voxel levels were strongly correlated with 2TCM-derived V_T (r = 0.96–0.99) and showed little bias (up to − 8%). SA was more weakly correlated to 2TCM-derived V_T (r = 0.93–0.98) and was more biased (~ 16%). IRF₁₂₀ showed a strong correlation with 2TCM-derived V_T (r = 0.96) but generated noisier parametric images. All techniques were comparable to 2TCM in terms of test–retest variability and reliability except IRF₁₂₀, which gave significantly worse results. Noninvasive SUVR values were not correlated with 2TCM-derived V_T, and arteriovenous equilibrium was never reached.

Conclusions

Compared to SA and IRF, Logan and MA1 are more suitable alternatives to 2TCM for quantifying [¹⁸F]LSN3316612 and generating good-quality parametric images.

Discovery of a Novel and Brain-Penetrant O-GlcNAcase Inhibitor via Virtual Screening, Structure-Based Analysis, and Rational Lead Optimization

O-GlcNAcase (OGA) has received increasing attention as an attractive therapeutic target for tau-mediated neurodegenerative disorders; however, its role in these pathologies remains unclear. Therefore, potent chemical tools with favorable pharmacokinetic profiles are desirable to characterize this enzyme. Herein, we report the discovery of a potent and novel OGA inhibitor, compound 5i, comprising an aminopyrimidine scaffold, identified by virtual screening based on multiple methodologies combining structure-based and ligand-based approaches, followed by sequential optimization with a focus on ligand lipophilicity efficiency. This compound was observed to increase the level of O-GlcNAcylated protein in cells and display suitable pharmacokinetic properties and brain permeability. Crystallographic analysis revealed that the chemical series bind to OGA via characteristic hydrophobic interactions, which resulted in a high affinity for OGA with moderate lipophilicity. Compound 5i could serve as a useful chemical probe to help establish a proof-of-concept of OGA inhibition as a therapeutic target for the treatment of tauopathies.