Exploration Of Chemical Space Research Articles

Traversing chemical space with active deep learning for low-data drug discovery.

Deep learning is accelerating drug discovery. However, current approaches are often affected by limitations in the available data, in terms of either size or molecular diversity. Active deep learning has high potential for low-data drug discovery, as it allows iterative model improvement during the screening process. However, there are several 'known unknowns' that limit the wider adoption of active deep learning in drug discovery: (1) what the best computational strategies are for chemical space exploration, (2) how active learning holds up to traditional, non-iterative, approaches and (3) how it should be used in the low-data scenarios typical of drug discovery. To provide answers, this study simulates a low-data drug discovery scenario, and systematically analyzes six active learning strategies combined with two deep learning architectures, on three large-scale molecular libraries. We identify the most important determinants of success in low-data regimes and show that active learning can achieve up to a sixfold improvement in hit discovery when compared with traditional screening methods.

Oxadiazolines as Photoreleasable Labels for Drug Target Identification.

Photoaffinity labeling is a widely used technique for studying ligand-protein and protein-protein interactions. Traditional photoaffinity labels utilize nonspecific C-H bond insertion reactions mediated by a highly reactive intermediate. Despite being the most widely used photoaffinity labels, diazirines exhibit limited compatibility with downstream organic reactions and suffer from storage stability concerns. This study introduces oxadiazolines as innovative and complementary photoactivatable labels for addition to the toolbox and demonstrates their application in vitro and through in cellulo labeling experiments. Oxadiazolines can be easily synthesized from ketone moieties and cleaved with 302-330 nm light to cleanly liberate a diazo reactive moiety that can covalently modify nucleophilic amino acid residues. Notably, oxadiazolines are compatible with various organic reaction conditions and functional groups, allowing for the exploration of a large chemical space. Several known inhibitors featuring the oxadiazoline functionality were prepared without affecting their binding affinity. Furthermore, we confirmed the ability of oxadiazolines to form covalent bonds with proteins upon UV-irradiation, both in vitro and in cellulo, yielding comparable results to those of the matched diazirine compounds.

Navigation through high-dimensional chemical space: discovery of Ba5Y13[SiO4]8O8.5 and Ba3Y2[Si2O7]2.

Two compounds were discovered in the well-studied BaO-Y2O3-SiO2 phase field. Two different experimental routines were used for the exploration of this system due to the differences of synthetic conditions and competition with a glass field. The first phase Ba5Y13[SiO4]8O8.5 was isolated through a combination of energy dispersive X-ray spectroscopy analysis and diffraction techniques which guided the exploration. The second phase Ba3Y2[Si2O7]2 was located using iterative algorithmic identification of target compositions. The structure solution of the new compounds was aided by continuous rotation electron diffraction, and the structures were refined against combined synchrotron and neutron time-of-flight powder diffraction. Ba5Y13[SiO4]8O8.5 crystallizes in I4̄2m, a = 18.92732(1), c = 5.357307(6) Å and represents its own structure type which combines elements of structures of known silicates embedded in columns of interconnected yttrium-centred polyhedra characteristic of high-pressure phases. Ba3Y2[Si2O7]2 has P21 symmetry with a pseudo-tetragonal cell (a = 16.47640(4), b = 9.04150(5), c = 9.04114(7) Å, β = 90.0122(9)°) and is a direct superstructure of the Ca3BaBi[P2O7]2 structure. Despite the lower symmetry, the structure of Ba3Y2[Si2O7]2 retains disorder in both Ba/Y sites and disilicate network, thus presenting a superposition of possible locally-ordered fragments. Ba5Y13[SiO4]8O8.5 has low thermal conductivity of 1.04(5) W m-1 K-1 at room temperature. The two discovered phases provide a rich structural platform for further functional material design. The interplay of automated unknown phase composition identification with multiple diffraction methods offers acceleration of the time-consuming exploration of high-dimensional chemical spaces for new structures.

Challenges Reconciling Theory and Experiments in the Prediction of Lattice Thermal Conductivity: The Case of Cu-Based Sulvanites.

