Enhancing the Small-Scale Screenable Biological Space beyond Known Chemogenomics Libraries with Gray Chemical Matter─Compounds with Novel Mechanisms from High-Throughput Screening Profiles.

Phenotypic assays have become an established approach to drug discovery. Greater disease relevance is often achieved through cellular models with increased complexity and more detailed readouts, such as gene expression or advanced imaging. However, the intricate nature and cost of these assays impose limitations on their screening capacity, often restricting screens to well-characterized small compound sets such as chemogenomics libraries. Here, we outline a cheminformatics approach to identify a small set of compounds with likely novel mechanisms of action (MoAs), expanding the MoA search space for throughput limited phenotypic assays. Our approach is based on mining existing large-scale, phenotypic high-throughput screening (HTS) data. It enables the identification of chemotypes that exhibit selectivity across multiple cell-based assays, which are characterized by persistent and broad structure activity relationships (SAR). We validate the effectiveness of our approach in broad cellular profiling assays (Cell Painting, DRUG-seq, and Promotor Signature Profiling) and chemical proteomics experiments. These experiments revealed that the compounds behave similarly to known chemogenetic libraries, but with a notable bias toward novel protein targets. To foster collaboration and advance research in this area, we have curated a public set of such compounds based on the PubChem BioAssay dataset and made it available for use by the scientific community.

Cheminformatics approach to identify andrographolide derivatives as dual inhibitors of methyltransferases (nsp14 and nsp16) of SARS-CoV-2.

The Covid-19 pandemic outbreak has accelerated tremendous efforts to discover a therapeutic strategy that targets severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to control viral infection. Various viral proteins have been identified as potential drug targets, however, to date, no specific therapeutic cure is available against the SARS-CoV-2. To address this issue, the present work reports a systematic cheminformatic approach to identify the potent andrographolide derivatives that can target methyltransferases of SARS-CoV-2, i.e. nsp14 and nsp16 which are crucial for the replication of the virus and host immune evasion. A consensus of cheminformatics methodologies including virtual screening, molecular docking, ADMET profiling, molecular dynamics simulations, free-energy landscape analysis, molecular mechanics generalized born surface area (MM-GBSA), and density functional theory (DFT) was utilized. Our study reveals two new andrographolide derivatives (PubChem CID: 2734589 and 138968421) as natural bioactive molecules that can form stable complexes with both proteins via hydrophobic interactions, hydrogen bonds and electrostatic interactions. The toxicity analysis predicts class four toxicity for both compounds with LD50 value in the range of 500-700 mg/kg. MD simulation reveals the stable formation of the complex for both the compounds and their average trajectory values were found to be lower than the control inhibitor and protein alone. MMGBSA analysis corroborates the MD simulation result and showed the lowest energy for the compounds 2734589 and 138968421. The DFT and MEP analysis also predicts the better reactivity and stability of both the hit compounds. Overall, both andrographolide derivatives exhibit good potential as potent inhibitors for both nsp14 and nsp16 proteins, however, in-vitro and in vivo assessment would be required to prove their efficacy and safety in clinical settings. Moreover, the drug discovery strategy aiming at the dual target approach might serve as a useful model for inventing novel drug molecules for various other diseases.

RDCanon: A Python Package for Canonicalizing the Order of Tokens in SMARTS Queries.

SMARTS is a widely used language in cheminformatics for defining substructural queries for database lookups, reaction templates for chemical transformations, and other applications. As an extension to SMILES, many SMARTS patterns can represent the same query. Despite this, no canonicalization algorithm invariant of the line notation sequence or atomic numbering is publicly available. Here, we introduce RDCanon, an open-source Python package that can be used to standardize SMARTS queries. RDCanon is designed to ensure that the sequence of atomic queries remains consistent for all graphs representing the same substructure query and to ensure a canonical sequence of primitives within each individual atom query; furthermore, the algorithm can be applied to canonicalize the order of reactants, agents, and products and their atom map numbers in reaction SMARTS templates. As part of its canonicalization algorithm, RDCanon provides a mechanism in which the canonicalized SMARTS is optimized for speed against specific molecular databases. Several case studies are provided to showcase improved efficiency in substructure matching and retrosynthetic analysis.

PaCh (Packed Chemicals): Computationally Effective Binary Format for Chemical Structure Encoding.

