Next generation 3D pharmacophore modeling

Simulation- and AI-directed optimization of 4,6-substituted 1,3,5-triazin-2(1H)-ones as inhibitors of human DNA topoisomerase IIα

The 4,6-substituted-1,3,5-triazin-2(1H)-ones are promising inhibitors of human DNA topoisomerase IIα. To further develop this chemical class targeting the enzyme´s ATP binding site, the triazin-2(1H)-one substitution position 6 was optimized. Inspired by binding of preclinical substituted 9H-purine derivative, bicyclic substituents were incorporated at position 6 and the utility of this modification was validated by a combination of molecular simulations, dynamic pharmacophores, and free energy calculations. Considering also predictions of Deepfrag, a software developed for structure-based lead optimization based on deep learning, compounds with both bicyclic and monocyclic substitutions were synthesized and investigated for their inhibitory activity. The SAR data showed that the bicyclic substituted compounds exhibited good inhibition of topo IIα, comparable to their mono-substituted counterparts. Further evaluation on a panel of human protein kinases showed selectivity for the inhibition of topo IIα. Mechanistic studies indicated that the compounds acted predominantly as catalytic inhibitors, with some exhibiting topo IIα poison effects at higher concentrations. Integration of STD NMR experiments and molecular simulations, provided insights into the binding model and highlighted the importance of the Asn120 interaction and hydrophobic interactions with substituents at positions 4 and 6. In addition, NCI-60 screening demonstrated cytotoxicity of the compounds with bicyclic substituents and identified sensitive human cancer cell lines, underlining the translational relevance of our findings for further preclinical development of this class of compounds. The study highlights the synergy between simulation and AI-based approaches in efficiently guiding molecular design for drug optimization, which has implications for further preclinical development of this class of compounds.

Structure-based pharmacophore modeling for precision inhibition of mutant ESR2 in breast cancer: A systematic computational approach

BACKGROUND

Breast cancer, a leading cause of female mortality, is closely linked to mutations in estrogen receptor beta (ESR2), particularly in the ligand-binding domain, which contributed to altered signaling pathways and uncontrolled cell growth.

OBJECTIVES/AIMS

This study investigates the molecular and structural aspects of ESR2 mutant proteins to identify shared pharmacophoric regions of ESR2 mutant proteins and potential therapeutic targets aligned within the pharmacophore model.

METHODS

This study was initiated by establishing a common pharmacophore model among three mutant ESR2 proteins (PDB ID: 2FSZ, 7XVZ, and 7XWR). The generated shared feature pharmacophore (SFP) includes four primary binding interactions: Hydrogen bond donors (HBD), hydrogen bond acceptors (HBA), hydrophobic interactions (HPho), and Aromatic interactions (Ar), along with halogen bond donors (XBD) and totalling 11 features (HBD: 2, HBA: 3, HPho: 3, Ar: 2, XBD: 1). By employing an in-house Python script, these 11 features distributed into 336 combinations, which were used as query to isolate a drug library of 41,248 compounds and subjected to virtual screening through the generated SFP.

RESULTS

The virtual screening demonstrated 33 hits showing potential pharmacophoric fit scores and low RMSD value. The top four compounds: ZINC94272748, ZINC79046938, ZINC05925939, and ZINC59928516 showed a fit score of more than 86% and satisfied the Lipinski rule of five. These four compounds and a control underwent molecular (XP Glide mode) docking analysis against wild-type ESR2 protein (PDB ID: 1QKM), resulting in binding affinity of -8.26, -5.73, -10.80, and -8.42 kcal/mol, respectively, along with the control -7.2 kcal/mol. Furthermore, the stability of the selected candidates was determined through molecular dynamics (MD) simulations of 200 ns and MM-GBSA analysis.

CONCLUSION

Based on MD simulations and MM-GBSA analysis, our study identified ZINC05925939 as a promising ESR2 inhibitor among the top four hits. However, it is essential to conduct further wet lab evaluation to assess its efficacy.

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

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

A Systematic Hierarchical Virtual Screening Model for RhlR Inhibitors Based on PCA, Pharmacophore, Docking, and Molecular Dynamics

RhlR plays a key role in the quorum sensing of Pseudomonas aeruginosa. The current structure-activity relationship (SAR) studies of RhlR inhibitors mainly focus on elucidating the functional groups. Based on a systematic review of previous research on RhlR inhibitors, this study aims to establish a systematic, hierarchical screening model for RhlR inhibitors. We initially established a database and utilized principal component analysis (PCA) to categorize the inhibitors into two classes. Based on the training set, pharmacophore models were established to elucidate the structural characteristics of ligands. Subsequently, molecular docking, molecular dynamics simulations, and the calculation of binding free energy and strain energy were performed to validate the crucial interactions between ligands and receptors. Then, the screening criteria for RhlR inhibitors were established hierarchically based on ligand structure characteristics, ligand-receptor interaction, and receptor affinity. Test sets were finally employed to validate the hierarchical virtual screening model by comparing it with the current SAR studies of RhlR inhibitors. The hierarchical screening model was confirmed to possess higher accuracy and a true positive rate, which holds promise for subsequent screening and the discovery of active RhlR inhibitors.

bbSelect - An Open-Source Tool for Performing a 3D Pharmacophore-Driven Diverse Selection of R-Groups

The design of compounds during hit-to-lead often seeks to explore a vector from a core scaffold to form additional interactions with the target protein. A rational approach to this is to probe the region of a protein accessed by a vector with a systematic placement of pharmacophore features in 3D, particularly when bound structures are not available. Herein, we present bbSelect, an open-source tool built to map the placements of pharmacophore features in 3D Euclidean space from a library of R-groups, employing partitioning to drive a diverse and systematic selection to a user-defined size. An evaluation of bbSelect against established methods exemplified the superiority of bbSelect in its ability to perform diverse selections, achieving high levels of pharmacophore feature placement coverage with selection sizes of a fraction of the total set and without the introduction of excess complexity. bbSelect also reports visualizations and rationale to enable users to understand and interrogate results. This provides a tool for the drug discovery community to guide their hit-to-lead activities.

FaissMolLib: An efficient and easy deployable tool for ligand-based virtual screening

Virtual screening-based molecular similarity and fingerprint are crucial in drug design, target prediction, and ADMET prediction, aiding in identifying potential hits and optimizing lead compounds. However, challenges such as lack of comprehensive open-source molecular fingerprint databases and efficient search methods for virtual screening are prevalent. To address these issues, we introduce FaissMolLib, an open-source virtual screening tool that integrates 2.8 million compounds from ChEMBL and ZINC databases. Notably, FaissMolLib employs the highly efficient Faiss search algorithm, outperforming the Tanimoto algorithm in identifying similar molecules with its tighter clustering in scatter plots and lower mean, standard deviation, and variance in key molecular properties. This feature enables FaissMolLib to screen 2.8 million compounds in just 0.05 seconds, offering researchers an efficient, easily deployable solution for virtual screening on laptops and building unique compound databases. This significant advancement holds great potential for accelerating drug discovery efforts and enhancing chemical data analysis. FaissMolLib is freely available at http://liuhaihan.gnway.cc:80. The code and dataset of FaissMolLib are freely available at https://github.com/Superhaihan/FiassMolLib.

Nature-inspired substituted 3-(imidazol-2-yl) morpholines targeting human topoisomerase IIα: Dynophore-derived discovery