The exploration of large chemical spaces in search of new thermoelectric materials requires the integration of experiments, theory, simulations, and data science. The development of high-throughput strategies that combine DFT calculations with machine learning has emerged as a powerful approach to discovering new materials. However, experimental validation is crucial to confirm the accuracy of these workflows. This validation becomes especially important in understanding the transport properties that govern the thermoelectric performance of materials since they are highly influenced by synthetic, processing, and operating conditions. In this work, we explore the thermal conductivity of Cu-based sulvanites by using a combination of theoretical and experimental methods. Previous discrepancies and significant variations in reported data for Cu3VS4 and Cu3VSe4 are explained using the Boltzmann Transport Equation for phonons and by synthesizing well-characterized defect-free samples. The use of machine learning approaches for extracting high-order force constants opens doors to charting the lattice thermal conductivity across the entire Cu-based sulvanite family-finding not only materials with κ l values below 2 W m-1 K-1 at moderate temperatures but also rationalizing their thermal transport properties based on chemical composition.

Conformational Space Profiling Enhances Generic Molecular Representation for AI-Powered Ligand-Based Drug Discovery.

The molecular representation model is a neural network that converts molecular representations (SMILES, Graph) into feature vectors, and is an essential module applied across a wide range of artificial intelligence-driven drug discovery scenarios. However, current molecular representation models rarely consider the three-dimensional conformational space of molecules, losing sight of the dynamic nature of small molecules as well as the essence of molecular conformational space that covers the heterogeneity of molecule properties, such as the multi-target mechanism of action, recognition of different biomolecules, dynamics in cytoplasm and membrane. In this study, a new model named GeminiMol is proposed to incorporate conformational space profiles into molecular representation learning, which extracts the feature of capturing the complicated interplay between the molecular structure and the conformational space. Although GeminiMol is pre-trained on a relatively small-scale molecular dataset (39290 molecules), it shows balanced and superior performance not only on 67 molecular properties predictions but also on 73 cellular activity predictions and 171 zero-shot tasks (including virtual screening and target identification). By capturing the molecular conformational space profile, the strategy paves the way for rapid exploration of chemical space and facilitates changing paradigms for drug design.

Identification of novel and potent inhibitors of SARS-CoV-2 main protease from DNA-encoded chemical libraries.

In vitro screening of large compound libraries with automated high-throughput screening is expensive and time-consuming and requires dedicated infrastructures. Conversely, the selection of DNA-encoded chemical libraries (DECLs) can be rapidly performed with routine equipment available in most laboratories. In this study, we identified novel inhibitors of SARS-CoV-2 main protease (Mpro) through the affinity-based selection of the DELopen library (open access for academics), containing 4.2 billion compounds. The identified inhibitors were peptide-like compounds containing an N-terminal electrophilic group able to form a covalent bond with the nucleophilic Cys145 of Mpro, as confirmed by x-ray crystallography. This DECL selection campaign enabled the discovery of the unoptimized compound SLL11 (IC50 = 30 nM), proving that the rapid exploration of large chemical spaces enabled by DECL technology allows for the direct identification of potent inhibitors avoiding several rounds of iterative medicinal chemistry. As demonstrated further by x-ray crystallography, SLL11 was found to adopt a highly unique U-shaped binding conformation, which allows the N-terminal electrophilic group to loop back to the S1' subsite while the C-terminal amino acid sits in the S1 subsite. MP1, a close analog of SLL11, showed antiviral activity against SARS-CoV-2 in the low micromolar range when tested in Caco-2 and Calu-3 (EC50 = 2.3 µM) cell lines. As peptide-like compounds can suffer from low cell permeability and metabolic stability, the cyclization of the compounds will be explored in the future to improve their antiviral activity.

Machine Learning Guided AQFEP: A Fast and Efficient Absolute Free Energy Perturbation Solution for Virtual Screening.

Structure-based methods in drug discovery have become an integral part of the modern drug discovery process. The power of virtual screening lies in its ability to rapidly and cost-effectively explore enormous chemical spaces to select promising ligands for further experimental investigation. Relative free energy perturbation (RFEP) and similar methods are the gold standard for binding affinity prediction in drug discovery hit-to-lead and lead optimization phases, but have high computational cost and the requirement of a structural analog with a known activity. Without a reference molecule requirement, absolute FEP (AFEP) has, in theory, better accuracy for hit ID, but in practice, the slow throughput is not compatible with VS, where fast docking and unreliable scoring functions are still the standard. Here, we present an integrated workflow to virtually screen large and diverse chemical libraries efficiently, combining active learning with a physics-based scoring function based on a fast absolute free energy perturbation method. We validated the performance of the approach in the ranking of structurally related ligands, virtual screening hit rate enrichment, and active learning chemical space exploration; disclosing the largest reported collection of free energy simulations to date.