In this work, we propose a versatile molecule and reaction encoding binary data format that aims to bridge the gap between the advantages of SMILES, like local stereo- and implicit hydrogen encoding, and block-structured MDL MOL with a 2D layout and explicit bond encoding, while addressing their respective limitations. Our new format introduces a balance between size efficiency, processing speed, and comprehensive representation, making it well-suited for various applications in cheminformatics, including deep learning, data storage, and searching. By offering an explicit approach to store atom connectivity (including implicit hydrogens), electronic state, stereochemistry, and other crucial molecular attributes, our proposal seeks to enhance data storage efficiency and promote interoperability among different software tools.

Cheminformatics-Guided Cell-Free Exploration of Peptide Natural Products.

There have been significant advances in the flexibility and power of in vitro cell-free translation systems. The increasing ability to incorporate noncanonical amino acids and complement translation with recombinant enzymes has enabled cell-free production of peptide-based natural products (NPs) and NP-like molecules. We anticipate that many more such compounds and analogs might be accessed in this way. To assess the peptide NP space that is directly accessible to current cell-free technologies, we developed a peptide parsing algorithm that breaks down peptide NPs into building blocks based on ribosomal translation logic. Using the resultant data set, we broadly analyze the biophysical properties of these privileged compounds and perform a retrobiosynthetic analysis to predict which peptide NPs could be directly synthesized in augmented cell-free translation reactions. We then tested these predictions by preparing a library of highly modified peptide NPs. Two macrocyclases, PatG and PCY1, were used to effect the head-to-tail macrocyclization of candidate NPs. This retrobiosynthetic analysis identified a collection of high-priority building blocks that are enriched throughout peptide NPs, yet they had not previously been tested in cell-free translation. To expand the cell-free toolbox into this space, we established, optimized, and characterized the flexizyme-enabled ribosomal incorporation of piperazic acids. Overall, these results demonstrate the feasibility of cell-free translation for peptide NP total synthesis while expanding the limits of the technology. This work provides a novel computational tool for exploration of peptide NP chemical space, that could be expanded in the future to allow design of ribosomal biosynthetic pathways for NPs and NP-like molecules.

Flexible scaffold-based cheminformatics approach for polypharmacological drug design.

Effective treatments for complex central nervous system (CNS) disorders require drugs with polypharmacology and multifunctionality, yet designing such drugs remains a challenge. Here, we present a flexible scaffold-based cheminformatics approach (FSCA) for the rational design of polypharmacological drugs. FSCA involves fitting a flexible scaffold to different receptors using different binding poses, as exemplified by IHCH-7179, which adopted a "bending-down" binding pose at 5-HT2AR to act as an antagonist and a "stretching-up" binding pose at 5-HT1AR to function as an agonist. IHCH-7179 demonstrated promising results in alleviating cognitive deficits and psychoactive symptoms in mice by blocking 5-HT2AR for psychoactive symptoms and activating 5-HT1AR to alleviate cognitive deficits. By analyzing aminergic receptor structures, we identified two featured motifs, the "agonist filter" and "conformation shaper," which determine ligand binding pose and predict activity at aminergic receptors. With these motifs, FSCA can be applied to the design of polypharmacological ligands at other receptors.

Computational resources and chemoinformatics for translational health research.

The integration of computational resources and chemoinformatics has revolutionized translational health research. It has offered a powerful set of tools for accelerating drug discovery. This chapter overviews the computational resources and chemoinformatics methods used in translational health research. The resources and methods can be used to analyze large datasets, identify potential drug candidates, predict drug-target interactions, and optimize treatment regimens. These resources have the potential to transform the drug discovery process and foster personalized medicine research. We discuss insights into their various applications in translational health and emphasize the need for addressing challenges, promoting collaboration, and advancing the field to fully realize the potential of these tools in transforming healthcare.

Exploiting Vector Pattern Diversity of Molecular Scaffolds for Cheminformatics Tasks in Drug Discovery.