The molecular nanomachine, human DNA topoisomerase IIα, plays a crucial role in replication, transcription, and recombination by catalyzing topological changes in the DNA, rendering it an optimal target for cancer chemotherapy. Current clinical topoisomerase II poisons often cause secondary tumors as side effects due to the accumulation of double-strand breaks in the DNA, spurring the development of catalytic inhibitors. Here, we used a dynamic pharmacophore approach to develop catalytic inhibitors targeting the ATP binding site of human DNA topoisomerase IIα. Our screening of a library of nature-inspired compounds led to the discovery of a class of 3-(imidazol-2-yl) morpholines as potent catalytic inhibitors that bind to the ATPase domain. Further experimental and computational studies identified hit compound 17, which exhibited selectivity against the human DNA topoisomerase IIα versus human protein kinases, cytotoxicity against several human cancer cells, and did not induce DNA double-strand breaks, making it distinct from clinical topoisomerase II poisons. This study integrates an innovative natural product-inspired chemistry and successful implementation of a molecular design strategy that incorporates a dynamic component of ligand-target molecular recognition, with comprehensive experimental characterization leading to hit compounds with potential impact on the development of more efficient chemotherapies.

PharmaCore: The Automatic Generation of 3D Structure-Based Pharmacophore Models from Protein/Ligand Complexes

In this work, we present PharmaCore: a new, completely automatic workflow aimed at generating three-dimensional (3D) structure-based pharmacophore models toward any target of interest. The proposed approach relies on using cocrystallized ligands to create the input files for generating the pharmacophore hypotheses, integrating not only the three-dimensional structural information on the ligand but also data concerning the binding mode of these molecules put in the protein cavity. We developed a Python library that, starting from the specific UniProt ID of the protein under investigation as the only element that requires user intervention, subsequently collects and aligns the corresponding structures bearing a known ligand in a fully automated fashion, bringing them all into the same coordinate system. The protocol includes a final phase in which the aligned small molecules are used to produce the pharmacophore hypotheses directly onto the protein structure using a specific software, e.g., Phase (Schrödinger LLC). To validate the entire procedure and highlight the possible applications in the field of drug discovery and repositioning, we first generated pharmacophores for soluble epoxide hydrolase (sEH) and compared with already-published ones. Then, we reproduced the binding profile of a reported selective binder of ATAD2 bromodomain (AM879), testing it against a panel of 1741 pharmacophores related to 16 epigenetic proteins and automatically generated with PharmaCore, finally disclosing putative unprecedented off-targets. The computational predictions were successfully validated with AlphaScreen assays, highlighting the applicability of the proposed workflow in drug discovery and repositioning. Finally, the process was also validated on tankyrase 2 and SARS-CoV-2 MPro, confirming the robustness of PharmaCore.

Identification of Druggable Binding Sites and Small Molecules as Modulators of TMC1

Our ability to hear and maintain balance relies on the proper functioning of inner ear sensory hair cells, which translate mechanical stimuli into electrical signals via mechano-electrical transducer (MET) channels, composed of TMC1/2 proteins. However, the therapeutic use of ototoxic drugs, such as aminoglycosides and cisplatin, which can enter hair cells through MET channels, often leads to profound auditory and vestibular dysfunction. Despite extensive research on otoprotective compounds targeting MET channels, our understanding of how small molecule modulators interact with these channels remains limited, hampering the discovery of novel compounds. Here, we propose a structure-based screening approach, integrating 3D-pharmacophore modeling, molecular simulations, and experimental validation. Our pipeline successfully identified several novel compounds and FDA-approved drugs that reduced dye uptake in cultured cochlear explants, indicating MET modulation activity. Molecular docking and free-energy estimations for binding allowed us to identify three potential drug binding sites within the channel pore, phospholipids, and key amino acids involved in modulator interactions. We also identified shared ligand-binding features between TMC and structurally related TMEM16 protein families, providing novel insights into their distinct inhibition, while potentially guiding the rational design of MET-channel-specific modulators. Our pipeline offers a broad application to discover small molecule modulators for a wide spectrum of mechanosensitive ion channels.

Targeted small molecule inhibitors blocking the cytolytic effects of pneumolysin and homologous toxins

Pneumolysin (PLY) is a cholesterol-dependent cytolysin (CDC) from Streptococcus pneumoniae, the main cause for bacterial pneumonia. Liberation of PLY during infection leads to compromised immune system and cytolytic cell death. Here, we report discovery, development, and validation of targeted small molecule inhibitors of PLY (pore-blockers, PB). PB-1 is a virtual screening hit inhibiting PLY-mediated hemolysis. Structural optimization provides PB-2 with improved efficacy. Cryo-electron tomography reveals that PB-2 blocks PLY-binding to cholesterol-containing membranes and subsequent pore formation. Scaffold-hopping delivers PB-3 with superior chemical stability and solubility. PB-3, formed in a protein-templated reaction, binds to Cys428 adjacent to the cholesterol recognition domain of PLY with a KD of 256 nM and a residence time of 2000 s. It acts as anti-virulence factor preventing human lung epithelial cells from PLY-mediated cytolysis and cell death during infection with Streptococcus pneumoniae and is active against the homologous Cys-containing CDC perfringolysin (PFO) as well.