Molecular similarity: Theory, applications, and perspectives

Molecular similarity pervades much of our understanding and rationalization of chemistry. This has become particularly evident in the current data-intensive era of chemical research, with similarity measures serving as the backbone of many Machine Learning (ML) supervised and unsupervised procedures. Here, we present a discussion on the role of molecular similarity in drug design, chemical space exploration, chemical “art” generation, molecular representations, and many more. We also discuss more recent topics in molecular similarity, like the ability to efficiently compare large molecular libraries.

Open Source Code Contributions to Global Health: The Case of Antimalarial Drug Discovery.

The discovery of treatments for infectious diseases that affect the poorest countries has been stagnant for decades. As long as expected returns on investment remain low, pharmaceutical companies' lack of interest in this disease area must be compensated for with collaborative efforts from the public sector. New approaches to drug discovery, inspired by the "open source" philosophy prevalent in software development, offer a platform for experts from diverse backgrounds to contribute their skills, enhancing reproducibility, progress tracking, and public discussion. Here, we present the first efforts of Ersilia, an initiative focused on attracting data scientists into contributing to global health, toward meeting the goals of Open Source Malaria, a consortium of medicinal chemists investigating antimalarial compounds using a purely open science approach. We showcase the chemical space exploration of a set of triazolopyrazine compounds with potent antiplasmodial activity and discuss how open source practices can serve as a common ground to make drug discovery more inclusive and participative.

Artificial Intelligence Methods and Models for Retro-Biosynthesis: A Scoping Review.

Retrosynthesis aims to efficiently plan the synthesis of desirable chemicals by strategically breaking down molecules into readily available building block compounds. Having a long history in chemistry, retro-biosynthesis has also been used in the fields of biocatalysis and synthetic biology. Artificial intelligence (AI) is driving us toward new frontiers in synthesis planning and the exploration of chemical spaces, arriving at an opportune moment for promoting bioproduction that would better align with green chemistry, enhancing environmental practices. In this review, we summarize the recent advancements in the application of AI methods and models for retrosynthetic and retro-biosynthetic pathway design. These techniques can be based either on reaction templates or generative models and require scoring functions and planning strategies to navigate through the retrosynthetic graph of possibilities. We finally discuss limitations and promising research directions in this field.

Designing of potential polymer donors with lower reorganization energy and chemical space exploration in search of materials for efficient organic solar cells

Machine learning with fascinating material predicting capabilities has emerged as a growing throughput artificial intelligence technology to accelerate the design, screening and prediction of material properties. In this study, chemical descriptors are employed to train machine learning (ML) models for the designing of potential polymer donor materials aimed at enhancing the efficiency of organic solar cells (OSCs). To achieve this, ten ML models are evaluated. Among all, the random forest regressor model has shown better predictive ability. In addition, to quantify the importance of features, the numerical scores of the significance of the features are determined through Feature importance analysis. Moreover, 10k new polymers are generated employing the python-based method and their reorganization energies are further predicted using a pre-trained ML model. The assemblage of polymers is done on the basis of the reorganization energy, those with lower reorganization energy being preferred. Furthermore, the synthetic accessibility and thus the structural diversity among the chosen polymers donors are examined, revealing a structural difference among the chosen polymers. This work paves a way for the designing and screening of efficient semiconductors with low reorganization energy to increase the efficiency of OSCs.

Exploration of Rhenium Bisquinoline Tricarbonyl Complexes for their Antibacterial Properties.

Metal complexes have emerged as a promising source for novel classes of antibacterial agents to combat the rise of antimicrobial resistance around the world. In the exploration of the transition metal chemical space for novel metalloantibiotics, the rhenium tricarbonyl moiety has been identified as a promising scaffold. Here we have prepared eight novel rhenium bisquinoline tricarbonyl complexes and explored their antibacterial properties. Significant activity against both Gram-positive and Gram-negative bacteria was observed. However, all complexes also showed significant toxicity against human cells, putting into question the prospects of this specific rhenium compound class as metalloantibiotics. To better understand their biological effects, we conduct the first mode of action studies on rhenium bisquinoline complexes and show that they are able to form pores through bacterial membranes. Their straight-forward synthesis and tuneability suggests that further optimisation of this compound class could lead to compounds with enhanced bacterial specificity.

User-friendly and industry-integrated AI for medicinal chemists and pharmaceuticals