Chemical diversity is challenging to describe objectively. Despite this, various notions of chemical diversity are used throughout the medicinal chemistry optimization process in drug discovery. In this work, we show the usefulness of considering exploited vectors during different phases of the drug design process to provide a quantitative and objective description of chemical diversity. We have developed a concise and fast approach to enumerate and analyze the exploited vector patterns (EVPs) of molecular compound series, which can then be used in archetypal compound selection tasks, from hit matter identification to hit expansion and lead optimization. We first show that EVPs can be used to assess the progressibility of compounds in a fragment library design exercise. By considering EVPs, we then show how a set of compounds can be prioritized for hit expansion using EVP-based, customizable diversity sampling approaches, reducing the time taken and mitigating human biases. We also show that EVPs are a useful tool to analyze SAR data, offering the chance to uncover correlations between different vectors without predetermining the molecular scaffold structures. The codes used to perform these tasks are presented as easy-to-use Jupyter notebooks, which can be readily adapted for further related tasks.

Synthesis and Chemoinformatic Analysis of Fluorinated Piperidines as 3D Fragments for Fragment-Based Drug Discovery.

The concise synthesis of a small library of fluorinated piperidines from readily available dihydropyridinone derivatives has been described. The effect of the fluorination on different positions has then been evaluated by chemoinformatic tools. In particular, the compounds' pKa's have been calculated, revealing that the fluorine atoms notably lowered their basicity, which is correlated to the affinity for hERG channels resulting in cardiac toxicity. The "lead-likeness" and three-dimensionality have also been evaluated to assess their ability as useful fragments for drug design. A random screening on a panel of representative proteolytic enzymes was then carried out and revealed that one scaffold is recognized by the catalytic pocket of 3CLPro (main protease of SARS-CoV-2 coronavirus).

Chemoinformatic approaches for navigating large chemical spaces.

INTRODUCTION: Large chemical spaces (CSs) include traditional large compound collections, combinatorial libraries covering billions to trillions of molecules, DNA-encoded chemical libraries comprising complete combinatorial CSs in a single mixture, and virtual CSs explored by generative models. The diverse nature of these types of CSs require different chemoinformatic approaches for navigation. AREAS COVERED: An overview of different types of large CSs is provided. Molecular representations and similarity metrics suitable for large CS exploration are discussed. A summary of navigation of CSs in generative models is provided. Methods for characterizing and comparing CSs are discussed. EXPERT OPINION: The size of large CSs might restrict navigation to specialized algorithms and limit it to considering neighborhoods of structurally similar molecules. Efficient navigation of large CSs not only requires methods that scale with size but also requires smart approaches that focus on better but not necessarily larger molecule selections. Deep generative models aim to provide such approaches by implicitly learning features relevant for targeted biological properties. It is unclear whether these models can fulfill this ideal as validation is difficult as long as the covered CSs remain mainly virtual without experimental verification.

Creation of Polymer Datasets with Targeted Backbones for Screening of High-Performance Membranes for Gas Separation.

A simple approach was developed to computationally construct a polymer dataset by combining simplified molecular-input line-entry system (SMILES) strings of a targeted polymer backbone and a variety of molecular fragments. This method was used to create 14 polymer datasets by combining seven polymer backbones and molecules from two large molecular datasets (MOSES and QM9). Polymer backbones that were studied include four polydimethylsiloxane (PDMS) based backbones, poly(ethylene oxide) (PEO), poly(allyl glycidyl ether) (PAGE), and polyphosphazene (PPZ). The generated polymer datasets can be used for various cheminformatics tasks, including high-throughput screening for gas permeability and selectivity. This study utilized machine learning (ML) models to screen the polymers for CO2/CH4 and CO2/N2 gas separation using membranes. Several polymers of interest were identified. The results highlight that employing an ML model fitted to polymer selectivities leads to higher accuracy in predicting polymer selectivity compared to using the ratio of predicted permeabilities.

Authentication of the species identity of squid rings using UHPLC-Q-Orbitrap MS/MS-based lipidome fingerprinting and chemoinformatics.

Species mislabeling of commercial loliginidae squid can undermine important conservation efforts and prevent consumers from making informed decisions. A comprehensive lipidomic fingerprint of Uroteuthis singhalensis, Uroteuthis edulis, and Uroteuthis duvauceli rings was established using high-resolution mass spectrometry-based lipidomics and chemoinformatics analysis. The principal component analysis showed a clear separation of sample groups, with R2X and Q2 values of 0.97 and 0.85 for ESI+ and 0.96 and 0.86 for ESI-, indicating a good model fit. The optimized OPLS-DA and PLS-DA models could discriminate the species identity of validation samples with 100 % accuracy. A total of 67 and 90 lipid molecules were putatively identified as biomarkers in ESI+ and ESI-, respectively. Identified lipids, including PC(40:6), C14 sphingomyelin, PS(O-36:0), and PE(41:4), played an important role in species discrimination. For the first time, this study provides a detailed lipidomics profile of commercially important loliginidae squid and establishes a faster workflow for species authentication.