Consensus Pharmacophore Strategy For Identifying Novel SARS-Cov-2 Mpro Inhibitors from Large Chemical Libraries

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main Protease (Mpro) is an enzyme that cleaves viral polyproteins translated from the viral genome and is critical for viral replication. Mpro is a target for anti-SARS-CoV-2 drug development, and multiple Mpro crystals complexed with competitive inhibitors have been reported. In this study, we aimed to develop an Mpro consensus pharmacophore as a tool to expand the search for inhibitors. We generated a consensus model by aligning and summarizing pharmacophoric points from 152 bioactive conformers of SARS-CoV-2 Mpro inhibitors. Validation against a library of conformers from a subset of ligands showed that our model retrieved poses that reproduced the crystal-binding mode in 77% of the cases. Using models derived from a consensus pharmacophore, we screened >340 million compounds. Pharmacophore-matching and chemoinformatics analyses identified new potential Mpro inhibitors. The candidate compounds were chemically dissimilar to the reference set, and among them, demonstrating the relevance of our model. We evaluated the effect of 16 candidates on Mpro enzymatic activity finding that seven have inhibitory activity. Three compounds (1, 4, and 5) had IC50 values in the midmicromolar range. The Mpro consensus pharmacophore reported herein can be used to identify compounds with improved activity and novel chemical scaffolds against Mpro. The method developed for its generation is provided as an open-access code (https://github.com/AngelRuizMoreno/ConcensusPharmacophore) and can be applied to other pharmacological targets.

Screening Techniques for Drug Discovery in Alzheimer's Disease

Alzheimer's disease (AD) is a neurodegenerative disorder characterized by progressive and irreversible impairment of memory and other cognitive functions of the aging brain. Pathways such as amyloid beta neurotoxicity, tau pathogenesis and neuroinflammatory have been used to understand AD, despite not knowing the definite molecular mechanism which causes this progressive disease. This review attempts to summarize the small molecules that target these pathways using various techniques involving high-throughput screening, molecular modeling, custom bioassays, and spectroscopic detection tools. Novel and evolving screening methods developed to advance drug discovery initiatives in AD research are also highlighted.

3D-QSAR pharmacophore modeling, virtual screening, molecular docking, MD simulations, in vitro and in vivo studies to identify potential anti-hyperplasia drugs

Psoriasis is a common immune-mediated skin condition characterized by aberrant keratinocytes and cell proliferation. The purpose of this study was to explore the FDA-approved drugs by 3D-QSAR pharmacophore model and evaluate their efficiency by in-silico, in vitro, and in vivo psoriasis animal model. A 3D-QSAR pharmacophore model was developed by utilizing HypoGen algorithm using the structural features of 48 diaryl derivatives with diverse molecular patterns. The model was validated by a test set of 27 compounds, by cost analysis method, and Fischer's randomization test. The correlation coefficient of the best model (Hypo2) was 0.9601 for the training set while it was 0.805 for the test set. The selected model was taken as a 3D query for the virtual screening of over 3000 FDA-approved drugs. Compounds mapped with the pharmacophore model were further screened through molecular docking. The hits that showed the best docking results were screened through in silico skin toxicity approach. Top five hits were selected for the MD simulation studies. Based on MD simulations results, the best two hit molecules, that is, ebastine (Ebs) and mebeverine (Mbv) were selected for in vitro and in vivo antioxidant studies performed in mice. TNF-α and COX pro-inflammatory mediators, biochemical assays, histopathological analyses, and immunohistochemistry observations confirmed the anti-inflammatory response of the selected drugs. Based on these findings, it appeared that Ebs can effectively treat psoriasis-like skin lesions and down-regulate inflammatory responses which was consistent with docking predictions and could potentially be employed for further research on inflammation-related skin illnesses such as psoriasis.

Identification of coumarin derivatives targeting acetylcholinesterase for Alzheimer's disease by field-based 3D-QSAR, pharmacophore model-based virtual screening, molecular docking, MM/GBSA, ADME and MD Simulation study

Alzheimer's disease (AD) leads to gradual memory loss including other compromised cognitive abilities. Acetylcholinesterase (AChE), an important biochemical enzyme from the cholinesterase (ChE) family, is recognized as primary pharmacological target for treating AD. Currently marketed drugs for AD treatment are primarily AChE inhibitors and coumarin derivatives comprising a wide variety of pharmacological activities have proved their efficacy towards AChE inhibition. Ensaculin (KA-672 HCl), a compound that belong to the coumarin family, is a clinical trial candidate for AD treatment. Therefore, a ligand library was prepared with 60 reported coumarin derivatives for field-based 3D-QSAR and pharmacophore modelling. The field-based 3D-QSAR model obtained at partial least square (PLS) factor 7, was the best validated model that predicted activity closer to original activity for each ligand introduced. The contour maps demonstrated spatial distribution of favourable and unfavorable steric, hydrophobic, electrostatic and H-bond donor and acceptor contours around coumarin nucleus. The best pharmacophore model, ADHRR_1 exhibited five essential pharmacophoric features of four different traits for optimum AChE inhibition. Virtual screening through ADHRR_1 accompanied with molecular docking and MM/GBSA identified 10 HITs from a 4,00,000 coumarin derivatives from PubChem database. HITs comprised docking scores ranging from -12.096 kcal/mol to -8.271 kcal/mol and compared with the reference drug Donepezil (-8.271 kcal/mol). ADME properties analysis led into detecting two leads (HIT 1 and HIT 2) among these 10 HITs. Molecular Dynamics Simulation indicated thermodynamic stability of the complex of lead compounds with AChE protein. Finally, thorough survey of the experimental results from 3D-QSAR modelling, pharmacophore modelling and molecular docking interactions led us to develop the lead formula I for future advancements in treating AD through AChE inhibitors.

Targeting ion channels with ultra-large library screening for hit discovery

Ion channels play a crucial role in a variety of physiological and pathological processes, making them attractive targets for drug development in diseases such as diabetes, epilepsy, hypertension, cancer, and chronic pain. Despite the importance of ion channels in drug discovery, the vastness of chemical space and the complexity of ion channels pose significant challenges for identifying drug candidates. The use of in silico methods in drug discovery has dramatically reduced the time and cost of drug development and has the potential to revolutionize the field of medicine. Recent advances in computer hardware and software have enabled the screening of ultra-large compound libraries. Integration of different methods at various scales and dimensions is becoming an inevitable trend in drug development. In this review, we provide an overview of current state-of-the-art computational chemistry methodologies for ultra-large compound library screening and their application to ion channel drug discovery research. We discuss the advantages and limitations of various in silico techniques, including virtual screening, molecular mechanics/dynamics simulations, and machine learning-based approaches. We also highlight several successful applications of computational chemistry methodologies in ion channel drug discovery and provide insights into future directions and challenges in this field.

In silico modeling revealed phytomolecules derived from Cymbopogon citratus (DC.) leaf extract as promising candidates for malaria therapy

The emergence of varying levels of resistance to currently available antimalarial drugs significantly threatens global health. This factor heightens the urgency to explore bioactive compounds from natural products with a view to discovering and developing newer antimalarial drugs with novel mode of actions. Therefore, we evaluated the inhibitory effects of sixteen phytocompounds from Cymbopogon citratus leaf extract against Plasmodium falciparum drug targets such as P. falciparum circumsporozoite protein (PfCSP), P. falciparum merozoite surface protein 1 (PfMSP1) and P. falciparum erythrocyte membrane protein 1 (PfEMP1). In silico approaches including molecular docking, pharmacophore modeling and 3D-QSAR were adopted to analyze the inhibitory activity of the compounds under consideration. The molecular docking results indicated that a compound swertiajaponin from C. citratus exhibited a higher binding affinity (-7.8 kcal/mol) to PfMSP1 as against the standard artesunate-amodiaquine (-6.6 kcal/mol). Swertiajaponin also formed strong hydrogen bond interactions with LYS29, CYS30, TYR34, ASN52, GLY55 and CYS28 amino acid residues. In addition, quercetin another compound from C. citratus exhibited significant binding energies -6.8 and -8.3 kcal/mol with PfCSP and PfEMP1, respectively but slightly lower than the standard artemether-lumefantrine with binding energies of -7.4 kcal/mol against PfCSP and -8.7 kcal/mol against PfEMP1. Overall, the present study provides evidence that swertiajaponin and other phytomolecules from C. citratus have modulatory properties toward P. falciparum drug targets and thus may warrant further exploration in early drug discovery efforts against malaria. Furthermore, these findings lend credence to the folkloric use of C. citratus for malaria treatment.Communicated by Ramaswamy H. Sarma.

Recent Deep Learning Applications to Structure-Based Drug Design

Identification and optimization of small molecules that bind to and modulate protein function is a crucial step in the early stages of drug development. For decades, this process has benefitted greatly from the use of computational models that can provide insights into molecular binding affinity and optimization. Over the past several years, various types of deep learning models have shown great potential in improving and enhancing the performance of traditional computational methods. In this chapter, we provide an overview of recent deep learning-based developments with applications in drug discovery. We classify these methods into four subcategories dependent on the task each method is aiming to solve. For each subcategory, we provide the general framework of the approach and discuss individual methods.

Design, synthesis, and pharmacological evaluation of indazole carboxamides of N-substituted pyrrole derivatives as soybean lipoxygenase inhibitors

In this paper, we attempted to develop a novel class of compounds against lipoxygenase, a key enzyme in the biosynthesis of leukotrienes implicated in a series of inflammatory diseases. Given the absence of appropriate human 5-lipoxygenase crystallographic data, solved soybean lipoxygenase-1 and -3 structures were used as a template to generate an accurate pharmacophore model which was further used for virtual screening purposes. Eight compounds (1-8) have been derived from the in-house library consisting of N-substituted pyrroles conjugated with 5- or 6-indazole moieties through a carboxamide linker. This study led to the discovery of hit molecule 8 bearing a naphthyl group with the IC50 value of 22 μM according to soybean lipoxygenase in vitro assay. Isosteric replacement of naphthyl ring with quinoline moieties and reduction of carbonyl carboxamide group resulted in compounds 9-12 and 13, respectively. Compound 12 demonstrated the most promising enzyme inhibition. In addition, compounds 8 and 12 were found to reduce the carrageenan-induced paw edema in vivo by 52.6 and 49.8%, respectively. In view of the encouraging outcomes concerning their notable in vitro and in vivo anti-inflammatory activities, compounds 8 and 12 could be further optimized for the discovery of novel 5-lipoxygenase inhibitors in future. A structure-based 3D pharmacophore model was used in the virtual screening of in-house library to discover novel potential 5-lipoxygenase inhibitors.

Applying atomistic neural networks to bias conformer ensembles towards bioactive-like conformations

Identifying bioactive conformations of small molecules is an essential process for virtual screening applications relying on three-dimensional structure such as molecular docking. For most small molecules, conformer generators retrieve at least one bioactive-like conformation, with an atomic root-mean-square deviation (ARMSD) lower than 1 Å, among the set of low-energy conformers generated. However, there is currently no general method to prioritise these likely target-bound conformations in the ensemble. In this work, we trained atomistic neural networks (AtNNs) on 3D information of generated conformers of a curated subset of PDBbind ligands to predict the ARMSD to their closest bioactive conformation, and evaluated the early enrichment of bioactive-like conformations when ranking conformers by AtNN prediction. AtNN ranking was compared with bioactivity-unaware baselines such as ascending Sage force field energy ranking, and a slower bioactivity-based baseline ranking by ascending Torsion Fingerprint Deviation to the Maximum Common Substructure to the most similar molecule in the training set (TFD2SimRefMCS). On test sets from random ligand splits of PDBbind, ranking conformers using ComENet, the AtNN encoding the most 3D information, leads to early enrichment of bioactive-like conformations with a median BEDROC of 0.29 ± 0.02, outperforming the best bioactivity-unaware Sage energy ranking baseline (median BEDROC of 0.18 ± 0.02), and performing on a par with the bioactivity-based TFD2SimRefMCS baseline (median BEDROC of 0.31 ± 0.02). The improved performance of the AtNN and TFD2SimRefMCS baseline is mostly observed on test set ligands that bind proteins similar to proteins observed in the training set. On a more challenging subset of flexible molecules, the bioactivity-unaware baselines showed median BEDROCs up to 0.02, while AtNNs and TFD2SimRefMCS showed median BEDROCs between 0.09 and 0.13. When performing rigid ligand re-docking of PDBbind ligands with GOLD using the 1% top-ranked conformers, ComENet ranked conformers showed a higher successful docking rate than bioactivity-unaware baselines, with a rate of 0.48 ± 0.02 compared to CSD probability baseline with a rate of 0.39 ± 0.02. Similarly, on a pharmacophore searching experiment, selecting the 20% top-ranked conformers ranked by ComENet showed higher hit rate compared to baselines. Hence, the approach presented here uses AtNNs successfully to focus conformer ensembles towards bioactive-like conformations, representing an opportunity to reduce computational expense in virtual screening applications on known targets that require input conformations.

Identification of potential inhibitors targeting yellow fever virus helicase through ligand and structure-based computational studies

Yellow fever is a flavivirus having plus-sensed RNA which encodes a single polyprotein. Host proteases cut this polyprotein into seven nonstructural proteins including a vital NS3 protein. The present study aims to identify the most effective inhibitor against the helicase (NS3) using different advanced ligand and structure-based computational studies. A set of 300 ligands was selected against helicase by chemical structural similarity model, which are similar to S-adenosyl-l-cysteine using infiniSee. This tool screens billions of compounds through a similarity search from in-built chemical spaces (CHEMriya, Galaxi, KnowledgeSpace and REALSpace). The pharmacophore was designed from ligands in the library that showed same features. According to the sequence of ligands, six compounds (29, 87, 99, 116, 148, and 208) were taken for pharmacophore designing against helicase protein. Subsequently, compounds from the library which showed the best pharmacophore shared-features were docked using FlexX functionality of SeeSAR and their optibrium properties were analyzed. Afterward, their ADME was improved by replacing the unfavorable fragments, which resulted in the generation of new compounds. The selected best compounds (301, 302, 303 and 304) were docked using SeeSAR and their pharmacokinetics and toxicological properties were evaluated using SwissADME. The optimal inhibitor for yellow fever helicase was 2-amino-N-(4-(dimethylamino)thiazol-2-yl)-4-methyloxazole-5-carboxamide (302), which exhibits promising potential for drug development.Communicated by Ramaswamy H. Sarma.

New 2-alkoxycyanopyridine derivatives as inhibitors of EGFR, HER2, and DHFR: Synthesis, anticancer evaluation, and molecular modeling studies

New series of substituted 2-alkoxycyanopyridine derivatives were synthesized and evaluated for their in vitro and in vivo anticancer activities. Comparing the evaluated activities against cancer cell lines to the broad-spectrum anticancer doxorubicin, and the kinase inhibitor sorafenib, compounds 3a, 4b, 4c, 7a, and 8d demonstrated superior anticancer efficacy with elevated safety profiles and selectivity indices, particularly against MCF7 breast cancer. For exploration of their mechanism of action, assays for inhibition of EGFR, HER2 kinase, and DHFR were performed. The promising synthesized compounds exhibited potent dual kinase EGFR/HER2 inhibitory activity with IC50values of 0.248/0.156 μM for 4b and 0.138/0.092 μM for 4c. Additionally, with IC50 values of 0.138 and 0.193 M, respectively, 4b and 4c had the greatest DHFR inhibitory activity that was comparable to methotrexate. In the MCF7 breast cancer cell line, they caused arrest at the S phase of the cell cycle and exhibited apoptosis induction activity. With restored caspase-3 immunoexpression, the anti-breast cancer assay performed in vivo of 4b and 4c demonstrated a substantial decrease in tumor volume. Results from molecular modeling were in agreement with biological assays proving the importance of the 3-caynopyridine, two substituted phenyl rings attached to central pyridine ring, and propoxy side chain moieties for binding with the receptors. As 4c works by inhibiting both EGFR/HER2 kinase, DHFR enzymes, in addition to cellular apoptosis, it could be viewed as a model of compounds possessing a multi-targeting anticancer activity. Collectively, compounds 4b and 4c might represent prototypes for further development as anticancer molecules.

Structural determinants of sphingosine-1-phosphate receptor selectivity

Fingolimod, the prodrug of fingolimod-1-phosphate (F1P), was the first sphingosine-1-phosphate receptor (S1PR) modulator approved for multiple sclerosis. F1P unselectively targets all five S1PR subtypes. While agonism (functional antagonism via receptor internalization) at S1PR1 leads to the desired immune modulatory effects, agonism at S1PR3 is associated with cardiac adverse effects. This motivated the development of S1PR3 -sparing compounds and led to a second generation of S1PR1,5 -selective ligands like siponimod and ozanimod. Our method combines molecular dynamics simulations and three-dimensional pharmacophores (dynophores) and enables the elucidation of S1PR subtype-specific binding site characteristics, visualizing also subtle differences in receptor-ligand interactions. F1P and the endogenous ligand sphingosine-1-phosphate bind to the orthosteric pocket of all S1PRs, but show different binding mode dynamics, uncovering potential starting points for the development of subtype-specific ligands. Our study contributes to the mechanistic understanding of the selectivity profile of approved drugs like ozanimod and siponimod and pharmaceutical tool compounds like CYM5541.

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

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

Development of Potential Multi-Target Inhibitors for Human Cholinesterases and Beta-Secretase 1: A Computational Approach

Alzheimer's disease causes chronic neurodegeneration and is the leading cause of dementia in the world. The causes of this disease are not fully understood but seem to involve two essential cerebral pathways: cholinergic and amyloid. The simultaneous inhibition of AChE, BuChE, and BACE-1, essential enzymes involved in those pathways, is a promising therapeutic approach to treat the symptoms and, hopefully, also halt the disease progression. This study sought to identify triple enzymatic inhibitors based on stereo-electronic requirements deduced from molecular modeling of AChE, BuChE, and BACE-1 active sites. A pharmacophore model was built, displaying four hydrophobic centers, three hydrogen bond acceptors, and one positively charged nitrogen, and used to prioritize molecules found in virtual libraries. Compounds showing adequate overlapping rates with the pharmacophore were subjected to molecular docking against the three enzymes and those with an adequate docking score (n = 12) were evaluated for physicochemical and toxicological parameters and commercial availability. The structure exhibiting the greatest inhibitory potential against all three enzymes was subjected to molecular dynamics simulations (100 ns) to assess the stability of the inhibitor-enzyme systems. The results of this in silico approach indicate ZINC1733 can be a potential multi-target inhibitor of AChE, BuChE, and BACE-1, and future enzymatic assays are planned to validate those results.

Clustering Protein Binding Pockets and Identifying Potential Drug Interactions: A Novel Ligand-Based Featurization Method

Protein-ligand interactions are essential to drug discovery and drug development efforts. Desirable on-target or multitarget interactions are the first step in finding an effective therapeutic, while undesirable off-target interactions are the first step in assessing safety. In this work, we introduce a novel ligand-based featurization and mapping of human protein pockets to identify closely related protein targets and to project novel drugs into a hybrid protein-ligand feature space to identify their likely protein interactions. Using structure-based template matches from PDB, protein pockets are featured by the ligands that bind to their best co-complex template matches. The simplicity and interpretability of this approach provide a granular characterization of the human proteome at the protein-pocket level instead of the traditional protein-level characterization by family, function, or pathway. We demonstrate the power of this featurization method by clustering a subset of the human proteome and evaluating the predicted cluster associations of over 7000 compounds.

Conformer Generation for Structure-Based Drug Design: How Many and How Good?

Conformer generation, the assignment of realistic 3D coordinates to a small molecule, is fundamental to structure-based drug design. Conformational ensembles are required for rigid-body matching algorithms, such as shape-based or pharmacophore approaches, and even methods that treat the ligand flexibly, such as docking, are dependent on the quality of the provided conformations due to not sampling all degrees of freedom (e.g., only sampling torsions). Here, we empirically elucidate some general principles about the size, diversity, and quality of the conformational ensembles needed to get the best performance in common structure-based drug discovery tasks. In many cases, our findings may parallel "common knowledge" well-known to practitioners of the field. Nonetheless, we feel that it is valuable to quantify these conformational effects while reproducing and expanding upon previous studies. Specifically, we investigate the performance of a state-of-the-art generative deep learning approach versus a more classical geometry-based approach, the effect of energy minimization as a postprocessing step, the effect of ensemble size (maximum number of conformers), and construction (filtering by root-mean-square deviation for diversity) and how these choices influence the ability to recapitulate bioactive conformations and perform pharmacophore screening and molecular docking.

Small-Molecule Conformer Generators: Evaluation of Traditional Methods and AI Models on High-Quality Data Sets

Small-molecule conformer generation (SMCG) is an extremely important task in both ligand- and structure-based computer-aided drug design, especially during the hit discovery phase. Recently, a multitude of artificial intelligence (AI) models tailored for SMCG have emerged. Despite developers typically furnishing performance evaluation data upon releasing their AI models, a comprehensive and equitable performance comparison between AI models and conventional methods is still lacking. In this study, we curated a new benchmarking data set comprising 3354 high-quality ligand bioactive conformations. Subsequently, we conducted a systematic assessment of the performance of four widely adopted traditional methods (i.e., ConfGenX, Conformator, OMEGA, and RDKit ETKDG) and five AI models (i.e., ConfGF, DMCG, GeoDiff, GeoMol, and torsional diffusion) in the tasks of reproducing bioactive and low-energy conformations of small molecules. In the former task, the AI models have no advantage, particularly with a maximum ensemble size of 1. Even the best-performing AI model GeoMol is still worse than any of the tested traditional methods. Conversely, in the latter task, the torsional diffusion model shows obvious advantages, surpassing the best-performing traditional method ConfGenX by 26.09 and 12.97% on the COV-R and COV-P metrics, respectively. Furthermore, the influence of force field-based fine-tuning on the quality of the generated conformers was also discussed. Finally, a user-friendly Web server called fastSMCG was developed to enable researchers to rapidly and flexibly generate small-molecule conformers using both traditional and AI methods. We anticipate that our work will offer valuable practical assistance to the scientific community in this field.

Expanding the Structural Diversity at the Phenylene Core of Ligands for the von Hippel-Lindau E3 Ubiquitin Ligase: Development of Highly Potent Hypoxia-Inducible Factor-1α Stabilizers

Hypoxia-inducible factor-1α (HIF-1α) constitutes the principal mediator of cellular adaptation to hypoxia in humans. The HIF-1α protein level and activity are tightly regulated by the ubiquitin E3 ligase von Hippel-Lindau (VHL). Here, we performed a structure-guided and bioactivity-driven design of new VHL inhibitors. Our iterative and combinatorial strategy focused on chemical variability at the phenylene unit and encompassed further points of diversity. The exploitation of tailored phenylene fragments and the stereoselective installation of the benzylic methyl group provided potent VHL ligands. Three high-resolution structures of VHL-ligand complexes were determined, and bioactive conformations of these ligands were explored. The most potent inhibitor (30) exhibited dissociation constants lower than 40 nM, independently determined by fluorescence polarization and surface plasmon resonance and an enhanced cellular potency, as evidenced by its superior ability to induce HIF-1α transcriptional activity. Our work is anticipated to inspire future efforts toward HIF-1α stabilizers and new ligands for proteolysis-targeting chimera (PROTAC) degraders.

Kinetics, Thermodynamics, and Structural Effects of Quinoline-2-Carboxylates, Zinc-Binding Inhibitors of New Delhi Metallo-β-lactamase-1 Re-sensitizing Multidrug-Resistant Bacteria for Carbapenems

Carbapenem resistance mediated by metallo-β-lactamases (MBL) such as New Delhi metallo-β-lactamase-1 (NDM-1) has become a major factor threatening the efficacy of essential β-lactam antibiotics. Starting from hit fragment dipicolinic acid (DPA), 8-hydroxy- and 8-sulfonamido-quinoline-2-carboxylic acids were developed as inhibitors of NDM-1 with highly improved inhibitory activity and binding affinity. The most active compounds formed reversibly inactive ternary protein-inhibitor complexes with two zinc ions as proven by native protein mass spectrometry and bio-layer interferometry. Modification of the NDM-1 structure with remarkable entropic gain was shown by isothermal titration calorimetry and NMR spectroscopy of isotopically labeled protein. The best compounds were potent inhibitors of NDM-1 and other representative MBL with no or little inhibition of human zinc-binding enzymes. These inhibitors significantly reduced the minimum inhibitory concentrations (MIC) of meropenem for multidrug-resistant bacteria recombinantly expressing blaNDM-1 as well as for several multidrug-resistant clinical strains at concentrations non-toxic to human cells.

Navigating bioactivity space in anti-tubercular drug discovery through the deployment of advanced machine learning models and cheminformatics tools: a molecular modeling based retrospective study

Mycobacterium tuberculosis is the bacterial strain that causes tuberculosis (TB). However, multidrug-resistant and extensively drug-resistant tuberculosis are significant obstacles to effective treatment. As a result, novel therapies against various strains of M. tuberculosis have been developed. Drug development is a lengthy procedure that includes identifying target protein and isolation, preclinical testing of the drug, and various phases of a clinical trial, etc., can take decades for a molecule to reach the market. Computational approaches such as QSAR, molecular docking techniques, and pharmacophore modeling have aided drug development. In this review article, we have discussed the various techniques in tuberculosis drug discovery by briefly introducing them and their importance. Also, the different databases, methods, approaches, and software used in conducting QSAR, pharmacophore modeling, and molecular docking have been discussed. The other targets targeted by these techniques in tuberculosis drug discovery have also been discussed, with important molecules discovered using these computational approaches. This review article also presents the list of drugs in a clinical trial for tuberculosis found drugs. Finally, we concluded with the challenges and future perspectives of these techniques in drug discovery.

TALAIA: a 3D visual dictionary for protein structures

MOTIVATION

Graphical analysis of the molecular structure of proteins can be very complex. Full-atom representations retain most geometric information but are generally crowded, and key structural patterns can be challenging to identify. Non-full-atom representations could be more instructive on physicochemical aspects but be insufficiently detailed regarding shapes (e.g. entity beans-like models in coarse grain approaches) or simple properties of amino acids (e.g. representation of superficial electrostatic properties). In this work, we present TALAIA a visual dictionary that aims to provide another layer of structural representations.TALAIA offers a visual grammar that combines simple representations of amino acids while retaining their general geometry and physicochemical properties. It uses unique objects, with differentiated shapes and colors to represent amino acids. It makes easier to spot crucial molecular information, including patches of amino acids or key interactions between side chains. Most conventions used in TALAIA are standard in chemistry and biochemistry, so experimentalists and modelers can rapidly grasp the meaning of any TALAIA depiction.

RESULTS

We propose TALAIA as a tool that renders protein structures and encodes structure and physicochemical aspects as a simple visual grammar. The approach is fast, highly informative, and intuitive, allowing the identification of possible interactions, hydrophobic patches, and other characteristic structural features at first glance. The first implementation of TALAIA can be found at https://github.com/insilichem/talaia.

Unveiling the interdomain dynamics of type II DNA topoisomerase through all-atom simulations: Implications for understanding its catalytic cycle

Type IIA DNA topoisomerases are complex molecular nanomachines that manage topological states of the DNA molecule in the cell and play a crucial role in cellular processes such as cell division and transcription. They are also established targets of cancer chemotherapy. Starting from the available crystal structure of a fully catalytic topoisomerase IIA homodimer from Saccharomyces cerevisiae, we constructed three states of this molecular motor primarily changing the configurations of the DNA segment bound in the DNA gate and performed μs-long all-atom molecular simulations. A comprehensive analysis revealed a sliding motion within the DNA gate and a teamwork between the N-gate and DNA gate that may be associated with the necessary molecular events that allow passage of the T-segment of DNA. The observed movement of the ATPase dimer relative to the DNA domain was reflected in different interaction patterns between the K-loops of the transducer domain and the B-A-B form of the bound DNA. Based on the obtained results, we mapped simulated configurations to the structures in the proposed catalytic cycle through which type IIA topoisomerases exert their function and discussed the possible transition events. The results extend our understanding of the mechanism of action of type IIA topoisomerases and provide an atomistic interpretation of some of the observed features of these molecular motors.

Structure-Based Drug Design of RdRp Inhibitors against SARS-CoV-2

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused a worldwide pandemic since 2019, spreading rapidly and posing a significant threat to human health and life. With over 6 billion confirmed cases of the virus, the need for effective therapeutic drugs has become more urgent than ever before. RNA-dependent RNA polymerase (RdRp) is crucial in viral replication and transcription, catalysing viral RNA synthesis and serving as a promising therapeutic target for developing antiviral drugs. In this article, we explore the inhibition of RdRp as a potential treatment for viral diseases, analysing the structural information of RdRp in virus proliferation and summarizing the reported inhibitors' pharmacophore features and structure-activity relationship profiles. We hope that the information provided by this review will aid in structure-based drug design and aid in the global fight against SARS-CoV-2 infection.

ProfhEX: AI-based platform for small molecules liability profiling

Off-target drug interactions are a major reason for candidate failure in the drug discovery process. Anticipating potential drug's adverse effects in the early stages is necessary to minimize health risks to patients, animal testing, and economical costs. With the constantly increasing size of virtual screening libraries, AI-driven methods can be exploited as first-tier screening tools to provide liability estimation for drug candidates. In this work we present ProfhEX, an AI-driven suite of 46 OECD-compliant machine learning models that can profile small molecules on 7 relevant liability groups: cardiovascular, central nervous system, gastrointestinal, endocrine, renal, pulmonary and immune system toxicities. Experimental affinity data was collected from public and commercial data sources. The entire chemical space comprised 289'202 activity data for a total of 210'116 unique compounds, spanning over 46 targets with dataset sizes ranging from 819 to 18896. Gradient boosting and random forest algorithms were initially employed and ensembled for the selection of a champion model. Models were validated according to the OECD principles, including robust internal (cross validation, bootstrap, y-scrambling) and external validation. Champion models achieved an average Pearson correlation coefficient of 0.84 (SD of 0.05), an R² determination coefficient of 0.68 (SD = 0.1) and a root mean squared error of 0.69 (SD of 0.08). All liability groups showed good hit-detection power with an average enrichment factor at 5% of 13.1 (SD of 4.5) and AUC of 0.92 (SD of 0.05). Benchmarking against already existing tools demonstrated the predictive power of ProfhEX models for large-scale liability profiling. This platform will be further expanded with the inclusion of new targets and through complementary modelling approaches, such as structure and pharmacophore-based models. ProfhEX is freely accessible at the following address: https://profhex.exscalate.eu/ .

Recent Advances in Computer-Aided Structure-Based Drug Design on Ion Channels

Ion channels play important roles in fundamental biological processes, such as electric signaling in cells, muscle contraction, hormone secretion, and regulation of the immune response. Targeting ion channels with drugs represents a treatment option for neurological and cardiovascular diseases, muscular degradation disorders, and pathologies related to disturbed pain sensation. While there are more than 300 different ion channels in the human organism, drugs have been developed only for some of them and currently available drugs lack selectivity. Computational approaches are an indispensable tool for drug discovery and can speed up, especially, the early development stages of lead identification and optimization. The number of molecular structures of ion channels has considerably increased over the last ten years, providing new opportunities for structure-based drug development. This review summarizes important knowledge about ion channel classification, structure, mechanisms, and pathology with the main focus on recent developments in the field of computer-aided, structure-based drug design on ion channels. We highlight studies that link structural data with modeling and chemoinformatic approaches for the identification and characterization of new molecules targeting ion channels. These approaches hold great potential to advance research on ion channel drugs in the future.

Combining the Classification and Pharmacophore Approaches to Understand Homogeneous Olfactory Perceptions at Peripheral Level: Focus on Two Aroma Mixtures

The mechanisms involved in the homogeneous perception of odorant mixtures remain largely unknown. With the aim of enhancing knowledge about blending and masking mixture perceptions, we focused on structure-odor relationships by combining the classification and pharmacophore approaches. We built a dataset of about 5000 molecules and their related odors and reduced the multidimensional space defined by 1014 fingerprints representing the structures to a tridimensional 3D space using uniform manifold approximation and projection (UMAP). The self-organizing map (SOM) classification was then performed using the 3D coordinates in the UMAP space that defined specific clusters. We explored the allocating in these clusters of the components of two aroma mixtures: a blended mixture (red cordial (RC) mixture, 6 molecules) and a masking binary mixture (isoamyl acetate/whiskey-lactone [IA/WL]). Focusing on clusters containing the components of the mixtures, we looked at the odor notes carried by the molecules belonging to these clusters and also at their structural features by pharmacophore modeling (PHASE). The obtained pharmacophore models suggest that WL and IA could have a common binding site(s) at the peripheral level, but that would be excluded for the components of RC. In vitro experiments will soon be carried out to assess these hypotheses.

Deep Generation Model Guided by the Docking Score for Active Molecular Design

A deep generation model, as a novel drug design and discovery tool, shows obvious advantages in generating compounds with novel backbones and has been applied successfully in the field of drug discovery. However, it is still a challenge to generate molecules with expected properties, especially high activity. Here, to obtain compounds both with novelty and high activity to a target, we proposed a conditional molecular generation model COMG by considering the docking score and 3D pharmacophore matching during molecular generation. The proposed model was based on the conditional variational autoencoder architecture constrained by the pharmacophore matching score. During Bayesian optimization, the docking score was applied to enhance the target relevance of generated compounds. Furthermore, to overcome the problem of high structural similarity caused by Bayesian optimization, the idea of the scaffold memory unit was also introduced. The evaluation results of COMG show that our model not only can improve the structural diversity of generated molecules but also can effectively improve the proportion of target-related drug-active molecules. The obtained results indicate that our proposed model COMG is a useful drug design tool.

Hybridization into a Bitopic Ligand Increased Muscarinic Receptor Activation for Isopilocarpine but Not for Pilocarpine Derivatives

Pilocarpine (1), a secondary metabolite of several Pilocarpus species, is a therapeutically used partial agonist of muscarinic acetylcholine receptors (mAChRs). The available pharmacological data and structure-activity relationships do not provide comparable data for all five receptor subtypes. In this study, pilocarpine (1), its epimer isopilocarpine (2), racemic analogues pilosinine (3) and desmethyl pilosinine (4), and the respective hybrid ligands with a naphmethonium fragment (5-C6 to 8-C6) were synthesized and analyzed in mini-G nano-BRET assays at the five mAChRs. In line with earlier studies, pilocarpine was the most active compound among the orthosteric ligands 1-4. Computational docking of pilocarpine and isopilocarpine to the active M2 receptor suggests that the trans-configuration of isopilocarpine leads to a loss of the hydrogen bond from the lactone carbonyl to N6.52, explaining the lower activity of isopilocarpine. Hybrid formation of pilocarpine (1) and isopilocarpine (2) led to an inverted activity rank, with the trans-configured isopilocarpine hybrid (6-C6) being more active. The hydrogen bond of interest is formed by the isopilocarpine hybrid (6-C6) but not by the pilocarpine hybrid (5-C6). Hybridization thus leads to a modified binding mode of the orthosteric moiety, as the binding mode of the hybrid is dominated by the high-affinity allosteric moiety.

Computational approaches streamlining drug discovery

Computer-aided drug discovery has been around for decades, although the past few years have seen a tectonic shift towards embracing computational technologies in both academia and pharma. This shift is largely defined by the flood of data on ligand properties and binding to therapeutic targets and their 3D structures, abundant computing capacities and the advent of on-demand virtual libraries of drug-like small molecules in their billions. Taking full advantage of these resources requires fast computational methods for effective ligand screening. This includes structure-based virtual screening of gigascale chemical spaces, further facilitated by fast iterative screening approaches. Highly synergistic are developments in deep learning predictions of ligand properties and target activities in lieu of receptor structure. Here we review recent advances in ligand discovery technologies, their potential for reshaping the whole process of drug discovery and development, as well as the challenges they encounter. We also discuss how the rapid identification of highly diverse, potent, target-selective and drug-like ligands to protein targets can democratize the drug discovery process, presenting new opportunities for the cost-effective development of safer and more effective small-molecule treatments.

Development of epitopephore-based rational hapten design strategy: A combination of theoretical evidence and experimental validation

Antibody is the key biomolecule that governing the sensitivity and specificity of an immunoassay for chemical compound, also named hapten molecule. Obviously, predication of hapten effectiveness before chemical synthesis is beneficial to boost success, save cost and improve controllability. Here, we proposed and evaluated an epitopephore based rational hapten design (ERHD) to assist antibody production to chemical compound, combining theoretical evidence and then experimental validation by using dinitrocarbanilide (DNC) as a model analyte. Briefly, epitopephores of DNC were firstly generated by HipHop algorithm after features mapping. A homemade drug database also containing reported fragment haptens (HFR) and new designed full hapten (HFU) were constructed, and then was virtually screened by using generated epitopephore followed by structural analysis and visual inspection. The DNC haptens based on the selected hits were further identified by Density Functional Theory before total synthesis. To prove and clarify the usability of the ERHD, two retrieved HFU haptens, one non-retrieved HFU hapten and three non-retrieved HFR haptens were all selected to produce monoclonal antibodies (mAbs) for comparison purpose. A maximal 6000-fold increased affinity of mAb from retrieved HFU than HFR was observed, while, non-retrieved HFU failed to produce antibody to DNC. More importantly, mAbs from HFU haptens provided highly specificity to DNC, while, mAbs from HFR haptens could recognize 15 others analogues. We then constructed antibody structure and investigated molecular recognition of the mAbs to DNC, well supporting the rationality of the ERHD. Lastly, an icELISA was developed for DNC with an IC₅₀ value as low as 0.19 ng mL^-1 with high specificity, which has never achieved before.

Phenotypic Discovery of Thiocarbohydrazone with Anticancer Properties and Catalytic Inhibition of Human DNA Topoisomerase IIα

Phenotypic screening of α-substituted thiocarbohydrazones revealed promising activity of 1,5-bis(salicylidene)thiocarbohydrazide against leukemia and breast cancer cells. Supplementary cell-based studies indicated an impairment of DNA replication via the ROS-independent pathway. The structural similarity of α-substituted thiocarbohydrazone to previously published thiosemicarbazone catalytic inhibitors targeting the ATP-binding site of human DNA topoisomerase IIα prompted us to investigate the inhibition activity on this target. Thiocarbohydrazone acted as a catalytic inhibitor and did not intercalate the DNA molecule, which validated their engagement with this cancer target. A comprehensive computational assessment of molecular recognition for a selected thiosemicarbazone and thiocarbohydrazone provided useful information for further optimization of this discovered lead compound for chemotherapeutic anticancer drug discovery.

Arabinofuranosyl Thymine Derivatives-Potential Candidates against Cowpox Virus: A Computational Screening Study

Cowpox is caused by a DNA virus known as the cowpox virus (CPXV) belonging to the Orthopoxvirus genus in the family Poxviridae. Cowpox is a zoonotic disease with the broadest host range among the known poxviruses. The natural reservoir hosts of CPXV are wild rodents. Recently, the cases of orthopoxviral infections have been increasing worldwide, and cowpox is considered the most common orthopoxviral infection in Europe. Cowpox is often a self-limiting disease, although cidofovir or anti-vaccinia gammaglobulin can be used in severe and disseminated cases of human cowpox. In this computational study, a molecular docking analysis of thymine- and arabinofuranosyl-thymine-related structures (1-21) on two cowpox-encoded proteins was performed with respect to the cidofovir standard and a 3D ligand-based pharmacophore model was generated. Three chemical structures (PubChem IDs: 123370001, 154137224, and 90413364) were identified as potential candidates for anti-cowpox agents. Further studies combining in vitro and in silico molecular dynamics simulations to test the stability of these promising compounds could effectively improve the future design of cowpox virus inhibitors, as molecular docking studies are not sufficient to consider a ligand a potential drug.

Revolutionizing antiretroviral therapy for human immunodeficiency virus/AIDS: A computational approach using molecular docking, virtual screening, and 3D pharmacophore building to address therapeutic failure and propose highly effective candidates

OBJECTIVES

In the context of human immunodeficiency virus (HIV) treatment, the emergence of therapeutic failures with existing antiretroviral drugs presents a significant challenge. This study aims to employ advanced molecular modeling techniques to identify potential alternatives to current antiretroviral agents.

METHODS

The study focuses on three essential classes of antiretroviral drugs: nucleoside reverse transcriptase inhibitors (NRTIs), non-nucleoside reverse transcriptase inhibitors (NNRTIs), and protease inhibitors (PIs). Computational analyses were performed on a database of 3,343,652 chemical molecules to evaluate their binding affinities, pharmacokinetic properties, and interactions with viral reverse transcriptase and protease enzymes. Molecular docking, virtual screening, and 3D pharmacophore modeling were utilized to identify promising candidates.

RESULTS

Molecular docking revealed compounds with high binding energies and strong interactions at the active sites of target enzymes. Virtual screening narrowed down potential candidates with favorable pharmacological profiles. 3D pharmacophore modeling identified crucial structural features for effective binding. Overall, two molecules for class 1, 7 molecules for class 2, and 2 molecules for class 3 were selected. These compounds exhibited robust binding affinities, interactions with target enzymes, and improved pharmacokinetic properties, showing promise for more effective HIV treatments in cases of therapeutic failures.

CONCLUSION

The combination of molecular docking, virtual screening, and 3D pharmacophore modeling yielded lead compounds that hold potential for addressing HIV therapeutic failures. Further experimental investigations are essential to validate the efficacy and safety of these compounds, with the ultimate goal of advancing toward clinical applications in HIV management.

Computer-Aided Drug Design: An Update

Computer-aided drug design (CADD) approaches are playing an increasingly important role in understanding the fundamentals of ligand-receptor interactions and helping medicinal chemists design therapeutics. About 5 years ago, we presented a chapter devoted to an overview of CADD methods and covered typical CADD protocols including structure-based drug design (SBDD) and ligand-based drug design (LBDD) approaches that were frequently used in the antibiotic drug design process. Advances in computational hardware and algorithms and emerging CADD methods are enhancing the accuracy and ability of CADD in drug design and development. In this chapter, an update to our previous chapter is provided with a focus on new CADD approaches from our laboratory and other peers that can be employed to facilitate the development of antibiotic therapeutics.

A Guide to In Silico Drug Design

The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.

Application of Computational Biology and Artificial Intelligence in Drug Design

Traditional drug design requires a great amount of research time and developmental expense. Booming computational approaches, including computational biology, computer-aided drug design, and artificial intelligence, have the potential to expedite the efficiency of drug discovery by minimizing the time and financial cost. In recent years, computational approaches are being widely used to improve the efficacy and effectiveness of drug discovery and pipeline, leading to the approval of plenty of new drugs for marketing. The present review emphasizes on the applications of these indispensable computational approaches in aiding target identification, lead discovery, and lead optimization. Some challenges of using these approaches for drug design are also discussed. Moreover, we propose a methodology for integrating various computational techniques into new drug discovery and design.

Towards an Enrichment Optimization Algorithm (EOA)-based Target Specific Docking Functions for Virtual Screening

Docking-based virtual screening (VS) is a common starting point in many drug discovery projects. While ligand-based approaches may sometimes provide better results, the advantage of docking lies in its ability to provide reliable ligand binding modes and approximated binding free energies, two factors that are important for hit selection and optimization. Most docking programs were developed to be as general as possible and consequently their performances on specific targets may be sub-optimal. With this in mind, in this work we present a method for the development of target-specific scoring functions using our recently reported Enrichment Optimization Algorithm (EOA). EOA derives QSAR models in the form of multiple linear regression (MLR) equations by optimizing an enrichment-like metric. Since EOA requires target-specific active and inactive (or decoy) compounds, we retrieved such data for six targets from the DUD-E database, and used them to re-derive the weights associated with the components that make up GOLD's ChemPLP scoring function yielding target-specific, modified functions. We then used the original ChemPLP function in small-scale VS experiments on the six targets and subsequently rescored the resulting poses with the modified functions. In addition, we used the modified functions for compounds re-docking. We found that in many although not all cases, either rescoring the original ChemPLP poses or repeating the entire docking process with the modified functions, yielded better results in terms of AUC and EF_1% , two metrics, common for the evaluation of VS performances. While work on additional datasets and docking tools is clearly required, we propose that the results obtained thus far hint to the potential benefits in using EOA-based optimization for the derivation of target-specific functions in the context of virtual screening. To this end, we discuss the downsides of the methods and how it could be improved.

Mind the Gap—Deciphering GPCR Pharmacology Using 3D Pharmacophores and Artificial Intelligence

G protein-coupled receptors (GPCRs) are amongst the most pharmaceutically relevant and well-studied protein targets, yet unanswered questions in the field leave significant gaps in our understanding of their nuanced structure and function. Three-dimensional pharmacophore models are powerful computational tools in in silico drug discovery, presenting myriad opportunities for the integration of GPCR structural biology and cheminformatics. This review highlights success stories in the application of 3D pharmacophore modeling to de novo drug design, the discovery of biased and allosteric ligands, scaffold hopping, QSAR analysis, hit-to-lead optimization, GPCR de-orphanization, mechanistic understanding of GPCR pharmacology and the elucidation of ligand–receptor interactions. Furthermore, advances in the incorporation of dynamics and machine learning are highlighted. The review will analyze challenges in the field of GPCR drug discovery, detailing how 3D pharmacophore modeling can be used to address them. Finally, we will present opportunities afforded by 3D pharmacophore modeling in the advancement of our understanding and targeting of GPCRs.

Pharmacophore study, molecular docking and molecular dynamic simulation of virgin coconut oil derivatives as anti-inflammatory agent against COX-2

Background

Virgin coconut oil is mostly made up of saturated fatty acids in which approximately 72% are medium chain triglycerides. Medium chain triglycerides can be digested into medium chain fatty acids and medium chain monoglycerides which are bioactive components. Therefore, it is very important to study the in-silico ability of some Virgin coconut oil derivatives, namely, medium chain fatty acids and medium chain monoglycerides to inhibit Cyclooxygenase 2 (COX-2) protein for prevention of excessive inflammatory response.

Results

Pharmacophore study displayed monolaurin with two hydrogen bond donor, three hydrogen bond acceptor and five hydrophobic interactions, while lauric acid presented two hydrogen bond acceptor, five hydrophobic interactions and a negative ion interaction. Molecular docking underlined the ability of monolaurin in the inhibition of COX-2 protein which causes inflammatory action with a decent result of energy binding affinity of - 7.58 kcal/mol and 15 interactions out of which 3 are strong hydrogen bond with TYR385 (3.00 Å), PHE529 (2.77 Å), and GLY533 (3.10 Å) residues of the protein. Monolaurin was employed as hydrogen bond acceptor to the side of residue TYR385 of COX-2 protein with an occupancy of 67.03% and was observed to be long-living during the entire 1000 frames of the molecular dynamic simulation. The analysis of RMSD score of the Monolaurin-COX-2 complex backbone was calculated to be low (1.137 ± 0.153 Å) and was in a stable range of 0.480 to 1.520 Å. Redocking of this complex still maintained a strong hydrogen bond (2.87 Å) with the main residue TYR385. AMDET results where promising for medium chain fatty acids and medium chain monoglycerides with good physicochemical drug scores.

Conclusions

This can be concluded from the results obtained that the monolaurin has strong interactions with COX-2 protein to disrupt its function due to significant hydrogen bonds and hydrophobic interactions with amino acid residues present in the target protein's active site. These results displayed a very significant anti-inflammatory potential of monolaurin and a new promising drug candidates as anti-inflammatory agent.

Graphical Abstract