Artificial intelligence has brought crucial changes to the whole field of natural sciences. Myriads of machine learning algorithms have been developed to facilitate the work of experimental scientists. Molecular property prediction and drug synthesis planning become routine tasks. Moreover, inverse design of compounds with tunable properties as well as on-the-fly autonomous process optimization and chemical space exploration became possible in silico. Affordable robotic platforms exist able to perform thousands of experiments every day, analyzing the results and tuning the protocols. Despite this, most of these developments get trapped at the stage of code or overlooked, limiting their use by experimental scientists. Meanwhile, visibility and the number of user-friendly tools and technologies available to date is too low to compensate for this fact, rendering the development of novel therapeutic compounds inefficient. In this Review, we set the goal to bridge the gap between modern technologies and experimental scientists to improve drug development efficacy. Here we survey advanced and easy-to-use technologies able to help medical chemists at every stage of their research, including those integrated in technological processes during COVID-19 pandemic motivated by the need for fast yet precise solutions. Moreover, we review how these technologies are integrated by industry and clinics to streamline drug development and production. These technologies already transform the current paradigm of scientific thinking and revolutionize not only medicinal chemistry, but the whole field of natural sciences.

Navigating pharmacophore space to identify activity discontinuities: A case study with BCR-ABL.

The exploration of chemical space is a fundamental aspect of chemoinformatics, particularly when one explores a large compound data set to relate chemical structures with molecular properties. In this study, we extend our previous work on chemical space visualization at the pharmacophoric level. Instead of using conventional binary classification of affinity (active vs inactive), we introduce a refined approach that categorizes compounds into four distinct classes based on their activity levels: super active, very active, active, and inactive. This classification enriches the color scheme applied to pharmacophore space, where the color representation of a pharmacophore hypothesis is driven by the associated compounds. Using the BCR-ABL tyrosine kinase as a case study, we identified intriguing regions corresponding to pharmacophore activity discontinuities, providing valuable insights for structure-activity relationships analysis.

Chemical space exploration with Molpher: Generating and assessing a glucocorticoid receptor ligand library.

Computational exploration of chemical space is crucial in modern cheminformatics research for accelerating the discovery of new biologically active compounds. In this study, we present a detailed analysis of the chemical library of potential glucocorticoid receptor (GR) ligands generated by the molecular generator, Molpher. To generate the targeted GR library and construct the classification models, structures from the ChEMBL database as well as from the internal IMG library, which was experimentally screened for biological activity in the primary luciferase reporter cell assay, were utilized. The composition of the targeted GR ligand library was compared with a reference library that randomly samples chemical space. A random forest model was used to determine the biological activity of ligands, incorporating its applicability domain using conformal prediction. It was demonstrated that the GR library is significantly enriched with GR ligands compared to the random library. Furthermore, a prospective analysis demonstrated that Molpher successfully designed compounds, which were subsequently experimentally confirmed to be active on the GR. A collection of 34 potential new GR ligands was also identified. Moreover, an important contribution of this study is the establishment of a comprehensive workflow for evaluating computationally generated ligands, particularly those with potential activity against targets that are challenging to dock.

Regioselective fluorination of allenes enabled by I(I)/I(III) catalysis

The prominence and versatility of propargylic fluorides in medicinal chemistry, coupled with the potency of F/H and F/OH bioisosterism, has created a powerful impetus to develop efficient methods to facilitate their construction. Motivated by the well-established conversion of propargylic alcohols to allenes, an operationally simple, organocatalysis-based strategy to process these abundant unsaturated precursors to propargylic fluorides would be highly enabling: this would consolidate the bioisosteric relationship that connects propargylic alcohols and fluorides. Herein, we describe a highly regioselective fluorination of unactivated allenes based on I(I)/I(III) catalysis in the presence of an inexpensive HF source that serves a dual role as both nucleophile and Brønsted acid activator. This strategy enables a variety of secondary and tertiary propargylic fluorides to be prepared: these motifs are prevalent across the bioactive small molecule spectrum. Facile product derivatisation, concise synthesis of multi-vicinal fluorinated products together with preliminary validation of enantioselective catalysis are disclosed. The expansive potential of this platform is also demonstrated through the highly regioselective organocatalytic oxidation, chlorination and arylation of allenes. It is envisaged that the transformation will find application in molecular design and accelerate the exploration of organofluorine chemical space.

Rapid Estimation of the Intermolecular Electronic Couplings and Charge-Carrier Mobilities of Crystalline Molecular Organic Semiconductors through a Machine Learning Pipeline.