Cheminformatics identification of phenolics as modulators of penicillin-binding protein-3 of Pseudomonas aeruginosa towards interventive antibacterial therapy.

Antibacterial resistance to ß-lactams in microorganisms has been attributed majorly to alterations in penicillin-binding proteins (PBPs) coupled with ß-lactams' inactivation by ß-lactamase. Consequently, the identification of a novel class of therapeutics with improved modulatory action on the PBPs is imperative and plant secondary metabolites, including phenolics, have found relevance in this regard. For the first time in this study, the over 10,000 phenolics currently known were computationally evaluated against PBP3 of Pseudomonas aeruginosa, a superbug implicated in several nosocomial infections. In doing this, a library of phenolics with an affinity for PBP3 of P. aeruginosa was screened using structure-activity relationship-based pharmacophore and molecular docking approaches. Subsequent thermodynamic screening of the top five phenolics with higher docking scores, more drug-likeness attributes, and feasible synthetic accessibility was achieved through a 120 ns molecular dynamic (MD) simulation. Four of the top five hits had higher binding free energy than cefotaxime (-18.72 kcal/mol), with catechin-3-rhamside having the highest affinity (-28.99 kcal/mol). All the hits were stable at the active site of the PBP3, with catechin-3-rhamside being the most stable (2.14 Å), and established important interactions with Ser294, implicated in the catalytic activity of PBP3. Also, PBP3 became more compact with less fluctuation of the active site amino acid residues following the binding of the hits. These observations are indicative of the potential of the test compounds as PBP3 inhibitors, with catechin-3-rhamside being the most prominent of the compounds that could be further improved for enhanced druggability against PBP3 in vitro and in vivo.Communicated by Ramaswamy H. Sarma.

GUIDEMOL: A Python graphical user interface for molecular descriptors based on RDKit.

GUIDEMOL is a Python computer program based on the RDKit software to process molecular structures and calculate molecular descriptors with a graphical user interface using the tkinter package. It can calculate descriptors already implemented in RDKit as well as grid representations of 3D molecular structures using the electrostatic potential or voxels. The GUIDEMOL app provides easy access to RDKit tools for chemoinformatics users with no programming skills and can be adapted to calculate other descriptors or to trigger other procedures. A command line interface (CLI) is also provided for the calculation of grid representations. The source code is available at https://github.com/jairesdesousa/guidemol.

When Yield Prediction Does Not Yield Prediction: An Overview of the Current Challenges.

Machine Learning (ML) techniques face significant challenges when predicting advanced chemical properties, such as yield, feasibility of chemical synthesis, and optimal reaction conditions. These challenges stem from the high-dimensional nature of the prediction task and the myriad essential variables involved, ranging from reactants and reagents to catalysts, temperature, and purification processes. Successfully developing a reliable predictive model not only holds the potential for optimizing high-throughput experiments but can also elevate existing retrosynthetic predictive approaches and bolster a plethora of applications within the field. In this review, we systematically evaluate the efficacy of current ML methodologies in chemoinformatics, shedding light on their milestones and inherent limitations. Additionally, a detailed examination of a representative case study provides insights into the prevailing issues related to data availability and transferability in the discipline.

Selection of Mexican Medicinal Plants by Identification of Potential Phytochemicals with Anti-Aging, Anti-Inflammatory, and Anti-Oxidant Properties through Network Analysis and Chemoinformatic Screening.

