Overlap of On-demand Ultra-large Combinatorial Spaces with On-the-shelf Drug-like Libraries

On-demand combinatorial spaces are shifting paradigms in early drug discovery, by considerably increasing the searchable chemical space to several billions of compounds while securing their synthetic accessibility. We here systematically compared the on-the-shelf available drug-like chemical space (9 million compounds) to three on-demand ultra-large (ODUL) combinatorial fragment spaces (REAL, CHEMriya, GalaXi) covering 32 billion of readily accessible molecules. Surprisingly, only one space (REAL) intersects almost entirely the currently available drug-like space, suggesting that it is the only ODUL widely suitable for in-stock hit expansion. Of course, expanding a preliminary ODUL hit in the same chemical space is the best possible strategy to rapidly generate structure-activity relationships. All three spaces remain well suited to early hit finding initiatives since they all provide numerous unique scaffolds that are not described by on-the shelf collections.

Unsupervised Representation Learning for Proteochemometric Modeling

In silico protein–ligand binding prediction is an ongoing area of research in computational chemistry and machine learning based drug discovery, as an accurate predictive model could greatly reduce the time and resources necessary for the detection and prioritization of possible drug candidates. Proteochemometric modeling (PCM) attempts to create an accurate model of the protein–ligand interaction space by combining explicit protein and ligand descriptors. This requires the creation of information-rich, uniform and computer interpretable representations of proteins and ligands. Previous studies in PCM modeling rely on pre-defined, handcrafted feature extraction methods, and many methods use protein descriptors that require alignment or are otherwise specific to a particular group of related proteins. However, recent advances in representation learning have shown that unsupervised machine learning can be used to generate embeddings that outperform complex, human-engineered representations. Several different embedding methods for proteins and molecules have been developed based on various language-modeling methods. Here, we demonstrate the utility of these unsupervised representations and compare three protein embeddings and two compound embeddings in a fair manner. We evaluate performance on various splits of a benchmark dataset, as well as on an internal dataset of protein–ligand binding activities and find that unsupervised-learned representations significantly outperform handcrafted representations.

Protein Ligand Interaction Fingerprints

A most challenging part in docking-based virtual screening is the scoring functions implemented in various docking programs in order to evaluate different poses of the ligands inside the binding cavity of the receptor. Precise and trustable measurement of ligand-protein affinity for Structure-Based Virtual Screening (SB-VS) is therefore, an outstanding problem in docking studies. Empirical post-docking filters can be helpful as a way to provide various types of structure-activity information. Different types of interaction have been presented between the ligands and the receptor so far. Based on the diversity and importance of PLIF methods, this chapter will focus on the comparison of different protocols. The advantages and disadvantages of all methods will be discussed explicitly in this chapter as well as future sights for further progress in this field. Different classifications approaches for the protein-ligand interaction fingerprints were also discussed in this chapter.

Target specific proteochemometric model development for BACE1 – protein flexibility and structural water are critical in virtual screening

Structural water and protein plasticity are important factors for BACE1 targeted ligand virtual screening.

Polypharmacology modelling using proteochemometrics (PCM): recent methodological developments, applications to target families, and future prospects

Proteochemometric (PCM) modelling is a computational method to model the bioactivity of multiple ligands against multiple related protein targets simultaneously.

Encoding protein-ligand interaction patterns in fingerprints and graphs

We herewith present a novel and universal method to convert protein-ligand coordinates into a simple fingerprint of 210 integers registering the corresponding molecular interaction pattern. Each interaction (hydrophobic, aromatic, hydrogen bond, ionic bond, metal complexation) is detected on the fly and physically described by a pseudoatom centered either on the interacting ligand atom, the interacting protein atom, or the geometric center of both interacting atoms. Counting all possible triplets of interaction pseudoatoms within six distance ranges, and pruning the full integer vector to keep the most frequent triplets enables the definition of a simple (210 integers) and coordinate frame-invariant interaction pattern descriptor (TIFP) that can be applied to compare any pair of protein-ligand complexes. TIFP fingerprints have been calculated for ca. 10,000 druggable protein-ligand complexes therefore enabling a wide comparison of relationships between interaction pattern similarity and ligand or binding site pairwise similarity. We notably show that interaction pattern similarity strongly depends on binding site similarity. In addition to the TIFP fingerprint which registers intermolecular interactions between a ligand and its target protein, we developed two tools (Ishape, Grim) to align protein-ligand complexes from their interaction patterns. Ishape is based on the overlap of interaction pseudoatoms using a smooth Gaussian function, whereas Grim utilizes a standard clique detection algorithm to match interaction pattern graphs. Both tools are complementary and enable protein-ligand complex alignments capitalizing on both global and local pattern similarities. The new fingerprint and companion alignment tools have been successfully used in three scenarios: (i) interaction-biased alignment of protein-ligand complexes, (ii) postprocessing docking poses according to known interaction patterns for a particular target, and (iii) virtual screening for bioisosteric scaffolds sharing similar interaction patterns.