Organic semiconductors (OSC) offer tremendous potential across a wide range of (opto)electronic applications. OSC development, however, is often limited by trial-and-error design, with computational modeling approaches deployed to evaluate and screen candidates through a suite of molecular and materials descriptors that generally require hours to days of computational time to accumulate. Such bottlenecks slow the pace and limit the exploration of the vast chemical space comprising OSC. When considering charge-carrier transport in OSC, a key parameter of interest is the intermolecular electronic coupling. Here, we introduce a machine learning (ML) model to predict intermolecular electronic couplings in organic crystalline materials from their three-dimensional (3D) molecular geometries. The ML predictions take only a few seconds of computing time compared to hours by density functional theory (DFT) methods. To demonstrate the utility of the ML predictions, we deploy the ML model in conjunction with mathematical formulations to rapidly screen the charge-carrier mobility anisotropy for more than 60,000 molecular crystal structures and compare the ML predictions to DFT benchmarks.

Asymmetric Carbene Transfer: Enhancing Chemical Diversity for Drug Discovery.

The quest to explore chemical space is vital for identifying novel disease targets, impacting both the effectiveness and safety profile of therapeutic agents. The tangible chemical space, currently estimated at a conservative 108 synthesized compounds, pales in comparison to the theoretically conceivable diversity of 1060 molecules. To bridge this vast gap, organic chemists are spearheading innovative methodologies that promise to broaden this limited chemical diversity. A beacon of this progressive wave is Asymmetric Carbene Transfer (ACT), a burgeoning strategy that significantly boosts molecular diversity with efficient bond-formation and precise chiral control. This review focuses on the capabilities of ACT in creating pharmaceutically significant molecules, encompassing an array of natural products and bioactive compounds. Through the lens of ACT, we discern its substantial influence on drug discovery, paving the way for novel therapeutic avenues by expanding the boundaries of molecular diversity. This review will shed light on prospective methodological developments of ACT and articulate their conceivable contributions to the medicinal chemistry arena.

Extensive targeting of chemical space at the prime side of ketoamide inhibitors of rhomboid proteases by branched substituents empowers their selectivity and potency

Rhomboid intramembrane serine proteases have been implicated in several pathologies, and emerge as attractive pharmacological target candidates. The most potent and selective rhomboid inhibitors available to date are peptidyl α-ketoamides, but their selectivity for diverse rhomboid proteases and strategies to modulate it in relevant contexts are poorly understood. This gap, together with the lack of suitable in vitro models, hinders ketoamide development for relevant eukaryotic rhomboid enzymes. Here we explore the structure-activity relationship principles of rhomboid inhibiting ketoamides by medicinal chemistry and enzymatic in vitro and in-cell assays with recombinant rhomboid proteases GlpG, human mitochondrial rhomboid PARL and human RHBDL2. We use X-ray crystallography in lipidic cubic phase to understand the binding mode of one of the best ketoamide inhibitors synthesized here containing a branched terminal substituent bound to GlpG. In addition, to extend the interpretation of the co-crystal structure, we use quantum mechanical calculations and quantify the relative importance of interactions along the inhibitor molecule. These combined experimental analyses implicates that more extensive exploration of chemical space at the prime side is unexpectedly powerful for the selectivity of rhomboid inhibiting ketoamides. Together with variations in the peptide sequence at the non-prime side, or its non-peptidic alternatives, this strategy enables targeted tailoring of potent and selective ketoamides towards diverse rhomboid proteases including disease-relevant ones such as PARL and RHBDL2.

Building Block-Centric Approach to DNA-Encoded Library Design.

DNA-encoded library technology grants access to nearly infinite opportunities to explore the chemical structure space for drug discovery. Successful navigation depends on the design and synthesis of libraries with appropriate physicochemical properties (PCPs) and structural diversity while aligning with practical considerations. To this end, we analyze combinatorial library design constraints including the number of chemistry cycles, bond construction strategies, and building block (BB) class selection in pursuit of ideal library designs. We compare two-cycle library designs (amino acid + carboxylic acid, primary amine + carboxylic acid) in the context of PCPs and chemical space coverage, given different BB selection strategies and constraints. We find that broad availability of amines and acids is essential for enabling the widest exploration of chemical space. Surprisingly, cost is not a driving factor, and virtually, the same chemical space can be explored with "budget" BBs.