Many natural products have been acquired from plants for their helpful properties. Medicinal plants are used for treating a variety of pathologies or symptoms. The axes of many pathological processes are inflammation, oxidative stress, and senescence. This work is focused on identifying Mexican medicinal plants with potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects through network analysis and chemoinformatic screening of their phytochemicals. We used computational methods to analyze drug-like phytochemicals in Mexican medicinal plants, multi-target compounds, and signaling pathways related to anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence mechanisms. A total of 1373 phytochemicals are found in 1025 Mexican medicinal plants, and 148 compounds showed no harmful functionalities. These compounds displayed comparable structures with reference molecules. Based on their capacity to interact with pharmacological targets, three clusters of Mexican medicinal plants have been established. Curatella americana, Ximenia americana, Malvastrum coromandelianum, and Manilkara zapota all have anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. Plumeria rubra, Lonchocarpus yucatanensis, and Salvia polystachya contained phytochemicals with anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence reported activity. Lonchocarpus guatemalensis, Vallesia glabra, Erythrina oaxacana, and Erythrina sousae have drug-like phytochemicals with potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. Between the drug-like phytochemicals, lonchocarpin, vallesine, and erysotrine exhibit potential anti-oxidant, anti-inflammatory, anti-aging, and anti-senescence effects. For the first time, we conducted an initial virtual screening of selected Mexican medicinal plants, which was subsequently confirmed in vivo, evaluating the anti-inflammatory activity of Lonchocarpus guatemalensis Benth in mice.

Vaccinomics-aided next-generation novel multi-epitope-based vaccine engineering against multidrug resistant Shigella Sonnei: Immunoinformatics and chemoinformatics approaches.

Shigella sonnei is a gram-negative bacterium and is the primary cause of shigellosis in advanced countries. An exceptional rise in the prevalence of the disease has been reported in Asia, the Middle East, and Latin America. To date, no preventive vaccine is available against S. sonnei infections. This pathogen has shown resistances towards both first- and second-line antibiotics. Therefore, an effective broad spectrum vaccine development against shigellosis is indispensable. In the present study, vaccinomics-aided immunoinformatics strategies were pursued to identify potential vaccine candidates from the S. sonnei whole proteome data. Pathogen essential proteins that are non-homologous to human and human gut microbiome proteome set, are feasible candidates for this purpose. Three antigenic outer membrane proteins were prioritized to predict lead epitopes based on reverse vaccinology approach. Multi-epitope-based chimeric vaccines was designed using lead B- and T-cell epitopes combined with suitable linker and adjuvant peptide sequences to enhance immune responses against the designed vaccine. The SS-MEVC construct was prioritized based on multiple physicochemical, immunological properties, and immune-receptors docking scores. Immune simulation analysis predicted strong immunogenic response capability of the designed vaccine construct. The Molecular dynamic simulations analysis ensured stable molecular interactions of lead vaccine construct with the host receptors. In silico restriction and cloning analysis predicted feasible cloning capability of the SS-MEVC construct within the E. coli expression system. The proposed vaccine construct is predicted to be more safe, effective and capable of inducing robust immune responses against S. sonnei infections and may be worthy of examination via in vitro/in vivo assays.

Anonymous Pattern Molecular Fingerprint and its Applications on Property Identification.

Molecular fingerprints are significant cheminformatics tools to map molecules into vectorial space according to their characteristics in diverse functional groups, atom sequences, and other topological structures. In this paper, we investigate a novel molecular fingerprint Anonymous-FP that possesses abundant perception about the underlying interactions shaped in small, medium, and large-scale atom chains. In detail, the possible atom chains from each molecule are sampled and extended as anonymous atom chains using an anonymous encoding manner. After that, the molecular fingerprint Anonymous-FP is embedded into vectorial space in virtue of the Natural Language Processing technique PV-DBOW. Anonymous-FP is studied on molecular property identification via molecule classification experiments on a series of molecule databases and has shown valuable advantages such as less dependence on prior knowledge, rich information content, full structural significance, and high experimental performance. During the experimental verification, the scale of the atom chain or its anonymous pattern is found significant to the overall representation ability of Anonymous-FP. Generally, the typical scale r = 8 could enhance the molecule classification performance, and specifically, Anonymous-FP gains the classification accuracy to above 93% on all NCI datasets.

Cheminformatics-driven discovery of hit compounds against Paracoccidioides spp.

Aims: The development of safe and effective therapies for treating paracoccidioidomycosis using computational strategies were employed to discover anti-Paracoccidioides compounds. Materials & methods: We 1) collected, curated and integrated the largest library of compounds tested against Paracoccidioides spp.; 2) employed a similarity search to virtually screen the ChemBridge database and select nine compounds for experimental evaluation; 3) performed an experimental evaluation to determine the minimum inhibitory concentration and minimum fungicidal concentration as well as cytotoxicity; and 4) employed computational tools to identify potential targets for the most active compounds. Seven compounds presented activity against Paracoccidioides spp. Conclusion: These compounds are new hits with a predicted mechanisms of action, making them potentially attractive to develop new compounds.